Статистическая физика

Размещено на http://www.allbest.ru/

Статистическая физика



1. Статистическое описание макросистем. Свойства функции распределения: условие нормировки, правило нахождения средних значений

канонический равновесие движение молекула

Для изучения систем многих частиц информация должна иметь обобщённый характер и относиться не к отдельным частицам, а к совокупности большого числа частиц. Соответствующие понятия также должны относиться не к отдельным частицам, а к большим совокупностям частиц. Новая форма информации и новые понятия требуют иного метода рассмотрения. Этот метод называется статистическим. Законы поведения совокупностей большого числа частиц, исследуемых статистическим методом, называются статистическими закономерностями. Статистические методы в физике имеют более широкое применение, чем динамические. Это связано с тем, что динамический метод эффективен только в применении к системам с небольшим числом степеней свободы. Большинство систем имеют громадное количество степеней свободы и могут изучаться лишь статистическими методами. Кроме того, кванто-механические закономерности по своей природе являются статистическими. Поэтому статистические методы используются так же и при изучении систем с небольшим числом степеней свободы, если в поведении этих систем существенны кванто-механические эффекты.

Рассмотрим математические понятия, которые понадобятся при изучении статистического метода. Пусть некоторая макроскопическая система находится в заданном состоянии. Предположим, что какая-то характерная для системы величина x может иметь дискретные значения:

x₁, x₂, x₃, …, x_i, …, x_k, …, x_s.

Осуществим над системой очень большое число N измерений величины x, приводя систему перед каждым измерением в одно и то же состояние. Вместо того чтобы производить повторные измерения над одной и той же системой, можно взять N одинаковых систем, находящихся в одном и том же состоянии, и осуществить однократное измерение величины х у всех этих систем. Набор одинаковых систем, находящихся в одинаковом состоянии, называется статистическим ансамблем.

Допустим, что N₁ измерений дали результат x₁, N₂ измерений - результат x₂,…, N_i измерений - результат x_i и т.д.(УN_i=N - числу систем в ансамбле.) Величина N_i/N называется относительной частотой появления результата x_i, а предел этой величины, получающийся при стремлении N к бесконечности, т.е.

(2.1)

называется вероятностью появления результата x_i. В дальнейшем для упрощения формул будем писать выражение для вероятности в виде N_i/N, производя предельный переход при N>?.

Поскольку УN_i = N, то

(2.2)

т.е. сумма вероятностей всех возможных результатов равна единице.

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

. (2.3)

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

Возьмём очень малую величину б (скажем б = 10^-6) и найдём число измерений ДN₀, при которых 0<x< б; ДN₁, при которых б<x<2б,…; ДN_x, при которых результат измерений находится в интервале от х до х+б, и т.д. Вероятность того, что результат измерений окажется в интервале от 0 до б равна ДP₀ = ДN₀/N, в интервале от б до 2б ДP₁ = ДN₁/N, …, в интервале от х до (х+б) ДP_х = ДN_х/N. Начертим ось х и отложим вверх от неё полоски шириной и высотой ДP_х/б (рис. 1). Полученная столбчатая диаграмма называется гистограммой. Площадь полоски, левый край которой имеет координату х, равна ДP_х, а площадь всей гистограммы единице.

Рис. 1

Рис. 2

Гистограмма наглядно характеризует вероятность получения результатов измерений, заключающихся в различных интервалах ширины б. Чем меньше ширина интервала б, тем детальнее будет охарактеризовано распределение вероятностей значений величины х. В пределе при б>0 ступенчатая линия, ограничивающая гистограмму, превратится в гладкую кривую (рис. 2). Функция ѓ(х), определяющая эту кривую, называется функцией распределения вероятностей (плотностью вероятности).

В соответствии со способом построения кривой распределения площадь столбика шириной dx (см. рис. 2) равна вероятности того, что результат измерения окажется в пределах от х до х+dx. Обозначив эту вероятность через dР_х, можно написать, что

(2.4)

Индекс «х» при dP указывает на то, что имеется в виду вероятность интервала, левый край которого лежит в точке с координатой х. Площадь, ограниченная кривой распределения, так же как и площадь гистограммы, равна единице. Это означает, что

(2.5)

Интегрирование производится по всему интервалу возможных значений величины х. Формула (2.5) является аналогом формулы (2.2).

