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

книги из ГПНТБ / Каули, Ч. Теория звездных спектров

.pdf
Скачиваний:
21
Добавлен:
21.10.2023
Размер:
9.13 Mб
Скачать

90

Г ЛАВА 3

и т. д. Ясно, что вероятность такой конфигурации пропорцио­ нальна числу способов, которыми может быть распределено Nx частиц с энергией Еі, N2 частиц с энергией Е2 и т. д. Мы огра­ ничены в этом выборе «краевыми» условиями:

а) общее число частиц должно быть N, т. е.

=

(3.2.1)

п

б) общая энергия системы должна быть постоянной:

=

(3.2.2)

П

 

Может существовать ряд различных элементарных состоя­ ний с энергией Еп. На языке квантовой механики энергетиче­ ское состояние Еп является ^„-кратно вырожденным, где gn — число элементарных состояний с энергией Еп.

Продолжим вывод статистики Ферми — Дирака или Бозе — Эйнштейна (в зависимости от того, проигнорируем мы или учтем принцип Паули). Здесь мы изложим вывод статистики Ферми — Дирака, а для статистики Бозе — Эйнштейна приведем лишь конечный результат, оставляя подробный вывод читателю.

Согласно принципу Паули, максимальное число частиц, ко­ торые могут занять п-е энергетическое состояние, есть gn, так что всегда Nn ^ gn. Рассмотрим теперь п-е энергетическое со­ стояние и спросим, сколькими способами Nn частиц можно рас­ ставить по gn состояниям. Первую частицу можно поместить в любое из gn состояний, вторую — в любое из gn — 1 оставшихся состояний и т. д. (Порядок, в котором помещаются частицы по состояниям, не важен.) Таким образом, мы насчитаем

g n ( g n - l ) ••• {gn — Nn+ \ ) = gn\l{gn — Nn)\

(3.2.3)

способов. В каждом наборе из Nn частиц можно произвести Nn\ перестановок, а в силу принципа неразличимости одинако­ вых частиц нельзя считать эти перестановки различными реа­ лизациями, поэтому общее число реализаций для п-го энерге­ тического состояния теперь равно

grMNn\ (gn Nn)\.

(3.2.4)

Для каждого возможного распределения частиц по т-му энерге­ тическому состоянию получим

gmW m\{gm- N m)\

(3.2.5)

различных реализаций, а общее число способов, которыми мож' но осуществить распределение по энергиям, составит

g n\

(3.2.6)

*-пп Nn\(gn -N ^\'

 

СТАТИСТИЧЕСКАЯ МЕХАНИКА

91

Это произведение берется по всем возможным энергетическим состояниям. Величина W называется термодинамической веро­ ятностью. Поскольку W ^ 1, это не истинная вероятность, но ясно, что эта величина пропорциональна вероятности.

Найдем теперь условный максимум ln W при выполнении двух краевых условий (3.2.1) и (3.2.2). Ясно, что W имеет ма­ ксимум, когда ln W максимален. Процесс максимизации ln W при двух краевых условиях проводится методом лагранжевых

множителей [90, р.

128]. Вводя неопределенные множители Я и р,

запишем

 

 

 

 

6 ln W + Я26УѴ„+ р 2£„6УѴ„ =

0.

(3.2.7)

 

п

п

 

 

Если использовать

формулу

Стирлинга для

In Х\

при боль­

ших X:

In Я! =

Я In X — X,

 

(3.2.8)

 

 

то легко получить

 

 

 

 

2 [— ln Nn4- ln (gnNn) Я p £ rt] бNn= 0.

(3.2.9)

Вариации Nn считаются произвольными, т. е. бУѴі может рав­ няться нулю, а бМ2 может быть конечной величиной или на­ оборот и т. д. Чтобы сумма всегда обращалась в нуль, необхо­ димо равенство нулю при всех п выражения в квадратных скоб­ ках. Следовательно, можно заключить, что

Xn/(gn — М„) = ехр(Я + р £ л),

(3.2.10)

или

________ Sn______ _

 

Nn

(3.2.11)

exp (— Я — цЕп) + 1

 

 

Для статистики Бозе — Эйнштейна выражениями, аналогичными

(3.2.6) и (3.2.11), являются:

хѵ/ ТТ (Nn + gn 1)!

1L Aalten-1)!

п

ДГ gn

