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

книги из ГПНТБ / Беккер Р. Теория теплоты

.pdf
Скачиваний:
15
Добавлен:
25.10.2023
Размер:
15.46 Mб
Скачать

 

Если g i , . . . ,

q?

означают

декартовы

координаты

Х\,

Ух, zu...,

 

X

j y

,

У1я, zN

 

N

частиц

(f=3N),

то, например, вы

ражение ~d§t\dx\

 

 

= Kx\

представляет

собой

 

проекцию

на ось х силы, действующей

на частицу № 1. Таким об­

разом,

последнее уравнение

означает:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

Й К/)

-

SNkT.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ - л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если к тому же принять,

что

подразумеваемое

здесь

усреднение по ансамблю частиц идентично

усреднению

по

времени

(32),

то

мы

получим

теорему

о

 

внриале

(28.2).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

Еще

раз

о максвелловском

распределении

 

скоростей

Мы

можем убедиться,

 

что

 

рассматривая

микроканонический

ан­

самбль с N свободными атомами газа, можно действительно полу­

чить распределение Максвелла для вероятности одной

компоненты

скорости | одного атома газа. Однако математически точное выра­

жение для этого

распределения

отсюда

получить

нельзя,

так

как

в распределении

Максвелла

е~^"

 

любое, даже весьма большое зна­

чение 4, может встретиться с конечной

вероятностью,

в

то

время

как

в микроканоническом ансамбле кинетическая энергия

т £ 2 / 2

никогда

не может стать больше постоянной энергии

Е

всего

 

ансамбля.

N

атомов газа

имеют

 

j = 3N

степеней свободы. При

числе

ЗЛ/

компо­

нент

импульсов pi,Pf

 

 

зависимая

от

импульсов часть

функции

Га­

мильтона равна

(р\

-f-...+p2 )/2/«.

Д л я

математического

 

описания

микроканонического

ансамбля используем б-функцию,

определяемую

с помощью выражений

6 ( s ) = 0

для s=^=0 и

j 6 ( s ) r f s = l

*. С

ее

по-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— оо

 

 

 

 

 

 

 

 

мощью микроканонический ансамбль описывается следующим об­

разом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(Pi

. . • •.

P\)dPx-

 

• • dpy-Ь

( р\

-|

 

VP)

2тЕ)

dpr..

 

 

dpf.

 

 

В данный момент нас интересует только одна компонента,

на­

пример

р/.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда с точностью до независимого от р/ множителя

 

вероят­

ность g(Pi)dpf

того,

что

pf

находится

в

интервале

dplt

 

будет

равна:

 

g(pf)dpj^dpf

 

 

+."°

 

 

 

 

 

 

 

 

 

 

dpf_,.

 

 

 

 

 

 

j

•••

j p ( p l , . . . , p f ) d p 1

, . . . ,

 

 

 

 

* Таким

образом,

для

любой

функции

U(s)

будет

выполняться

 

 

ь

 

 

 

 

 

 

 

 

 

ь

 

 

 

 

 

 

 

 

 

 

 

условие

j U(s)6s

ds — U(0)

или

\U(s)6sds

=

0

в зависимости

от

 

 

а

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

того, расположена ли нулевая точка в пределах интервала от а до Ь или нет.

180

Используя вышеприведенное выражение для р, получаем:

e(pf)=

J

•••^{Р1+•••+p2f-l-{2mE-p2f))dp1...dpl_1.

Подынтегральное выражение

зависит

только от значения векто­

ра г (/ 1)-мерного пространства:

r2=p\

. В этом простран­

стве объем шара радиусом г пропорционален величине г*~1, объем «шарового слоя» dr равен Gr>-2dr. Если мы обозначим временно

и — 2тЕ~р2,

то, опять-таки

с

точностью

до

численного

множителя

 

 

 

g (Pf) =

JS

2

-

a)

r^-dr.

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

Теперь

введем

в

качестве

независимой

переменной

s аргумент

5-функции,

т. е. примем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г 2 a =

s;

г

=

(s +

 

1/

 

dr—

 

ds

 

 

 

 

a ) / 2

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2V'i

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g(P/)

=

 

- j

"

\

 