Зная функцию распределения ѓ(х), можно найти среднее значение результатов измерения величины х. В случаях получается результат, равный х. Сумма таких результатов определяется выражением . Сумма всех возможных результатов равна . Разделив эту сумму на число измерений N, получим среднее значение величины х

(2.6)



Эта формула является аналогом формулы (2.3).

Для произвольной функции ц(x) можно записать правило нахождения среднего значения

.

2. Микроканоническое и каноническое распределения Гиббса

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

, (2.7.)

где E - заданное значение полной энергии системы; H - гамильтониан, который представляет собой сумму полных механических энергий всех частиц системы; - дельта-функция, обладающая следующим свойством: она равна бесконечности, когда ее аргумент равен нулю, и равна нулю во всех остальных случаях; Щ - статистический вес (нормировочная константа); X - совокупность координат и импульсов всех частиц системы; f_N - функция распределения системы из N частиц.

Формула (2.7) называется микроканоническим распределением Гиббса и означает, что в изолированной системе нет таких частиц, которые не удовлетворяли бы закону сохранения полной энергии.

Для случая, когда система находится в термостате, Гиббс доказал, что из микроканонического распределения следует каноническое распределение

,

где Z - нормировочная константа.

3. Распределение Максвелла

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

Максвелл предположил, что вероятность того, что скорости частицы в равновесном идеальном газе находятся в интервале (dv_xdv_ydv_z), может быть записана в виде:

, (2.8)



где C - константа, которую можно найти из условия нормировки:

. (2.9)

Интеграл (3.2) можно разбить на произведение трех независимых интегралов по каждой компоненте скорости, каждый из которых является интегралом Пуассона:

. (2.10)

Тогда получим функцию распределения Максвелла по компонентам скорости:

. (2.11)

Эта функция распределения распадается на произведение трех функций распределения по каждой из компонент скорости. Например, для x-компоненты получим

. (2.12)

Заметим, что функция распределения по компоненте скорости является симметричной относительно замены x на -x. Отсюда следует вывод о том, что среднее значение компоненты скорости x равно нулю (это можно так же подтвердить прямым взятием интеграла). Поскольку все компоненты скорости равноправны, то отсюда следует, что и средней вектор скорости равен нулю

.

Это равенство верно только в условиях равновесия.

Заметим, что в уравнении (2.9) интегрирование ведется от - до + , хотя очевидно, что для любой системы с конечным числом частиц скорость каждой частицы ограничена, кроме того, она не может быть больше скорости света.

Найдем теперь распределение Максвелла по модулю скорости. Для этого нужно перейти от компонент скорости к ее модулю по формуле

,

а также определить физически бесконечно малый объем в пространстве модулей скорости. Такой объем соответствует объему тонкой сферы радиуса v и равен 4рv²v. Проводя соответствующие замены переменных, получим функцию распределения Максвелла по модулю скорости:

. (2.13)

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

Исследуем свойства функции распределения по модулю скорости.

Во-первых, беря производную от функции F(v), и приравнивая ее к нулю, можно найти максимум функции распределения, который соответствует наиболее вероятной скорости частиц:

Два других экстремума v = 0 и v = ? являются минимумами функции распределения.

Найдем среднее значение модуля скорости. По правилу нахождения средних получим

. (2.14)

Занося одну из скоростей под дифференциал, получим интеграл, который можно взять по частям. В результате получим:

. (2.15)

Найдем так же средний квадрат скорости:

(2.16)

Этот интеграл является табличным, и может быть легко получен из интеграла Пуассона



.

Полученная величина позволяет найти среднее значение кинетической энергии молекул

.

Этот результат без вывода был получен в элементарной кинетической теории.

Размещено на http://www.allbest.ru/

Рис. 3

Вид распределения Максвелла для молекул азота показан на рис. 3. Температуры на рисунке равны: T₁ = 50 К, T₂ = 100 К, T₃ = 300 К. С увеличением скорости максимум распределения смещается в сторону больших скоростей, а высота кривой в максимуме несколько понижается. Наличие максимума объясняется тем, что кривая отражает результат двух противоборствующих тенденций: вероятность состояний с ростом скорости падает, а плотность состояний увеличивается. При малых скоростях преобладает тенденция роста плотности состояний. При скоростях после максимума кривой преобладает тенденция уменьшения вероятности состояний.

При комнатной температуре характерные скорости молекул кислорода и азота в воздухе равны примерно 400-500 м/с. Скорости молекул водорода при этом примерно в четыре раза больше. С повышением температуры скорости молекул растут как