пexp (—Я (х£„) 1 '

(3.2.12)

(3.2.13)

Промежуточные выкладки приводятся в [151, § 87, 88] и в боль­ шинстве книг по статистической механике.

Теперь нужно определить постоянные Я и р. Мы сделаем это для статистики Больцмана, в которую переходят при gn >> Nn и статистика Ферми — Дирака, и статистика Бозе — Эйнштейна.

92 ГЛАВА 3

Производя упрощение (3.2.10) и используя формулу (3.2.1),

получим

V

ехр (Я)

(3.2.14)

2

Sn ехр (\іЕп)

П

 

Сумма в уравнении (3.2.14) называется суммой по состояниям.

Вобщем случае обозначим ее через Q.

Вследующем разделе мы покажем, что множитель р опре­ деляется соотношением

ц = — 1 ЦкТ),

(3.2.15)

где k — постоянная Больцмана. Окончательно формула Больц­ мана приобретает вид

 

-jf- = Щ- ехр (— EJkT).

(3.2.16)

В астрономии

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

Ѳ= 5040/7".

В этом случае,

если энергия Еп выражается в электрон-вольтах

и обозначается символом %п, то

 

 

=

(3.2.17)

3.3. ОПРЕДЕЛЕНИЕ ПОСТОЯННОЙ р

Самый простой вывод постоянной р в уравнении (3.2.14) делается на основе принципа Больцмана, который выражает связь между энтропией S системы и термодинамической вероят­ ностью W. В соответствии с вышеизложенной точкой зрения мы будем говорить, что все физические системы стремятся к наибо­ лее вероятному состоянию. Формулу Больцмана (3.2.16) мы вы­ водили, имея в виду это предположение.

Но, согласно термодинамике, все физические системы стре­ мятся к состоянию с максимальной энтропией. Известно, что энтропия системы обладает экстенсивным свойством, т. е. за­ висит от количества вещества в системе. Энтропия системы, составленной из двух частей с энтропиями Si и S2, равна их сумме

S = S, + S2.

(3.3.1)

С другой стороны, из теории вероятности известно, что если W1 есть термодинамическая вероятность одной части системы, а W2— другой, то термодинамическая вероятность W всей систе­ мы выражается произведением

W = W lW2.

(3.3.2)

СТАТИ СТИЧЕСКАЯ МЕХАНИКА

93

Сравнение формул (3.3.1) и (3.3.2) наводит на мысль, что энтропия и термодинамическая вероятность должны быть свя­ заны логарифмически, скажем

S = £ l n r ,

(3.3.3)

где k — постоянная. Уравнение (3 3.3)

выражает принцип

Больцмана, а постоянная k есть постоянная Больцмана *). Согласно термодинамике, введение элементарного количе­

ства тепла dq в систему приводит к росту энтропии по закону **)

dS = dq/T.

(3.3.4)

В соответствии с первым законом термодинамики это тепло dq может повысить внутреннюю энергию на dE, а также совер­ шить работу PdV, где Р — давление, а V— объем. Таким об­ разом,

dS = Y dE + T dV-

(3.3.5)

Из того факта, чтоэнтропия является полным дифференциа­ лом, и из (3.3.5) следует

[dS/dE)v = \ / T .

(3.3.6)

Теперь оценим эту частную производную исходя из вели­ чины S, полученной из принципа Больцмана. Если воспользо­ ваться приближением Стирлинга при оценке натуральных лога­ рифмов всех множителей уравнения (3.2.6), то можно получить

1п1Г = - 2

( Ь 1 п 1 ^ - +

^ 1 п

(3.3.7)

 

п

 

 

Введем упрощающее допущение: gn

Nn. Первое слагаемое в

правой части (3.3.7)

обратится в нуль,

если мы пренебрежем Nn

вчислителе выражения под знаком логарифма. Но мы запишем аргумент логарифмической функции в виде 1 — {NJën) и исполь­ зуем формулу 1п(1 — х) X —л: при малых х. Таким образом, первое слагаемое в (3.3.7) дает полное число частиц N, если просуммировать по п.

Представляем читателю убедиться, что при отбрасывании Nn

взнаменатель второго слагаемого (3.3.7) вносится относительно малая ошибка. Если при этом используется соотношение

== ехр (ц£п),

(3.3.8)

*) Согласно Зоммерфельду [137], как уравнение (3.3.3), гак и постоян­ ная Больцмана k были введены Планком.

**) См. ниже разд. 3.4. Чтобы выполнялось (3.3.4), количество тепла dq должно вводиться в систему медленно (или обратимо).

94

ГЛАВА 3

 

ТО

 

 

 

S =

k [N (l -fln -g-) — ц£ ] .

(3.3.9)

Таким образом,

 

 

(3.3.10)

(dS/dE)v =

— i ik = 1IT,

отсюда

ц =

— 1 IkT.