6 ( s ) ( s + a )

 

ds

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g (P/)

=

a

2

при

a

>

0

и g

(pf)

= 0

 

 

 

 

 

 

 

 

 

 

при

a

< 0.

 

 

 

 

 

 

 

Подставив значения

и,

выделив

не

зависящий

от

pf

множитель

£-_3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2тЕ)

 

и опустив

индекс f, получим

окончательно:

 

 

 

 

 

g(p)dp

=

C^l—

 

 

 

 

2

Ф

для

 

 

 

 

 

 

 

 

 

 

 

р 2 < 2 т £ .

 

 

 

 

 

 

 

 

Данная

формула

строго

выполняется

для любого

числа

f > l .

J. Если

теперь

/

бесконечно

 

велико,

то g(p)

заметно

отличается

от нуля

лишь тогда,

когда

величина

р2/2тЕ

очень

мала. Но

в этом

случае

 

 

 

 

 

 

 

 

 

 

 

 

f-з

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р 2

 

 

 

 

 

 

 

 

 

 

 

 

g(p)xe

 

 

•1тЕ

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р!

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если пренеоречь

числом 3 по сравнению с /, то g = e

2m£

2

 

 

Но из закона равнораспределения мы знаем, что E=fkTj2.

Таким об­

разом

при

больших

/ мы

действительно

получили

распределение

Максвелла.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

181

34.Э Н Т Р О П И Я

а) «Параметр» в ffl-функции.

Пусть функция Гамильтона, кроме координат и импуль­ сов, зависит от одного или нескольких «параметров» а. Под параметром мы понимаем величину, которая в нор­

мальном случае не изменяется

при движении, однако

ее численное значение можно

выбирать произвольным

путем воздействия извне. Примерами таких величин яв­ ляются, например, объем сосуда, в котором заключена наша система или налагаемое извне магнитное поле.

Если для функции Ж(а, t, а) параметр изменяется во времени, то это изменение не влияет на уравнения дви­ жения (31.1). Из этого следует, что для изменения Ж, а следовательно, и энергии во времени в каждое мгновение справедливо соотношение

dWe

дЖ

da

 

дЗ€

,

dt

да

dt

т.е. аЖ=

да

da,

Если за время х а изменилось на очень малую вели­ чину da, то для приращения Е за время т имеет место соотношение

t+z

 

 

 

dE

 

=

\

да

adx'.

 

 

 

 

 

 

 

 

 

 

.1

 

 

 

 

 

 

г-,

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

da

и исполь-

 

Подставляя

для малого изменения а =

 

 

 

 

 

 

 

д§?т

1

Г

 

дше ,

,

зуя

усреднение

по

в р е м е н и - ^ -

= —

] - ^ - а т

,получа­

ем

для изменения

Е

 

 

 

 

 

 

 

 

 

 

 

 

,

с

дЖт

,

 

 

 

 

 

 

 

 

dt

 

да

da.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дЖх

Наконец, при достаточно

медленном

изменении для

 

 

 

 

 

 

 

 

 

 

 

 

 

можно принять усреднение по микроканоническому

ансамблю:

 

 

dE

= —da.

 

 

 

 

(34.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

да

 

 

 

 

'

 

Если, например, подставить а— V, то, следовательно,

выражение р=

Шт

 

должно

иметь значение

давления.

182

 

 

dV

 

 

 

 

 

 

 

 

 

Поясним ситуацию на одном специальном примере (рис. 63).

Возьмем цилиндр с поршнем, имеющий единичное поперечное сече­

ние. Координата х отсчитывается от днища

сосуда

(х =

0)

до

порш­

ня (х — 1). В этом

случае / численно

равно

объему. Пусть Xj, у , , г,

будут координаты

/-го атома. Для того чтобы этот атом

действи­

тельно отражался

от поршня,

в функцию

должна

входить

соот­

ветствующая потенциальная

энергия

[см.

С т е н к и

в

(31.2)].

Она

может, например, иметь форму

Рис. 63. Позиция I

поршня как «пара­ метр» В еЖ-фуНКЦИИ.

с положительными коэффициентами с и А. Если А очень велико, то эта энергия при X j < l практически равна нулю, а при X j > l бесконеч­ но велика. Мы можем тогда в функции <Я? сумму

 

N

 

W(xv...,xN;

1) = с Л

e-A('~xi)

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

дЖ

dW

al

dl

Она, естественно, равна по величине и противоположна по знаку сумме всех сил Ki, которые испытывают отдельные молекулы со сто­ роны стенки:

N N

/=1 s=l

Приведенную в (34.2), усредненную во времени силу мы ощу­ щаем как давление. В каждое мгновение суммирование производит­ ся лишь для тех атомов, которые находятся вблизи стенки (XJ почти равно / ) .

б) Адиабатная инвариантность Ф *

Одновременно с Жот а будет зависеть и фазовый объем Ф*. Ф* становится тем самым функцией двух перемен­ ных Е и а и определяется с помощью выражения

Ф*(Е,а)=

J---

j dqv...dpf.

(34.3)

 

$f(<7;

Р)-а)<Е

 