Частота ударов молекулы о стенку. Направим ось перпендикулярно стенке (рис. 2.4) и обозначим концентрацию молекул. Тогда плотность потока молекул в направлении стенки со скоростями между и равна

, (2.17)

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

Рис. 4

(2.18)



Принимая во внимание формулу (2.15), окончательно напишем

. (2.19)

Число молекул в различных участках распределения Максвелла. Если плотность молекул , то число молекул, скорости которых распределены между и , равно

(2.20)

При вычислении (2.18) учтено, что в можно представить как , что более наглядно при анализе формы кривой распределения Максвелла. Имеются таблицы интеграла

. (2.21)

С их помощью величина (2.20) вычисляется по следующей очевидной формуле:

(2.22)

Из таблицы, в частности, находим:



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

Пример 2.1. Найти число молекул кислорода О₂, скорости которых заключены в пределах от 195 до 205 м/с при . Масса кислорода 0,1 кг.

Поскольку интервал скоростей от 195 до 205 м/с достаточно мал, можно воспользоваться теоремой о среднем и по формуле (2.13) написать:

(2.23)

где м/с, м/с.

Относительная молекулярная масса кислорода и, следовательно, масса молекулы кг кг. Молярная масса кислорода кг/моль, поэтому в 0,1 кг кислорода имеется молекул. Далее учтем, что Дж Дж, поэтому

Экспериментальная проверка распределения Максвелла. Ввиду принципиальной важности распределения Максвелла для статической физики оно было много раз подвергнуто тщательной экспериментальной проверке. Принципиальная схема наиболее типичной экспериментальной установки состоит в следующем. Через отверстие выходит пучок исследуемых молекул. Чтобы в процессе движения пучка распределение молекул в нем не изменялось, они должны двигаться практически без взаимодействия друг с другом. Поэтому на пути движения пучка создается высокий вакуум, а в сосуде газ находится под низким давлением. Для того чтобы в выходящем из сосуда пучке молекулы имели такое же распределение скоростей, как в сосуде, необходимо обеспечить истечение газа через отверстие без гидродинамического напора. Это возможно, если в области вблизи отверстия молекулы не успевают сталкиваться друг с другом. Тогда молекула, попадающая в отверстие, вылетает из сосуда, не возмущая состояния всех остальных молекул в сосуде. В результате число молекул в сосуде медленно уменьшается, но их равновесное состояние не изменяется. Чтобы обеспечить такое «бесстолкновительное» покидание сосуда молекулами, отверстие по размерам должно быть много меньше длины свободного пробега молекул, т.е. среднего расстояния между последовательными столкновениями.

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

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

Для сортировки молекул наиболее часто используют метод вращающихся дисков с щелями вдоль радиуса (рис. 2.5, а). С помощью этого метода Физо в прошлом столетии измерял скорость света в земных условиях. Если щели повернуты на угол друг относительно друга, а диски располагаются на расстоянии друг от друга, то при угловой скорости диски повернутся на угол в течение времени и (где ), соответствующем нескольким поворотам дисков за время прохождения молекулами расстояния . Регистрация молекул производится различными методами в зависимости от их свойств. В простейшем случае они осаждаются на экран и по толщине осажденного слоя можно судить об их числе. Так поступают, например, когда в качестве объекта исследования берется пучок атомов серебра, газ из которых образуется при нагревании в результате испарения.

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

Рис. 5

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

Проведенные эксперименты подтвердили справедливость распределения Максвелла.

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

Возьмем, например, некоторые четыре элемента объема 1-4 (см. рис. 6, а) и представим себе обмен частицами между ними. Каждая из стрелок изображает определенное число частиц, которые покидают рассматриваемый объем или проходят в него в единицу времени. Например, на диаграмме б (рис. 6) из объема 1 частицы уходят в объем 2, но такое же количество приходит в него из объема 4 и т.д. На диаграмме в (рис. 6) в объем 3 приходят частицы из объемов 2 и 4, но зато равное их сумме число частиц уходит в объем 1. В результате осуществления указанных обменов частицами обеспечивается постоянство частиц во всех объемах.

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

а

Размещено на http://www.allbest.ru/

б

Размещено на http://www.allbest.ru/

в

Рис. 6



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

4. Распределение Больцмана

