Файл: Смирнова Н.А. Методы статистической термодинамики в физической химии учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 183

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

где e(i) — энергия i-ü частицы. Одноатомный идеальный газ представ­ ляем как систему N практически невзаимодействующих материальных точек, движущихся в некотором объеме V. Будем полагать, что внеш­ нее поле отсутствует и, следовательно, энергия газа есть кинетическая энергия поступательного движения частиц.

В j A - n р о с т р а н с т в е , имеющем шесть измерений, осями являются X, у, z, рх, ру, рг. Можно выделить подпространство координат с осями ж, у, z и подпространство импульсов с осями рх, ру, рг (рис. 7).

Элемент

фазового объема есть

 

 

 

 

 

dr\ = dx dy

dz dpx dpy

dpz

= d^v d-\pl

 

где

d^v

= dxdydz = dV — элемент

в

подпространстве

координат;

dyp

— dpxdpydpz — элемент

объема

в

подпространстве

импульсов.

Координаты частиц могут принимать любые значения в пределах объе­

ма V. Энергия частицы е есть кинетическая энергия

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

движения

 

 

 

2

2

2

 

- Р х

+

А . .

(11.63)

 

Величина

s зависит только от модуля импульса Р = \р\ —

= ѴрІ +

РІ + PÎ ;

(11.64)

 

Определим, какой объем в фазовом пространстве отвечает состоя­ ниям с энергией частицы в интервале е, е - f de. Этот объем назовем объемом энергетического слоя и обозначим di(s). Величину dy(e) можно найти, взяв интеграл от d-; по всем состояниям, совместимым с данной энергией и условием, что частица находится в объеме V:

