химический каталог




СТАТИСТИЧЕСКАЯ ТЕРМОДИНАМИКА

Автор Химическая энциклопедия г.р. Н.С.Зефиров

СТАТИСТИЧЕСКАЯ ТЕРМОДИНАМИКА, раздел стати-стич. физики, посвященный обоснованию законов термодинамики на основе законов взаимодействие и движения составляющих систему частиц. Для систем в равновесном состоянии СТАТИСТИЧЕСКАЯ ТЕРМОДИНАМИКА т. позволяет вычислять термодинамические потенциалы, записывать уравнения состояния, условия фазовых и химический равновесий. Неравновесная СТАТИСТИЧЕСКАЯ ТЕРМОДИНАМИКА т. дает обоснование соотношений термодинамики необратимых процессов (уравений переноса энергии, импульса, массы и их граничных условий) и позволяет вычислять входящие в уравения переноса кинетическая коэффициенты. СТАТИСТИЧЕСКАЯ ТЕРМОДИНАМИКА т. устанавливает количественное связь между микро- и макросвойствами физических и химический систем. Расчетные методы СТАТИСТИЧЕСКАЯ ТЕРМОДИНАМИКА т. используются во всех направлениях современной теоретич. химии.

Основные понятия. Для статистич. описания макроскопич. систем Дж. Гиббсом (1901) предложено использовать понятия статистич. ансамбля и фазового пространства, что позволяет применять к решению задач методы теории вероятности. Статистич. ансамбль-совокупность очень большого числа одинаковых систем многие частиц (т. е. "копий" рассматриваемой системы), находящихся в одном и том же макросостоянии, которое определяется параметрами состояния; микросостояния системы при этом могут различаться. Осн. статистич. ансамбли-микроканонич., канонич., большой канонич. и изобарно-изотермический.

Микроканонич. ансамбль Гиббса используетя при рассмотрении изолированных систем (не обменивающихся энергией E с окружающей средой), имеющих постоянный объем V и число одинаковых частиц N (Е, V и N- параметры состояния системы). Канонич. ансамбль Гиббса используется для описания систем постоянного объема, находящихся в тепловом равновесии с окружающей средой (абс. температура Т) при постоянном числе частиц N (параметры состояния V, Т, N). Большой канонич. ансамбль Гиббса используется для описания открытых систем, находящихся в тепловом равновесии с окружающей средой (температура Т) и материальном равновесии с резервуаром частиц (осуществляется обмен частицами всех сортов через "стенки", окружающие систему объемом V). Параметры состояния такой системы-V, Т и m—химический потенциал частиц. Изобарно-изотермодинамически ансамбль Гиббса используется для описания систем, находящихся в тепловом и механические равновесии с окружающей средой при постоянном давлении P (параметры состояния Т, P, N).

Фазовое пространство в статистич. механике-многомерное пространство, осями которого служат все обобщенные координаты qi и сопряженные им импульсы pi(i =1,2,..., М) системы с М степенями свободы. Для системы, состоящей из N атомов, qi и pi соответствуют декартовой координатеи компоненте импульса (a = х, у, z) некоторого атома j и М = 3N. Совокупность координат и импульсов обозначаются q и p соответственно. Состояние системы изображается точкой в фазовом пространстве размерности 2М, а изменение состояния системы во времени-движением точки вдоль линии, называют фазовой траекторией. Для статистич. описания состояния системы вводятся понятия фазового объема (элемента объема фазового пространства) и функции распределения f(p, q), которая характеризует плотность вероятности нахождения точки, изображающей состояние системы, в элементе фазового пространства вблизи точки с координатами р, q. В квантовой механике вместо фазового объема используют понятие дискретного энергетич. спектра системы конечного объема, так как состояние отдельной частицы определяется не импульсом и координатами, а волновой функцией, которой в стационарном динамич. состоянии системы соответствует энергетич. спектр квантовых состояний.

Функция распределения классич. системы f(p, q)характеризует плотность вероятности реализации данного микро состояния (р, q) в элементе объема dГ фазового пространства. Вероятность пребывания N частиц в бесконечно малом объеме фазового пространства равна:


где dГN - элемент фазового объема системы в единицах h3N, h-постоянная Планка; делитель N! учитывает тот факт, что перестановка тождеств. частиц не меняет состояния системы. Функция распределения удовлетворяет условию нормировки тf(p, q)dГN= 1, так как система достоверно находится в к.-л. состоянии. Для квантовых систем функция распределения определяет вероятность wi, N нахождения системы из N частиц в квантовом состоянии, задаваемом набором квантовых чисел i, с энергией Ei,N при условии нормировки


