Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Метода термодинамика(лекции)

.pdf
Скачиваний:
37
Добавлен:
23.02.2015
Размер:
2.65 Mб
Скачать
= mg(z z0 ).

δ = 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 ekT 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 = mgzdw(z) = mg

zekT 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 1026 кг,

h = 6, 6 1034 Дж·с,

kT = 4 1021

Дж) N /G 106 . Такие малые

значения критерия

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

Больцмана. Электронный газ в металлах (n = 5 1028 м–3, m = 9 1031 кг)

оказывается

вырожденным вплоть до температур T 105

К

(N /G 0, 5 ) и

в реальных

условиях подчиняется квантовой

статистике

Ферми-Дирака.

Невырожденное состояние электронного газа

может

быть достигнуто

уменьшением концентрации n . Уже при n = 1022

м–3

критерий

N /G 103 . Такая и меньшая концентрация электронов

имеет

место в

невырожденных полупроводниках.

 

 

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 = eE /kTd(E).

(2.39)

Поскольку все частицы газа являются независимыми, то

 

s

dq dp

dΩ = Π′di = Π′

i

i

,

3s

 

i

i

h

 

 

 

 

 

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

s

dp dq

i

 

Z = ∫′eE /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

dΓвр = dM1dM2 sin θdθdϕ;

 

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