Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

WATER VAPOR CONDENSATION. CLASSICAL AND QUANTUM APPROACH

Khlopkov A.Y. 1
1 Moscow Institute of Physics and Technology (State University)
1063 KB
Formation process of condensation’s initial nucleuses is one of the most difficult and incompletely determinate problems in new phase formation theory. These nucleuses are able to have a small size and consist only of some molecules, called “clusters” or micro-cluster. Micro-cluster is the smallest units included from two to hundred and more atoms. Modern quantum-chemical methods let to find an inertia moment, frequencies of intermolecular and intermolecular oscillations in the clusters and their dissociation energy. On the basis of these calculations it’s possible to calculate all the thermodynamic and also kinetic characteristics of the clusters. The clustered method underlies in thermodynamic of natural gas components. In this paper compared the two theories in order to find the nucleation rate of water- the classical and quantum mechanical.
clusters
nucleation
condensation
thermodynamics of natural gas
classical and quantum views
molecular collisions
1. Artyukhin A.S., Egorov B.V., Zababurin E.A., Koshcheev A.V., Markachev Yu.E., Plekhanov E.A., Ufimtsev I.S., Khlopkov A.Yu., Chugreev A.L. Khimicheskaya fizika. 2004. vol. 23, no. 4. pр. 28–46.
2. Gordiets B.F., Shelepin L.A., Shmotkin Yu.S. Trudy FIAN. 1984. vol. 145. pр. 189–219.
3. Poltorak O.N. Termodinamika v fizicheskoi khimii. Moscow, Vysshaya shkola, 1991.
4. Robinson P., Kholbruk K.K. Monomolekulyarnye reaktsii. Moscow, Mir, 1975.
5. David J. Di Mattio. In partial Fulfillment of the Requirements for the Degree, 1999.
6. Tompson S.M., Gubbins K.E., Walton J.P.R.B., Chantry R.A.R., Rowlinson J.S. Journal of Chem. Phys. 1984. vol. 81. no. 1. pр. 530–542.
7. Volmer M. Kinetik der Phasenbildung. Dresden: Steinkopff, 1939.

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

Вычисление скорости зародышеобразования кластеров воды по классической теории

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

hlopkov01.wmf (1)

где n1 – концентрация мономеров (предполагается, что пар мономеров не истощается в процессе конденсации); βi – скорость присоединения молекулы к кластеру; αi – скорость испарения; g* – критический размер зародыша. Используя далее капельную модель, т.е. предполагая, что кластеры сферические и однородные, получим для скорости нуклеации выражение

hlopkov02.wmf

где a – коэффициент прилипания; v – объем одной молекулы жидкости; γ – поверхностное натяжение; hlopkov03.wmf – пересыщение; B – поправка Варнарда, которая возникает при последовательном кинетическом подходе и которая отсутствует при расчете из классических термодинамических соображений, т.е. при вычислении работы, необходимой для создания зародыша размера g. Поверхностное натяжение для кластеров воды малых размеров было вычислено с большой точностью в работе [1] (рис. 1). Данные результаты были сделаны по расчету размера кластера по точным координатам атомов в кластере. Отсюда видно, что реальное поверхностное натяжение кластеров малых размеров на 30–50 % меньше, чем для обычных капель жидкости. g* – в этом случае для рассматриваемых температур T = 230–250 K получается порядка 3, в то время как по классической теории g* ~ 5, для данных температур.

pic_98.tif

Рис. 1. Зависимость поверхностного натяжения от размера кластера

Выражение (1) переходит в

hlopkov04.wmf,

где

hlopkov05.wmf

Для функции f(g) была взята зависимость, полученная методом молекулярной динамики в эксперименте [6]. В классической теории нуклеации скорость испарения кластера, содержащего g молекул,

hlopkov06.wmf

где hlopkov07.wmf – площадь поверхности кластера, состоящего из g – молекул; v – объем одной молекулы (мономера), и тогда

hlopkov08.wmf – изменение поверхности за счет испарения одной молекулы.

hlopkov09.wmf – скорость испарения с единицы площади плоской границы раздела; E – теплота испарения с плоскости.

