(28)

Его решения имеют вид:

                                (29)

или                            (30).

Величина 4μ2/(M1M2) при любых M1, M2 не превосходит единицы, поэтому подкоренное выражение всегда неотрицательно.

1.3. Трехмерный кристалл

Мы рассмотрели колебания в одномерной цепочке. Подобным образом могут быть описаны и колебания решетки трехмерного кристалла.

Предположим, что примитивная ячейка кристалла состоит из l атомов. Каждый атом ячейки будем обозначать индексом s, этот индекс принимает l различных значений. Любой атом кристалла однозначно определяется радиус-вектором , задающим положение ячейки (соответствующего узла решетки Браве), и индексом s, характеризующим положение атома внутри ячейки (тип атома).

Смещение атомов при колебаниях решетки является линейной комбинацией плоских гармонических волн (точнее, их вещественных частей):    (40).

Частота колебаний одинакова для всех атомов кристалла. Амплитуда колебаний зависит от типа атома (индекса s), то есть одинакова для всех однотипных атомов. Направление вектора амплитуды может, вообще говоря, быть каким угодно.

Индекс j обозначает ветвь колебаний. Волновой вектор  и ветвь j однозначно определяют частоту и относительные амплитуды атомов всех типов. Для каждой ветви зависимости  и являются непрерывными функциями.

Если примитивная ячейка кристалла содержит l атомов, то число ветвей равно 3l. Таким образом, каждому значению волнового вектора соответствуют 3l разных колебаний.

Три из этих ветвей – акустические, в предельном случае длинных волн их частота пропорциональна длине волнового вектора ω = |k|. Однако скорость звука  зависит от направления распространения волны, то есть от направления вектора . В случае длинноволновых акустических колебаний амплитуды всех атомов примитивной ячейки примерно одинаковы.

Остальные 3l–3 ветвей – оптические , при их частота отлична от нуля.

По направлению амплитуды относительно волнового вектора акустические колебания можно разделить на продольное и два поперечных. Как правило, скорость звука у продольного колебания больше, чем у поперечных.

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

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

Количество разрешенных волновых векторов в зоне Бриллюэна равно N = V/v0 – числу примитивных ячеек в нормировочном объеме кристалла V = L3 (v0 – объем примитивной ячейки). Действительно, плотность разрешенных волновых векторов в обратном пространстве равна V/(2π)3, т. е. в объеме обратного пространства Δ3k содержится Δ3k· V/(2π)3 разрешенных волновых векторов. Объем зоны Бриллюэна – объем примитивной ячейки обратной решетки — равен (2π)3/v0, и для числа разрешенных состояний получаем (2π)3/v0· V/(2π)3 = V/v0 = N.

Число ветвей — 3l, поэтому полное число колебаний равно 3lN — утроенному числу атомов кристалла в объеме L3, т. е. числу степеней свободы механической системы.

Глава 2. Фононы. Фононный газ

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

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

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

Р. Паерлс в 1929 году ввел в теорию Дебая квантовые  (фононные) явления и показал, что тепловое сопротивление решетки обусловлено взаимодействием фононов. Фонон, в отличии от обычных частиц, может существовать лишь в некоторой среде, которая пребывает в состоянии теплового возбуждения. Нельзя вообразить фонон, который распространялся бы в вакууме, поскольку он описывает квантовый характер тепловых колебаний решетки и навечно замкнут в кристалле. Корпускулярный аспект малых колебаний атомов решетки кристалла приводит к понятию фонона, и распространение упругих тепловых волн в кристалле можно рассматривать как перенесение фононов.

Фононы являются бозе-частицами: число фононов, соответствующих определенному колебанию (число фононов одном состоянии), может быть сколь угодно большим. В состоянии термодинамического равновесия среднее число фононов njk ветви j с волновым вектором  зависит только от энергии фонона (частоты колебания):

                (31).

Здесь k – постоянная Больцмана. С точки зрения квантовой (да и классической) механики, нормальные колебания решетки ведут себя как набор независимых гармонических осцилляторов. Роль координаты осциллятора играет при этом амплитуда колебания, число фононов является уровнем энергии осциллятора.

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

При высоких температурах, kb T >> ħω, число фононов пропорционально температуре:                          (32).

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