Барометрическая формула. Рассмотрим газ, находящийся в равновесии в поле силы тяжести. В этом случае сумма действующих сил на каждый элемент объема газа равна нулю. Выделим малый объем газа на высоте h (рис. 7) и рассмотрим действующие на него силы:

Рис. 7

На выделенный объем действует сила давления газа снизу, сила давления газа сверху и сила тяжести. Тогда баланс сил запишется в виде



,

где dm - масса выделенного объема. Для этого объема можно записать уравнение Менделеева-Клапейрона

.

Выражая величину dm, можно получить уравнение

.

Разделяя переменные, получим

.

Проинтегрируем полученное уравнение, учтя, что температура постоянна,

.

Пусть давление на поверхности равно p₀, тогда полученное уравнение легко преобразовать к виду



. (2.24)

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

Распределение Больцмана. Запишем барометрическую формулу (2.24) через концентрацию частиц, воспользовавшись тем, что p = nkT:

, (2.25)

где m₀ - масса молекулы газа.

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

. (2.26)

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

Найдем число частиц газа, координаты которых находятся в элементе объема dV = dxdydz



.

Полное число частиц в системе может быть записано в виде

.

Здесь интеграл формально записан по всему пространству, но надо иметь в виду, что объем системы конечен, что приведет к тому, что интегрирование будет вестись по всему объему системы. Тогда отношение

как раз и даст вероятность того, что частица попадет в элемент объема dV. Тогда для этой вероятности запишем

,

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



. (2.27)

Это и есть функция распределения Больцмана по координатам частиц (или по потенциальным энергиям, имея в виду, что потенциальная энергия зависит от координат). Легко показать, что полученная функция нормирована на единицу.

Связь распределений Максвелла и Больцмана. Распределения Максвелла и Больцмана являются составными частями распределения Гиббса. Температура определяется средней кинетической энергией. Поэтому возникает вопрос, почему в потенциальном поле температура постоянная, хотя по закону сохранения энергии при изменении потенциальной энергии частиц должна также изменяться их кинетическая энергия, а следовательно, как кажется на первый взгляд, и их температура. Другими словами, почему в поле тяжести при движении частиц вверх у всех них кинетическая энергия уменьшается, а температура остается постоянной, т.е. остается постоянной их средняя кинетическая энергия, а при движении частиц вниз энергия всех частиц увеличивается, а средняя энергия остается постоянной?

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

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

Атмосфера планет. Потенциальная энергия частицы массой в поле тяготения шарообразного небесного тела равна

, (2.28)

где - масса тела; - расстояние от центра тела до частицы; - гравитационная постоянная. Атмосфера планет, в том числе и Земли, не находится в равновесном состоянии. Например, вследствие того, что атмосфера Земли находится в неравновесном состоянии, ее температура не постоянна, как это должно было быть, а изменяется с высотой (уменьшается с увеличением высоты). Покажем, что равновесное состояние атмосферы планеты в принципе невозможно. Если бы оно было возможно, то плотность атмосферы должна была бы изменяться с высотой по формуле (2.26), которая принимает вид

(2.29)



где учтено выражение (2.28) для потенциальной энергии, - радиус планеты. Формула (2.29) показывает, что при плотность стремится к конечному пределу

(2.30)

Это означает, что если в атмосфере имеется конечное число молекул, то они должны быть распределены по всему бесконечному пространству, т.е. атмосфера рассеяна.

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

Для количественного рассмотрения вопроса о потере атмосферы планетами необходимо принять во внимание распределение молекул по скоростям. Силу земного притяжения могут преодолеть лишь молекулы, скорость которых превосходит вторую космическую. Эти молекулы находятся в «хвосте» распределения Максвелла и их относительное число незначительно. Тем не менее, за значительные промежутки времени потеря молекул является чувствительной. Поскольку вторая космическая скорость у тяжелых планет больше, чем у легких, интенсивность потери атмосферы у массивных небесных тел меньше, чем у легких, т.е. легкие планеты теряют атмосферу быстрее, чем тяжелые. Время потери атмосферы зависит также от радиуса планеты, состава атмосферы и т.д. Полный количественный анализ этого вопроса является сложной задачей.

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

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

(2.31)

Поэтому распределение концентраций этих частиц по высоте дается формулой

(2.32)