В классической теории нуклеации скорость роста зародышей βg выражают через число столкновений с мономерами в единицу времени:

hlopkov10.wmf

где a – коэффициент прилипания, величина неопределенная и может лежать от 0,01 до 1.

На основе вышеприведенных данных была вычислена скорость зародышеобразования для жидкокапельной модели с поправкой Варнарда (2) (рис. 2). Полученные данные были сравнены с последними экспериментальными данными работы [5]. Как видно, классический подход к теории зародышеобразования не дает удовлетворительных результатов. Далее будет использован и продемонстрирован квантово-механический подход для расчета скорости зародышеобразования.

Вычисление скорости зародышеобразования кластеров воды по квантово-механической теории

Рассмотрим систему, состоящую из газофазных молекул А1 конденсирующегося вещества (мономеров) и молекулярных агрегатов An (кластеров), состоящих из двух и более молекул: n ≥ 2. Кластер An размерности n обладает энергией диссоциации (разрыва связи) En, тремя моментами инерции IA(n), IB(n), IC(n), а также 3nl – 6 колебательными степенями свободы, из которых (3l – 6)n относятся к внутримолекулярным, а (6n – 6) – к межмолекулярным колебаниям (l – число атомов в мономере).

Рассмотрим сначала упрощенную модель кластеризации, предполагая, что рост кластеров происходит за счет присоединения к ним одиночного мономера молекулы-ассоциации, а их разрушение происходит через промежуточную стадию – образование и распад активированного комплекса [2]. Опуская стадию активированного комплекса, учитывая только конечный продукт реакции, цепочка элементарных процессов образования и гибели кластеров An может быть записана в суммарном виде:

Ka(n), An + A1 + M ó An + 1 + M, Kd(n + 1). (2)

pic_99.tif

Рис. 2. Скорость зародышеобразования по классической теории (жидкокапельная модель)

Здесь M – любая молекула или кластер газовой фазы; Ka(n) и Kd(n + 1) – константы скоростей прямых и обратных реакций; Ka(n) – константа скорости присоединения мономера к кластеру An, а Kd(n + 1) – константа скорости диссоциации кластера An + 1 на кластер An и мономер.

Для концентрации Nn кластеров вида An на основе (2) можно записать уравнение изменения их со временем в виде [2]:

hlopkov11.wmf (3)

где hlopkov12.wmf

Здесь Jn – поток кластеров в пространстве числа мономеров в них, т.е. число кластеров, переходящих в единицу времени из точки n в точку n + 1.

В условиях термодинамического равновесия скорость прямых и обратных реакций Ka(n), Kd(n + 1) связаны через константу равновесия:

hlopkov13.wmf

Здесь Pi = aiP, hlopkov14.wmf – парциальное давление и мольная доля компонента i; hlopkov15.wmf – константа равновесия реакции (2), выраженная в терминах парциального давления компонентов ([Kp] = атм–1). Расчет константы скорости диссоциации кластера основан на теории мономолекулярных реакций РРКМ [4], здесь приведен только окончательный результат:

hlopkov16.wmf

где E* = E+ + E0;

hlopkov17.wmf

здесь L – число путей, эквивалентных в смысле энергетики (например, L = 2 для диссоциации H2O на H и OH); Q1 и hlopkov18.wmf – статистические суммы для вращений в молекуле и активированном комплексе соответственно (в расчетах принималось hlopkov19.wmf); Q2 – полная (колебательная и вращательная) статистическая сумма для диссоциирующей молекулы; hlopkov20.wmf – количество колебательно-вращательных состояний активированного комплекса в интервале энергий [0, Е+]; N* – плотность колебательно-вращательных уровней реагента. Число состояний активизированных комплексов и плотность энергии молекул реагента малых размеров подсчитывались путем прямого пересчета (для этого все частоты нормальных колебаний были разбиты на 9 групп). Для больших кластеров использовалось аналитическое приближение Виттена – Рабиновича [4]. Константа образования кластера Ка(n) определялась через константу диссоциации Кd(n + 1) и константу равновесия Kn.. В квазиравновесных условиях dNn/dt = 0, следовательно, J не зависит от n, так как согласно (3) hlopkov21.wmf для любого n. В этом случае