Они играют существенную роль при низких температурах. Если k T << ħω, то среднее число фононов экспоненциально мало:

                                  (33).

Можно сказать, что колебания, частота которых превосходит величину kT/ħ, ''вымерзают''. Энергия колебания не может быть меньше энергии одного фонона ħωjk а энергия фонона много больше характерной тепловой энергии kT, поэтому такие колебания практически не возбуждаются.

Глава 3. Акустическая и оптическая ветки колебаний.

Итак, для каждого волнового вектора k, согласно уравнения (30) существуют две частоты ω, удовлетворяющие дисперсионному уравнению. Точнее, есть две непрерывные функции ω(k), которые отличаются знаком перед корнем. Говорят, что существуют две ветви колебаний.

Исследуем обе ветви.

Напомним, что волновые вектора, отличающиеся на вектор обратной решки, описывают одно и то же колебания. (Вследствие этого функция ω(k) периодична с периодом обратной решетки 2π/a, а в трехмерном случае обладает трансляционной симметрией обратной решетки). Поэтому мы считаем, что волновой вектор лежит в пределах первой зоны Бриллюэна: –π/a<k<π/a.

Решение со знаком ''минус''

В точке k = 0:

                                       (34).

На границе зоны Бриллюэна:

                               (35)

При ka<< 1 (длинные волны):

           (36)

Или другими словами:

        (37)

Мы видим, что в длинноволновом пределе закон дисперсии этой ветви линеен, то есть, как и в случае цепочки с одним атомом в примитивной ячейке, описывает акустические колебания. По этой причине вся ветвь (решение со знаком ''–'') называется акустической (рис.3.1).

Выражение для скорости звука имеет такой же вид, что и соответствующее выражение для цепочки с одним атомом в ячейке (20) и зависит от тех же макроскопических характеристик: линейной плотности и упругой постоянной цепочки:

                                       (38).

Линейная плотность двухатомной цепочки равна (M1+M2)/a, а упругая постоянная — γ· a/2 (так как длина одной пружинки в наших обозначениях равна a/2).

Это и неудивительно. Мы уже видели, изучая цепочку с одним атомом в примитивной ячейке, что длинноволновые акустические колебания можно получить, рассматривая цепочку, как непрерывную упругую среду. Атомы ячейки при таких колебаниях движутся вместе, как единое целое, поэтому структура примитивной ячейки не играет роли, а важны лишь макроскопические, усредненные характеристики цепочки.

То, что атомы ячейки при длинноволновых акустических колебаниях движутся вместе, можно получить и непосредственно, решив систему (25). Эта система разрешима, когда ее определитель равен нулю, а определитель равен нулю, когда ω и k связаны законом дисперсии. При этом уравнения системы уже не являются независимыми, и мы можем взять любое из них, чтобы найти отношение амплитуд A и B.

Из первого уравнения системы (25) получаем:

                                            (39),

откуда в пределе длинноволновых акустических колебаний (k→ 0, ω = s |k|→ 0) следует B/A→ 1, т. е. A = B: атомы движутся в фазе с одинаковыми амплитудами.

Отметим также, что на границе зоны Бриллюэна групповая скорость ∂ω/∂ k равна нулю. Это утверждение справедливо для обеих ветвей колебаний.

Решение со знаком ''плюс''.

В точке k = 0:

                                       (40)

На границе зоны Бриллюэна:

                                (41)

Групповая скорость этой ветви ∂ω/∂ k равна нулю как на границе зоны Бриллюэна, так и при k = 0.

Эта ветвь целиком лежит выше акустической ветви: ее минимальная частота  больше максимальной частоты акустических колебаний . Таким образом, в цепочке могут распространяться волны в частотами от 0 до и от до . Интервал частот является ''запрещенной зоной'': волн с такими частотами не существует. Относительная ширина этого интервала тем больше, чем больше отношение масс M2/M1.

Чтобы понять, что представляют собой длинноволновые колебания этой ветви, найдем отношение амплитуд колебаний B/A при k = 0 с помощью (36):

                         (42)