(3.3.11)

 

Этот метод вывода постоянной р, не очень удобен для тех, кто не знаком с принципом Больцмана или с понятием энтро­ пии. Поэтому рассмотрим другой способ получения соотношения (3.3.11) из уравнения состояния идеального газа

 

PV = nNakT = nRT,

(3.3.12)

где V — объем газа,

Na — число Авогадро, равное

6,024-ІО23,

а п — число молей в

объеме V. Сначала получим

максвеллов­

ское распределение скоростей из формулы Больцмана.

Левую часть выражения (3.2.16) можно интерпретировать как вероятность того, что частица будет иметь энергию Еп. Пе­ рейдем от этого соотношения, в котором постулируются дискрет­ ные энергетические состояния, к схеме, в которой возможна непрерывная область энергий, как в классическом случае моле­ кулярного газа. Если выразить энергию через импульсы, то по­ лучим вероятность, определенную как функция импульсов, а не энергии. Переход от дискретного к непрерывному распреде­ лению осуществляется по следующей схеме:

N n ^ N

(р х , р у , р г ) d p x d p y d p z

(а),

 

d p x d p y d p z T7

 

V d p

gn->

h 3

V,

или

h 3

(6),

(3.3.13)

И exp[~ {P l + Pl + Pl)l2 mkT]n jr- (e ) .

Здесь V— выделенный объем, в котором содержится N частиц; через d p обозначено произведение d p x d p y d p z , а не вектор. Мы использовали концепцию ячейки фазового пространства объема h3 (разд. 3.1).

Выражение для суммы по состояниям легко проинтегриро­ вать. Поскольку импульсы по направлениям х , у и г взаимоне­ зависимы, имеем произведение трех интегралов ошибок (см., например, [105, § 13.29]). Следовательно,

Q __ ( 2 n m k T ) ^ у

(3.3.14)

**

h 3

 

СТАТИСТИЧЕСКАЯ МЕХАНИКА

95

В разд. 3.4 мы подробно выведем это соотношение для частицы, заключенной в ящик.

Величина gn была взята пропорциональной объему, посколь­ ку Nn— число частиц в этом объеме. Ясно, что когда V растет, должно пропорционально увеличиваться число возможных ячеек объема /г3.

Вероятность Р(рх, ру, pz)dp того, что импульсы примут зна­

чения от рх, Ру, рг до рх + dpx, ру +

dpy, pz + dpz, дается соотно­

шением N(px, ру, pz)dp/N — NJN, т.

е.

 

 

Р(рх, Ру, pz)dp = (2nmkT) 3/2exp

Рж +

Р2у + р\

dp. (3.3.15) *)

 

 

2mkT

 

Величина V/h3 одинаково входила в суммы по состояниям Q

и статистические веса g n и потому

в (3.3.15) сократилась.

Средний квадрат импульса, скажем, в направлении х полу­

чается умножением выражения (3.3.15)

на р2х

и интегрирова­

нием по всем импульсам. Вследствие независимости направле­ ний X, у и z мы снова получаем три интеграла. Интегралы по у и z по-прежнему равны интегралу ошибок и легко отыски­ ваются. Таким образом,

 

РІ = (2nmkT)~'u f РІ ыр { — 2ш ) йРх-

(3.3.16)

Оставшийся интеграл легко берется,

так что получаем

 

 

р| == mkT.

(3.3.17)

Если

учесть изотропию,

полагая,

что полная

скорость

ѵ2=

Зѵ2х, то

_

 

 

 

 

тгР/З = kT.

(3.3.18)

Умножая обе части соотношения на число частиц N, равное про­ изведению числа молей п на число частиц в моле Na, найдем

mNv2!3 — nRT = PV.

(3.3.19)

Уравнение (3.3.19) следует из элементарной кинетической теории газов (см., например, [133]). Этот результат согласуется с определением р как — \/ІгТ = Na/RT и, следовательно, под­ тверждает определение р, даваемое соотношением (3.3.11).

*) Это выражение получается подстановкой (3.3.136) и (3.3.14) в (3.3.8) с учетом того, что рЕп = — (р \ + р2 + pf) j2mkT. Прим, перев.

96

ГЛАВА 3

3.4.

ПОНЯТИЕ РАВНОВЕСИЯ

Прежде чем выводить уравнения ионизации и диссоциации, нужно кратко рассмотреть вопрос о химическом равновесии с точки зрения термодинамики. Хотя читатель, по-видимому, не­ много знаком с понятиями энтропии и свободной энергии, мы сначала кратко напомним эти понятия. Мы будем пользоваться обозначениями, принятыми в химии, чтобы облегчить обраще­ ние к таким книгам, как учебник физической химии [58].

