Фононы являются бозе-частицами: число фононов, соответствующих определенному колебанию (число фононов одном состоянии), может быть сколь угодно большим. В состоянии термодинамического равновесия среднее число фононов 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).
Рис. 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: атомы движутся в фазе с одинаковыми амплитудами.
Рис. 3.2. Амплитуды атомов цепочки в случае длинноволновых акустических колебаний. |
Отметим также, что на границе зоны Бриллюэна групповая скорость ∂ω/∂ k равна нулю. Это утверждение справедливо для обеих ветвей колебаний.
Решение со знаком ''плюс''.
В точке k = 0:
(40)
На границе зоны Бриллюэна:
(41)
Групповая скорость этой ветви ∂ω/∂ k равна нулю как на границе зоны Бриллюэна, так и при k = 0.
Эта ветвь целиком лежит выше акустической ветви: ее минимальная частота больше максимальной частоты акустических колебаний . Таким образом, в цепочке могут распространяться волны в частотами от 0 до и от до . Интервал частот является ''запрещенной зоной'': волн с такими частотами не существует. Относительная ширина этого интервала тем больше, чем больше отношение масс M2/M1.
Чтобы понять, что представляют собой длинноволновые колебания этой ветви, найдем отношение амплитуд колебаний B/A при k = 0 с помощью (36):
(42)
Мы видим, что атомы в каждой ячейке движутся в противофазе, то сближаясь, то удаляясь друг от друга, причем одновременно во всех ячейках (если k = 0). Амплитуда движения легкого атома больше амплитуды тяжелого в M2/M1 раз, т. е. центр тяжести ячейки остается на месте.
Рис. 3.3. Амплитуды атомов цепочки в случае длинноволновых оптических колебаний. |
Если атомы заряжены, то при колебаниях такого типа каждая ячейка представляет собой переменный дипольный момент. Дипольные моменты взаимодействуют с электромагнитным полем, и колебания легко возбуждаются электромагнитными волнами соответствующих частот. В связи с этим, вся ветвь колебаний называется оптической.
При длинноволновых акустических колебаниях атомы ячейки движутся в фазе и никакого дипольного момента не возникает. Поэтому акустические колебания с электромагнитным полем взаимодействуют слабо.
Энергия длинноволнового оптического фонона имеет тот же порядок величины, что и энергия фонона акустического колебания с максимальной частотой, которую мы оценили в 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 соответствует одно и то же колебание и поэтому считаем эти два значения за одно). Мы уже упоминали выше об этом свойстве зоны Бриллюэна.
Так как колебание однозначно определяется волновым вектором и ветвью, то различных колебаний столько, сколько атомов содержит цепочка. Это общее свойство линейных колебательных систем: количество независимых колебаний (нормальных мод) равно числу степеней свободы системы.