Мы видим, что атомы в каждой ячейке движутся в противофазе, то сближаясь, то удаляясь друг от друга, причем одновременно во всех ячейках (если k = 0). Амплитуда движения легкого атома больше амплитуды тяжелого в M2/M1 раз, т. е. центр тяжести ячейки остается на месте.

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

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

Энергия длинноволнового оптического фонона имеет тот же порядок величины, что и энергия фонона акустического колебания с максимальной частотой, которую мы оценили в 0,05 эВ. Энергии оптических фононов большинства полупроводниковых кристаллов лежат в диапазоне 0,03÷ 0,1 эВ.

Посмотрим теперь, как колеблются атомы, когда длина волны минимальна, т. е. когда волновой вектор лежит на границе зоны Бриллюэна.

В случае акустических колебаний ω2 = 2γ/M2, коэффициент при B во втором уравнении системы (25) обращается в ноль, откуда следует что A = 0.

В случае оптических колебаний ω2 = 2γ/M1, и из первого уравнения (25) следует что B = 0.

Таким образом, при k = π/a в случае акустических волн колеблются тяжелые атомы, а легкие неподвижны, в случае оптических, наоборот: колеблются легкие, тяжелые стоят на месте.

Обобщим теперь полученные результаты. Нетрудно показать, что если примитивная ячейка одномерной цепочки содержит l атомов, то спектр колебаний состоит из l ветвей, одна из которых акустическая, а остальные – оптические.

Мы рассматривали бесконечную цепочку, не накладывая никаких ограничений на длины волн упругих колебаний. В результате, мы пришли к выводу, что в цепочке могут распространяться колебания с любыми волновыми векторами, лежащими в первой зоне Бриллюэна. (Было показано, что из-за дискретности цепочки волновые вектора, отличающиеся на произвольный вектор обратной решетки, описывают одни и те же колебания. Поэтому можно брать волновой вектор из любой зоны Бриллюэна. Естественней всего описывать колебание наименьшим волновым вектором, т.е. вектором из первой зоны Бриллюна.)

Чтобы иметь дело не с непрерывным, а с дискретным набором волновых векторов, можно потребовать, чтобы отклонение атомов от равновесия было периодической функцией: u(xn) = u(xn+L). Иными словами — поставить граничные условия Борна-Кармана. Период L должен быть кратен постоянной решетки цепочки.

Условиям Борна-Кармана удовлетворяют только гармонические колебания с ''разрешенными'' волновыми векторами kn = 2π n/L. Нетрудно подсчитать, что в зоне Бриллюэна размещается L/a разрешенных волновых векторов, т. е. ровно столько, сколько примитивных ячеек укладывается на длине L. (Волновым векторам –π/a и π/a соответствует одно и то же колебание и поэтому считаем эти два значения за одно). Мы уже упоминали выше об этом свойстве зоны Бриллюэна.

Так как колебание однозначно определяется волновым вектором и ветвью, то различных колебаний столько, сколько атомов содержит цепочка. Это общее свойство линейных колебательных систем: количество независимых колебаний (нормальных мод) равно числу степеней свободы системы.

Глава 4. Энергия колебаний и теплоемкость кристаллической решетки

Энергию колебаний и теплоемкость решетки будем рассчитывать для единичного объема кристалла, т. е. положим нормировочный объем равным единице: V = L3 = 1.

Чтобы вычислить среднюю энергию колебаний кристаллической решетки, нужно просуммировать среднюю энергию всех типов колебаний (всех состояний фононов):

                    (43).

Проще всего это сделать при высоких температурах, когда для частот всех колебаний выполняется неравенство ħωjk<< kT (классический предел). Тогда средняя энергия, приходящаяся на каждое колебание, равна kBT, всего колебаний 3lN = 3lN, для полной энергии E получаем:

                                       (44).

Так как N – число примитивных ячеек кристалла в единице объема, то N = 1/v0, где v0 – объем примитивной ячейки.

Теплоемкость решетки при высоких температурах постоянна (закон Дюлонга и Пти): CV = 3lNk                                      (45).

При невысоких температурах все сложнее. Чтобы точно вычислить энергию решетки, то есть сосчитать сумму (45), необходимо знать дисперсионные зависимости для всех ветвей колебаний. И даже при условии, что зависимости эти известны, аналитическое выражение для энергии получить практически невозможно.

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