Среднее значение в момент времени т (т.е. по бесконечно малому интервалу времени от т до т + dт)любой физических величины А(р, q), являющейся функцией координат и импульсов всех частиц системы, с помощью функции распределения вычисляется по правилу (в т.ч. и для неравновесных процессов):


Интегрирование по координатам проводится по всему объему системы, а интегрирование по импульсам от — , до +,. Состояние термодинамическое равновесия системы следует рассматривать как предел т:,. Для равновесных состояний функции распределения определяются без решения уравения движения составляющих систему частиц. Вид этих функций (одинаковый для классич. и квантовых систем) был установлен Дж. Гиббсом (1901).

В микроканонич. ансамбле Гиббса все микросостояния с данной энергией Е равновероятны и функция распределения для классич. систем имеет вид:

f(p,q) = Ad[H(p,q)-E],

где d-дельта-функция Дирака, Н(р,q)-функция Гамильтона, представляющая собой сумму кинетическая и потенц. энергий всех частиц; постоянная А определяется из условия нормировки функции f(p, q). Для квантовых систем при точности задания квантового состояния, равной величине DE, в соответствии с соотношением неопределенностей между энергией и временем (между импульсом и координатой частицы), функция w(Ek) = [g(E, N, V)]-1, если ЕEkE + DE, и w(Ek) = 0, если Ek < Е и Ek > E + DE. Величина g(E, N, V)-так называемой статистич. вес, равный числу квантовых состояний в энергетич. слое DE. Важное соотношение СТАТИСТИЧЕСКАЯ ТЕРМОДИНАМИКА т.-связь энтропии системы со статистич. весом:

S(E, N, V) = klng(E, N, V), где k-Больцмана постоянная.

В канонич. ансамбле Гиббса вероятность нахождения системы в микросостоянии, определяемом координатами и импульсами всех N частиц или значениями Ei,N, имеет вид: f(p, q) = exp {[F - H(p, q)]/kT}; wi,N = exp[(F - Ei,N)/kT], где F-свободный энергия (энергия Гельмгольца), зависящая от значений V, Т, N:

F = -kTlnZN,

где ZN-статистич. сумма (в случае квантовой системы) или статистич. интеграл (в случае классич. системы), определяемые из условия нормировки функций wi,N или f(p, q):


ZN = тexp[-H(р, q)/kT]dpdq/(N!h3N)

(сумма по г берется по всем квантовым состояниям системы, а интегрирование проводится по всему фазовому пространству).

В большом канонич. ансамбле Гиббса функция распределения f(p, q) и статистич. сумма X, определяемая из условия нормировки, имеют вид:


где W-термодинамическое потенциал, зависящий от переменных V, Т, m (суммирование ведется по всем целым положит. N). В изобарно-изотермодинамически ансамбле Гиббса функция распределения и статистич. сумма Q, определяемая из условия нормировки, имеют вид:


где G- энергия Гиббса системы (изобарно-изотермодинамически потенциал, свободный энтальпия).

Для вычисления термодинамическое функции можно использовать любое распределение: они эквивалентны друг другу и соответствуют разным физических условиям. Микроканонич. распределение Гиббса применяется главным образом в теоретич. исследованиях. Для решения конкретных задач рассматривают ансамбли, в которых есть обмен энергией со средой (канонич. и изобарно-изотермодинамически) или обмен энергией и частицами (большой канонич. ансамбль). Последний особенно удобен для изучения фазового и химический равновесий. Статистич. суммы ZN и Q позволяют определить энергию Гельмгольца F, энергию Гиббса G, а также термодинамическое свойства системы, получаемые дифференцированием статистич. суммы по соответствующим параметрам (в расчете на 1 моль вещества): внутр. энергию U = RT2(9lnZN/9T)V, энтальпию H = RT2(9lnQ/9T)P, энтропию S = RlnZN + RT(9lnZN/9T)V = = R ln Q + RT(9ln Q/9T)P, теплоемкость при постоянном объеме СV = 2RT(9lnZN/9T)V + RT2(92 lnZN/9T2)V, теплоемкость при постоянном давлении СР = 2RT (9lnZN/9T)P + + RT2(92lnZN/9T2)P и т.д. Соотв. все эти величины приобретают и статистич. смысл. Так, внутренняя энергия отождествляется со средней энергией системы, что позволяет рассматривать первое начало термодинамики как закон сохранения энергии при движении составляющих систему частиц; свободный энергия связана со статистич. суммой системы, энтропия-с числом микросостояний g в данном макросостоянии, или статистич. весом макросостояния, и, следовательно, с его вероятностью. Смысл энтропии как меры вероятности состояния сохраняется по отношению к произвольным (неравновесным) состояниям. В состоянии равновесия энтропия изолир. системы имеет максимально возможное значение при заданных внешний условиях (Е, V, N), т. е. равновесное состояние является наиболее вероятным состоянием (с макс. статистич. весом). Поэтому переход из неравновесного состояния в равновесное есть процесс перехода из менее вероятных состояний в более вероятное. В этом заключается статистич. смысл закона возрастания энтропии, согласно которому энтропия замкнутой системы может только увеличиваться (см. Второе начало термодинамики). При температуре абс. нуля любая система находится в основные состоянии, в котором w0 = 1 и S = 0. Это утверждение представляет собой третье начало термодинамики (см. Тепловая теорема). Существенно, что для однозначного определения энтропии нужно пользоваться квантовым описанием, так как в классич. статистике энтропия может быть определена только с точностью до произвольного слагаемого.