Мы уже говорили об энтропии системы с точки зрения при­ ближения системы к ее наиболее вероятному состоянию. Столь же важным является понятие энтропии системы как меры ко­ личества энергии, которая не способна производить работу. Это понятие наиболее наглядно иллюстрируется при изучении цикла Карно, на примере которого показано, что к. п. д. тепловой машины всегда меньше 100%.

Рассмотрение цикла Карно приводит к определению свобод­ ной энергии системы Л, которая равна полной энергии, способ­

ной совершить работу:

 

А = Е — TS.

(3.4.1)

Даже эту энергию, называемую свободной энергией Гельм­ гольца, нельзя полностью превратить в полезную работу, по­ скольку часть энергии всегда идет на изотермическое сжатие.

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

К = Л + РК

(3.4.2)

и рассмотрим газ, который расширяется, передвигая поршень. Пусть давление Р постоянно, а изменение объема равно АѴ =

— Ѵ2— Vi, тогда

AF = AA + PdV.

(3.4.3)

Мы видим, что Р dV положительно. Поскольку газ производит работу, свободная энергия Гельмгольца уменьшается и АЛ от­ рицательно. Таким образом, левая часть соотношения (3.4.3) представляет собой разность между изменением свободной энер­ гии Гельмгольца АЛ и работой, затраченной на расширение. Эта разность АF может расходоваться не только на изотерми­ ческое сжатие. Из (3.4.1) и (3.4.2) имеем

К = Я — TS,

(3.4.4)

где Н — Е + РѴ называют

теплосодержанием,

или энтальпией,

газа.

 

 

Приращение энтропии системы dS можно записать [ср.

уравнение (3.3.5)] в форме

 

 

dS =

(dE + PdV)/T,

(3.4.5)

СТАТИСТИЧЕСКАЯ МЕХАНИКА

97

где числитель соответствует приращению тепла, которое ме­ няет внутреннюю энергию dE, и производит работу против дав­ ления Р. Для реакций, которые происходят при постоянных дав­ лении и температуре, из (3.4.4) и (3.4.5) следует

(dF)T'P = dE + PdV — TdS = 0.

(3.4.6)

Поскольку обратимые изменения достигаются посредством смены ряда равновесных состояний, можно считать уравнение (3.4.6) уравнением, определяющим равновесные процессы. Рас­ смотрим, например, равновесный процесс

А + В ^ С ,

(3.4.7)

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

нент А, В и С являются Fa, Fb и F c ,

т о условие

равновесия

(3.4.7) есть

 

 

й{РА + Рв -Рс)т.Р =

Ъ-

(3-4.8)

Уравнение (3.4.8) представляет условие минимума свобод­ ной энергии Гиббса при постоянных температуре и давлении. Это только одна из трех альтернатив, которые можно выбрать при выражении условия равновесия в зависимости от того, ка­ кие переменные фиксированы. Другими возможностями являет­ ся максимизация S при фиксированных V и Е или минимиза­ ция А при фиксированных V и Т.

Таким образом, наше исследование химического равновесия приводит к вычислению свободных энергий FA, FB и Fc-

Сначала найдем выражение для энтропии (отметим, что нам нужна именно абсолютная энтропия). Имеется несколько спо­ собов, которыми можно это выполнить. Здесь мы воспользуемся принципом Больцмана [уравнение (3.3.3)] и выражением (3.2.6) для термодинамической вероятности. Возможны и другие ме­ тоды, например в [58, р. 871] всюду, за исключением абсолют­ ного нуля, применяется статистика Больцмана и выводится выражение для энтропии интегрированием уравнения (3.3.4). Читатель может убедиться сам, что окончательный результат

будет одним и тем же, если для W вместо формулы

Ферми —

Дирака взять формулу Бозе — Эйнштейна (3.2.12).

 

Из уравнения (3.3.9)

при р = — \/кТ имеем

 

lnW =

N + Nln(Q/N) + (E/kT),

(3.4.9)

где полная энергия Е дается соотношением (3.2.2).

 

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

 

F = £ + Py — Ш пІГ .

(3.4. ІО)

4 Ч . К а у л и

98

ГЛАВА 3

 

Применяя уравнение состояния идеального газа

(3.3.12), упро­

стим это выражение:

 

 

F = — nRT ln (Q/N),

(3.4.11)

где n =

NINa— число молей.

 

Вернемся к уравнению (3.4.8). Дифференциал выражения в

круглых

скобках отвечает случаю, когда малое