4.1. Модель Эйнштейна

В модели Эйнштейна предполагается, что частоты всех фононов одинаковы:         ωjk = ω1                                         (46).

Тогда для энергии получаем:

                             (47).

При высоких температурах, kBT>>ħω1, эта зависимость приводит к выражению (45) для энергии и закону Дюлонга и Пти (46) для теплоемкости.

При низких температурах, kT<<ħω1, энергия колебаний и теплоемкость экспоненциально уменьшаются:

Модель Эйнштейна хорошо описывает вклад в энергию и теплоемкость оптических ветвей фононов, у которых частота слабо зависит от волнового вектора и ее можно считать постоянной. Чтобы учесть только оптические ветви, частоту которых мы полагаем равной ω1, нужно вместо 3l писать число этих ветвей. В общем случае, частоты разных оптических ветвей могут сильно отличаться друг от друга и их вклад в энергию и теплоемкость нужно учитывать отдельно.

4.2. Модель Дебая

Опыт показывает, что теплоемкость действительно падает с уменьшением температуры, но не экспоненциально, а пропорционально T3. Дело в том, что при любых, сколь угодно низких температурах в кристалле найдутся колебания, энергия фонона которых меньше kBT. Это – длинноволновые акустические колебания. Именно такие колебания, точнее те из них, частота которых меньше kBT/ħ, вносят основной вклад в энергию при низких температурах. Колебания с большими частотами (оптические и более коротковолновые акустические) ''заморожены'': фононов этих колебаний экспоненциально мало.

Сделаем простую оценку. Вклад в энергию вносят фононы, энергия которых меньше kT. Пусть скорость звука j-й акустической ветви равна j и не зависит от направления волнового вектора: ω = j|k|. Тогда вклад в энергию дают колебания с волновыми векторами, меньшими kmax = kBT/(ħ j). Плотность разрешенных значений волновых векторов в k-пространстве кристалла равна V/(2π)3, поэтому внутри сферы радиуса kmax содержится  разрешенных значений волновых векторов. Это число колебаний одной акустической ветви, вносящих существенный вклад в энергию. На каждое такое колебание приходится энергия порядка kT. Для энергии колебаний одной акустической ветви получаем:

                                  (50).

Так как мы вычисляем энергию и теплоемкость единицы объема кристалла, то в (50) мы положили V = 1.

Таким образом, вклад одной акустической ветви в теплоемкость пропорционален T3:

                                      (51).

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

                                         (52),

где через j обозначена скорости звука j-й акустической ветви.

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

Посчитаем теперь энергию решетки при низких температурах более аккуратно.

Формула (44) имеет вид суммы по различным колебаниям (различным состояниям фононов) определенной величины, которая зависит только от энергии фонона:

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

                                    (54).

Здесь - плотность состояний фононов. Напомним, что  - это число состояний квазичастиц (фононов) в единице объема с энергиями от до , то есть число различных колебаний с такими энергиями.

Суммарная плотность состояний складывается из плотности состояний разных ветвей: ; плотность состояний ветви определяется ее законом дисперсии . Аналитически получить законы дисперсии и плотности состояний фононов реальных кристаллов практически невозможно.

Однако при низких температурах энергия и теплоемкость определяются длинноволновыми акустическими фононами. Плотность состояний акустических фононов нам известна, мы получили ее в качестве примера, когда вводили само понятие плотности состояний . Если для j-й акустической ветви ω = j|k|, то:

                                           (55).

Плотность состояний длинноволновых колебаний всех акустических ветвей получается суммированием по трем акустическим ветвям:

 (56), где  - ''усредненная'' скорость звука:

                                                (57).

Линейный закон дисперсии ω = |k| и соответствующая плотность состояний верны только для малых k. При больших значениях волнового вектора закон дисперсии и плотность состояний имеют более сложный вид.