Чтобы эффект был достаточно хорошо заметен, частицы должны быть достаточно малыми. Число таких частиц на разных высотах в сосуде считают с помощью микроскопа. Эксперименты такого рода впервые были выполнены начиная с 1906 г. Ж.Б. Перреном (1870-1942).

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

Таким путем Перрен измерил и получил результат, весьма близкий к современному. Другим независимым способом значение было получено Перреном из опытов с броуновским движением.

В последующем были проведены также эксперименты другого типа, полностью подтвердившие распределение Больцмана. Из экспериментов другого типа можно указать, например, на проверку зависимости поляризации полярных диэлектриков от температуры, рассмотренную выше.

Пример 2.2. Перрен использовал распределение гуммигутовых зерен в воде для измерения постоянной Авогадро. Плотность частиц гуммигута составляла = 1,2110³кг/м³, их объем = 1,0310^-19м³. Температура, при которой проводился эксперимент, была равна . Найти высоту , на которой плотность распределения гуммигутовых зерен уменьшилась в два раза.

Принимая во внимание, что, по условию задачи, ( - ₀) = 0,2210^-16 кг, получаем на основе формулы (2.32) h = kTln2/[( - ₀)g] = 12,310^-6 м.

Пример 2.3. В воздухе при температуре и давлении Па взвешены шарообразные частицы радиусом 10^-7 м. Найти массу взвешенной частицы.

По формуле (2.32) находим ( - ₀) = kTln2/gh = 1,0610^-23 кг.

Учитывая, что = 4,1910^-21 м³, находим ( - ₀) = 2,5310^-3 кг/м³. Поскольку ₀ = 1,293 кг/м³, получаем = 1,296 кг/м³ и, следовательно, масса частицы

5. Статистический смысл энтропии

Рассмотрим первое начало термодинамики с точки зрения статистической физики. Как известно , откуда

. (2.33)

Внутренняя энергия находится по правилу нахождения средних от гамильтониана:

, (2.34)

где a - внешний параметр.

Возьмём дифференциал по а и Т от (2.34)

(2.35)

Теперь рассмотрим микроскопические силы. Они определяются выражением а макроскопическая сила находится по правилу средних от микроскопических сил: т.е. работа макроскопических сил . Сравнивая это выражение с (2.35) и (2.33), можно прийти к выводу, что

. (2.36)

Умножим и разделим (2.36) на kT и вспомним каноническое распределение Гиббса:

,

и представим его в виде:

, (2.37)

т.е. (2.37а)

Из условия нормировки возьмём дифференциал по а и Т

(2.38)

Учитывая это, представим (2.36) в виде

. (2.39)



Или, учитывая (2.37)

. (2.40)

Количество теплоты в свою очередь можно представить, как , откуда энтропия равна

. (2.41)

И, в заключение, несколько слов о величине F(a,T). Рассмотрим формулу для свободной энергии, полученную ранее . Подставим в эту формулу (2.41) и (2.34)

. (2.42)

Сравнив (2.42) и (2.37), получим

, (2.43)

откуда следует, что введённая нами в (2.37) величина имеет смысл свободной энергии.

Для дискретного случая энтропию можно представить в виде

,

где p_i - вероятность i-го микросостояния системы. Легко показать, что распределение вероятностей, когда одно из значений p_i равно единице, а остальные нулю, приводит к минимальному значению энтропии, равному нулю.

С другой стороны, максимально значение энтропии соответствует равновероятному распределению. В этом случае энтропия равна

,

где Щ - число микросостояний (ранее введенная величина - статистический вес). Таким образом, энтропия может выступать в качестве меры беспорядка системы.

6. Термодинамические функции идеального газа

Размещено на http://www.allbest.ru/

Рис. 8

Пусть ИГ заключён в сосуд с объёмом V (рис. 8). Из формулы (2.37а) . Так как мы взяли за модель идеальный газ, то частицы не взаимодействуют друг с другом, т.е. сумма всех энергий это сумма кинетических энергий.

.

С учётом вышесказанного представим это выражение как произведение интегралов по отдельным частицам



Подставим это выражение в (2.37а):

. (2.44)

А отсюда следует:

или ,

Пусть , а (термодинамический предел), то:

(2.45)

Но S и F должны быть аддитивными величинами, т.е. пропорциональными N, а не lnN. Это значит, что классические распределения не удовлетворительны, т.е. нужны квантовые распределения.



7. Распределение Гиббса для квантовых систем

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

