2 |
Второй классический труд Гиббса посвящен основам статистической механики. Несмотря на то, что развитие этой науки в дальнейшем пошло по пути, указанному Гиббсом, многое из
этого труда еще ждет своего дальнейшего использования и развития.
Пожалуй, наиболее поучительным примером является теория образования новой фазы. Это — по существу проблема кинетики. В то же время ее решение^очетает как термодинамические, так и молекулярно-статистические расчеты. Первые нужны, как показал Гиббс, для вычисления работы образования критического зародыша — понятия, введенного им же [113. Вторые — для расчета вероятности достижения и перехода, через активационный барьер, отвечающий критическому зародышу. Основную трудность представляет расчет этой вероятности. Так, в случае конденсации, например, полный статистический расчет процессов дорастания молекулярных комплексов вплоть до критического зародыша представляет в общем случае невероятно сложную задачу как в силу математических трудностей, так и необходимости знания многочисленных и трудноопределимых констант, характеризующих различные стадии процессов агрегации и дезагрегации молекулярных комплексов. По сути процессы роста агрегатов представляют собой сложно разветвленную цепную реакцию.
Еще труднее провести строгий расчет процесса зарождения новой фазы внутри жидкости, например, при ее вскипании. Начальная стадия этого процесса даже в качественном отношении служит, скорее, предметом догадок. Во всяком случае, она определяется коллективными взаимодействиями, расчет которых представляет одну из сложнейших задач статистической физики. Поэтому разумно, прежде чем пытаться давать несомненно мало точное решение задачи для общего случая, вначале рассмотреть более доступный строгому решению частный случай.
Существенное облегчение задачи расчета вероятности ну — клеации возникает в предельном случае малых пересыщений. Ему соответствует критический зародыш существенно макроскопических размеров. В силу этого определение его размеров и работы образования может быть проведено наиболее строго на основе термодинамики Гиббса. Второе упрощение заключается в возможности в этом случае использования метода перевала. Суть последнего, как известно, заключается в том, что, беря дифференциальное уравнение роста докритических зародышей в форме Крамерса [12] или Зельдовича [13], применяют его только к области размеров, близких к размеру критического зародыша. Это может быть оправдано, только если размер критического зародыша велик.
Здесь, однако, возникают две трудности. Во-первых, при интегрировании уравнения, например Крамерса для перевадь — ной области, появляется неопределенный сомножитель, определимый однако, как показал Френкель [14, гл. 7 3, для разреженных паров. Для более общих случаев, в особенности при нуклеа — ции в жидкой фазе, имеющиеся попытки определения этого множителя [13, 15] как было показано [в 16, 17], являются необоснованными и ведут к ошибочным результатам. В то же время указаны случаи [18], когда важнее знать этот множитель, чем экспоненциальный.
Во-вторых, рассматривают, как правило, стационарное решение кинетических уравнений роста зародышей и их перехода через барьер. Для его получения налагают предельное условие, при котором концентрации «закритических» зародышей равны нулю. Это означает, что выросшие закритические зародыши, несмотря на гарантированную возможность их дальнейшего роста, изымаются «из игры», а их вещество превращается в исходную фазу. Это поддерживает пересыщение последней постоянным. Одновременно при интегрировании кинетического уравнения рассматривается равновесная функция распределения зародышей по размерам. Законность такого рассмотрения трудно оправдать, в особенности поскольку рассмотрение относится, как правило, к среде неограниченного объема.
Все эти трудности и нестрогости, однако, исключаются, если применить к задаче нуклеации построение большого ансамбля Гиббса [16—18]. Для этого достаточно рассматривать нуклеа — цию в ограниченном объеме V, отделенном от остальной части системы заданного объема W и энергии Е перегородкой, допускающей обмен молекулами, но не «околокритическими» зародышами. Ввиду ограниченности объема V, в котором рассматривается возможность образования новой фазы, строя большой ансамбль Гиббса для объема, мы можем рассматривать во всей строгости равновесное статистическое распределение состояний ансамбля, отличающихся, в частности, по размерам зародышей новой фазы. При этом считается возможным возврат, хотя бы и по прошествии чрезвычайно большого времени, за счет маловероятной флуктуации от состояния новой фазы к первоначальному состоянию.
Для нахождения вероятности прямого перехода от начальной, пересыщенной фазы через критический зародыш достаточно разделить поток зародышей через перевал, выраженный на основе подсчета состояний большого ансамбля Гиббса, на число состояний этого ансамбля для исходного метастабильного состояния. При таком подсчете происходит сокращение аналогичного числа состояний в числителе и получается классическое выражение вероятности нуклеации через работу образования критического зародыша. Однако, в отличие от ранее примененных методов, не использующих большой ансамбль Гиббса, полученная формула содержит строго и однозначно определенный пред — экспоненциальный множитель. Для его вывода вовсе не требуется знание и учет начальных стадий роста зародышей. Достаточно рассмотреть его кинетику в околокритической области.
В особенности незаменим разработанный метод для расчета вероятности нуклеации в жидкости, например вероятности вскипания жидкости произвольной летучести и вязкости [1]. В этом последнем случае задача осложняется тем, что флуктуа- ционный рост зародыша необходимо характеризовать двумя параметрами — размером и числом молекул пара в нем или же объемом и давлением.
Не останавливаясь на способе и результатах решения этой задачи, опубликованной в [18], ограничимся следующим. Применение большого ансамбля Гиббса в сочетании с обобщенным на два измерения уравнением Крамерса—Зельдовича позволяет для случая, когда критический зародыш велик, получить строгую формулу для вычисления v — вероятности образования за критического пузырька, не рассматривая начальную, микроскопическую стадию роста его. В общую формулу входит безразмерный параметр
To = pcRJvПР
Где рс и Rc — давление и радиус критического зародыша; V — тепловая скорость молекул; т] — вязкость жидкости; Р — коэффициент конденсации.
Если (а > 1, то при 2O/PcRc > 3, где а — поверхностное натяжение жидкости, формула принимает вид:
А
, = (15)
2л
Здесь V" — объем среды;
Р7 — плотность жидкости;
[зс — работа образования критического зародыша.
При 2O/PcRc <3
Р/ V‘Py т/ ° ~ТГ V = — — У WE (16)
Г»
Где р£ — плотность пара в критическом зародыше.
Эта формула отличается от формулы (36) работы [15] на
Множитель р/р£.
При ю < 1 общая формула принимает вид:
■Ф
V.J^^STW— ,17,
Формула (17) совпадает с формулой, выведенной ранее 117], и с формулой (15).