Файл: Смирнова Н.А. Методы статистической термодинамики в физической химии учеб. пособие.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) |
2т |
2т |
2т |
|
Величина |
s зависит только от модуля импульса Р = \р\ — |
= ѴрІ + |
РІ + PÎ ; |
2т |
(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 |
|
|
|
|
|
|
|||
dï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т |
2т |
2т j |
2т |
к=\ |
|
|
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]