Метода термодинамика(лекции)
.pdfδ = kTmg (для газов δ (1 ÷ 30) 104 м). На всех высотах имеет место
максвелловское распределение по скоростям с постоянной температурой T . На первый взгляд может показаться, что температура с высотой меняется. Когда молекула поднимается вверх, ее потенциальная энергия
mg(z −z0 ) увеличивается и должна меняться кинетическая энергия
mv2 02 − mv2 2
Но кинетическая энергия связана со скоростью молекулы. Поэтому на большей высоте молекула будет иметь меньшую скорость и энергию. А температура определяется средней квадратичной скоростью (см. п. 2.4.3)
v 2 = 3kT m . Следовательно, температура должна уменьшаться с высотой.
Ошибочность этого рассуждения коренится в том, что в нем фигурирует только одна молекула, рассматриваемая без учета ее столкновений с другими молекулами. Максвелловское распределение скоростей устанавливается благодаря столкновениям. Здесь столкновительный процесс полностью игнорируется и рассматривается не имеющая смысла «температура молекулы». Фактически подниматься вверх будут преимущественно те молекулы, которые имеют большую скорость, и максвелловское распределение устанавливается автоматически на всех высотах.
Определим распределение по высоте молекул, обладающих различной массой. Чем больше масса молекулы, тем быстрее убывает их число с
высотой (2.27). Если число молекул с массами m1 |
и m2 на высоте z = 0 |
||||||
одинаково n1(m1, 0) = n2(m2, 0) = n0 , то |
отношение |
числа молекул обоих |
|||||
сортов на высоте h равно |
|
|
|
|
|
|
|
|
n1(m1,h) |
=e− |
(m1 −m2 ) |
gh . |
(2.28) |
||
|
kT |
||||||
|
|
||||||
|
n |
(m |
,h) |
|
|
|
|
2 |
2 |
|
|
|
|
|
Отношение означает, что доля более тяжелых молекул в равновесном газе увеличивается по мере приближения к центру гравитации.
82
2.7.3. Вес и потенциальная энергия . Найдем вес газа, который находится в сосуде – цилиндре высотой h и площадью основания S . Можно поступить двояким образом. Во-первых, считать вес столба газа равным весу всех входящих в него молекул, во-вторых, найти разность давлений, оказываемых газом на дно (z = 0 ) и крышку (z = h ) сосуда. В первом случае имеем
|
|
|
|
|
h |
mgz |
|
F = mg ∫dN(z,x,y) = mgn0 ∫e− kT dz ∫dxdy = |
|
||||||
|
|
|
|
|
0 |
S |
(2.29) |
|
|
kT |
|
−mgh |
|
|
|
|
0 |
|
( |
) |
|
||
|
|
|
|
|
|||
= Smgn |
|
|
−e |
kT |
= SkT n(0) −n(h) , |
|
|
|
1 |
|
|
||||
|
|
mg |
|
|
|
|
|
|
|
|
|
|
|
|
|
где S |
площадь сечения сосуда, |
а n(0) и n(h) – плотности газа на высоте |
||||
z = 0 |
и z = h , соответственно. |
Во втором случае, используя уравнение |
||||
состояния газа p(z) = n(z) kT , можно записать |
||||||
|
F = S |
( |
|
) |
( |
) |
|
|
p(0) −p(h) |
= S n(0) |
−n(h) kT , |
где p(0) и p(h) – давления на соответствующих высотах. Последние выражения совпадают, что означает: вес газа суть разность давлений на дно и крышку сосуда. В этом заключается различие между весом газа и весом тела, лежащем на тарелке весов.
Молекулы газа в поле тяжести обладают некоторой средней потенциальной энергией, избыточной по сравнению с энергией газа в отсутствие внешних сил. Средняя потенциальная энергия газовой молекулы в поле тяжести равна
|
|
|
|
mg |
|
|
|
|
= mgz = mg∫ zdw(z) = mg |
∫ ze−kT zdz |
|
||
|
u |
∫e |
−mg z |
. |
(2.30) |
|
|
|
|
kT |
dz |
|
Результат вычисления существенно зависит от высоты столба газа. Если высота значительно больше характеристической длины распределения
83
z >> δ = mgkT , то пределы интегрирования по распространяются от z = 0 до
z → ∞ и
u = kT .
Средняя потенциальная энергия одной молекулы в бесконечном столбе газа оказывается пропорциональна температуре. Средняя потенциальная энергия грамм-моля газа в бесконечном столбе равна
U = Nu = NkT ,
т.е. сравнима со средней кинетической энергией ( 3NkT2 ).
Иной результат получается для сосуда с малой высотой h << δ . Показатель экспоненты в исходной формуле очень мал mgz / kT = z / δ << 1, и подынтегральную функцию можно разложить в ряд, ограничившись первым членом разложения:
u = mgh2 ; U = Nu = Nmgh2 .
В этом приближении средняя потенциальная энергия столба газа вообще не зависит от температуры.
Сопоставление свойств идеального газа и земной атмосферы показали на их существенное несоответствие: состав атмосферы не меняется с высотой, как это следует из формулы (2.28). Общеизвестно понижение температуры с увеличением высоты, что противоречит требованиям постоянства температуры.
Эти противоречия ставят два вопроса: 1) правильны ли были основные положения теории; 2) если да, то применимы ли полученные закономерности к явлениям, происходящем в атмосфере. На первый вопрос должен, без сомнения, следовать положительный ответ. Поскольку изложенная теория равновесного идеального газа получила подтверждение в газовых законах, при экспериментальной проверке распределения Максвелла, длины свободного пробега и т.п. Что касается второго вопроса, то ответ –
84
отрицательный. Земная атмосфера не является адекватной моделью равновесного газа. Земная атмосфера явно система неравновесная; явления, происходящие в ней, весьма сложные: это и конвекционные потоки, солнечное излучение, механическое перемешивание слоев, связанное с движением Земли, и другие осложняющие факторы.
2.8.Число и функция состояний идеального газа
2.8.1.Число состояний и критерий невырожденности. Для расчета числа состояний моделируем идеальный газ квазиклассической системой. Каждому состоянию свободной частицы в фазовом пространстве
(три степени свободы, f = 3 ) отвечает элементарная ячейка объемом h3 (h –
постоянная Планка). Объем фазового пространства, которому соответствует в координатном пространстве объем V , где может находиться частица, а
импульс изменяется в интервале от p |
до p +dp , равен ∆Γ = 4πp2dp V . |
||||
Тогда число состояний определяется формулой |
|
|
|||
Ω(p)dp = 4πp2dp V |
, |
(2.31.а) |
|||
|
|
h3 |
|
|
|
или, учитывая связь между энергией частицы и импульсом ( ε = p2 |
2m ), |
||||
|
|
|
3/2 |
|
|
|
2m |
εdε. |
(2.31.б) |
||
Ω(ε)dε = 2πV |
|
|
|
||
|
|
2 |
|
|
|
h |
|
|
|
|
|
Плотность числа состояний свободной частицы пропорциональна |
ε . Если в |
качестве частицы рассматривать электрон, который имеет в каждой ячейке два состояния, отличающиеся направлением спина, число его состояний удваивается, Ωэ(ε) = 2Ω(ε) . Формулы легко обобщаются на системы с произвольным числом степеней свободы
Ω(ε)dε = ∆Γhf = h1f ∂Γ∂ε dε.
85
Общее число |
состояний |
свободной |
частицы в интервале энергий от |
|||||||||||||
ε = 0 до ε = 3kT 2 определяется интегралом |
|
|
|
|
|
|||||||||||
|
ε |
|
|
|
|
3/2 |
|
|
|
|
|
|
|
3/2 |
|
|
|
|
|
|
|
|
|
2 |
|
|
3/2 |
|
|
|
|
||
G = |
∫ |
|
|
2m |
|
ε |
|
2πmkT |
(2.32) |
|||||||
Ω(ε)dε = 2πV |
|
|
|
|
|
≈V |
|
|
. |
|||||||
|
|
|
|
2 |
|
|
3 |
|
|
|
|
h |
2 |
|
|
|
|
0 |
|
h |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
При числе частиц N в объемеV критерий невырожденности принимает вид |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
3/2 |
|
|
|
|
|
|
|
N |
|
|
|
h2 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.33) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
= n |
|
|
|
|
. |
|
|
|
|||||
|
|
G |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
2πmkT |
|
|
|
|
|
Для одноатомного молекулярного газа, например, азота при нормальных
условиях |
(n = 1026 м–3, |
m = 4, 5 10−26 кг, |
h = 6, 6 10−34 Дж·с, |
kT = 4 10−21 |
Дж) N /G 10−6 . Такие малые |
значения критерия |
свидетельствуют, что обычные молекулярные газы в нормальных условиях являются невырожденными и описываются статистикой Максвелла-
Больцмана. Электронный газ в металлах (n = 5 1028 м–3, m = 9 10−31 кг)
оказывается |
вырожденным вплоть до температур T ≤105 |
К |
(N /G 0, 5 ) и |
|
в реальных |
условиях подчиняется квантовой |
статистике |
Ферми-Дирака. |
|
Невырожденное состояние электронного газа |
может |
быть достигнуто |
уменьшением концентрации n . Уже при n = 1022 |
м–3 |
критерий |
N /G 10−3 . Такая и меньшая концентрация электронов |
имеет |
место в |
невырожденных полупроводниках. |
|
|
2.8.2. Функция распределения идеального газа суть плотность частиц в данном состоянии. Связь между термодинамическими параметрами системы и характеристиками микрочастиц задается полной статистической функции распределения F(ε)dε, имеющей смысл числа частиц с энергией от ε до ε +dε. Её можно представить произведением числа состояний Ω(ε)dε на вероятность их заполнения: F(ε)dε = f (ε)Ω(ε)dεdε.
Функцию f (ε) называют функцией распределения. Кроме вероятности
86
заполнения состояний, она имеет смысл среднего числа частиц, приходящегося на одно состояние. Из известных распределений Максвелла dw(p) и dw(ε) для импульсов и энергий одноатомной молекулы (2.5) и (2.7)
нетрудно получить число соответствующих свободных молекул в объеме V (nV dw), т.е. полную функцию распределения:
|
|
|
|
|
|
|
|
|
|
|
|
|
3/2 |
|
p |
2 |
|
|
|
|
|
||
|
|
|
N |
|
|
|
h2 |
|
|
|
|
|
|
− |
|
|
4πp2dpV |
||||||
|
|
|
|
|
|
|
|
|
|
2mkT |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
dn(p) |
= |
|
|
|
|
|
|
|
|
e |
|
|
|
|
|
|
|
|
3 = |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
V |
|
2πmkT |
|
|
|
|
|
|
|
|
|
h |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
3/2 |
|
|
ε |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
N |
|
h2 |
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
= |
|
|
|
|
|
|
e |
|
|
kT |
Ω(ε)dε = f (ε)Ω(ε)dε. |
||||||||||||
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
2πmkT |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Функция распределения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3/2 |
|
|
ε |
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
h2 |
|
|
|
|
− |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
fM (ε) = |
|
|
|
|
|
|
|
|
|
|
e |
|
kT |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
2πmkT |
|
|
|
|
(2.34)
(2.35)
получила название функция распределения Максвелла-Больцмана. Её можно записать и в другом виде – в форме распределения Гиббса для системы с переменным числом частиц, для чего воспользуемся тождеством
|
|
|
N |
|
|
|
3/2 |
|
|
|
|
µ/kT |
|
|
|
h2 |
|
|
|
||
e |
|
≡ |
|
|
|
|
|
|
(2.36) |
|
|
|
|
|
|
||||||
|
|
|
|
|
. |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
2πmkT |
|
|
||||
Распределение Максвелла-Больцмана принимает вид |
|
|||||||||
|
|
|
|
|
|
µ−ε |
|
|
|
|
|
|
f |
(ε) =e kT . |
|
|
|
(2.37) |
|||
|
|
M |
|
|
|
|
|
|
|
|
Введенный здесь параметр |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
h2 |
3/2 |
|
||||
|
|
|
N |
|
|
|
|
|||
µ = kT ln |
|
|
|
|
|
, |
(2.38) |
|||
|
|
|
||||||||
|
|
|
|
|||||||
|
|
|
V |
|
|
|
|
|
|
|
|
|
|
2πmkT |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
как будет показана в |
следующих |
разделах, является |
химическим |
|||||||
(парциальным) потенциалом. |
|
|
|
|
|
|
|
|
|
2.8.3. Функция ( интеграл ) состояний . Согласно определению, функция состояний идеального одноатомного газа в
квазиклассическом приближении имеет вид
87
Z = ∫e−E /kTdΩ(E). |
(2.39) |
|||
Поскольку все частицы газа являются независимыми, то |
||||
|
s |
dq dp |
||
dΩ = Π′dΩi = Π′ |
i |
i |
, |
|
3s |
|
|||
i |
i |
h |
|
|
|
|
|
где произведение берется по всем координатам и импульсам частиц. Штрих в произведении означает, что при его образовании следует включать в него только такие члены, которые отвечают различным состояниям системы в целом. Последующее интегрирование ведется по соответствующим областям фазового пространства
s |
dp dq |
i |
|
|
Z = ∫′e−E /kT Π′ |
i |
|
. |
|
h |
3s |
|
||
i |
|
|
|
|
|
|
|
|
Если бы молекулы в газе отличались друг от друга, то определенным значениям импульса и координаты частиц соответствовало одно состояние
газа, т.е. Π′ = Π. Но частицы однородного газа – молекулы, атомы ничем не
i i
отличаются друг от друга.
Для определения области интегрирования рассмотрим случай двух частиц. Пусть первая частица имеет импульс p1 и координату q1 , а вторая – p2 и q2 . В силу тождественности частиц состояние системы не изменится,
если значения импульса и координаты второй молекулы станут равными первой и наоборот. В такой системе из двух частиц следует интегрировать по всем состояниям – координатам и импульсам двух частиц ( p1 , q1 , p2 , q2 ), но полученный результат уменьшить вдвое. При этом автоматически будет учтено, что «перестановка» в фазовом пространстве изобразительных точек (состояний частиц) не приводит к разным состояниям системы.
Полученный результат обобщается на произвольные состояния газа, содержащего N тождественных молекул. Можно сказать, что совокупность всех состояний, которые получаются перестановкой между собой изобразительных точек в фазовом пространстве, являются тождественными
88
физическими состояниями – отвечают только одному физическому состоянию. Поэтому при выполнении интегрирования по всем возможным значениям координат и импульсов молекул нужно разделить результат на число перестановок изобразительных точек между собой, т.е. на N !. В частном случае для двух молекул N ! = 2 . Итак можно записать
dΩ = Π′dΩ = |
( |
ΠdΩ |
N ! = |
1 |
Πdq dp . |
|||
i |
i |
i |
i ) |
|
h3N N ! i |
i i |
Соответственно для функции состояний имеем:
|
1 |
|
e− |
E |
|
|
||
Z = |
|
|
ΠdΩ , |
(2.40) |
||||
|
kT |
|||||||
N ! ∫ |
||||||||
|
|
|
i |
i |
|
где интегрирование ведется по всему фазовому пространству. Учитывая, что полная энергия идеального газа равна сумме энергий составляющих его
молекул E = ∑εi , находим
i
|
|
|
Z = |
1 |
|
|
1 |
|
∏ |
∫ |
e− εi /kTdpidqi |
= zN / N ! |
(2.41) |
||||||||
|
|
|
|
|
3N |
||||||||||||||||
|
|
|
|
|
N ! h |
|
i |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
3/2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где |
2πmkT |
|
– |
|
ранее |
определенный |
|
интеграл состояний |
одной |
||||||||||||
z =V |
|
|
|
|
|
||||||||||||||||
|
|
h |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
молекулы, поэтому |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3N /2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
2πmkT |
|
|
N |
. |
(2.42) |
|||||
|
|
|
|
|
|
Z = |
|
|
|
|
|
|
|
|
V |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
h |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
! |
|
|
|
|
|
|
|
|
Через эту функцию состояний находятся все термодинамические потенциалы, а следовательно, свойства идеального газа как системы, подчиняющейся классической статистике.
2.9. Классическая теория теплоемкости газа двухатомных молекул
Большая часть веществ в газообразном состоянии существует в виде неодноатомных молекул, например, H2 , O2 , HCl , CO , CO2 , H2O , NH4 .
89
Наша задача обобщение полученных в предыдущих разделах результатов на газ многоатомных молекул и, в особенности, на самую простую структуру – двухатомный газ. Многоатомная молекула суть устойчивое образование из двух, трех и большего числа одинаковых или разных атомов. Их особенность – наличие у них вращательных и колебательных степеней свободы. Считаем,
что молекула представляет систему, подчиняющуюся законам классической механики. Поскольку внутреннее движение электронов в атоме нас не будет интересовать, то каждый атом отождествляем материальной точкой.
Моделью двухатомной молекулы может служить миниатюрная гантель,
на концах которой находятся бесконечно малые шарики с массами m1 и m2
(различными массами атомов). Связь между атомами в молекуле не является абсолютно жесткой. По этой причине возможны следующие виды движения: поступательное (три степени свободы); вращение вокруг двух осей, перпендикулярных оси, соединяющей оба атома (две степени свободы); колебания атомов вдоль линии, их соединяющей. Общая энергия молекулы состоит из суммы энергий соответствующих движений
ε = εпост + εвр + εколеб .
Состояние молекулы определяется заданием шести пространственных координат и шести импульсов. Координаты (x,y,z) относятся к центру инерции, углы θ и ϕ определяют положение оси молекулы в пространстве,
координата q описывает отклонение атомов от равновесного расстояния a .
Импульсы, соответствующие этим координатам, равны
px , py , pz , M1 = I ω1 , M2 = I ω2 , pq = µq ,
где I = µa2 – момент инерции, µ = m1m2 (m1 +m2 ) – приведенная масса,
ω – угловая скорость молекулы. Соответствующее фазовое пространство имеет 2 (3N ) = 12 измерений. Считая молекулу квазизамкнутой системой,
воспользуемся распределением Гиббса
90
|
1 |
|
ε |
+ ε |
+ ε |
|
|
|
|
пост |
вр |
колеб |
dΓ. |
||
dw = |
|
exp − |
|
|
|
|
|
6 |
|
|
|
||||
|
|
|
kT |
|
|
||
|
h z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Элемент фазового пространства двухатомной молекулы может быть представлен в виде
dΓ = dΓпост dΓвр dΓколеб ,
где введены обозначения
dΓпост = dpxdpydpzdxdydz ;
dΓколеб = dpqdq .
Очевидно, что распределение Гиббса распадается на три независимых множителя, соответствующих поступательному, вращательному и колебательному движениям. Каждый из этих видов движения не влияет на другие и может рассматриваться независимо. Применяя теорему о
распределении |
средней энергии по степеням свободы, получим |
|||
εпост = 3 kT 2 , |
εвр = 2 kT 2 , |
εколеб = 1 kT . Этот |
результат |
можно |
получить через соответствующее распределение Гиббса |
|
|
||
|
εj = ∫ εjdw(εj ), |
|
|
|
учитывая, что, εвр |
=(Iω12 +Iω22 ) |
2 , εколеб = pq2 2µ+µω2q2 |
2 , где ω = |
χ µ – |
частота колебаний атомов, связанная с постоянной квазиупругой силы χ.
Поступательное движение двухатомной молекулы не отличается от движения одноатомной, поскольку оно сводится к движению центра тяжести системы. Таким образом, суммарное значение средней энергии двухатомной молекулы оказывается равным ε = 3, 5kT .
Алгоритм такого расчета может быть перенесен на случай многоатомных молекул. Рассмотрим, например, линейную трехатомную молекулу CO2 с девятью степенями свободы (рис. 2.7). Здесь имеют место следующие виды движения: поступательное (три степени свободы);
вращение вокруг двух осей, перпендикулярных оси молекулы (две степени
91