Во-вторых, в квантовой механике нельзя абсолютно точно одновременно определить скорость и координату частицы. Существует соотношение неопределённостей Гейзенберга: , где h - постоянная Планка, показывающее предел точности определения скорости и координаты частицы. Следовательно, в фазовом пространстве существует элементарная ячейка, объем которой равен h³. Положение частицы в фазовом пространстве можно определить лишь с точностью до элементарной ячейки. Более точное определение невозможно.

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

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



,

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

Каноническое распределение можно представить в следующем виде:

,

где z - статистическая сумма, которая находится по формуле

.

8. Связь квантовых и классических распределений Гиббса

Чтобы от квантовых распределений перейти к классическим функциям, надо перейти от суммирования к интегрированию . Для этого выделим в фазовом пространстве (пространстве координат-импульсов всех частиц) элементарную ячейку (а в классической механике можно сделать объём сколь угодно малым). Объём этой ячейки, учитывая соотношение неопределённостей Гейзенберга, (ДxДp_x?h³, V=ДxДp_x). h^3N - объём всех частиц.

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

.

9. Парадокс Гиббса

Рис. 9

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

. (2.46)

Уберём перегородку: . Изменяется и значение энтропии:

, (2.47)

, (2.48)



т.е. энтропия возрастает, хотя не должна. Ответ на это противоречие заключается в том, что энтропия определена с точностью до константы (). Гиббс предложил следующее выражение для S₀:

. (2.49)

Но, если использовать квантовое распределение Гиббса - парадокса можно избежать. Здесь выражение для S₀ выглядит по-другому:

, (2.50)

, т.е. если , то формула совпадает с (2.49).

10. Третье начало термодинамики с точки зрения статистической физики

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

. (2.51)

Это утверждение представляет собой содержание так называемой теоремы Нернста. Иногда утверждение (2.51) называют третьим началом термодинамики.

С точки зрения статистической физики это утверждение может быть интерпретировано следующим образом. Рассмотрим систему, находящуюся в равновесии при температуре, стремящейся к нулю кельвинов. При этом система находится в состоянии с наименьшей энергией (основном состоянии). Тогда в каноническом распределении Гиббса наибольшим будет слагаемое, отвечающее наименьшей энергии, а остальные будут стремиться к нулю (в силу условия нормировки). Как было показано ранее, такое распределение соответствует энтропии системы и равно нулю. Таким образом, третье начало термодинамики с точки зрения статистической физики можно сформулировать так: при нуле кельвинов система находится, в основном, в полностью определенном состоянии, что соответствует нулевой энтропии (как мере неопределенности).

11. Флуктуации

Пусть ѓ -- любая физическая величина, испытывающая флуктуации. Флуктуацией величины ѓ называется отклонение мгновенного значения этой величины от ее среднего значения. Очевидно, что <Дѓ> = 0. Поэтому обычно пользуются средним квадратом флуктуации, т.е. величиной (<Дѓ>)². Квадратный корень из этой величины называется -- среднеквадратичной относительной флуктуацией.

Усредняя выражение

,

получим

.

Но ѓ - величина постоянная, а потому. Следовательно,



. (2.52)

Усредним теперь произведение двух флуктуирующих величин:

fg = (<f>+ Дѓ ) (<g>+ Дg)=,

Так как f и g -- величины постоянные, а (<Дѓ>)= <Дg >= 0, то

. (2.53)

Формула (2.52) содержится здесь как частный случай, который получается при g =f.

Величины f и g называются статистически независимыми, если

<ДfДg> = 0.

Для таких величин:

<fg>=<f><g>. (2.54)

Рассмотрим теперь любую физическую систему, состоящую из N независимых одинаковых частей. Примером такой системы может служить идеальный газ, а составных частей -- отдельные молекулы. Пусть f_t -- произвольная аддитивная величина, характеризующая i-ю подсистему, например, в приведенном примере - кинетическая энергия r'-й молекулы. Тогда в силу предполагаемой аддитивности соответствующая величина для всей системы будет F = Уѓ_i. Выразим средний квадрат флуктуации величины F через аналогичный квадрат для величины ѓ_i. Очевидно, <F> = У<ѓ>_i = N<f>, где индекс i опущен, так как предполагается, что все составные части системы тождественны. Далее,

.

А так как эти части независимы, то <ѓ_i ѓ_j >= <ѓ_i >< ѓ_j >=(<ѓ>)². Следовательно, .

