228 |
Г л а в а 4. Статистическая термодинамика |
Таким образом, зная сумму по состояниям, можно определить внутреннюю энергию как функцию температуры и объема. Дифференцируя (15.19) по температуре, находим выражение для изохорной теплоемкости:
CV |
|
∂U |
|
∂ ln Z |
+ kT |
2 |
|
∂ 2 ln Z |
= |
|
= 2kT |
∂T |
|
|
|
∂T |
2 |
. |
|
|
∂T V |
|
V |
|
|
|
|
V |
Энергия Гельмгольца связана с внутренней энергией уравнением Гиббса–Гельмгольца:
|
∂(F / T ) |
= − |
U |
|
∂T |
|
|
|
. |
T |
2 |
|
V |
|
|
|
Подставляя в это уравнение выражение (15.19) и выполняя интегрирование, получаем формулу, связывающую сумму по состояниям и энергию Гельмгольца:
F − U 0 = −kT ln Z ,
где U0 – постоянная интегрирования, равная энергии при абсолютном нуле.
Дифференцируя функцию F по температуре и объему, находим энтропию и давление:
|
|
∂F |
∂ ln |
Z |
|
U − U 0 |
|
|
S = − |
|
= k ln Z + kT |
|
|
= k ln Z + |
|
, |
|
∂T |
T |
|
|
∂T V |
|
V |
|
|
|
∂F |
|
∂ ln Z |
p = − |
|
= kT |
∂V |
. |
|
∂V T |
|
T |
Последнее соотношение есть не что иное, как термическое уравнение состояния, т.е. зависимость давления от объема и температуры.
Интересный результат получается при попытке представить энтропию в виде среднего значения по каноническому ансамблю подобно тому, как это сделано для внутренней энергии в (15.17). Обозначим вероятность нахождения частицы на i-м уровне через pi, тогда из распределения Больцмана (14.29) и определения статистической суммы
(15.11) следует:
|
(15.25) |
− ln pi |
= |
Ei |
+ ln Z . |
|
kT |
|
|
|
|
|
Усредним это тождество по всем уровням, то есть умножим его на pi и просуммируем с учетом ∑pi = 1:
|
|
i |
|
|
|
|
(15.26) |
−∑pi ln pi |
= |
U − U 0 |
+ ln Z . |
|
kT |
|
|
i |
|
|
230 |
Г л а в а 4. Статистическая термодинамика |
(15.30)
(15.31)
(15.32)
(15.33)
щего из N одинаковых частиц, выражается через произведение сумм по состояниям одной частицы Q:
Z ид = Q N ,
N !
где множитель 1/N! учитывает квантовый принцип неразличимости частиц.
Энергия молекулы идеального газа складывается из поступательной и внутренней энергии:
E= Eпост + Eвнутр.
Всвою очередь, внутреннюю энергию молекулы в хорошем приближении можно представить как сумму ядерной, электронной, колебательной и вращательной энергии:
Eвнутр = Eяд + Eэл + Eкол + Eвращ,
поэтому молекулярная сумма по состояниям факторизуется, то есть представляется в виде произведения сумм по состояниям, соответствующих отдельным видам движения:
Q = Qпост Qвнутр = Qпост Qяд Qэл Qкол Qвращ.
Энергии, соответствующие различным внутренним видам движения, заметно отличаются друг от друга по порядку величины. В табл. 15.1 приведены типичные значения энергии переходов между соседними уровнями и указан диапазон электромагнитного излучения, соответствующий этим переходам1. Более подробная классификация различных видов электромагнитного излучения содержится в табл. П-14 в Приложе-
нии III.
Таблица 15.11 |
Свойства различных видов внутренней энергии |
|
|
|
|
|
|
Вид энергии |
Ядерная |
Электронная |
Колебательная |
Вращательная |
|
|
|
|
|
|
Энергия перехода, см–1 |
~ 1010 |
~ 104 ÷ 105 |
~ 103 ÷ 104 |
~ 100 ÷ 101 |
Длина волны |
~ 10–3 нм |
~ 102 ÷ 103 нм |
~ 103 ÷ 104 нм |
100 ÷ 101 мм |
Диапазон излучения |
γ-излучение |
УФ, видимое |
ИК |
микроволновое |
Рассмотрим способы расчета отдельных сумм по состояниям.
а) Поступательную сумму по состояниям рассчитывают в классическом приближении по формуле (15.13) с функцией Гамильтона H(p,q) = p2 / 2m (m – масса молекулы). Интегрирование по трем координатам и трем проекциям импульса производится раздельно и дает:
1 Длина волны перехода λ (см) связана с разностью энергий уровней ∆E (см–1) соотношением: λ = 1 / ∆E.
Г л а в а 4. Статистическая |
термодинамика |
231 |
2πmkT |
3 / 2 |
|
|
Qпост = |
|
|
|
V , |
(15.34) |
h |
2 |
|
|
|
|
|
|
где V – объем, в котором движется молекула.
б) Внутренние суммы по состояниям рассчитывают с использованием квантового определения (15.11). Значительный вклад в сумму (15.11) вносят только те уровни энергии, для которых Ei ≤ kT. Если энергия измеряется в см–1, это условие приобретает вид: hcEi ≤ kT, где c = 3 1010 см с–1 – скорость света. Подставляя в это неравенство значения постоянных Планка и Больцмана, находим:
Ei (см–1) < T (К). |
(15.35) |
Только такие уровни энергии необходимо учитывать в сумме по состояниям, вкладом от остальных уровней обычно пренебрегают. Если теперь обратиться к табл. 15.1, то можно увидеть, какие именно температуры нужны для возбуждения отдельных видов движения. При комнатной температуре из внутренних видов движения большинства молекул надо учитывать только вращательное.
в) Вращательная сумма по состояниям зависит от симметрии молекулы. В простейшем случае, в модели жесткого ротатора, которая описывает линейные молекулы с постоянными межъядерными расстояниями, уровни энергии зависят только от вращательного квантового числа J:
EJ = hcBJ(J + 1),
где B – вращательная постоянная (размерность – см–1), которая определяется моментом инерции I молекулы:
Каждый вращательный уровень имеет статистический вес gJ = 2J + 1.
Если ввести эффективную вращательную температуру
Tвращ = hcBk = 1.44B (см−1 ) ,
то вращательная сумма по состояниям приобретает вид:
|
∞ |
|
|
Tвращ |
|
|
Qвращ = ∑(2J + 1)exp − |
|
J (J + 1) |
|
T |
|
J = |
0 |
|
|
(15.36)
(15.37)
(15.38)
(15.39)
При не очень низких температурах (T >> Tвращ) суммирование в (15.39) можно заменить интегрированием по J, что дает:
Qвращ = |
T |
= |
kT |
. |
(15.40) |
Tвращ |
|
|
|
hcB |
|
232 |
Г л а в а 4. Статистическая термодинамика |
Для симметричных молекул это значение надо поделить на число симметрии (для двухатомных гомоядерных молекул оно равно 2), так как не все вращательные уровни энергии доступны.
При низких температурах (T ≤ Tвращ) вращательную сумму по состояниям находят суммированием по нескольким низшим значениям J.
г) Простейшая модель, описывающая колебания ядер, – гармонический осциллятор. Уровни энергии в этой модели линейно зависят от колебательного квантового числа:
En = hcωn,
где ω – частота колебаний (в см–1); энергия состояния с n = 0 принята за точку отсчета. Колебательные уровни энергии невырождены, статистический вес равен 1. Если определить эффективную колебательную температуру для данного колебания:
Tкол = hckω = 1.44ω (см−1 ) ,
то сумма по состояниям гармонического осциллятора с частотой ω приобретает вид:
|
|
∞ |
|
|
Tкол |
|
|
|
|
|
Tкол |
−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(15.43) |
Qω = ∑exp |
− |
|
n |
= |
1 |
− exp |
− |
|
. |
|
T |
T |
|
|
n=0 |
|
|
|
|
|
|
|
|
|
Эта сумма заметно отличается от 1 только тогда, когда дробь в показателе экспоненты меньше 1, т.е. для температур T > Tкол. При высоких температурах, когда Tкол/T << 1, колебательная сумма по состояниям прямо пропорциональна T:
|
|
T >>Tкол |
T |
|
|
(15.44) |
Qω → |
|
. |
|
Tкол |
|
|
|
|
Если же температура ниже эффективной колебательной температуры, то сумма по состояниям почти равна 1, в таких случаях говорят, что «колебание заморожено».
В молекуле, состоящей из n атомов, происходит 3n – 6 (в линейной молекуле 3n – 5) разных колебаний, каждое со своей частотой ωi, поэтому колебательная сумма по состояниям молекулы равна произведению сумм по состояниям для каждого из этих колебаний:
(15.45) |
Qкол = |
|
|
− exp |
|
− |
Tкол,i |
−1 |
|
1 |
|
|
. |
|
|
|
|
∏ |
|
|
T |
|
|
|
|
i |
|
|
|
|
|
|
д) Электронные и ядерные уровни энергии в молекуле обычно отстоят очень далеко друг от друга (см. табл. 15.1), и при не слишком высоких температурах вклад в соответствующую сумму по состояниям вносит только основной уровень, энергия которого принимается равной 0. Элек-
234 |
Г л а в а 4. Статистическая термодинамика |
поведения. Давление системы определяется только конфигурационным интегралом:
|
|
|
∂ln Z |
конф |
|
|
(15.51) |
p = kT |
|
. |
|
∂V |
|
|
|
T |
Точный теоретический расчет статистической суммы газов или жидкостей с произвольным гамильтонианом (15.48) – задача, которая лежит далеко за пределами возможностей современной статистической теории. Тем не менее, можно сделать ряд разумных и достаточно хороших приближений, которые позволяют оценить статистическую сумму (15.49) и конфигурационный интеграл (15.50) для реальных газов, состоящих из валентно насыщенных молекул.
Главное из этих приближений состоит в том, что общая потенциальная энергия взаимодействия частиц может быть представлена в виде суммы парных потенциалов, каждый из которых зависит только от расстояния между двумя частицами rij = |qi – qj|:
|
N |
N −1 |
(15.52) |
V (q1 ,..., qN ) = ∑∑u(rij ) . |
|
i> j |
j=1 |
Для того, чтобы конфигурационный интеграл с потенциалом (15.52) сходился, необходимо, чтобы парный потенциал при больших расстояниях стремился к нулю быстрее, чем r–3. Фактически это означает, что существует некоторая конечная область действия межмолекулярного потенциала, за пределами которой потенциал можно считать равным 0. Это – довольно мягкое требование, которое удовлетворяется для большинства реальных газов. Некоторые модельные парные потенциалы, которые используются при расчете свойств реальных газов, перечислены в табл. 15.2 и изображены на рис. 15.1.
Таблица 15.2 Наиболее распространенные межмолекулярные потенциалы
Название |
Формула |
|
|
|
|
Потенциал твердых сфер |
|
|
∞, |
r ≤ r0 |
|
|
|
u(r) = |
0, |
r > r0 |
|
|
|
|
|
|
|
|
Потенциал прямоугольной ямы |
∞, |
r ≤ r0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u(r) = −u0 , r0 < r < r1 |
|
|
0, r ≥ r |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
Потенциал Леннард-Джонса |
|
r0 |
|
12 |
|
r0 |
|
6 |
|
u(r) = u0 |
|
|
|
|
|
− 2 |
|
|
|
|
|
|
|
|
|
|
|
r |
|
|
|
|
r |
|
|
Потенциал Сазерленда |
∞, r < r |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
u(r) = |
|
|
−m |
, r |
≥ r |
|
|
|
−cr |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
exp – 6 потенциал |
u(r) = b exp (−ar )− cr −6 |
|
236 |
Г л а в а 4. Статистическая термодинамика |
При больших r каждый сомножитель стремится к 1, поэтому довольно удобно выделить единицу из экспоненты и определить функцию Майера1 следующим образом:
f (rij ) = |
fij |
|
|
u(rij |
) |
|
= exp |
− |
|
|
|
− 1 . |
kT |
|
|
|
|
|
|
|
|
При малых r потенциал стремится к бесконечности, а функция Майера к –1; при больших r и потенциал, и функция Майера довольно быстро стремятся к нулю. На рис. 15.1 изображены функции Майера для некоторых межмолекулярных потенциалов.
Подставляя определение (15.54) и разложение (15.53) в определение конфигурационного интеграла, находим:
Z конф = ∫ ∏ (1 + |
|
∑ |
fij |
+ ∑∑ fij fi ' j ' |
|
d 3N q . |
fij )d 3N q =∫ 1+ |
+ ... |
N ≥i> j≥1 |
|
N ≥i> j≥1 |
|
|
|
|
Скобка под знаком интеграла содержит 2N(N–1)/2 слагаемых, каждое из которых представляет собой либо функцию Майера для двух конкретных частиц, либо произведение нескольких (вплоть до N) таких функций.
Введем еще одно, довольно грубое приближение, которое справедливо только для сильно разреженных газов, а именно: в интеграле (15.55) пренебрежем всеми произведениями функций Майера:
Z конф |
|
∑ |
|
d 3N q . |
≈ ∫ 1+ |
fij |
|
|
N ≥i> j≥1 |
|
|
|
Это приближение означает, что в конфигурационном интеграле мы |
|
учитываем только те конфигурации, при которых две частицы находят- |
|
ся рядом, а остальные удалены друг от друга на значительные расстоя- |
|
ния. Это соответствует состоянию сильно разреженного газа. |
|
Интеграл в (15.56) разбивается на сумму интегралов. Интеграл от |
|
1 по координатам всех частиц равен VN. Каждая функция Майера fij за- |
|
висит от координат только двух частиц: i-й и j-й. Интеграл по коорди- |
|
натам оставшихся (N – 2) частиц дает множитель VN – 2. В оставшемся |
|
интеграле по d 3qi d 3q j |
можно перейти к сферическим координатам. |
|
В результате получим: |
|
(15.57) |
∫ fij d 3N q = V N − 2 ∫ f (qi |
− q j )d 3 qi d 3 q j = V N −1 ∫4πr 2 f (r) dr = V N −1β |
1 Эту функцию ввел Г. Урселл, а Дж. Майер и М. Гепперт-Майер использовали ее для построения общей статистической теории реальных газов. Иногда эту функцию называют функцией Урселла-Майера.