Идеальные системы. Расчет статистич. сумм большинства систем представляет сложную задачу. Она существенно упрощается в случае газов, если вкладом потенц. энергии в полную энергию системы можно пренебречь. В этом случае полная функция распределения f(p, q) для N частиц идеальной системы выражается через произведение одно-частичных функций распределения f1(p, q):


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

Статистика Ферми-Дирака описывает распределение в системе тождеств. частиц с полуцелым спином 1/2, 3/2,... в единицах ђ = h/2p. Частица (или квазичастица), подчиняющаяся указанной статистике, называют фермионом. К фер-мионам относятся электроны в атомах, металлах и полупроводниках, атомные ядра с нечетным атомным номером, атомы с нечетной разностью атомного номера и числа электронов, квазичастицы (например, электроны и дырки в твердых телах) и т.д. Данная статистика была предложена Э.Ферми в 1926; в том же году П.Дирак выяснил ее квантовомеханические смысл. Волновая функция системы фермионов антисимметрична, т.е. меняет свой знак при перестановке координат и спинов любой пары тождеств. частиц. В каждом квантовом состоянии может находиться не более одной частицы (см. Паули принцип). Среднее число частиц ni идеального газа фермионов, находящихся в состоянии с энергией Ei, определяется функцией распределения Ферми-Дирака:

ni={1+exp[(Ei-m)/kT]}-1,

где i-набор квантовых чисел, характеризующих состояние частицы.

Статистика Бозе-Эйнштейна описывает системы тождеств. частиц с нулевым или целочисленным спином (0, ђ, 2ђ, ...). Частица или квазичастица, подчиняющаяся указанной статистике, называют бозоном. Данная статистика была предложена Ш. Бозе (1924) для фотонов и развита А. Эйнштейном (1924) применительно к молекулам идеального газа, рассматриваемым как составные частицы из четного числа фермионов, например атомные ядра с четным суммарным числом протонов и нейтронов (дейтрон, ядро 4Не и т.д.). К бозонам относятся также фононы в твердом теле и жидком 4Не, экситоны в полупроводниках и диэлектриках. Волновая функция системы симметрична относительно перестановки любой пары тождеств. частиц. Числа заполнения квантовых состояний ничем не ограничены, т.е. в одном состоянии может находиться любое число частиц. Среднее число частиц ni идеального газа бозонов, находящихся в состоянии с энергией Еi описывается функцией распределения Бозе-Эйнштейна:

ni={exp[(Ei-m)/kT]-1}-1.

Статистика Больцмана представляет собой частный случай квантовой статистики, когда можно пренебречь квантовыми эффектами (высокие температуры). В ней рассматривается распределение частиц идеального газа по импульсам и координатам в фазовом пространстве одной частицы, а не в фазовом пространстве всех частиц, как в распределениях Гиббса. В качестве миним. единицы объема фазового пространства, имеющего шесть измерений (три координаты и три проекции импульса частицы), в соответствии с квантовомеханические соотношением неопределенностей, нельзя выбрать объем меньший, чем h3. Среднее число частиц ni идеального газа, находящихся в состоянии с энергией Ei, описывается функцией распределения Больцмана:

ni=exp[(m-Ei)/kT].

Для частиц, которые движутся по законам классич. механики во внешний потенц. поле U(r), статистически равновесная функция распределения f1(p,r) по импульсам p и координатам r частиц идеального газа имеет вид: f1(p,r) = A ехр{ - [р2/2m + U(r)]/kT}. Здесь р2/2т-кинетическая энергия молекул массой ш, постоянная А определяется из условия нормировки. Данное выражение часто называют распределением Максвелла-Больцмана, а распределением Больцмана называют функцию

