Физическая кинетика

Эта статья находится на начальном уровне проработки, в одной из её версий выборочно используется текст из источника, распространяемого под свободной лицензией
Материал из энциклопедии Руниверсалис

Физи́ческая кине́тика (др.-греч. κίνησις — движение) — микроскопическая теория процессов в неравновесных средах. В кинетике методами квантовой или классической статистической физики изучают процессы переноса энергии, импульса, заряда и вещества в различных физических системах (газах, плазме, жидкостях, твёрдых телах) и влияние на них внешних полей[⇨]. В отличие от термодинамики неравновесных процессов и электродинамики сплошных сред, кинетика исходит из представления о молекулярном строении рассматриваемых сред, что позволяет вычислить из первых принципов кинетические коэффициенты, диэлектрические и магнитные проницаемости и другие характеристики сплошных сред. Физическая кинетика включает в себя кинетическую теорию газов из нейтральных атомов или молекул, статистическую теорию неравновесных процессов в плазме[⇨], теорию явлений переноса в твёрдых телах (диэлектриках, металлах и полупроводниках) и жидкостях, кинетику магнитных процессов и теорию кинетических явлений, связанных с прохождением быстрых частиц через вещество. К ней же относятся теория процессов переноса в квантовых жидкостях и сверхпроводниках и кинетика фазовых переходов[⇨].

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

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

Кинетическое уравнение

Основной метод физической кинетики — решение кинетического уравнения Больцмана для одночастичной функции распределения [math]\displaystyle{ f(x,\;p,\;t) }[/math] молекул в фазовом пространстве их координат [math]\displaystyle{ x }[/math] и импульсов [math]\displaystyle{ p }[/math]. Это уравнение ввёл Больцман в 1872 году[1]. Функция распределения удовлетворяет кинетическому уравнению[2]:

[math]\displaystyle{ \frac{d f}{d t}=\frac{\partial f}{\partial t}+\frac{\vec{p}}{m}\frac{\partial f}{\partial\vec{x}}+\vec{F}\frac{\partial f}{\partial\vec{p}}=\mathrm{St}\,f, }[/math]

где [math]\displaystyle{ \mathrm{St} }[/math] — интеграл столкновений, определяющий разность числа частиц, приходящих в элемент объёма вследствие прямых столкновений и убывающих из него вследствие обратных столкновений. Для одноатомных молекул или для многоатомных, но без учёта их внутренних степеней свободы[3]

[math]\displaystyle{ \mathrm{St}\,f=\int\omega\cdot(f'f'_1-ff_1)\,dp_1dp'dp'_1, }[/math]

где [math]\displaystyle{ \omega }[/math] — вероятность столкновения, связанная с дифференциальным эффективным сечением рассеяния.

[math]\displaystyle{ \omega\,dp'dp'_1=|v-v_1|\,d\sigma, }[/math]

где [math]\displaystyle{ p }[/math], [math]\displaystyle{ p_1 }[/math] — импульсы молекул до столкновения, [math]\displaystyle{ v }[/math], [math]\displaystyle{ v_1 }[/math] — соответственно скорости, [math]\displaystyle{ p' }[/math], [math]\displaystyle{ p'_1 }[/math] — их импульсы после столкновения, [math]\displaystyle{ f }[/math], [math]\displaystyle{ f_1 }[/math] — функции распределения молекул до столкновения, [math]\displaystyle{ f' }[/math], [math]\displaystyle{ f'_1 }[/math] — их функции распределения после столкновения.

Для газа из сложных молекул, обладающих внутренними степенями свободы, их следует учитывать в функции распределения. Например, для двухатомных молекул с собственным моментом вращения M функции распределения будут зависеть также от [math]\displaystyle{ M }[/math].

Из кинетического уравнения следует теорема Больцмана — убывание со временем [math]\displaystyle{ H }[/math]-функции Больцмана (среднего логарифма функции распределения) или возрастание энтропии, так как она равна [math]\displaystyle{ H }[/math]-функции Больцмана с обратным знаком[4].

Уравнения переноса