m

Выясним теперь значения частных производных этой функции. Производную дФ*/дЕ (д* (Е, а) мы уже про­ анализировали в (31.10). Для расчета дФ*1да учтем, что выражение

^ 1 б а = Ф* (Е,а+ 8а) — Ф* (Е,а)

да

определяет ограниченный двумя гиперповерхностями

(I) Ж (q, р;а) = Е

и (II) Щч, р;а + 8а) = Е

объем фазового пространства

(рис. 64).

Поверхность (II) равнозначна с

Ж(Я,р-а)

=

Е-^-8а.

 

 

да

Если мы, как и ранее, обозначим через dO элемент поверхности Ж=Е и через 6s расстояние по нормали между двумя поверхностями, то будем иметь:

-^8a=[8sdO. (34.4)

да J

Величина 6s следует из условия, что 6s должно иметь направление grad Ж и что при удалении на расстояние

8эЖ также должна уменьшаться на дда$€ 6а. Следовательно,

 

Igrad^l

6 s = - ^ 6 a

 

 

 

 

да

и поэтому

с.

 

 

 

дФ*

с Г

дШб dO

да

оа = — 8а \

да

 

 

J

|GRAD^|

Но согласно

определению

(32.5) и выражениям

(32.3а) и (32.4) для микроканонического среднего спра­ ведливо

 

* дШ

dO

 

dfe

J, да

|gradgg|

 

да

со* (Е)

 

следовательно,

 

 

 

дФ*

* , с. . Ш

,

= — (о* (Е, а)

да

 

да

 

Тем самым дли любого изменения ( t) : , : (£) будет иметь

место:

 

с1Ф* = (О* dE-

da\

 

да

Это чрезвычайно важный

результат.

Из (34.1) мы знаем, что связанные с медленным из­ менением а затраты работы увеличивают энергию на ве-

Рис. 64. I и

I I представляют

собой гиперпо­

верхности Ж

(р, q\ а) —Е кЖ

(р, q\ а-\-Ьа) — Е.

личину (дЖ1да) da. Назовем такое воздействие на систе­ му «адиабатным» процессом.

Из (34.4) следует, что при адиабатном1 изменении параметра (или нескольких параметров) фазовый объем Ф* остается неизменным. Это и есть адиабатная инва­ риантность фазового объема.

в) k In Ф* как энтропия

Способ записи уравнения (34.4)

dE = — ЛФ* + ~ da

(34.5)

ш* да

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

!_аф* = 8<2.

(34.6)

со*

Оба способа увеличения энергии различаются весьма характерным образом. Во время совершения работы мы изменяем механическое состояние системы, харак­ теризуемое изменением а, при этом движение, описывае­ мое уравнениями (31.1), не нарушается. В противопо­ ложность этому при подводе тепла 6Q механическое сос-

1 Адиабатный здесь имеет значение бесконечно медленного, од­ новременно данный процесс является адиабатным в том смысле, что согласно (34.6) он не сопровождается переносом тепла.

185

тояние не изменяется; вместо этого мы произвольно воз­ мущаем характер хода движения, изменяя, например, значения отдельных импульсов pi не предусмотренным уравнениями (31.1) образом. Если наша система вообще

имеет

свойства нагретого

тела,

то уравнения

(34.4) и

(34.6)

заставляют нас рассматривать

величину /г1пФ*

Как энтропию. В самом деле, после умножения

на k/Ф*

из уравнения

(34.4)

(с учетом a* = d<t>*/dE)

 

следует:

 

 

 

 

 

 

 

 

 

d(k

1пФ*)

 

д (k In Ф + )

dE-

da

 

 

 

 

 

дЕ

 

 

 

 

 

 

 

 

да

 

в полном соответствии с выражением

 

 

 

 

dS=

\dE — 8А],

 

 

если согласно

(33.4)

представить

 

 

 

 

д ( £ 1 п Ф * ) _ 1

 

/ г ! п Ф *

+ S 0 .

(34.7)

 

 

дЕ

Т

И о :

 

 

 

 

 

 

 

 

Здесь S0 означает вначале неизвестную, независимую от Е и а величину. Таким образом, введенный уравнением (31.9) фазовый объем оказывается параметром, опреде­ ляющим термическое поведение нашей системы.

Поясним адиабатную инвариантность Ф* на очень простом при­ мере (рис. 65). Рассмотрим маятник на нити массой т, длиной / с координатой а и энергией

1

•£кин = — ml2 а 2 ; ЕП0Т =mgl cosa. m

Соответствующий импульс

будет

р

=dE^Jda

— tnl2a, следова­

тельно, ЕКШ2х/2т12.

Ограничимся

небольшим

отклонением

а, за­

меним

 

 

 

 

 

 

 

 

cos a « 1

 

 

 

 

в связи с чем имеем функцию Гамильтона:

 

 

Ща,

р; I) = •

