Материал: Неравновесные системы

Причина «странного» поведения ионов - это не учтённая нами непрерывная бомбардировка ионов электронами. Попробуем разобраться в существе дела. Выясним, какое влияние оказывают взаимные соударения электронов и ионов на движение электронов и ионов при наличии внешнего электрического поля.

Для наглядности будем считать ион массивным упругим шаром, на который налетает пучок лёгких частиц - электронов (масса иона в несколько тысяч раз больше массы электрона). Пусть до столкновения пучок частиц имеет среднюю направленную скорость u и соответственно ненулевой импульс. После столкновения пучок рассеивается равномерно во все стороны (рис. 1) (так же равномерно во все стороны зеркальный шар рассеивает направленный пучок света). То есть сразу после соударения все направления движения равновероятны, а значит, суммарный импульс электронов равен нулю.

Раз импульс пучка из-за соударений с ионами меняется, то это означает, что на электроны со стороны ионов действует сила. Если u - средняя скорость электронов в пучке, то mu - средний импульс электрона. Пусть время между столкновениями электрона с ионами - , тогда можно считать, что именно за это время теряется средний импульс mu каждого электрона. Так что средняя сила, отнесённая к одному электрону, равна

.

Потеря импульса электронами в целом из-за столкновений с ионами и описывается средней тормозящей силой f. Конечно, на отдельный электрон действует не обязательно такая сила, но для движения пучка в целом важна именно эта средняя сила.

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

При наличии внешнего электрического поля напряжённости Е «рассеянные» электроны «подхватываются» полем, которое вновь упорядочивает их движение, сообщая им некоторый направленный импульс, растрачиваемый при следующем соударении, и т.д. Таким образом, можно сказать, что устанавливается некоторая постоянная средняя скорость движения электронного пучка. Иначе говоря, сумма сил, действующих на электронный пучок, равна нулю. В расчёте на один электрон получаем

еЕ + f = 0.

Итак, при движении с установившейся средней скоростью, когда полная сила равна нулю, тормозящая сила со стороны ионов в расчёте на один электрон равна

f = - eE. (*)

Какие же силы действуют на ионы? Рассмотрим чистый жидкий металл, у которого все ионы одинаковы, причём каждый атом отдал в общее пользование по Z электронов проводимости. Электрическое поле действует на каждый ион с силой где  - заряд иона (заряд электрона взят по модулю, ибо ион положителен). Со стороны каждого электрона действует в среднем сила f₁ = - f (по третьему закону Ньютона). Если число ионов N, то на каждый ион со стороны одного электрона приходится сила . Тогда сила, с которой на один ион действуют все электроны (их число ZN), равна Zf₁. Но сумма сил

Чтобы в этом убедиться, подставим значение силы f₁, отличающейся от f только знаком, из выражения (*):

F= Z|e|E+ZeE=0.

Итак, силы, действующие на ионы в чистом металле со стороны электрического поля, компенсируются силами, действующими со стороны «электронного ветра» (движущихся электронов). Равенство полной силы нулю означает либо покой, либо равномерное движение. Но при течении жидкости в сосуде возникает вязкое трение (в конечном счёте из-за взаимодействия со стенками сосуда), оно и обеспечивает покой жидкости относительно сосуда в случае компенсации остальных сил, действующих на жидкость.

А если в металле имеется небольшая примесь «чужих» ионов? Пусть заряд чужака равен . Из-за того, что размеры примесного иона другие, изменится и число ударов электронов об него. Если площадь сечения иона примеси , а «своего» иона , то чужак подвергнется в  раз большему числу ударов, чем свой ион. Во столько же раз изменится и сила, действующая на примесный ион со стороны электронов.

Если , то ионы примеси будут двигаться в направлении этой силы, их скорость будет направлена по полю. Если , то ионы будут двигаться против поля.

Как же объяснить тот факт, что в чистой ртути ионы движутся по направлению электронного ветра? Мы ведь показали, что в чистом металле ионы должны быть неподвижны! Дело в том, что предполагалась полная одинаковость всех ионов, но это не совсем так. Хотя большинство ионов, действительно, находится в одном и том же «нормальном» энергетическом состоянии, некоторые ионы всегда имеют энергию больше нормальной. Такие ионы называют «активированными». Вероятность столкновений электронов с ионами увеличивается с ростом энергии ионов. Можно сказать, что сечения активированных ионов как бы увеличились по сравнению с нормальными. Таким образом, активированные ионы можно считать чужаками с тем же зарядом, что и нормальный ион, но с бóльшим сечением. Такие ионы увлекаются электронным ветром и сносятся к положительному полюсу источника. На этом основано разделение изотопов ртути в электрическом поле. У ионов разных изотопов заряд и сечение практически одинаковы, но активированные ионы лёгких изотопов сносятся электронным ветром быстрее, чем активированные ионы тяжёлых изотопов, из-за их меньшей массы. Это приводит к увеличению их концентрации у «плюса», то есть, если около положительного полюса отсасывать ртуть, то она будет обогащённой лёгким изотопом.