Физическая кинетика позволяет получить уравнения баланса для средней плотности вещества, импульса и энергии. Например, для простого газа плотность [math]\displaystyle{ \rho }[/math], гидродинамическая скорость [math]\displaystyle{ \textbf{V} }[/math] и средняя энергия [math]\displaystyle{ \bar{E} }[/math] удовлетворяют уравнениям баланса[5]:

[math]\displaystyle{ \frac{\partial\rho}{\partial t}+\mathrm{div}(\rho \textbf{V})=0, }[/math] — также известное как уравнение непрерывности
[math]\displaystyle{ \frac{\partial}{\partial t}(\rho V_\alpha)+\sum_\beta{\frac{\partial\Pi_{\alpha\beta}}{\partial x_\beta}}=0, }[/math]
[math]\displaystyle{ \frac{\partial}{\partial t}n\bar{E}+\mathrm{div}(\textbf{q})=0, }[/math]
[math]\displaystyle{ \Pi_{\alpha\beta}=\int mV_\alpha V_\beta f\,dp, }[/math]

где [math]\displaystyle{ \Pi_{\alpha\beta} }[/math] — тензор плотности потока импульса, [math]\displaystyle{ m }[/math] — масса частиц, [math]\displaystyle{ n }[/math] — плотность числа частиц, [math]\displaystyle{ \textbf{q}=\int E\textbf{V}f\,dp }[/math] — плотность потока энергии.

Если состояние газа мало отличается от равновесного, то в малых элементах объёма устанавливается распределение, близкое к локально равновесному распределению Максвелла, с температурой, плотностью и гидродинамической скоростью, соответствующими рассматриваемой точке газа. В этом случае неравновесная функция распределения мало отличается от локально равновесной, и решение кинетического уравнения даёт малую поправку к последней, пропорциональную градиентам температуры [math]\displaystyle{ \nabla T }[/math] и гидродинамической скорости [math]\displaystyle{ \nabla \textbf{V} }[/math], так как [math]\displaystyle{ \mathrm{St}\,f_0=0 }[/math].

С помощью неравновесной функции распределения можно найти поток энергии (в неподвижной жидкости) [math]\displaystyle{ \textbf{q}=-\lambda\nabla T }[/math], где [math]\displaystyle{ \lambda }[/math] — коэффициент теплопроводности, и тензор плотности потока импульса[6]