Pi

mgl

• mgl = Е

(34.8)

 

2

 

 

2 т / 2

 

 

 

 

Фазовый объем Ф*( £ , /) =

^dadpa

представляет собой

эллипс

в плоскости а, ра

с полуосями

 

 

 

 

 

I Vim

V E+mgl

и Vi V E + 2mgl/

V mgl .

 

186

Следовательно,

Ф* ( £ , I) = 2я

- j - ( £ + m g / ) .

(34.9)

При небольшом изменении £ и / будет выполняться

/

/ . „

, £ + 3mg/

d/

(34.10)

dQ*

d £

2/

 

 

 

 

Рассмотрим далее длину нити / как параметр, который мы изме­ няем адиабатно. С этой целью, как показано на рисунке, нить введе-

Рис. 65. Маятник с длиной нити / в качестве «пара­ метра». Пример адиабатной инва­ риантности Ф*.

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

 

 

К=—

дЖ

Pa

 

mga,2

 

 

 

 

 

 

 

- г — = — — — — — + m g .

 

 

 

 

 

 

 

dl

ml3

 

 

2

 

 

 

 

 

Значение этих трех слагаемых легко уяснить. Первое представ­

ляет собой центробежную силу, два

других

равны

умноженной

на

cos а

силе тяжести mg.

 

 

 

 

 

 

 

 

 

 

 

 

Усреднение

по

микроканоническому

ансамблю

означает

здесь

усреднение во времени по одному колебанию

(при

постоянной

/ ) .

Для

гармонического

колебания

средние

значения

потенциальной

и кинетической энергий равны между собой

 

(потенциальную

энер­

гию в

состоянии

покоя

mgl

в данном

случае

учитывать не следует).

Тогда

из (34.8)

следует:

 

 

 

 

 

 

 

 

 

 

 

 

 

-

2

 

mgl

 

1

,

 

 

 

 

 

 

 

 

Pa

=

 

 

mgl).

 

 

 

 

 

2m/2

2

а 2

= —

£ +

 

 

 

 

 

 

 

 

2

 

 

ё

;

 

 

 

 

 

 

 

 

 

 

— о

 

— о

 

 

 

 

 

 

С этими средними значениями ра

и а

получаем:

 

 

 

 

 

К =—

дЖ

=

1

(E +

3mgl).

 

 

 

 

 

 

 

 

dl

 

21

 

 

ё

 

'

 

 

 

При адиабатном изменении / работа, совершенная над системой, равна —Kdl (при укорочении нити dl отрицательно). Таким образом,