Однако при низких температурах вклад в энергию и теплоемкость вносят как раз только длинноволновые фононы, а при высоких температурах вид плотности состояний не важен, так как в этом случае на каждое колебание приходится энергия kT. Чтобы получить выражение, которое давало бы правильные предельные зависимости при низких и высоких температурах, Дебай предложил считать, что закон дисперсии ω = |k| выполняется и при больших k. Максимальное значение волнового вектора kD при этом выбирается так, чтобы в шаре радиуса kD содержалось столько разрешенных значений волновых векторов, сколько их содержится в зоне Бриллюэна, N = 1/v0. Иными словами, объем этого шара должен быть равен объему зоны Бриллюэна (2π)3/v0, откуда

                                            (58).

Таким образом, сохраняя число акустических колебаний, мы заменяем первую зону Бриллюэна сферой, а реальный закон дисперсии – линейным. Фонон с волновым вектором kD имеет энергию . Соответствующая  температура:

 (59), называется температурой Дебая.

В таком приближении мы можем вычислить вклад акустических ветвей в энергию и теплоемкость решетки:

При низких температурах, T<<θ, верхний предел интеграла много больше единицы. Благодаря экспоненте в знаменателе интеграл сходится очень быстро, что позволяет положить верхний предел равным бесконечности. Значение такого интеграла известно:         (61).

Для энергии акустических колебаний при низких температурах получаем:

                                   (62)

Откуда следует, что теплоемкость решетки при низких температурах пропорциональна T3:

                              (63).

При высоких температурах, T>>θ, верхний предел интегрирования мал, поэтому можно считать, что exp(x)–1≈ x, таким образом:

                                             (64).

Тогда: E = 3NkT и CV = 3Nk.

Это закон Дюлонга и Пти, только вместо полного числа колебаний 3lN стоит число колебаний акустических ветвей 3N. (При высоких температурах на каждое колебание приходится средняя энергия kT, полное число акустических колебаний равно 3N, поэтому вклад акустических ветвей в энергию равен 3NkT).

В пределе низких и высоких температур модель Дебая дает точные значения для вклада акустических ветвей в энергию и теплоемкость. В области же промежуточных температур, T~θ, эта модель лишь аппроксимирует реальную зависимость энергии и теплоемкости от температуры.

Температура Дебая разделяет две температурные области. В области низких температур на энергию и теплоемкость решетки сильное влияние оказывают квантовые эффекты (''вымерзание'' высокочастотных колебаний). В области высоких температур эти эффекты не существенны, и теплоемкость может быть вычислена в классическом приближении. Для большинства кристаллов температура Дебая лежит в интервале от 100 до 300K.

Чтобы получить полную энергию и теплоемкость кристаллической решетки, надо к вкладу акустических колебаний прибавить вклад оптических ветвей, для которого хорошим приближением является модель Эйнштейна. Этот вклад пренебрежимо мал при низких температурах. При высоких температурах вклады всех ветвей в энергию и теплоемкость равны.

Выводы

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

Акустические колебания играют основную роль в определении тепловых свойств кристаллов – теплоемкости, теплопроводности, теплового расширения.

В качестве примера рассмотрено нахождение теплоемкости тел при различных температурах и вклад оптических и акустических колебаний в величину теплоемкости при различных температурах.

СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ

1.                Ансельм А. И. Введение в теорию полупроводников. – М.: Мир, 1965. – 588 с.

2.                 Басс Ф. Г. Электроны и фононы в ограниченных полупроводниках. – М.: Наука, 1984. – 287 с.

3.                 Дущенко В. П., Кучерук И. М. Общая физика. – К.: Высшая школа, 1995. – 430 с.

4.                Епифанов Г. И. Физические основы микроэлектроники. М.:
Советское радио, 1971, 374 с.

5.                 Зисман Г. А., Тодес О. М. Курс общей физики. В 3 т. – М.: Наука, 1995. – 343 с.  

6.                 Кухлинг Х. Справочник по физике: Пер. с нем. – М.: Мир, 1983. – 520 с.

7.                Случинская И. А. Основы материаловедения и технологи
полупроводников. М.: Либрус, 2002, 376 с.

8.                Харрисон У. Теория твёрдого тела. – М. :Мир. – 1978. – 616 с.

9.                Шалимова К. В. Физика полупроводников. М.: Энергия, 1976, 417 с.

10.             Яворский Б. М., Детлаф А. А. Справочник по физике. – М.: Наука, 1982. – 846 с.

Страницы: 1, 2