n(r) = n0 ехр[-U(r)]/kT],

где n(r) = тf1(p, r)dp - плотность числа частиц в точке r (n0-плотность числа частиц в отсутствие внешний поля). Распределение Больцмана описывает распределение моле кул в поле тяготения (барометрич. формула), молекул и высокодисперсных частиц в поле центробежных сил, электронов в невырожденных полупроводниках, а также используется для расчета распределения ионов в разбавл. растворах электролитов (в объеме и на границе с электродом) и т. п. При U(r) = 0 из распределения Максвелла - Больц-мана следует распределение Максвелла, описывающее распределение по скоростям частиц, находящихся в ста-тистич. равновесии (Дж. Максвелл, 1859). Согласно этому распределению, вероятное число молекул в единице объема компоненты скоростей которых лежат в интервалах от ui до ui + dui (i = x, у, z), определяется функцией:


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

Сумма по состояниям молекулы. Статистич. сумма идеального газа в канонич. ансамбле Гиббса выражается через сумму по состояниям одной молекулы Q1:

где Еi - энергияi-го квантового уровня молекулы (i = О соответствует нулевому уровню молекулы), gi-статистич. вес i-го уровня. В общем случае отдельные виды движения электронов, атомов и групп атомов в молекуле, а также движение молекулы как целого взаимосвязаны, однако приближенно их можно рассматривать как независимые. Тогда сумма по состояниям молекулы может быть представлена в виде произведения отдельных составляющих, связанных с по-ступат. движением (Qпост) и с внутримол. движениями (Qвн):

Q1 = Qпост•Qвн, Qпост= l(V/N),

где l = (2pmkТ/h2)3/2. Для атомов Qвн представляет собой сумму по электронным и ядерным состояниям атома; для молекул Qвн- сумма по электронным, ядерным, колебательное и вращательное состояниям. В области температур от 10 до 103 К обычно используют приближенное описание, в котором каждый из указанных типов движения рассматривается независимо: Qвн = Qэл•Qяд•Qвращ•Qкол/g, где g - число симметрии, равное числу тождество. конфигураций, возникающих при вращении молекулы, состоящей из одинаковых атомов или групп атомов.

Сумма по состояниям электронного движения Qэл равна статистич. весу Рт основные электронного состояния молекулы. Во многие случаях основные уровень невырожден и отделен от ближайшего возбужденного уровня значительной энергией: (Рт = 1). Однако в ряде случаев, например для молекулы О2, Рт = з, в основные состоянии момент количества движения молекулы отличен от нуля и имеет место вырождение энергетических уровней, а энергии возбужденных состояний может быть достаточно низкими. Сумма по состояниям Qяд, обусловленная вырождением ядерных спинов, равна:

где si-спин ядра атома i, произведение беретсяпо всем атомам молекулы. Сумма по состояниям колебательное движения молекулы где vi-частоты нор мальных колебаний, n-число атомов в молекуле. Сумму по состояниям вращательное движений многоатомной молекулы с большими моментами инерции можно рассматривать классически [высокотемпературное приближение, T/qi1, где qi = h2/8p2kIi (i = x, у, z), It-главный момент инерции вращения вокруг оси i]: Qвр = (pT3/qxqyqz)1/2. Для линейных молекул с моментом инерции I статистич. сумма Qвр = T/q, где q = h2/8p2*kI.

При расчетах при температурах выше 103К необходимо учитывать ангармонизм колебаний атомов, эффекты взаимодействие колебательное и вращательное степеней свободы (см. Нежесткие молекулы), а также мультиплетности электронных состояний, заселенности возбужденных уровней и т. д. При низких температурах (ниже 10 К) необходимо учитывать квантовые эффекты (особенно для двухатомных молекул). Так, вращательное движение гетеро-ядерной молекулы АВ описывается по формуле:

l-номервращательное состояния, а для гомоядерных молекул А2(особенно для молекул водорода Н2, дейтерия D2, трития Т2) ядерные и вращательное степени свободы взаимодействие друг с другом: Qяд.вращ. Qяд•Qвращ.

Знание суммы по состояниям молекулы позволяет рассчитать термодинамическое свойства идеального газа и смеси идеальных газов, в т.ч. константы химический равновесия, равновесную степень ионизации и т.п. Важное значение в теории абс. скоростей реакций имеет возможность расчета константы равновесия процесса образования активир. комплекса (переходного состояния), которое представляется как модифицир. частица, одна из колебательное степеней свободы которой заменена степенью свободы поступат. движения.