* Т ( 0 =

j

d T = = j j J

а х а У а г

j

dtp = v

j

dy.

(s,

e + rfs, V)

(V)

(«,

E+rfe)

(e, s +

de)

 

40


Так как согласно (11.64) всем состояниям с заданным значением модуля импульса р отвечает одна и та же энергия, то в подпростран­ стве импульсов состояния с энергией е изображаются точками, лежа­ щими на сфере радиуса р = Уітг. Сфера в подпространстве ур является поверхностью постоянной энергии. Состояниям с энергией от е до s + ds в подпространстве ур отвечает сферический слой радиуса р и толщины dp; объем этого слоя

dtP(p)=*A*pdp.

(11.65)

Указанная величина dyp(p) является результатом интегрирования

элементарного объема

 

dfp = dpx dpy dpz = p2 sin Ѳ dp d9 dtp

(11.66)

(p, Ѳ, ф —сферические координаты в подпространстве импульсов) по

всем значениям Ѳ (от 0 до it) и ф (от 0 до 2іг).

Таким

образом, фа­

зовый объем, отвечающий значениям модуля

импульса

частицы от

р до р + dp, есть

 

 

df (р) = V 4г.р2 dp.

 

(11.67)

Чтобы получить эту величину как функцию энергии, сделаем замену

переменных,

учитывая равенство

(11.64):

 

 

J .

 

 

1

JL _ - L

 

 

p = (2me)2 ;

dp =

y ( 2 m ) 2 e 2 dz.

( I L 6 8 )

В силу соотношений (11.67)

и (11.68)

 

 

 

!

J

 

!_

- L

 

d 7

(в) = V 4* 2m е — (2т)2

s

2

de = 4пт V (2ms)2 de.

(11.69)

Формула (11.69) определяет объем энергетического слоя в р.-простран- стве как функцию энергии частицы е. Введем понятие энергетической плотности состояний

dT (e)

g ( s ) = = d T 1

( I L 7 0 )

g(e) — это объем энергетического слоя при заданных е и de, от­ несенный к единичному интервалу изменения энергии. Для рассмат­ риваемого случая

j

_

 

g (s) = Ar.rnV (2ms)2

.

(11.71)

Фазовый

объем, в

котором находятся

изображающие

точки частиц

с энергией, равной или меньше заданной (0 ^ Я < ; е), обозначим Y(S).

Величину

его находим интегрированием выражения

(11.69):

 

«

_і_ в j _

_i_ JL

 

7

(е) =5 Г g(e)ds = 4ктѴ (2m)2 Г е2

de = - | _ rcmV (2m)2

s2 . (11.72)

 

о

о

 

 

В подпространстве

импульсов изображающие точки частиц с энерги-

41


ей, равной или меньше заданной

е = р2І2т, лежат

внутри сферы

радиуса р, так что

 

 

7 P ( P ) = - J * P 3 и

T ( P ) = V - f */>»•

( п - 73)

После замены р на s согласно соотношениям (11.68) получаем равен­

ство (11.72).

 

 

 

 

 

содержащего N

 

Г-пространство

одноатомного

газа,

частиц,

6ІѴ-мерное. Элемент

объема

пространства

есть

 

 

 

з.ѵ

 

 

 

 

 

dT =

П

dpi

dqt

=dTpdTv,

 

(H.74)

 

 

i'=i

 

 

 

 

где

3N

 

N

 

 

N

 

 

 

 

 

 

v

= ïldqt

 

= П

dxK dyK dzK = П dVK

(11.75)

 

/=1

 

K = l

 

K=\

 

(к — номер частицы, i — номер обобщенной координаты или импуль­ са);

ЗЛГ

N

 

dVp = П

dpt = П dpXK dpyK dpZK.

(11.76)

{»1

*—1

 

Кинетическая энергия газа равна

я = £N (А+А+А\=У;-3N?1. (11.77)

\ 2т

j

к=\

 

 

1=1

Учитывая, что для системы величины N и V фиксированы, опре­ делим фазовый объем, в котором находятся изображающие точки системы с энергией, равной или меньше Е:

 

 

Т(Е)=

j

dT =

J j

dTvdYp=TvTp(E)

.

(11.78)

 

 

 

{0<H<E,

V)

{0<H<E,

V)

 

 

 

 

 

Величина

ГѴ есть

результат

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

по координатам:

 

 

 

 

Г ѵ

= j

j . . .

f dVi...

dVN

= V";

 

(11.79)

Tp (Е) — результат

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

по

импульсам

при

условии

О ^

H <^ Е.

Оценить

эту величину

можно

следующим

способом.

Уравнение (11.77) есть уравнение ЗЛ/-мерной сферы

радиуса " | / 2 т £

(уравнение n-мерной сферы радиуса г

в общем виде записывается

как

2 х /

= г 2

) .

В подпространстве

импульсов

изображающие

точки

для

і—І

 

 

 

 

 

 

 

 

 

 

 

 

всех состояний с энергией меньше или равной Е лежат внутри сферы радиуса V2mEt так что величина Тр(Е) равна объему этой сферы. Воспользуемся асимптотической формулой для объема п-мер-

42


ного шара радиуса г при п > 1 (см. [5], стр. 50):

 

 

 

п

 

2кег2

 

 

 

lnVn--jln——,

 

 

(11-80)

где Ѵп — объем

я-мерного

шара;

е — основание

натуральных лога­

рифмов.

 

 

 

 

 

 

 

Положив п =

3/Ѵ и г = У 2тЕ,

найдем

 

 

 

. „

 

ЗЛ^ ,

2ue2m £

 

 

1 п Г

р ( £ ) ~ — !п-

 

 

 

 

 

2

 

З.Ѵ

 

 

гак что

 

 

 

 

 

 

 

 

 

 

 

 

ЗуѴ

3N

 

 

Г , ( £ ) Ц і і ! 2 * - ) а £

2 .

(11.81)

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

 

 

 

 

 

 

 

 

 

 

 

З.Ѵ

3N

 

 

Г (£) = ^ 4

^ е

j 2

£

2 .

(11.82)

Объем энергетического слоя в фазовом пространстве есть дифферен­ циал от величины Т{Е):

З.Ѵ

3N

 

 

dT (£) = -у- ( J ^ L J "*

£ " 2

' d £ .

(Il.83)

откуда найдем энергетическую плотность состояний

З Ѵ

3N

 

g (£) = *Ш « H . ( i g L ) - ^

£ - - 1 .

(II .84)

Формулы (11.82) и (11.84) показывают, что фазовый объем Г(£) при большом числе частиц N является чрезвычайно быстро возрастающей

зѵ

функцией энергии системы, поскольку Т(Е) ~ Е 2 .


III. ОСНОВНЫЕ ПОЛОЖЕНИЯ КЛАССИЧЕСКОЙ СТАТИСТИЧЕСКОЙ ТЕРМОДИНАМИКИ. МИКРОКАНОНИЧЕСКОЕ

И КАНОНИЧЕСКОЕ РАСПРЕДЕЛЕНИЯ

§ 1. Метод ансамблей Гиббса

Задача, стоящая перед молекулярной теорией макроскопических процессов, состоит в объяснении на базе законов микроявлений наблю­ даемого на опыте поведения системы. Таким образом, объектом изу­ чения должно быть изменение состояния данной системы во времени, что связано, если описывать движение частиц законами классической механики, с изучением фазовой траектории системы. Пусть М(р, q) — некоторая однозначная функция обобщенных импульсов и координат. Среднее значение ее за время опыта і может быть рассчитано соглас­ но формуле

 

 

( П і . і )

где dt(p,

q) — время, в течениеокоторого фазовая точка системы

нахо­

дилась в

элементе объема dpdq около точки с координатами

р и с .

Величина. Мх, измеряемая на опыте, является, следовательно,

сред­

ней по фазовой траектории, которую изображающая точка системы описала за время т.

Чтобы по формуле (III.1) рассчитать среднее М т , надо решить механическую задачу о движении системы, т. е. определить фазовую траекторию системы при заданных начальных условиях. Как уже отмечалось во введении, подобный путь решения в применении к макроскопической системе наталкивается на огромные практические трудности, хотя развитие вычислительной техники открывает здесь широкие перспективы*. Путем решения уравнений движения (если практически такое решение доступно) можно получить наиболее полные, в рамках классической теории, сведения о поведении кон­ кретной рассматриваемой системы. Однако чисто механический под­ ход имеет ограничения принципиального характера, о которых гово­ рилось ранее, и не достаточен для анализа общих закономерностей наблюдаемого на опыте поведения макроскопических систем (тер­ модинамических закономерностей). Такие фундаментальные термо­ динамические параметры, как температура, энтропия, химический

потенциал, не являются средними значениями механических

величин

и по формуле ( I I I . 1)

рассчитать эти

параметры нельзя (в

формуле

( I I I . I) интересующие

нас параметры

просто отсутствуют).

 

Чтобы вскрыть смысл термодинамических параметров «немеха­ нического характера», необходимо перейти к вероятностному описа-

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

44]