Слайд 2
Структура раздела
Общие замечания
Описание движения частиц в т.т.
Гармоническое приближение
Выражение
для смещений как функция времени и координат
Закон дисперсии
Зоны Бриллюэна
Взаимодействие
с Э.М. полем
Теплоёмкость кристаллической решётки
Слайд 3
Структура раздела
Теплоёмкость кристаллической решётки
Модели Дебая и Эйнштейна
Плотность колебательных
состояний и фактор Дебая-Уоллера
Температура плавления. Формула Линдемана
Тепловое расширение
Теплопроводность
Слайд 4
Общие замечания
Следующий шаг в изучении механических свойств т.т.
Учёт
дискретной структуры вещества
Учёт квантования энергии колебаний
Существующие теоретические подходы имеют
свои ограничения (гармоническое приближение, взаимодействие между ближайшими соседями и т.д.)
Теряется информация о непосредственной связи между механическим воздействием и откликом системы
Слайд 5
Общие замечания
Дискретный характер строения вещества оказывает влияние на
свойства деформационных колебаний в кристалле
Когда длина волны становится сравнимой
с межатомным расстоянием, изменяется зависимость ω(k) (закон дисперсии)
Скорость распространения колебаний становится функцией волнового вектора
Слайд 6
Общие замечания
Квантование колебаний приводит к тому, что теплоёмкость
т.т. Стремится к нулю при Т0
Оно так же приводит
к особенностям взаимодействия фононов с материальными частицами (нейтроны, электроны) и электромагнитными волнами
Эти особенности заключаются в существовании неупругого рассеяния, когда происходит рождение или уничтожение кванта колебаний среды. При этом наблюдается скачкообразное изменение характеристик потока частиц, взаимодействующих с твёрдым телом
Слайд 7
Описание движения частиц в т.т.
Функцию, описывающую колебания решётки
можно получить как решение уравнений движения:
Fs=Müs, где Fs -
сила действующая на атом плоскости s, M – масса атома, üs его ускорение
Слайд 8
Описание движения частиц в т.т.
В приближении закона Гука:
Можно
показать, что для системы из двух атомов силовая постоянная
связана с потенциалом взаимодействия U:
a
Слайд 9
Описание движения частиц в т.т.
Использование закона Гука соответствует
гармоническому приближению
Существует и другой подход к составлению уравнений движения:
Ws,p
– тензорная величина. Имеет смысл силы, действующей на частицу s при смещении частицы р на up.
Слайд 10
Описание движения частиц в т.т.
Решение уравнения движения ищем
в виде:
Подставив это выражение в уравнение движения и учитывая,
что Cp=C-p, получим закон дисперсии:
Слайд 11
Описание движения частиц в т.т.
Анализ закона дисперсии показывает,
что при малых K, ω≈const·K=vsK, где vs – скорость
звука
При K→±π/a, ω→const
Если учитывать только взаимодействие между соседними атомами, то можно получить:
Слайд 12
Описание движения частиц в т.т.
Область независимых значений волнового
вектора K:
Эта область называется (первой) зоной Бриллюэна
Значения K, лежащие
за её пределами, можно привести к значениям, лежащим в первой зоне, прибавляя (вычитая) nπ, где n – целое число. Эти значения являются физически идентичными
Слайд 13
Описание движения частиц в т.т.
Закон дисперсии фононов можно
определить экспериментально по рассеянию нейтронов
Зная закон дисперсии, можно вычислить
силовые постоянные Cp:
Установлено, что в металлах межатомные силы могут быть достаточно дальнодействующими (р~20)
Слайд 14
Описание движения частиц в т.т.
Аналогичным образом можно провести
анализ для поперечных колебаний
Во всех формулах будут отличаться только
значения силовых постоянных и подразумеваться смещение в направлении перпендикулярном волновому вектору
Слайд 15
Описание движения частиц в т.т.
Если в кристаллической решётке
содержится больше одного атома, то в спектре колебаний возникает
новая особенность
Появляются оптические ветви колебаний
Слайд 16
Взаимодействие с Э.М. полем
Оптические колебания имеют иной закон
дисперсии, чем акустические. У них ω(0)≠0
Если атомы, входящие в
элементарную ячейку несут избыточный заряд, то при их колебаниях возникают колебания дипольного момента. Это приводит к излучению электромагнитных волн
С другой стороны, электромагнитное излучение может приводить к возбуждению колебаний решётки
Слайд 17
Взаимодействие с Э.М. полем
Частоты фотонов, взаимодействующих с колебаниями
решётки лежат в инфракрасной области
ИК спектроскопия является важным методом
исследования вещества
Слайд 18
Взаимодействие с Э.М. полем
Если в элементарной ячейке содержится
n атомов, то возникает 3n ветвей колебаний. 3 из
них акустические. Остальные – оптические
В кристаллах содержащих дефекты могут возникать дополнительные (локальные) колебания. Они могут так же проявляться в оптических спектрах т.т.
Слайд 19
Описание движения частиц в т.т.
Использованное выше гармоническое приближение
подразумевало разложение потенциальной энергии как функции координат атомов в
ряд по малым смещениям этих атомов из положений равновесия
Слайд 20
Описание движения частиц в т.т.
Нулевой член ряда от
смещений не зависит и на результаты не влияет
Первый член
ряда, линейный по смещениям, в точности равен нулю, т.к. рассматривается состояние вблизи равновесия
Разложение ограничивается квадратичным слагаемым
Слайд 21
Описание движения частиц в т.т.
Важной особенностью гармонического приближения
является представление колебаний кристаллической решётки в виде суперпозиции невзаимодействующих
между собой колебательных мод
Математически этот результат следует из того факта, что функция Гамильтона, описывающая колебания, является положительно определённой квадратичной формой
Слайд 22
Описание движения частиц в т.т.
С помощью преобразований переменных
такую форму можно привести к сумме слагаемых, не содержащих
перекрёстных членов, а только квадраты смещений и импульсов (диагонализация)
Уравнения движения можно получить из функции Гамильтона. Если она приведена к диагональному виду, то получается несколько уравнений движения, зависящих каждое от одной координаты
Слайд 23
Описание движения частиц в т.т.
Такие координаты называются нормальными
Недостатки
этого подхода:
Отсутствует механизм установления теплового равновесия
Исчезает эффект теплового расширения
Нельзя
описать процесс теплопроводности
Теплоёмкость не зависит от типа термодинамического процесса
Слайд 24
Теплоёмкость кристаллической решётки
Различают теплоёмкости CP и CV
В экспериментах
определяют CP, в теоретических расчётах – CV
разница между ними
невелика: CP-CV=9α2BVΔT, где α – температурный коэффициент линейного расширения, V – объём, В – модуль всестороннего сжатия
Слайд 25
Теплоёмкость кристаллической решётки
Основные экспериментальные факты:
При комнатной температуре теплоёмкости
твёрдых тел близки к 3NkB, т.е. 25 Дж/(моль·К)
Вблизи Т=0
теплоёмкость диэлектриков пропорциональна Т3, а металлов – Т
Слайд 26
Теплоёмкость кристаллической решётки
В состоянии теплового равновесия число фононов
с частотой ω определяется с помощью формулы Планка:
Энергия колебаний
с частотой ω: Eω=ħω
Слайд 27
Теплоёмкость кристаллической решётки. Модель Эйнштейна
Модель Эйнштейна: энергия Е
системы N осцилляторов с частотой ω равна сумме их
энергий
Теплоёмкость:
Слайд 28
Теплоёмкость кристаллической решётки. Модель Эйнштейна
При высоких температурах CV→3NkB
– закон Дюлонга и Пти
При низких температурах: CV~exp(-ħω/kBT)
Слайд 29
Теплоёмкость кристаллической решётки
Более сложная модель:
Имеются осцилляторы с различными
частотами ω(k):
Слайд 30
Теплоёмкость кристаллической решётки
Теплоёмкость находится дифференцированием энергии по температуре
Таким
образом, надо знать функцию плотности (колебательных) состояний D(ω)
Слайд 31
Теплоёмкость кристаллической решётки
Нахождение D(ω)
Представим одномерный кристалл как ограниченную
цепочку атомов длины L
Потребуем, чтобы в его объёме укладывалось
целое число волн. Тогда, допустимые значения k=n2π/L, где n=0,1,….
Из-за дискретности структуры вещества существует верхнее ограничение на k и, следовательно, на n. k≤π/a=(N-1)π/L, где a – постоянная решётки, а N – число атомов
Слайд 32
Теплоёмкость кристаллической решётки
Имеем (N-1) колебаний приходящихся на интервал
0≤k≤(N-1)π/L с равномерной плотностью dNk/dk=L/π, и некоторый закон дисперсии
ω(k)
Слайд 33
Теплоёмкость кристаллической решётки
Слайд 34
Теплоёмкость кристаллической решётки. Приближение Дебая
Слайд 35
Теплоёмкость кристаллической решётки. Приближение Дебая
Мы рассмотрели линейную цепочку
Для
трёхмерного кристалла выкладки проводятся аналогично
Для каждой моды звуковых колебаний
получим:
Слайд 36
Теплоёмкость кристаллической решётки
В приближении Эйнштейна: D(ω)=Nδ(ωωE)
Слайд 37
Теплоёмкость кристаллической решётки
Для практических целей выбирают некоторую дебаевскую
частоту ωD, которая для данного конкретного вещества позволяет наилучшим
образом согласовать теоретическую зависимость с экспериментальной зависимостью теплоёмкости от температуры
Эти значения приводятся в справочниках
Температура Дебая определяется из соотношения: ħωD=kBTD
Слайд 38
Теплоёмкость кристаллической решётки
Приближение Дебая относительно хорошо работает для
структур не обладающих оптическими колебаниями
Для оптических колебаний лучше работает
модель Эйнштейна
Слайд 39
Дифракция на кристалле
Рассеяние частиц или рентгеновского излучения на
периодическом потенциале описывается матричными элементами переходов
В случае идеальной решётки
матричные элементы пропорциональны фурье-образу потенциала
Рассеяние идёт в дискретных направлениях
Слайд 40
Дифракция на кристалле
Рассмотрим случай колеблющейся решётки
Матричный элемент рассеяния
можно представить в виде произведения фурье-образа атомного потенциала и
структурного фактора
Слайд 41
Дифракция на кристалле
Положения атомов задаются векторами Rℓ
Это выражение
подставляется в структурный фактор, который затем раскладывается в ряд
по малым смещениям из положений равновесия
Показывается, что происходит рассеяние в любом направлении. Его интенсивность определяется амплитудой колебаний с волновыми векторами, определённым образом связанными с волновыми векторами падающего и рассеянного излучения
Слайд 42
Фактор Дебая-Уоллера
Можно показать, что матричные элементы для упругого
и неупругого рассеяния содержат множитель e-2W, называемый фактором Дебая-Уоллера
Для
его расчёта используется модель Дебая. При высоких температурах:
Слайд 43
Фактор Дебая-Уоллера
Аналогичные рассуждения используются при объяснении температурной зависимости
эффекта Мёссбауэра и люминесценции в твёрдом теле
При расчёте фактора
Дебая-Уоллера можно так же найти величину среднего квадрата смещения атома из положения равновесия:
Слайд 44
Формула Линдемана
Можно предположить, что плавление твёрдого тела происходит,
когда амплитуда колебаний атомов начинает составлять некоторую долю xm
от среднего значения параметра элементарной ячейки rs. Тогда, температуру плавления можно связать, с характеристическими постоянными
xm=0,2-0,25
Слайд 45
Тепловое расширение
Рассмотрим двухатомную молекулу с потенциалом взаимодействия U(x)
Разложим
потенциал в ряд Тейлора вблизи положения равновесия с точностью
до членов четвёртого порядка: U(x)=U(0)+cx2-gx3-fx4
Используя распределение Больцмана, можно показать, что:
Слайд 47
Тепловое расширение
Коэффициент линейного теплового расширения α определяется как
относительное изменение межатомного расстояния в расчёте на единицу изменения
температуры
Таким, образом, эта модель даёт линейную зависимость изменения длины от температуры и показывает связь константы линейного расширения с коэффициентом ангармоничности
Слайд 48
Изменение частот колебаний
Представления о нормальных колебаниях являются следствием
решения уравнений движения в гармоническом приближении
При учёте слагаемого третьего
порядка в разложении потенциала изменится вид уравнений движения
Слайд 49
Изменение частот колебаний
Решение уравнений движения можно искать методом
последовательных приближений
При этом появятся дополнительные решения в виде колебаний
с комбинационными частотами: ωα±ωβ
Амплитуды комбинационных частот пропорциональны произведениям амплитуд соответствующих нормальных колебаний aαaβ
Слайд 50
Изменение частот колебаний
При учёте членов разложения потенциала более
высокого порядка появятся частоты, являющиеся комбинацией большего числа частот
нормальных колебаний
Ещё одним эффектом, обусловленным ангармоничностью, будет смещение частот колебаний осцилляторов
Слайд 51
Теплопроводность
Экспериментально можно установить зависимость, связывающую поток тепла j
с градиентом температуры
В одномерном случае:
j=KT/x, где K – коэффициент
теплопроводности (с точностью до знака)
В трёхмерном:
Слайд 52
Теплопроводность
Явление теплопроводности не согласуется с представлениями о невзаимодействующих
между собой колебаниях решётки (фононах)
Можно сохранить понятие фононов дополнив
его представлениями об их взаимодействии (рассеянии)
Это соответствует учёту ангармоничности в уравнениях движения
Кроме того, механизм взаимодействия фононов необходим для установления теплового равновесия между колебательными состояниями
Слайд 53
Теплопроводность
В кинетической теории газов можно получить выражение: K=1/3Cvℓ,
где C – теплоёмкость единицы объёма, v – средняя
скорость частиц, ℓ – длина свободного пробега
Эту формулу можно применить к твёрдым диэлектрикам, подразумевая под частицами фононный газ
Слайд 54
Теплопроводность
Задача рассмотрения теплопроводности кристаллической решётки – сложная
Установлено, что
теплопроводность обусловлена такими взаимодействиями, в которых импульс фононов изменяется
на вектор обратной решётки (процессы переброса)
Слайд 55
Заключение
Дискретная структура вещества и квантование колебательной энергии приводят
к ряду особенностей в свойствах твёрдого тела, обусловленных колебаниями
кристаллической решётки
Область независимых значений волнового вектора колебаний решётки называется зоной Бриллюэна
Слайд 56
Заключение
Существуют оптические и акустические колебания, отличающиеся законом дисперсии
Использование
гармонического приближения приводит к выводу о существовании невзаимодействующих «нормальных»
колебаний – фононов
Гармоническое приближение не описывает многие важные эффекты
Слайд 57
Заключение
При высоких температурах теплоёмкости твёрдых тел близки к
25 Дж/(моль·К) – закон Дюлонга и Пти
При Т→0 теплоёмкость→0
В
условиях теплового равновесия число фононов с определённой частотой описывается формулой Бозе-Эйнштейна
Слайд 58
Заключение
В модели теплоёмкости Эйнштейна учитывается лишь одна мода
колебаний
В модели Дебая учитываются различные колебательные моды с линейным
законом дисперсии
В общем случае для определения теплоёмкости т.т. надо знать функцию плотности состояний
Слайд 59
Заключение
Важным параметром, использующимся при описании различных свойств твёрдого
тела, является температура Дебая
При рассеянии излучения на кристалле возникает
фон, обусловленный тепловым движением атомов решётки
С увеличением температуры уменьшается интенсивность брэгговского рассеяния и резонансного поглощения/излучения
Слайд 60
Заключение
Тепловое расширение и теплопроводность обусловлены ангармоничностью колебаний частиц
в т.т.
Учёт ангармоничности приводит к изменению частот и конечному
времени жизни колебаний
В гармоническом приближении невозможно установление теплового равновесия между колебательными состояниями
Слайд 61
Контрольные задания
Какое влияние оказывает дискретная структура вещества на
механические колебания распространяющиеся в нём?
В каких эффектах проявляется квантовый
характер колебаний атомов в твёрдом теле?
Чем отличаются оптические и акустические колебания решётки т.т.?
Слайд 62
Контрольные задания
Сколько имеется акустических ветвей колебаний кристаллической решётки?
Сколько
имеется оптических ветвей колебаний кристаллической решётки?
На чём основан метод
ИК спектроскопии вещества?
Слайд 63
Контрольные задания
В чём состоит гармоническое приближение?
В чём заключается
особенность результатов, получаемых при гармоническом приближении?
Что такое нормальные
колебания?
Каковы недостатки гармонического приближения?
Слайд 64
Контрольные задания
Какая теплоёмкость больше, CP или CV, почему?
Почему
пренебрегают разностью теплоёмкостей твёрдого тела при постоянном давлении и
постоянном объёме?
Как зависит теплоёмкость твёрдого тела от температуры при нормальных условиях?
Слайд 65
Контрольные задания
Как ведёт себя теплоёмкость твёрдого тела при
низких температурах?
Что описывает функция распределения Бозе-Эйнштейна? Как она выглядит?
В
чём заключается модель теплоёмкости Эйнштейна?
Какую температурную зависимость теплоёмкости предсказывает модель Эйнштейна?
Слайд 66
Контрольные задания
Для чего используется функция плотности (колебательных) состояний?
Как
находится функция плотности (колебательных) состояний?
В чём заключается модель теплоёмкости
Дебая?
Какой вид имеет функция плотности (колебательных) состояний в модели Дебая?
Слайд 67
Контрольные задания
Что описывает фактор Дебая-Уоллера?
Что описывает формула Линдемана?
Какой
вид имеет функция плотности (колебательных) состояний в модели Эйнштейна?
Как определяются частота и температура Дебая?
Какие эффекты возникают при учёте ангармоничности колебаний?