Неидеальные системы. В реальных газах молекулы взаимодействие друг с другом. В этом случае сумма по состояниям ансамбля не сводится к произведению сумм по состояниям отдельных молекул. Если считать, что межмол. взаимодействие не влияют на внутр. состояния молекул, статистич. сумма системы в классич. приближении для газа, состоящего из N тождеств. частиц, имеет вид:

где

Здесь <2N-конфигурац. интеграл, учитывающий взаимодействие молекул. Наиб, часто потенц. энергия молекул U рассматривается в виде суммы парных потенциалов: U = =где U(rij)- потенциал центральное сил, зависящий от расстояния rij между молекулами i и j. Учитывают также многочастичные вклады в потенц. энергию, эффекты ориентации молекул и т.д. Необходимость расчета конфигурац. интеграла возникает при рассмотрении любых конденсир. фаз и границ раздела фаз. Точное решение задачи многие тел практически невозможно, поэтому для расчета статистич. суммы и всех термодинамическое свойств, получаемых из статистич. суммы дифференцированием по соответствующим параметрам, используют различные приближенные методы.

Согласно так называемой методу групповых разложений, состояние системы рассматривается в виде совокупности комплексов (групп), состоящих из разного числа молекул, и конфигурац. интеграл распадается на совокупность групповых интегралов. Такой подход позволяет представить любую термодинамическое функцию реального газа в виде ряда по степеням плотности. Наиб. важное соотношение такого рода - вириальное уравение состояния.

Для теоретич. описания свойств плотных газов, жидкостей и твердых тел, растворов неэлектролитов и электролитов и границ раздела в этих системах более удобным, чем прямой расчет статистич. суммы, является метод n-частичных функций распределения. В нем вместо подсчета статистич. веса каждого состояния с фиксир. энергией используют соотношения между функциями распределения fn, которые характеризуют вероятность нахождения частиц одновременно в точках пространства с координатами r1,..., rn; при n = N fN = b тf(p, r)dp (здесь и ниже qi = ri). Одночастичная функция f1(r1) (n = 1) характеризует распределение плотности вещества. Для твердого тела это периодической функция с максимумами в узлах кристаллич. структуры; для газов или жидкостей в отсутствие внешний поля это постоянная величина, равная макроскопич. плотности вещества р. Двухчастичная функция распределения (п = 2) характеризует вероятность нахождения двух частиц в точках 1 и 2, она определяет так называемой корреляционную функцию g(|r1 — r2|) = f2(r1, r2)/r2, характеризующую взаимную корреляцию в распределении частиц. Соответствующую опытную информацию дает рентгеновский структурный анализ.

Функции распределения размерности n и n + 1 связаны бесконечной системой зацепляющихся интегродифференц. уравений Боголюбова-Борна-Грина-Кирквуда-Ивона, решение которых чрезвычайно сложно, поэтому эффекты корреляции между частицами учитывают введением различные аппроксимаций, которые определяют, каким бразом функция fn выражается через функции меньшей размерности. Соотв. разработано несколько приближенных методов расчета функций fn, а через них-всех термодинамическое характеристик рассматриваемой системы. Наиб. применение имеют приближения Перкус-Иевика и гиперцепное.

Решеточные модели конденсир. состояния нашли широкое применение при термодинамическое рассмотрении практически всех физических-химический задач. Весь объем системы разбивается на локальные области с характерным размером порядка размера молекулы u0. В общем случае в разных моделях размер локальной области может быть как больше, так и меньше u0; в большинстве случаев они совпадают. Переход к дискретному распределению молекул в пространстве существенно упрощает подсчет различные конфигураций молекул. Решеточные модели учитывают взаимодействие молекул друг с другом; энергия взаимодействие описывается энергетич. параметрами. В ряде случаев решеточные модели допускают точные решения, что позволяет оценить характер используемых приближений. С их помощью возможно рассмотрение многочастичных и специфический взаимодействие, ориентац. эффектов и т. п. Решеточные модели являются основными при изучении и проведении прикладных расчетов растворов неэлектролитов и полимеров, фазовых переходов, критических явлений, сильно неоднородных систем.

Численные методы определения термодинамическое свойств приобретают все большее значение по мере развития вычислит. техники. В методе Монте-Карло осуществляется прямой расчет многомерных интегралов, что позволяет непосредственно получить статистич. среднее наблюдаемой величины А(r1.....rN) по любому из статистич. ансамблей (например, А - энергия системы). Так, в канонич. ансамбле термодинамическое среднее имеет вид:


Данный метод применим практически ко всем системам; получаемые с его помощью средние величины для ограниченных объемов (N = 102-105) служат хорошим приближением для описания макроскопич. объектов и могут рассматриваться как точные результаты.