Подставляя эти значения в формулу (2.52), получим

. (2.55)

Отсюда на основании (2.52)

. (2.56)

С увеличением N относительная флуктуация величины F убывает обратно пропорционально При больших N относительные флуктуации ничтожны. Этот вывод качественно верен и для неаддитивных величин. С ним связана достоверность термодинамических результатов для больших макроскопических систем.

В соответствии со сказанным видим, что в объемах с большим средним числом частиц N относительные флуктуации малы и трудно доступны наблюдению. Наоборот, при малых N относительные флуктуации велики. Более общий метод вычисления флуктуации плотности, применимый также к жидкостям и твердым телам, основан на теореме о равномерном распределении кинетической энергии по степеням свободы. Рассмотрим малую часть жидкости или газа, окруженную такой же жидкой или газообразной средой, температура Т которой поддерживается постоянной (термостатом). С целью упрощения и наглядности вычислений предположим, что малая часть жидкости или газа заключена в цилиндр с поршнем. Стенки цилиндра идеально проводят тепло, а поршень может ходить в нем без трения. Тогда наличие стенок цилиндра и поршня не будет препятствовать обмену энергией и выравниванию давлений между веществом в цилиндре и термостатом. Тепловое движение молекул вещества вызовет броуновское движение поршня. К этому движению поршня мы и применим теорему о равномерном распределении кинетической энергии по степеням свободы. Поршень можно рассматривать как гармонический осциллятор, совершающий беспорядочные тепловые колебания. Среднее значение его потенциальной энергии при смещении на х из положения равновесия х = 0 равно (1/2)кх² = (1/2)кТ, где к -- жесткость, соответствующая такому смещению. Если S -- площадь поршня, a AV -- изменение объема системы, то V = Sx. Таким образом, <(ДV)² > = <S²x²>= S²kT/K. Сила, возвращающая поршень в положение равновесия, будет , где Р -- давление газа или жидкости. Поэтому к = -- SdP/dх =S²dP/dV. В результате получим

. (2.57)

Знак Т указывает, что в выводе предполагалось постоянство температуры окружающей среды (термостата). Если бы вещество внутри объема V было адиабатически изолировано, то индекс Т следовало бы заменить на S (постоянство энтропии), т.е.

<(ДV ²)>_s = - kT(dV/dP)_s. (2.58)



Формулы (2.57) и (2.58) выражают флуктуации объема одной и той же массы вещества, находящейся в термодинамическом равновесии с окружающей средой.

Для идеального газа при постоянной температуре PV = const, так что (dV/dP) _т = -V/P. А так как PV - NкТ, где N - число молекул в объеме V, то из формулы (2.57) получаем (Д V ²)=V²/N.

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

, (2.59)

где б = 1/кT, dГ - элементы объема пространства скоростей, а Z определяется условием нормировки:

. (2.60)

Дифференцируя это соотношение по параметру б, получим

и формула (2.59) перейдет в

. (2.61)



Отсюда

.

Аналогично

. (2.62)

Сравнение этой формулы с предыдущей дает

.

Или после подстановки параметра б = 1/кТ

. (2.63)

Так как для идеального одноатомного газа = (3/2)кТ, то в этом случае

. (2.64)

Приведенный метод вычисления <Де²> может быть распространен практически без изменений на случай величин внешнего силового поля. Надо только максвелловское распределение заменить распределением Больцмана и производить интегрирование не только по скоростям, но и по координатам обычного пространства. В результате снова получится формула (2.63), в которой, однако, под е следует понимать уже полную энергию молекулы, состоящую из кинетической и потенциальной.

Важно отметить расширение области применимости формулы (2.63) для флуктуации полной энергии е. Выделим малую часть (подсистему) изотропной среды (жидкости или газа), находящуюся в статистическом равновесии со всей средой, температура Т которой поддерживается постоянной. Подсистемой может быть и отдельная молекула. По отношению к выделенной подсистеме окружающая среда играет роль термостата. Из-за обмена энергией между термостатом и подсистемой энергия последней будет непрерывно флуктуировать. Флуктуации полной энергии е подсистемы определяются уравнением (2.63). Не имеет значения, меняется ли энергия е непрерывно (классическая система), или принимает дискретный ряд значений (квантовая система). Доказательство этого утверждения дается в статистической физике. Оно основано на каноническом распределении Гиббса, частными случаями которого являются распределения Максвелла и Больцмана.