Впервые «парадоксальный» перенос ионов в амальгамах ртути был обнаружен в 1907 году. Тогда же было и введено представление о взаимном трении ионов и электронов. В последующих теориях, просуществовавших до 1959 года, был сделан шаг назад: движения ионов и электронов считались независимыми. Отчасти это объяснялось неубедительностью результатов экспериментов. В 1953 году был открыт эффект разделения изотопов ртути постоянным током. После этого появился целый ряд работ по исследованию электропереноса ионов. Первые теоретические работы, где был вскрыт физический механизм действия электронного ветра, относятся к началу 1959 года.

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

12. Уравнение Мастера. Диффузия тяжелого газа в легком

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

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

.(*)

Первое слагаемое в (*) характеризует прямые переходы , а второе обратные . Уравнение (*) называется основным кинетическим уравнением (уравнение Мастера). Основное кинетическое уравнение представляет собой уравнение баланса для плотности вероятности каждого состояния . Первый член в правой части уравнения (*) соответствует возрастанию вероятности из-за переходов из других состояний , а второй - уменьшению вероятности из-за переходов в другие состояния.

Для дискретных величин аналогом уравнения мастера является уравнение Паули:

.

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

.

Это уравнение еще называется уравнением детального баланса. Оно является следствием законов сохранения.

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

Вывод уравнения Фоккера-Планка из уравнения мастера.

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

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

Будем для определенности говорить именно об этом примере и выведем кинетическое уравнение, которому удовлетворяет в таком случае функция распределения частиц примеси по импульсам, f(t, p).

Обозначим через ω(р, q) d³q отнесенную к единице времени вероятность изменения импульса р → р-q тяжелой частицы при элементарном акте - ее столкновении с легкой частицей. Тогда кинетическое уравнение для функции f(t, p) запишется в виде

,

где справа стоит разность между числом частиц, поступающих (в 1 с) в заданный элемент импульсного пространства d³p и покидающих его за то же время. Согласно сделанным предположениям, функция ω(р, q) быстро убывает с увеличением q, так что основную роль в интеграле играют значения q, малые по сравнению со средним импульсом частиц. Это обстоятельство позволяет произвести в подынтегральном выражении разложение

,

, .

Выражение в правой части (21.2) имеет вид дивергенции в импульсном пространстве, , от вектора

.

Если коэффициент постоянен или почти постоянен, то это откуда-то должно следовать! Превращение из тензора в скаляр тоже не так очевидно, поскольку не совсем ясно, почему при  эта величина должна обращаться в ноль.

Другими словами, полученное уравнение имеет, как и следовало, вид уравнения непрерывности в импульсном пространстве; тем самым автоматически соблюдается сохранение числа частиц при процессе. Вектор же S является плотностью потока частиц в импульсном пространстве.

Согласно приведенным формулам коэффициенты в кинетическом уравнении выражаются через средние характеристики столкновений, и в этом смысле их вычисление представляет собой механическую задачу. Фактически, однако, нет необходимости в раздельном вычислении коэффициентов A_α и B_αβ; они могут быть выражены друг через друга из условия обращения потока в нуль в статистическом равновесии. В данном случае равновесная функция распределения есть

,

где М - масса частиц тяжелого газа, а Т - температура основного (легкого) газа. Подстановка этого выражения в уравнение S = 0 дает

.

Таким образом, кинетическое уравнение принимает вид

. (*)

Единственный вектор, от которого могут зависеть коэффициенты B_αβ, - импульс тяжелых частиц р. Но если скорости этих частиц, р/М, в среднем малы по сравнению со скоростями легких частиц (при чем тут скорости?), то при столкновениях их можно считать неподвижными; в этом приближении величины В_αβ не будут зависеть от р. Другими словами, тензор В_αβ сведется к постоянному скаляру В:

,

. (**)

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

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

. Причины необратимости в макросистемах

Рассмотрим вопрос о необратимом характере процессов релаксации в макроскопических системах. Как известно, законы, определяющие эволюцию микроскопической системы во времени, инвариантны по отношению к изменению знака времени как в классической, так и в квантовой механике. Это значит, что все микропроцессы принципиально обратимы, т.е. при надлежащем изменении условий (например, при изменении направлений начальных скоростей) любой микропроцесс может протекать, как в прямом, так и в обратном направлении, проходя через одни и те же промежуточные состояния. Возникает вопрос: каким образом статистическая физика, основанная на обратимых во времени законах микропроцессов, может приводить к необратимым законам макроскопических процессов, в частности, к описанию процессов релаксации и к закону возрастания энтропии в замкнутых системах?