В методе мол. динамики эволюция состояния системы рассматривается с помощью численного интегрирования уравений Ньютона для движения каждой частицы (N = = 102-105) при заданных потенциалах межчастичного взаимодействия. Равновесные характеристики системы получаются при усреднении по фазовым траекториям (по скоростям и координатам) на больших временах, после установления максвелловского распределения частиц по скоростям (так называемой период термализации).

Ограничения в использовании численных методов в основные определяются возможностями ЭВМ. Спец. вычислит. приемы позволяют обходить сложности, связанные с тем, что рассматривается не реальная система, а небольшой объем; это особенно важно при учете дальнодействующих потенциалов взаимодействие, анализе фазовых переходов и т.п.

Физическая кинетика - раздел статистич. физики, который дает обоснование соотношениям термодинамики необратимых процессов, описывающим перенос энергии, импульса и массы, а также влияние на эти процессы внешний полей. Кинетич. коэффициенты-макроскопич. характеристики сплошной среды, определяющие зависимости потоков физических величин (теплоты, импульса, массы компонентов и др.) от вызывающих эти потоки градиентов температуры, концентрации, гидродинамич. скорости и др. Необходимо различать коэффициенты Онсагера, входящие в уравения, связывающие потоки с термодинамическое силами (термодинамическое уравения движения), и коэффициенты переноса (диффузии, теплопроводности, вязкости и т. п.), входящие в уравения переноса. Первые может быть выражены через вторые с помощью соотношений между макроскопич. характеристиками системы, поэтому в дальнейшем будут рассматриваться лишь коэффициент переноса.

Для расчета макроскопич. коэффициент переноса необходимо провести усреднение по вероятностям реализаций элементарных актов переноса с помощью неравновесной функции распределения. Главная сложность заключается в том, что аналит. вид функции распределения f(р, q, т) (т-время) неизвестен (в отличие от равновесного состояния системы, которое описывается с помощью функций распределения Гиббса, получаемых при т : ,). Рассматривают n-частичные функции распределения fn(r, q, т), которые получают из функций f(р, q, т) усреднением по координатам и импульсам остальных (N — п) частиц:


Для них может быть составлена система уравений, позволяющая описать произвольные неравновесные состояния. Решение этой системы уравений очень сложно. Как правило, в кинетическая теории газов и газообразных квазичастиц в твердом теле (фермионов и бозонов) используется лишь уравение для одно-частичной функции распределения f1. В предположении об отсутствии корреляции между состояниями любых частиц (гипотеза мол. хаоса) получено т. называют кинетическая уравение БоЛьцмана (Л. Больцман, 1872). Это уравение учитывает изменение функции распределения частиц под действием внешний силы F(r, т) и парных столкновений между частицами:


где f1(u, r, т) и -функции распределения частиц до столкновения, f"1(u", r, т) и-функции распределения после столкновения; u и-скорости частиц до столкновения, u" и -скорости тех же частиц после столкновения, и = |u —|-модуль относит. скорости сталкивающихся частиц, q - угол между относит. скоростью u — сталкивающихся частиц и линией, соединяющей их центры, s(u,q)dW-дифференц. эффективное сечение рассеяния частиц на телесный угол dW в лабораторная системе координат, зависящее от закона взаимодействие частиц. Для модели молекул в виде упругих жестких сфер, имеющих радиус R, принимается s = 4R2cosq. В рамках классич. механики дифференц. сечение выражается через параметры столкновения b и e (соответственно прицельное расстояние и азимутальный угол линии центров): sdW = bdbde, а молекулы рассматриваются как центры сил с потенциалом, зависящим от расстояния. Для квантовых газов выражение для дифференц. эффективного сечения получают на основе квантовой механики, с учетом влияния эффектов симметрии на вероятность столкновения.

Если система находится в статистич. равновесии, интеграл столкновений Stf равен нулю и решением кинетическая уравения Больцмана будет распределение Максвелла. Для неравновесных состояний решения кинетическая уравнения Больцмана обычно ищут в виде разложения в ряд функции f1(u, r, т) по малым параметрам относительно функции распределения Максвелла. В простейшем (релаксационном) приближении интеграл столкновений аппроксимируется как Stf = = -(f1 - )/т*, где т*-среднее время релаксации. Зная решение уравения Больцмана, можно определить плотность числа частиц газа в точке r в момент времени т: n = = тf1(u, r, т)du и др. средние значения величин, характеризующих состояние газа: средняя (по абс. значениям) скорость молекул, средний квадрат их скорости и т.д., а через них рассчитать гидродинамич. характеристики процессов пере носа энергии, импульса и числа частиц, которые могут в нем происходить.