(по сравнению

с полным числом) число молей атомов А реагирует с малым числом молей атомов В с образованием того же самого числа молей молекул С. Поэтому в выражениях свободных энергий коэффициенты nRT при логарифмических членах равны для А, В и С и сокращаются в (3.4.8). Тогда мы можем приравнять аргумент алгебраической суммы логарифмов единице и получим

NaNв/^с — QaQb/Qc-

(3.4.12)

Частицы А, В и С должны свободно перемещаться, следова­ тельно, нужно считать, что сумма по состояниям

Q = h g n e x p ( - E J k T )

(3.4.13)

П

 

содержит и внутреннюю энергию, и кинетическую энергию поступательного движения, т. е. Еп — (£п)вн+ (Дп)пост. Поэтому правая часть уравнения (3.4.13) записывается в виде произве­ дения

Q = QBItQnocT.

(3.4.14)

Теперь выведем выражение для суммы по состояниям посту­ пательного движения квантовомеханической частицы, помещен­ ной в прямоугольный ящик со сторонами a, b и с. Забегая впе­ ред, скажем, что мы придем к результату, который согласуется с уравнением (3.4.14), полученным в предположении, что объем ячейки фазового пространства равен /і3, и такое совпадение убе­ дит нас в пригодности постулата об объеме фазовой ячейки.

Волновое уравнение для частицы в прямоугольном ящике записывается в виде

 

 

h2 ( дгЧ>

I

дгФ

,

дгЧ>\

р

 

(3.4.15)

 

 

8л2т \ дх2

^

ду2

^

dz2 ) ~

w’

 

 

 

 

 

где

Е — энергия

частицы.

Подстановкой

ф(х,

у, z) =

=

X(х) Y(у)Z(z)

можно разделить

переменные

этого

уравне­

ния. А если затем принять, что Ех есть кинетическая энергия

поступательного

движения в

направлении х, Еѵ— в направле­

нии у, а Ег — в

направлении

z, то уравнение (3.4.15) распа­

дается на два уравнения, одно из которых содержит только переменную х, а другое — только у и г. Рассматривая лишь дви­ жение в направлении х, получим уравнение

— (fi2/8n2m) (Х"ІХ) = Ех,

(3.4.16)

СТАТИСТИЧЕСКАЯ МЕХАНИКА

99

решение которого

 

Х = А sin(Y8я2тЕх//і2х) + В.

(3.4.17)

Требование, чтобы X имело узловые точки при х — 0 и х — а, приводит к равенству нулю постоянной интегрирования В и квантовомеханическому условию

Ех = tiji2/8а2т, пх = 1, 2, 3, . .. .

(3.4.18)

Решения для у и z аналогичны. Сумма посостояниям для ки­ нетической энергии поступательного движения Quoct есть

Q n O C T

/ 2

2

2

 

/ nY

п„

пу

(3.4.19)

V ехр [— \ о 2

^ b2

^ с2) h2/8mkT

_ J L

J _______L

j _______ 1

 

Поскольку расстояния между энергетическими уровнями

малы, перейдем от суммы (3.4.19)

к интегралу

оо

2

9

2

 

QnoCT === ■§" J СJ SXp

Пх

I пу

, пг

hV8mkT dnx dn,, dtiy

а2

"Г й2

“Г С2

(3.4.20)

где множитель Vs появился вследствие того, что пределы п взяты от —оо до оо, а не от 0 до оо. После интегрирования окончательно получаем

Qnocx = (2jmikT)sl1 V/h3,

(3.4.2 1)

где V = abc — объем ящика.

Рассмотрим произведение (l/s)dnxdny dnz. При помощи урав­ нения (3.4.18), подобных уравнений для у, z и выражений типа

Ех = рЦ2т получим соотношение

 

dnx dny dnj8 = dpx dpy dpzV/h3,

(3.4.22)

которое находится в согласии с постулатом (3.3.136).

Прежде чем вернуться к уравнению (3.4.12), сделаем не­ сколько замечаний относительно сумм по состояниям для каждой из частиц А, В и С. Начало отсчета энергии можно выбрать до некоторой степени произвольно, но нужно убедиться в том, что нуль-пункты сумм по состояниям атомов А, В и молекулы С са­ мосогласованны. Если мы примем за нуль-пункт атомов А и В основные состояния, то основное состояние молекулы будет на D (эВ) ниже этого нуль-пункта. Здесь D — энергия диссоциации молекулы. Обычно сумма по состояниям для молекулы отсчиты­

вается от основного состояния вверх.

Если Е°п— энергетиче­

ское состояние молекулы, отнесенное

к основным состояниям

4

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