|
|
E |
n |
|
|
Z = ∑exp |
− |
|
|
(3.43) |
|
|
|
||||
|
|
kT |
|
||
называется статистической суммой (статсуммой). Отсюда, с учетом (3.42), связь свободной энергии со статистической сумой имеет вид
F = −kT ln Z . |
(3.44) |
Таким образом, статистическое определение энтропии и распределение Гиббса позволяют получить статистические определения свободной энергии и статистической суммы (статсуммы).
Обсудим некоторые аспекты нормировки функции распределения. С физической точки зрения бывает удобно рассматривать систему в фазовом пространстве, как совокупности обычного трехмерного пространства координат и соответствующего ему пространства обобщенных импульсов. При этом функция распределения имеет смысл средней локальной концентрации частиц в фазовом пространстве.
Если в системе N частиц, то по определению объем фазового пространства всего ансамбля частиц равен
Φ = Nh s , |
(3.45) |
где h – постоянная Планка, s – число степеней свободы (в трехмерном декартовом пространстве s = 3).
При этом если нормировку (3.40) записать в интегральном виде (как интеграл от концентрации точек ансамбля по фазовому про-
странству) |
|
|
|
|
∫ |
f (p, q) |
dpdq |
= 1 . |
(3.46) |
|
||||
Φ |
|
Φ |
|
|
|
|
|
|
|
Если функция распределения везде в фазовом пространстве равна единице, получим тождество
∫ |
Φ |
|
|
|
dpdq |
= 1. |
(3.47) |
|
|
||
Φ
С учетом специфики интересующих нас приложений, когда представляет интерес рассмотрение системы не как точки в 6N-мерном
56
фазовом пространстве, а как ансамбля N частиц в 6-мерном фазовом пространстве, объединяющем трехмерное евклидово конфигурационное пространство с сопряженным ему пространством импульсов, перепишем (3.46) в виде
∫ |
f (p, q) |
dpdq |
= N . |
(3.48) |
||
|
s |
|||||
Φ |
|
|
h |
|
||
|
|
|
|
|
|
|
Обозначая |
|
|
|
|
|
|
dΓ = |
dpdq |
, |
(3.49) |
|||
|
||||||
|
|
h s |
|
|||
получим |
|
|
|
|
|
|
∫ f (p, q)dΓ = N . |
(3.50) |
|||||
Γ |
|
|
|
|
|
|
Данным видом нормировки мы будем пользоваться в дальнейшем.
57
Глава 4. Термодинамические функции реального газа
4.1. Статистическая сумма
Рассмотрим реальный газ, представляющий собой совокупность частиц, возможно, имеющих внутреннюю структуру, с функцией распределения в фазовом пространстве, соответствующей каноническому распределению Гиббса:
|
|
H (p, q) |
|
|
f (p, q) = A exp |
− |
|
, |
(4.1) |
|
||||
|
|
kT |
|
|
удовлетворяющей условию нормировки |
|
|
|
|
∫ f (p, q)dΓ = N , |
|
(4.2) |
||
Γ |
|
|
|
|
где dΓ = h −s dpdq – элемент объема фазового пространства сис-
темы, по которому производится интегрирование; s – размерность координатного или сопряженного ему импульсного пространств; N – число частиц в системе; h = 2πH – постоянная Планка. Отсюда для нормировочной константы A получаем
A = |
|
|
|
Nh s |
|
||
|
|
|
|
, |
(4.3) |
||
|
|
|
|
||||
|
|
|
H (p, q) |
|
|||
|
∫ exp |
− |
|
dpdq |
|
||
|
|
||||||
|
Γ |
|
|
kT |
|
||
где H – гамильтониан системы (энергия частицы с координатами q и сопряженными импульсами p). Заметим, что представление (4.1) соответствует каноническому распределению Гиббса, а функция (4.1) тождественно удовлетворяет уравнению Лиувилля для равновесной системы.
Если гамильтониан может быть представлен в виде
H = K (p) + u(q) , |
(4.4) |
58
т.е., как обычно, в виде суммы кинетической и потенциальной энергий частицы или микросистемы, то интеграл (4.2) разбивается на произведение двух интегралов:
|
|
− |
Ah −s ∫ exp |
||
P |
|
|
K (p) |
|
|
|
|
|
dp∫ exp |
− |
||
|
||||
kT |
Q |
|
|
|
u(q) |
|
|
|
dq = N , |
(4.5) |
|
||
kT |
|
|
где интегрирование производится по координатному и импульсному пространствам раздельно.
Если имеется некоторая физическая величина g(p, q) , характе-
ризующая отдельную частицу (микросистему) с соответствующими координатами в фазовом пространстве, то ее значение для всей системы, очевидно, можно найти по формуле
G = h− s ∫ g(p, q) f (p, q)dpdq . |
(4.6) |
Γ |
|
Функцию распределения (4.1) можно представить в виде произведения функций распределения по неполным пространствам, на-
пример, по импульсному пространству – n p и по координатному пространству nq :
f (p, q) = n p (p)nq (q) . |
(4.7) |
Рассмотрим фазовый интеграл для поступательных степеней свободы молекул в трехмерном евклидовом пространстве.
Если распределение частиц в пространстве импульсов симметрично относительно нулевого значения импульса, то с учетом того, что
|
|
|
|
|
K ( p) = |
p 2 |
|
|
|
|
(4.8) |
||
|
|
|
|
∞ |
2m |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
K (p) |
|
p 2 |
|
2 |
3 |
|
|||||
|
|
|
|
|
|||||||||
∫ exp − |
|
|
|
|
− |
|
|
|
|
|
2 |
, (4.9) |
|
|
dp = 4π ∫ exp |
|
|
|
p dp = (2πmkT ) |
|
|
||||||
P |
|
kT |
0 |
|
|
2mkT |
|
|
|
||||
что соответствует распределению Максвелла.
После интегрирования (4.1) по импульсному пространству, с учетом (4.7) получаем
59
n(x, y, z) = (2πmkT ) |
3 |
h−3 A exp |
− |
u(x, y, z) |
|
|
|
||||||||
2 |
, |
(4.10) |
|||||||||||||
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
kT |
|
|
|||
подставляя (4.9) в (4.5), находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
− |
3 |
|
|
|
|
|
|
|
|
A = |
nh 3 (2πmkT ) 2 |
. |
|
|
|
(4.11) |
|||||||||
|
|
|
|
u(q) |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
||||||
|
∫ exp |
− |
|
dq |
|
|
|||||||||
|
|
|
|
||||||||||||
|
Q |
|
|
|
kT |
|
|
|
|
|
|
||||
Заметим, что если потенциальная энергия частицы постоянна по всему объему, то из (4.11) получаем
− 3
A = nh3 (2πmkT ) 2
|
u |
|
|
exp |
|
, |
(4.12) |
|
|||
kT |
|
||
где n – объемная концентрация в обычном трехмерном пространстве. С учетом (4.7) при этом концентрация в импульсном пространстве тоже постоянна.
Если система представляет собой две фазы с потенциальными энергиями частиц U1 и U2 и объемамиV1 и V2, то
|
|
|
|
|
|
|
|
− |
3 |
|
|
|
|
|
|
A = |
|
Nh 3 (2πmkT ) |
2 |
|
|
|
|
. |
(4.13) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
− |
u1 |
+ V2 |
|
|
|
|
− |
u 2 |
|
|
|||
|
|
|
|
|
|
||||||||||
V1 |
exp |
|
|
exp |
|
|
|
||||||||
|
|
|
|||||||||||||
|
|
|
|
kT |
|
|
|
|
|
|
kT |
|
|||
Энергия системы в соответствии с (4.6) может быть найдена по формуле
E = ∫ H (p, q) f (p, q)dpdq . |
(4.14) |
Γ |
|
Кинетическая энергия для случая постоянной потенциальной энергии определяется формулой
K = h−3 ∫ K (p) f (p, q)dpdq =
Γ |
|
− u ∞ |
|
|
|
|
(4.15) |
= 4πAVh−3 |
|
|
|
K ( p) |
|||
|
|
||||||
exp |
∫ K ( p) exp |
− |
|
p 2 dp, |
|
||
|
|
||||||
|
kT 0 |
|
|
kT |
|
||
|
|
|
60 |
|
|
|
|