Данный подход реализуется при исследовании процессов в газовых смесях, в многоатомных газах с учетом внутр. степеней свободы молекул (колебательное, вращательное и т.д.), в плотных газах, при изучении влияния стенок сосудов на распределения молекул газа в приповерхностной области и многие др. задачах. Анализ решений кинетическая уравения Больц-мана позволяет обосновать область применимости условия локального термодинамическое равновесия и определить вклады в поток, обусловленные неравновесностью потока. Неравновесный поток импульса дает сдвиговую вязкость; для газов с внутр. степенями свободы молекул он дополнительно содержит член, обусловленный объемной вязкостью. Плотность потока энергии пропорциональна градиенту температуры (обычная теплопроводность), а в случае смеси газов она содержит член, пропорциональный градиенту концентраций (эффект Дюфура). Поток вещества в смеси газов содержит член, пропорциональный градиенту концентрации (обычная диффузия), и член, пропорциональный градиенту температуры (термодиффузия). Физ. кинетика дает для этих коэффициент пропорциональности выражения через эффективные сечения столкновения, следовательно через потенциалы межмол. взаимодействий. Коэф. переноса удовлетворяют принципу симметрии, выражающему симметрию уравений механики относительно изменения знака времени (теорема Онсагера).

Для сильно разреженных газов, когда длина свободный пробега молекулы l между двумя последоват. соударениями значительно превосходит характерный размер системы L (1/L1), гидродинамич. уравения неприменимы и необходимо решать кинетическая уравение Больцмана с определенными граничными условиями на поверхности. Эти условия определяются функцией распределения молекул, рассеянных из-за взаимодействие со стенкой.

Кинетич. уравение Больцмана для фононов кристаллич. решетки позволяет исследовать теплопроводность и поглощение звука в диэлектриках, а кинетическая уравение Больцмана для электронов с учетом их взаимодействие с фононами-электрич. сопротивление и объясняет гальваномагн., термоэлектрич., термомагн. и др. явления в металлах, а также особенности поведения сверхпроводников в высокочастотных полях.

Для описания неравновесных процессов в жидкостях од-ночастичная функция распределения f1 не раскрывает специфики явлений и требуется рассмотрение двухчастичной функции распределения f2. Однако для достаточно медленных процессов и в случаях, когда масштабы пространств. неодно-родностей значительно меньше масштаба корреляции между частицами жидкости, можно использовать локально равновесную одночастичную функцию распределения с температурой, химический потенциалами и гидродинамич. скоростью, которые соответствуют рассматриваемому малому объему жидкости. К ней можно найти поправку, пропорциональную градиентам температуры, гидродинамич. скорости и химический потенциалам компонентов, и вычислить потоки импульсов, энергии и вещества, а также обосновать уравения Навье-Стокса, теплопроводности и диффузии. В этом случае коэффициент переноса оказываются пропорциональными пространственно-временным корреляц. функциям потоков энергии, импульса и вещества каждого компонента.

Для описания процессов переноса вещества в твердых телах и на границах раздела с твердым телом широко используется решеточная модель конденсир. фазы. Эволюция состояния системы описывается основные кинетическая уравением (master equation) относительно функции распределения P(q, т):