d £ = ~ K d / = - ^ £ ' d / .

187

Следовательно, согласно (34.10) при адиабатном изменении / действительно ЙФ* = 0.

Применительно к квантовой теории внесем некоторые дополне­

ния к

(34.9). Если мы назовем E-\-mgl

= Екол

энергией колеоэния

нашего

1

.

представляет собой ча-

маятника и учтем, 4 T O V = — —

| g/l

 

 

 

 

стоту его колебаний, то из уравнения (34.9)

следует:

v

т. е. отношение энергии к частоте при адиабатном изменении остает­ ся постоянным.

35.

Д Е Л Е Н И Е НА N1 И Д Е Й С Т В И Т Е Л Ь Н Ы Й Ф А З О В Ы Й О Б Ъ Е М

Ф

(Я, V, N)

а) Расчет Ф * для идеального газа

Начнем с расчета Ф* для газа, состоящего из N молекул и заключенного в объеме V. Функция Гамильтона для него имеет вид:

3N

2

/ I ~7Г Г- ^ с т е н к и -

2m

i=l

В интеграле

Ф* = f... \dxv..dx3Ndpr..dp3N

интегрирование для любой тройки координат dxjdtjjdzj дает значение V. В целом, следовательно, появляется множитель VN. Интеграл по импульсам распространяет­ ся на ту часть ЗЛ^-мерного пространства импульсов, для которой

ЗЛ/

£р2<2шЕ.

/-==1

Это будет ЗМ-мерный

шар с

радиусом

V2mE. Его

объем равен A3N(2mE)?N^.Зяесъ

A3N — объем

ЗМ-мерно-

го единичного шара. Следовательно, используя

приведен­

ное в § 35 значение AZN,

имеем

 

 

 

 

ЗЛ/

3/V

 

 

Ф* = A3N (2mE) 2 VN =

e"j 2

VN

(35.1)

188

&1пФ*

— klnE + NklnV

+

f

3Nk

ln(2m)+ kin A3

(35.2)

 

2

 

 

От E и К зависят только два первых слагаемых. Ес­ ли мы предварительно заменим S = £lrU'D*, то для одно­ атомного газа будем иметь уже известный результат

dS

\дЕ

Т

dV

V

Т

Следовательно, в отношении зависимости энтропии от

Е и V значение

энтропии

в виде

klnO*

могло бы нас

удовлетворить.

В то же время это определение оказывает­

ся неприемлемым, если мы захотим также

воспроизвести

правильную зависимость 5

от N.

Зависимость энтропии

от числа частиц мы уже обсуждали в § 10, руководству­ ясь тем фактом, что энтропия газа, заполняющего объ­ ем V, не изменяется, если с помощью ввода разделитель­ ной стенки разделить этот объем на два других, напри­ мер V\ и V2. Этот факт обусловил то требование, что в

выражении (10.3) 5 (Т, V)=N(cvlnT+klnV+C)

ве­

личина С должна зависеть от Л' в виде

С——klnN-\-e,

где а теперь более не зависит от N.

 

Рассмотрим теперь изменение, которое претерпевает фазовый объем Ф* после удаления упомянутой стенки.

Пусть ранее в объеме в Vi находилось

Ni определенных

молекул, а в V2 N2

молекул. После

удаления

стенки

какая-либо из N\ молекул поменялась

местами

с одной

из N2 молекул. Каждый раз это дает новую точку в Г-про-

странстве. Так как в

целом имеется

Nl/N\IN2\

различ­

ных возможностей выбрать из N атомов Ni, то величина Ф* увеличивается точно в такое же число раз. Это неже­ лательное увеличение вызывается тем, что мы получаем в Г-пространстве новую точку в каждом случае, когда мы меняем координаты и импульсы каких-либо двух оди­ наковых Частиц. Следовательно, если мы хотим избе­ жать раЗдуйания фазового объема при устранении раз­ делительной стенки, то нужно договориться приписывать одному и тому же состоянию все то точки, которые обус­ ловлены взаимной переменой места двух одинаковых ато­ мов. Но за счет произвольной перестановки N атомов

189

Соседние файлы в папке книги из ГПНТБ