В случае макроскопической подсистемы, объем которой поддерживается постоянным, имеет смысл внутренней - энергии подсистемы, a dе/dT -ее теплоемкости С_н при постоянном объеме. Тогда из формулы (2.63) получается

<(Де²)>_V = kT ²C_v. (2.65)

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

Рассмотрим теперь флуктуации энтальпии H подсистемы. Для этого воспользуемся следующим искусственным приемом. Предположим, что подсистема заключена в оболочку с идеально проводящими подвижными стенками, так что объем подсистемы не сохраняется постоянным. Пусть оболочка снаружи подвергается действию постоянных внешних сил, поддерживающих внешнее давление Р постоянным. Эти силы увеличивают
потенциальную энергию подсистемы на величину PV. Если под понимать ту же энергию, что и в предыдущем выводе, то с учетом дополнительной потенциальной энергии PV среднее значение полной энергии подсистемы будет . Но это есть энтальпия подсистемы H. Все предыдущие рассуждения можно повторить без изменений, заменив на + PV. В результате вместо (2.63) получится

. (2.56)

Но при Р = const производная dH/dT есть теплоемкость Ср подсистемы при постоянном давлении, а потому

. (2.67)

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

12. Связь флуктуаций с энтропией системы. Распределение Гаусса

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

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

В работах Эйнштейна и Смолуховского было показано, что целый ряд давно известных физических процессов обусловлен явлением флуктуаций, и была развита количественная теория этих процессов, оказавшаяся в прекрасном согласии с экспериментами. Рассмотрим флуктуации в замкнутой системе. Пусть система находится в состоянии статического равновесия и имеет энтропию S₀. Предположим теперь, что состояние системы изменяется так, что она переходит в неравновесное состояние с энтропией S. Будем считать, что изменение состояния системы можно характеризовать изменением некоторого внутреннего параметра x, значение которого зависит от состояния всей системы. В равновесии x = x₀. Пример x - плотность с газа, находящегося в замкнутом теплоизолированном сосуде в состоянии равновесия плотность постоянна по всему объему сосуда x₀ = с₀ = const.

В результате флуктуации система может самопроизвольно перейти в неравновесное состояние с переменной плотностью x = (r).

Энтропия системы будет некоторой функцией параметра x, S = S(x). При этом в состоянии равновесия S₀ = S(x₀).

Вероятность застать систему в интервале значений x, x+dx должна быть пропорциональна статистическому весу (числу микросостояний) и величине интервала dx:

dP(x) (x)dx.

Для изолированной системы верна формула Больцмана:

S = kln

dP(x) = const exp(S(x)/k)dx.

Константа определяется из условия нормировки



1 = ?dP(x) = const?exp(S(x)/k)dx.

Или для функции распределения

F(x) = const exp(S(x)/k). (2.68)

Выберем x так, что энтропия имеет максимум при x = = 0. Поэтому:

S/x_x=0 = 0, ²S/x²_x=0 < 0.

Величина x при флуктуациях очень мала. Разлагая S(x) в ряд по степеням x и ограничиваясь членом второго порядка, получим

S(x) = S(0) (1/2)x²,

где > 0. Тогда получим распределение вероятностей в виде:

f(x)dx = Aexp(x²/2k)dx.

Нормировочная константа A определяется условием

.

Хотя выражение для f(x) относится к малым x, но ввиду быстрого убывания подынтегральной функции с увеличением x область интегрирования можно распространить на все значения от до . Произведя интегрирование, получим:



A = (интеграл Пуассона).

Таким образом, распределение вероятностей для различных значений флуктуаций x определяется формулой

f(x)dx = exp(x²/2k)dx.

Распределение такого вида называется распределением Гаусса.

Оно имеет максимум при x = 0 и быстро спадает с увеличениемx симметрично в обе стороны.

Средний квадрат флуктуаций равен

x² = = k/.

Поэтому распределение Гаусса можно записать в виде:

f(x)dx = dx.

Функция f(x) имеет тем более острый максимум, чем меньше x².

Отметим, что по известному x² можно найти аналогичную величину для любой функции (x). Ввиду малости x имеем:

()² = ,



(подразумевается, что функция (x) мало меняется при значениях x x² и что производная d/dx отлична от нуля при x = 0).

Если , то распределение Гаусса запишется в виде:

dP = f(x)dx = dx.

Размещено на Allbest.ru

...