где P(q,т)= тf(p,q,т)du- функция распределения, усредненная по импульсам (скоростям) всех N частиц, описывающая распределение частиц по узлам решеточной структуры (их число равно Ny, N < Ny), q- номер узла или его координата. В модели "решеточного газа" частица может находиться в узле (узел занят) или отсутствовать (узел свободен); W(q : q")-вероятность перехода системы в единицу времени из состояния q, описываемого полным набором координат частиц, в др. состояние q". Первая сумма описывает вклад всех процессов, в которых осуществляется переход в данное состояние q, вторая сумма-выход из этого состояния. В случае равновесного распределения частиц (т : ,) P(q) = exp[-H(q)/kT]/Q, где Q-статистич. сумма, H(q)-энергия системы в состоянии q. Вероятности перехода удовлетворяют детального равновесия принципу: W(q":q)exp[-H(q")/kT] = W(q : q")ехр[-H(q)/kТ]. На базе уравений для функций P(q,т) строят кинетическая уравения для n-частичных функций распределения, которые получают путем усреднения по расположениям всех остальных (N — п) частиц. Для малых h кинетическая уравения может быть решены аналитически или численно и с их помощью может быть получены коэффициент диффузии, самодиффузии, сдвиговой вязкости, подвижности и т.п. Такой подход применим к процессам переноса вещества в моноатомных кристаллах, сплавах, оксидных кристаллах, ионитах и т.д., к процессам переноса вещества через границу с твердым телом, роста кристаллов, фазовым превращениям и т. п. Для межфазного переноса, из-за различий в характерных временах протекания элементарных процессов миграции частиц, важную роль играет вид граничных условий на границах раздела фаз.

Для малых систем (число узлов Ny = 102 — 105) система уравений относительно функции P(q,т) может быть решена численно методом Монте-Карло. Этап релаксации системы к равновесному состоянию позволяет рассмотреть различные переходные процессы при исследовании кинетики фазовых превращений, роста кристаллов, кинетики поверхностных реакций и т.д. и определить их динамич. характеристики, в том числе и коэффициент переноса.

Для расчета коэффициент переноса в газообразных, жидких и твердых фазах, а также на границах раздела фаз активно используются разнообразные варианты метода мол. динамики, который позволяет детально проследить за эволюцией системы от времен ~10-15с до ~10-10с (на временах порядка 10-10— 10-9с и более используются так называемой уравения Ланжевена, это уравения Ньютона, содержащие в правой части стохастич. слагаемое).

Для систем с химический реакциями на характер распределения частиц большое влияние оказывает соотношение между характерными временами переноса реагентов и их химический превращения. Если скорость химический превращения мала, распределение частиц не сильно отличается от случая, когда реакция отсутствует. Если скорость реакции велика, ее влияние на характер распределения частиц велико и использовать средние концентрации частиц (т.е. функции распределения с n = 1), как это делается при использовании закона действующих масс, нельзя. Необходимо более детально описывать распределение реагентов с помощью функций распределения fn с n > 1. Важное значение при описании реакционное потоков частиц на поверхности и скоростей диффузионно-контролируемых реакций имеют граничные условия (см. Макрокинетика).

Литература: Кубо Р., Статистическая механика, пер. с англ., М.,1967; Зубарев Д. Н., Неравновесная статистическая термодинамика, М., 1971; Исихара А., Статистическая физика, пер. с англ., М., 1973; Ландау Л. Д., Лифшиц E. МL Статистическая физика, 3 изд., ч. 1, М., 1976; Ферцигер Дж., Капер Г., Математическая теория процессов переноса в газах, пер. с англ., М., 1976; Румер Ю. Б., Рывкин М. Ш., Термодинамика. Статистическая физика и кинетика, 2 изд., М., 1977; Смирнова H. А., Методы статистической термодинамики в физической химии, 2 изд., М., 1982; Берклеевский курс физики, пер. с англ., 3 изд., т. 5-Рейф Ф., Статистическая физика, М., 1986; Товбин Ю.К., Теория физико-химических процессов на границе газ-твердое тело, М., 1990. Ю.К. Товбин.


Химическая энциклопедия. Том 4 >> К списку статей


[каталог]  [статьи]  [доска объявлений]    [обратная связь]

 

 

Реклама
Huawei S5720-36C-EI-28S-AC
рекламно-производственные компании наружной рекламы в москве
вентилятор с монтажной панелью
автоматизация спб

Рекомендуемые книги

Введение в химию окружающей среды.

Книга известных английских ученых раскрывает основные принципы химии окружающей среды и их действие в локальных и глобальных масштабах. Важный аспект книги заключается в раскрытии механизма действия природных геохимических процессов в разных масштабах времени и влияния на них человеческой деятельности. Показываются химический состав, происхождение и эволюция земной коры, океанов и атмосферы. Детально рассматриваются процессы выветривания и их влияние на химический состав осадочных образований, почв и поверхностных вод на континентах. Для студентов и преподавателей факультетов биологии, географии и химии университетов и преподавателей средних школ, а также для широкого круга читателей.

Химия и технология редких и рассеянных элементов.

Книга представляет собой учебное пособие по специальным курсам для студентов химико-технологических вузов. В первой части изложены основы химии и технологии лития, рубидия, цезия, бериллия, галлия, индия, таллия. Во второй части книги изложены основы химии и технологии скандия, натрия, лантана, лантаноидов, германия, титана, циркония, гафния. В третьей части книги изложены основы химии и технологии ванадия, ниобия, тантала, селена, теллура, молибдена, вольфрама, рения. Наибольшее внимание уделено свойствам соединений элементов, имеющих значение в технологии. В технологии каждого элемента описаны важнейшие области применения, характеристика рудного сырья и его обогащение, получение соединений из концентратов и отходов производства, современные методы разделения и очистки элементов. Пособие составлено по материалам, опубликованным из советской и зарубежной печати по 1972 год включительно.

 

 



Рейтинг@Mail.ru Rambler's Top100

Copyright © 2001-2012
(22.10.2017)