книги из ГПНТБ / Каули, Ч. Теория звездных спектров
.pdf90 |
Г ЛАВА 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 (gn— Nn) Я 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