hlopkov22.wmf

В равновесных условиях, когда скорости любых прямых и обратных процессов уравновешиваются

hlopkov23.wmf

Здесь hlopkov24.wmf – квазиравновесная функция распределения в отсутствии потока (J0 = 0) и решение hlopkov25.wmf hlopkov26.wmf – пересыщение. Здесь hlopkov27.wmf – числовая плотность мономеров; N1eq – равновесная числовая плотность мономеров над кластерами предельно большого размера.

hlopkov28.wmf

hlopkov29.wmf

где Edis(∞) – энергия отрыва мономера от кластера бесконечно большого размера, «среднестатистическая» характеристическая колебательная температура кластера размера n – hlopkov30.wmf определяется из следующего соотношения:

hlopkov31.wmf

Квазиравновесное распределение принимает вид

hlopkov32.wmf

Здесь hlopkov33.wmf Edis(∞) – энергия отрыва мономера от кластера бесконечно большого размера; Edis(2) – энергия диссоциации димера; α ∈ (0, 1) – параметр [2]. В работе [2] было получено выражение для квазиравновесного распределения кластеров, содержащих n мономеров вида

hlopkov34.wmf

полностью совпадающее с выражением, полученным в приближении жидко-капельной модели классической теории конденсации, если

hlopkov35.wmf

и α = 2/3, где rn, σ – радиус и поверхностное натяжение.

Выражение (2) отличается от классического и выражения работы [2] явным учетом структуры и энергетических характеристик межмолекулярных квантовых движений в кластере размерности n.

Статистическая сумма идеального газа с внутренними степенями свободы, состоящего из N неразличимых идентичных частиц – кластеров размерности n, равна [3]:

hlopkov36.wmf

Таким образом, для одного моля кластеров размерности n

hlopkov37.wmf (4)

где NA – число Авогадро; M1 – масса мономера; ge – вырождение электронного уровня кластера; S1 – ядерный спин мономера.

Знание Zn позволяет вычислить все необходимые термодинамические функции для одного моля кластеров размерности n [3]. В квазиравновесном приближении приращение свободной энергии Гиббса может быть приведено к виду

hlopkov38.wmf

pic_100.tif

Рис. 3. Скорость зародышеобразования по квантово-механической теории для различных T (217, 239, 244, K)

Решение неоднородной системы в квазистационарном приближении ищется в виде [7]:

hlopkov39.wmf

где hlopkov40.wmf – квазиравновесное распределение кластеров по размерам при (J0 = 0).

Решение неоднородной системы в случае (J0 ≠ 0) имеет вид

hlopkov41.wmf

где hlopkov42.wmf

При известных константах ассоциации и диссоциации кластеров находится квазиравновесное решение hlopkov43.wmf и вычисляется число зародышей в единице объема, переходящих в единицу времени из точки n (кластера размерности n) в точку n + 1 (кластера размерности n + 1) для системы с плотностью мономеров N1.

Числовая плотность кластеров в стационарном состоянии меньше соответствующей числовой плотности в квазиравновесном состоянии. Уменьшение числовой плотности кластеров в стационарном состоянии по сравнению с квазиравновесным тем больше, чем больше его размер n.

Для подсчета J0 была использована компьютерная программа работы [1]. Как видно из рис. 3, при использовании данного квантово-механического подхода получается результат близкий к экспериментальному, что подтверждает необходимость использования данного подхода. Далее будет продемонстрирован квантово-механический подход для вычисления термодинамических функций кластеров воды.

Выводы

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

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

Работа выполнена при поддержке РФФИ (Грант № 14-07-00564-а).

Рецензенты:

Николаев В.С., д.т.н., профессор, кафедра математики, Московский физико-технический институт (государственный университет), г. Жуковский;

Чернышев С.Л., д.ф.-м.н., профессор, ФГУП «Центральный аэрогидродинамический институт им. профессора Жуковского», г. Жуковский.