[math]\displaystyle{ \Pi_{\alpha\beta}=\rho V_\alpha V_\beta+\delta_{\alpha\beta}P-\sigma'_{\alpha\beta}, }[/math]

где

[math]\displaystyle{ \sigma'_{\alpha\beta}=\eta\left[\left(\frac{\partial V_\alpha}{\partial x_\beta}+\frac{\partial V_\beta}{\partial x_\alpha}\right)-\frac{2}{3}\delta_{\alpha\beta}\,\mathrm{div}\,\textbf{V}\right]+\zeta\delta_{\alpha\beta}\mathrm{div}\,\textbf{V} }[/math]

— тензор вязких напряжении, [math]\displaystyle{ \eta }[/math] — коэффициент сдвиговой вязкости, [math]\displaystyle{ P }[/math] — давление. Эти два соотношения известны в механике сплошных сред как закон теплопроводности Фурье и закон вязкости Ньютона. Последнее слагаемое в [math]\displaystyle{ \sigma'_{\alpha\beta} }[/math] для газов с внутренними степенями свободы, где [math]\displaystyle{ \zeta }[/math] — коэффициент «второй», объёмной вязкости, проявляющейся лишь при движениях, в которых [math]\displaystyle{ \mathrm{div}\,\textbf{V}\ne 0 }[/math]. Для кинетических коэффициентов [math]\displaystyle{ \lambda }[/math], [math]\displaystyle{ \eta }[/math], [math]\displaystyle{ \zeta }[/math] получаются выражения через эффективные сечения столкновений, которые, в свою очередь, рассчитываются через константы молекулярных взаимодействий. В многокомпонентной смеси поток какого-либо компонента включает в себя диффузионный поток, пропорциональный градиенту концентрации вещества в смеси с коэффициентом диффузии, и поток за счет термодиффузии (эффект Соре), пропорциональный градиенту температуры с коэффициентом термодиффузии. Поток тепла включает помимо обычного потока за счёт теплопроводности, пропорционального градиенту температуры, дополнительную составляющую, пропорциональную градиентам концентраций компонентов и описывающую диффузионную теплопроводность (эффект Дюфура). Кинетическая теория даёт выражения для этих кинетических коэффициентов через эффективные сечения столкновений, при этом кинетические коэффициенты для перекрёстных явлений вследствие теоремы Онсагера оказываются равными. Эти соотношения являются следствием микроскопической обратимости уравнений движения частиц системы, то есть инвариантности их относительно обращения времени.

Уравнение баланса импульса с учётом выражения для плотности потока импульса через градиент скорости даёт уравнения Навье — Стокса, уравнение баланса энергии с учётом выражения для плотности потока тепла даёт уравнение теплопроводности, уравнение баланса числа частиц определённого сорта с учётом выражения для диффузионного потока даёт уравнение диффузии. Такой гидродинамический подход справедлив, если длина свободного пробега [math]\displaystyle{ \lambda }[/math] значительно меньше характерных размеров областей неоднородности.

Газы и плазма

Физическая кинетика позволяет исследовать явления переноса в разреженных газах, когда отношение длины свободного пробега [math]\displaystyle{ \lambda }[/math] к характерным размерам задачи [math]\displaystyle{ L }[/math] (то есть число Кнудсена [math]\displaystyle{ \lambda/L }[/math]) уже не очень мало́ и имеет смысл рассматривать поправки порядка [math]\displaystyle{ \lambda/L }[/math] (слабо разреженные газы)[7]. В этом случае кинетика объясняет явления температурного скачка и течения газов вблизи твёрдых поверхностей[8].

Для сильно разреженных газов, когда [math]\displaystyle{ \lambda/L\gg 1 }[/math], гидродинамические уравнения и обычное уравнение теплопроводности уже не применимы и для исследования процессов переноса необходимо решать кинетическое уравнение с определёнными граничными условиями на поверхностях, ограничивающих газ. Эти условия выражаются через функцию распределения молекул, рассеянных из-за взаимодействия со стенкой. Рассеянный поток частиц может приходить в тепловое равновесие со стенкой, но в реальных случаях это не достигается. Для сильно разреженных газов роль коэффициента теплопроводности играют коэффициенты теплопередачи[9]. Например, количество тепла [math]\displaystyle{ Q }[/math], отнесённое к единице площади параллельных пластинок, между которыми находится разреженный газ, равно [math]\displaystyle{ Q=\varkappa(T_2-T_1)/L }[/math], где [math]\displaystyle{ T_1 }[/math] и [math]\displaystyle{ T_2 }[/math] — температуры пластинок, [math]\displaystyle{ L }[/math] — расстояние между ними, [math]\displaystyle{ \varkappa }[/math] — коэффициент теплопередачи.

Теория явлений переноса в плотных газах и жидкостях значительно сложнее, так как для описания неравновесного состояния уже недостаточно одночастичной функции распределения, а нужно учитывать функции распределения более высокого порядка. Частичные функции распределения удовлетворяют цепочке зацепляющихся уравнений (так называемых уравнений Боголюбова или цепочке ББГКИ, то есть уравнений Боголюбова — Борна — Грина — Кирквуда — Ивона). С помощью этих уравнений можно уточнить кинетическое уравнение для газов средней плотности и исследовать для них явления переноса.

Физическая кинетика двухкомпонентной плазмы описывается двумя функциями распределения (для электронов [math]\displaystyle{ f_e }[/math], для ионов [math]\displaystyle{ f_i }[/math]), удовлетворяющими системе двух кинетических уравнений (уравнений Власова). На частицы плазмы действуют силы

[math]\displaystyle{ F_e=-e\left(E+\frac{v\times B}{c}\right),\quad F_i=-Z_eF_e, }[/math]

где [math]\displaystyle{ Z_e }[/math] — заряд иона, [math]\displaystyle{ E }[/math] — напряжённость электрического поля, [math]\displaystyle{ B }[/math] — магнитная индукция, удовлетворяющие уравнениям Максвелла. Уравнения Максвелла содержат средние плотности тока [math]\displaystyle{ j }[/math] и заряда [math]\displaystyle{ \rho }[/math], определяемые с помощью функций распределения[10]:

[math]\displaystyle{ \textbf{j}=e\int \textbf{v}(Zf_i-f_e)\,d\textbf{p},\quad p=e\int (Zf_i-f_e)\,d\textbf{p}. }[/math]

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

Конденсированные среды

Физическая кинетика неравновесных процессов в диэлектриках основана на решении кинетического уравнения Больцмана для фононов решётки[11]. Взаимодействие между фононами вызвано ангармоническими членами гамильтониана решётки относительно смещения атомов из положения равновесия. При простейших столкновениях один фонон распадается на два или происходит слияние двух фононов в один, причём сумма их квазиимпульсов либо сохраняется (нормальные процессы столкновений), либо меняется на вектор обратной решётки (процессы переброса). Конечная теплопроводность возникает при учёте процессов переброса. При низких температурах, когда длина свободного пробега больше размеров образца [math]\displaystyle{ L }[/math], роль длины свободного пробега играет [math]\displaystyle{ L }[/math]. Кинетическое уравнение для фононов позволяет исследовать теплопроводность[12] и поглощение звука в диэлектриках[13]. Если длина свободного пробега для нормальных процессов значительно меньше длины свободного пробега для процессов переброса, то система фононов в кристалле при низких температурах подобна обычному газу. Нормальные столкновения устанавливают внутреннее равновесие в каждом элементе объёма газа, которьй может двигаться со скоростью [math]\displaystyle{ V }[/math], мало меняющейся на длине свободного пробега для нормальных столкновений. Поэтому можно построить уравнения гидродинамики фононного газа в диэлектрике[14].

Физическая кинетика металлов основана на решении кинетического уравнения для электронов, взаимодействующих с колебаниями кристаллической решётки. Электроны рассеиваются на колебаниях атомов решётки[15], примесях и дефектах, нарушающих её периодичность, причём возможны как нормальные столкновения, так и процессы переброса[16]. Электрическое сопротивление возникает в результате этих столкновений. Физическая кинетика объясняет термоэлектрические, гальваномагнически и термомагнинтные явления[17], аномальный скин-эффект[18], циклотронный резонанс в высокочастотных полях и другие кинетические эффекты в металлах. Для сверхпроводников она объясняет особенности их высокочастотного поведения.

Физическая кинетика магнитных явлений основана на решении кинетического уравнения для магнонов. Она позволяет вычислить динамическии восприимчивости магнитных систем в переменных полях, изучить кинетику процессов намагничивания.

Физическая кинетика явлений при прохождении быстрых частиц через вещество основана на решении системы кинетических уравнений для быстрых частиц и вторичных частиц, возникающих при столкновениях, например для [math]\displaystyle{ \gamma }[/math]-лучей (фотонов) с учётом различных процессов в среде (фотоэффекта, комптоновского рассеяния, образования пар). В этом случае кинетика позволяет вычислить коэффициенты поглощения и рассеяния быстрых частиц.

Фазовые переходы

Физическая кинетика фазовых переходов первого рода, то есть со скачком энтропии, связана с образованием и ростом зародышей новой фазы. Функция распределения зародышей по их размерам (если зародыши считать макроскопическими образованиями, а процесс роста — медленным) удовлетворяет уравнению Фоккера — Планка[19]:

[math]\displaystyle{ \frac{\partial f}{\partial t}=\frac{\partial}{\partial \alpha}\left(D\frac{\partial f}{\partial\alpha}-Af\right), }[/math]

где [math]\displaystyle{ \alpha }[/math] — радиус зародыша, [math]\displaystyle{ D }[/math] — «коэффициент диффузии зародышей по размерам», [math]\displaystyle{ A }[/math] — пропорционально минимальной работе, которую нужно затратить на создание зародыша данного размера. Кинетика фазовых переходов второго рода в наиболее простом приближении основана на уравнении релаксации параметра порядка [math]\displaystyle{ \eta }[/math], характеризующего степень упорядоченности, возникающей при фазовом переходе (уравнение Ландау — Халатникова)[20]:

[math]\displaystyle{ \frac{\partial\eta}{\partial t}=-\gamma\frac{\partial\Omega}{\partial\eta}, }[/math]

где [math]\displaystyle{ \gamma }[/math] — постоянный коэффициент, [math]\displaystyle{ \Omega }[/math] — термодинамический потенциал в переменных [math]\displaystyle{ T }[/math] и [math]\displaystyle{ \eta }[/math], вблизи точки фазового перехода зависящий от [math]\displaystyle{ \eta }[/math]. Для этой зависимости используется разложение по степеням [math]\displaystyle{ \eta }[/math] и [math]\displaystyle{ T-T_c }[/math], где [math]\displaystyle{ T_c }[/math] — температура фазового перехода.

Явления переноса в жидкостях

Теорию явлений переноса в жидкостях также можно отнести к физической кинетике. Xотя для жидкостей метод кинетических уравнений непригоден, для них возможен более общий подход, основанный на иерархии времён релаксации. Для жидкости время установления равновесия в макроскопически малых (но содержащих ещё большое число молекул) элементарных объёмах значительно меньше, чем время релаксации во всей системе, вследствие чего в малых элементах объёма приближённо устанавливается статистическое равновесие. Поэтому в качестве исходного приближения при решении уравнения Лиувилля можно принять локально равновесное распределение Гиббса с температурой [math]\displaystyle{ T(x,\;t) }[/math], химическим потенциалом [math]\displaystyle{ \mu(x,\;t) }[/math] и гидродинамической скоростью [math]\displaystyle{ V(x,\;t) }[/math], соответствующими рассматриваемой точке жидкости. Например, для однокомпонентной жидкости локально равновесная функция распределения (или матрица плотности) имеет вид

[math]\displaystyle{ f=\frac{1}{Z}\exp\left(-\int\beta(x,\;t)[\tilde{H}(x)-\mu(x,\;t)n(x)]\,dx\right), }[/math]

где

  • [math]\displaystyle{ \beta(x,\;t)=\frac{1}{kT(x,\;t)}, }[/math]
  • [math]\displaystyle{ \tilde{H}(x)= H(x)-p(x)V(x,\;t)+\frac{1}{2}mn(x)V^2(x,\;t) }[/math] — плотность энергии в системе координат, движущейся вместе с элементом жидкости,
  • [math]\displaystyle{ H(x) }[/math] — плотность энергии в неподвижной системе координат,
  • [math]\displaystyle{ p(x) }[/math] — плотность импульса,
  • [math]\displaystyle{ n(x) }[/math] — плотность числа частиц, рассматриваемые как фазовые функции, то есть функции от координат и импульсов всех частиц, например [math]\displaystyle{ n(x)=\sum_j^N\delta(x-x_j) }[/math].

Приближённое решение уравнения Лиувилля для состояний, близких к статистически равновесному, позволяет вывести уравнения теплопроводности и Навье — Стокса для жидкости и получить микроскопические выражения для кинетических коэффициентов теплопроводности и вязкости через пространственно-временные корреляционные функции плотностей потоков энергии и импульсов всех частиц системы. Этот же подход возможен и для смеси жидкостей. Подобное решение уравнения Лиувилля есть его частное решение, зависящее от времени лишь через параметры [math]\displaystyle{ \beta(x,\;t) }[/math], [math]\displaystyle{ \mu(x,\;t) }[/math], [math]\displaystyle{ V(x,\;t) }[/math], соответствующие сокращённому гидродинамическому описанию неравновесного состояния системы, которое справедливо, когда все гидродинамические параметры мало меняются на расстояниях порядка длины свободного пробега (для газов) или длины корреляций потоков энергии или импульса (для жидкостей).

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

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

См. также

Примечания

Литература