Поверхность потенциальной энергии и термодинамические функции молекулярного комплекса H2O·O3 Текст научной статьи по специальности «Физика»

УДК 544.13

ПОВЕРХНОСТЬ ПОТЕНЦИАЛЬНОЙ ЭНЕРГИИ И ТЕРМОДИНАМИЧЕСКИЕ ФУНКЦИИ МОЛЕКУЛЯРНОГО КОМПЛЕКСА H2OO3

С.К. ИГНАТОВ* , А.Г. РАЗУВАЕВ*, И.Н. АНДРОНОВ**

*Нижегородский государственный университет им. Н.И. Лобачевского, г. Нижний Новгород

**Ухтинский государственный технический университет, г. Ухта ignatov@ichem.unn.ru

Квантовохимическим методом исследована поверхность потенциальной энергии (ППЭ) слабосвязанного бинарного молекулярного комплекса озона и воды H2OO3, включающая 315 900 точек, среди которых 14 016 уникальных по симметрии. На основе оптимизации локальных минимумов, найденных на дискретной ППЭ, определена структура глобального минимума, вычислены его энергия связывания и колебательные частоты. С использованием приближения, учитывающего ангармонизм ППЭ, определены термодинамические функции данного комплекса в состоянии идеального газа, оценены его константа равновесия и ее температурная зависимость.

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

S.K. IGNATOV, A.G. RAZUVAEV, I.N. ANDRONOV. SURFACE OF POTENTIAL ENERGY AND THERMODYNAMIC FUNCTIONS OF A MOLECULAR COMPLEX H2OO3

The surface of potential energy of loosely coupled binary molecular complex of ozone and water H2O-O3, including 315 900 points, among which 14 016 are unique in symmetry, is investigated by means of quantum-chemical method. On the basis of optimization of the local minima found on a discrete surface of potential energy, the structure of a global minimum is defined, its energy of linkage and oscillatory frequencies are calculated. With use of the approach considering anharmonicity of surfaces of potential energy, thermodynamic functions of the given complex in condition of ideal gas are defined, its constant of balance and its temperature dependence is estimated.

Key words: loosely coupled molecular complexes, nonrigid molecules, ozone, water, thermodynamic functions

Комплексы и кластеры воды с кислородсодержащими соединениями и их влияние на различные физические и химические процессы, протекающие в газовой фазе, вызывают в настоящее время особый интерес в связи с проблемами атмосферной химии, климатических изменений и обеднения озонового слоя [1, 2]. Среди сложных физико-химических процессов с участием следовых газов атмосферы, протекающих в тропосфере (перистые облака) и стратосфере (полярные стратосферные и серебристые облака), а также на поверхности снега и льда в полярных областях Земли, особенно важно образование комплексов и кластеров воды с активными кислородсодержащими молекулами и, в частности, с озоном. Физическим аспектам этих процессов были посвящены обзоры [3, 4].

Вода во всех ее состояниях играет главную роль в поглощении солнечной и земной радиации и является важным парниковым газом [5]. Поскольку образование водных комплексов и кластеров вызывает сдвиг и уширение спектральных полос моно-

мера в спектрах поглощения, а также ведет к появлению новых полос поглощения и повышению интенсивности запрещенных электронных переходов, водные комплексы способны изменять радиационный баланс атмосферы [6]. Это может привести, например, к искажению результатов при определении концентрации атмосферных примесей лазерными методами. Кластеры воды, участвующие в фотохимических процессах, могут также в значительной степени изменять скорости диссоциации мономера и приводить к совершенно новым каналам фотодиссоциации. Некоторые аспекты роли малых водных кластеров (в том числе, комплексов 1:1) в атмосферных процессах были недавно рассмотрены в работах [7, 8].

Химические свойства озона представляют большой интерес для многих областей химии, однако существует только незначительное число исследований, посвященных изучению его комплексов с малыми молекулами, в том числе с Н20. Результаты исследований структуры комплекса Н2ОО3,

полученные методами микроволновой спектроскопии [9] и низкотемпературной матричной изоляции

[10], были опубликованы практически одновременно. Найдено [9], что структура этого комплекса в газовой фазе (в импульсной сверхзвуковой струе) имеет симметрию Cs: молекула H2O и центральный атом кислорода молекулы озона расположены в плоскости симметрии комплекса с атомом водорода воды, направленным к концевому атому кислорода озона. Такая структура характеризуется противоположными направлениями дипольных моментов двух мономеров, на основании чего можно заключить, что диполь-дипольное электростатическое взаимодействие дает главный вклад в энергию связи. Этот результат был подтвержден квантово-химическим расчетом на уровне MP2/6-31G(d,p). Симметрия Cs этой структуры определяется куло-новским взаимодействием между положительно заряженным водородным атомом и двумя отрицательно-заряженными атомами кислорода молекулы озона. Энергия связи, оцененная на уровне MP2/6-31G(d), составляет 14.9 кДж моль-1 (AECP = 2.9 кДж моль-1), неполная оптимизация на уровне MP4/6-31G(d) (с геометрией мономера, взятой из расчета MP2/6-31G(d)) дает AE=10.0 кДж моль-1, AECP=4.2 кДж моль-1).

Иная структура (открытая конформация OOO-...H-OH с большим расстоянием O-O) была найдена впервые [11] на основе оптимизации MP4/4-31G(d,p). Хотя уровень теории, очевидно, не позволяет сделать надежное заключение о тонких энергетических различиях, авторы привлекли внимание к важному факту, что самая распространенная структура в атмосфере не та, у которой наиболее низкая энергия связи, а структура, имеющая максимально низкую свободную энергию Гиббса при данной температуре. Если различие энергии между конформациями - мало, то термодинамически благоприятнее конформация с большим расстоянием O.O.

Оптимизация геометрии комплекса 1:1 между водой и озоном на уровне QCISD/6-311++G(d,p) была выполнена в работе [12]. Расчет показывает, что три структуры abc - устойчивые формы комплекса H2OO3. Наиболее выгодная структура H2OO3 a) имеет симметрию Cs, причем центральный кислород O3 и всех атомов H2O расположен в плоскости симметрии (что близко к структуре, найденной в работе [9]). Диполь H2O ориентирован по направлению к центральному атому кислорода O3. Другие две конформации - cis и trans формы комплекса - сходны со структурами, рассмотренными в

[11]. В этих конформациях все атомы расположены в плоскости молекулы, и водород H2O связан с одним из концевых атомов кислорода O3 водородной связью. Энергии связи O3 с H2O для дипольной, cis-, и trans-форм составляют 13.5, 8.4, и 7.9 кДж моль-1, соответственно, на уровне QCISD/6-311++G(d,p) и 10.0, 9.5, 9.6 кДж моль-1 на уровне QCISD(T)/6-311 ++G(3df,3pd)//QCISD/6-311 ++G(d,p). Таким образом, диполь-дипольная форма энергетически более выгодна, чем cis- и trans- конформации, которые по стабильности очень близки. К сожалению, авторы

Тачикава и Абэ не привели величину BSSE для данного расчета, которая может значительно понизить указанную энергию комплекса. Вращательные константы для дипольной формы составляют А = 11.897, В = 4.177 и С = 3.318 ГГц, что находится в хорошем согласии с экспериментальными значениями (А = 11.961, В = 4.174 и С = 3.265 ГГц). С другой стороны, вращательные константы для ав-и trans-форм гораздо хуже согласуются с экспериментом. Эти результаты показывают, что диполь-дипольная конформация комплекса - наиболее вероятный кандидат на структуру Н2003, наблюдаемую экспериментально в газовой фазе [12].

Несмотря на успехи в установлении структуры комплекса озона и воды, его термодинамические функции остаются практически неизученными. Немногочисленные работы, в которых сообщалось о расчете колебательных частот и термодинамических параметрах, основаны на приближении «жесткий ротатор - гармонический осциллятор» (ЖРГО), которое, очевидно, не применимо к такой нежесткой и существенно ангармонической системе, как слабосвязанный комплекс Н2О О3. Так, концентрационный профиль комплекса Н2ОО3 в атмосфере в зависимости от высоты был рассчитан в работе [7].

Таким образом, достигнутые в последнее время результаты исследования структуры, энергетики и термодинамических функций бинарных комплексов Н2003, образующихся в газовой фазе при низких температурах и давлениях, нельзя считать исчерпывающими. В настоящей работе мы предпринимаем исследование структуры, энергии образования бинарных комплексов Н20 03 и полной межмолекулярной поверхности потенциальной энергии (ППЭ) квантовохимическим методом МР2, а также оцениваем термодинамические функции на основе нового оригинального метода расчета термодинамических функций, базирующегося на более точных приближениях, чем использованное ранее гармоническое приближение.

Теория и метод расчетов

Для описания термодинамических функций молекулярного комплекса Н2003 в данной работе используется приближение классических ангармонических ровибрационных движений / квантовых колебаний (КЛАРВ/КГК), специально разработанное для описания термодинамики слабосвязанных комплексов. В этом приближении бинарный комплекс АВ рассматривается как система из двух жестких нелинейных молекул (мономеров) А и В. Набор из шести координат (а,р^,Я,6,ф) позволяет однозначно описать положение одного мономера относительно другого. Все возможные движения атомов разделяются на межмолекулярные движения мономеров и внутримолекулярные колебания. Межмолекулярные ровибрационные движения, спаривание между которыми не позволяет выделить отдельные колебания или вращения, представлены классическими уравнениями движения. Внутримолекулярные движения полностью независимы друг от друга и от межмолекулярных движений и описываются уравнениями движения квантового гармо-

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

^ = 1

J * (A)

2 3

-z

J * (B)

2 3

P2

к=1 21 к к=1 21В к=1 21 к 2тАВ где 1А,1В, I0 - мгновенные моменты инерции мономеров А, В и комплекса в целом, тАВ - приведенная масса комплекса. Потенциальная энергия есть ангармоническая функция шести межмолекулярных координат и=Ща,р,у,Я,6,ф), определяемая путем сканирования полной шестимерной ППЭ на дискретной сетке. В узлах этой сетки энергия комплекса определяется квантовохимическим методом, в межузельных положениях - линейной интерполяцией энергии ближайших узлов.

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

2 = 2 2 2 2

0 е\ес 1гап угЪ гоу1Ь '

Статсуммы 2е1ес, 2^, 2уЛ определяются по стандартным формулам статистической термодинамики, аналогичным приближению ЖРГО. Стат-сумма ровибрационнго движения мономеров вычисляется путем суммирования по всем возможным энергиям фазового объема системы П(Е) :

z=i

exp(-E IdE.

dE ч kT '

Классический фазовый объем при данной энергии определяется как интеграл всех возможных фазовых состояний:

0.(Е) = | dpdq

Е< Н (р,д)< Е+ йЕ

и рассчитывается методом интегрирования Монте-Карло с простой случайной выборкой. После того, как статсумма комплекса 2 рассчитана, термодинамические функции комплекса в состоянии идеального газа при данной температуре и давлении определяются как:

F = -k ln Z

S = R ln Z + RT

d (ln Z) dt

U = F + TS

H = U + RT (для трансляционных степеней свободы) H = U (для остальных)

G = H - TS .

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

AbG0 (T) = AbE + G0 (T; complex) - Z G° (T; monomer)

К0(Т) = ехр(-Д^0(Т)/ ЯТ).

Здесь ДЬЕ - электронно-ядерная энергия

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

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

2„

7 = F 7 , F =

corr 1 Q classic * Q

J quantHO

Z„

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

Описанный подход реализован в виде оригинальной компьютерной программы ADANIMEHS и применен к исследованию ППЭ комплекса Н2003 с целью определения его термодинамических параметров. Геометрические параметры, описывающие ППЭ:

1. Вращение Н20 на угол Эйлера а (0° < а < 360°; да = 45°)

2. Вращение Н20 на угол Эйлера р (0° < р < 180°; др = 45°)

3. Вращение Н20 на угол Эйлера Y (0° < Y < 360°; dY = 45°)

4. Трансляция Н20 по оси OZ на [ (2.75 А < [ < 5.50 А; = 0.25 А)

5. Вращение 03 вокруг оси OY на угол 6 (0° < 8 < 180°; д6 = 45°)

6. Вращение 03 вокруг OZ на угол ф (0° < ф < 360°; дф = 30°).

Дискретная ППЭ, соответствующая данной процедуре, состоит из 315 900 точек, среди которых 14 016 уникальных. Энергии уникальных точек ППЭ рассчитывались методом MP2/6-311++G(2d,2p). Для расчета энергии использовалась программа РС GAMESS [15], работающая под управлением программы ADANIMEHS. Расчет проводился в полностью автоматическом режиме, с автоматической генерацией исходных файлов и считыванием результатов. Для расчетов использовались четыре компьютера Репйит4/3 ГГц и два компьютера Реп-йит3/800 Мгц. Время расчета на компьютере Реп-йит4 составляло около 21 час (в среднем 64 с на уникальную точку, 2.28 с на 1 точку ППЭ). Расчет полной ППЭ продолжался около семи дней.

Для определения геометрии мономеров проводилась оптимизация структур Н20 и 03 на том же уровне расчета, что и расчет комплекса (МР2/6-311++G(2d,2p)). Геометрические параметры мономера Н20 составляют: г(0Н)=0.9583 А, а(н-О-Н) = 104.3°, мономера 03: г(О-О) = 1.2899 А, а(0-0-0) =

116.8°. Полные энергии мономеров: Е^(Н20) = -76.2988008673 а.е., Е^(03) = -225.0549293698 а.е.

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

Результаты и обсуждение

На дискретной ППЭ комплекса Н2003, рассчитанной методом МР2/ 6-311++G(2d,2p), обнаруживается 110 локальных минимумов (включая точки глобального минимума), которые классифицируются по семи типам. Их энергии и геометрические параметры приведены в табл. 1, структурные параметры показаны на рис.1. Глобальный минимум (ГМ) соответствует асимметричной структуре, близкой к той, что была найдена в [12]. Глобальный минимум соответствует межмолекулярному расстоянию R = 3.00 А, а его структура показана на рис. 2. Полная энергия глобаль-

ного минимума составляет -301.3578124128 а.е., что соответствует энергии связывания 10.72 кДж моль-1.

Полностью оптимизированная (МР2/6-311++ G(2d,2p)) структура ГМ сходна со структурой глобального минимума ППЭ (рис.2).

Для уточнения энергии комплекса были оценены величины BSSE и нулевой энергии колебаний для всех найденных локальных минимумов. Оценка BSSE проводилась стандартным методом Бойса.

Таблица 1

Структурные и энергетические параметры оптимизированных структур, соответствующих локальным минимумам (ЛМ) ППЭ

системы H2O•O3.

Локальный минимум ППЭ Группа симметрии V т, см-1 а.е. Еь (Еь+ВББЕ), кДж/моль Eb+ZPE (Eb+ZPE+BSSE), кДж/моль

ЛМ 1 (ГМ) ЛМ 2 с5 51/ -301.44638263 11.83 (8.59) (перегруппировывается в ЛМ 1) 8.06а (4.82)

ЛМ 3 (перегруппировывается в ЛМ 1)

ЛМ 4 Сs - -301.44542875 9.33 4.66

ЛМ 5 Сs -301.44541893 9.30

ЛМ 6 С2м 134/, 37/ -301.44496112 8.10 4.33а

ЛМ 7 С2м со со -301.44358527 4.49 3.18а

а - ТД вклад колебаний с мнимой частотой не учитывается

(*> 119.3 2.981

с^сВЧрР06-5) (2-929) 1 Ж)

------^у^ч

---------(К)

3.390 ^

(2.951)

0.9бв з 747 ~ ~ ~ ~ ЧО)

Ф)

(С)

165.7

(172.5) 2.547

(2.238) __

100.1 (111.9)

ГС

О)'

2.526 .(2.553)

(Ф

92.4 (95.7)

153.2 (144.1)

1.290 (1.282)

2.984 (3.068)

Рис. 1. Структура ГМ (а) и ЛМ (Ь^) дискретной ППЭ комплекса Н20-03 (MP2/6-311++G(2d,2p)). Значения без скобок - геометрические параметры ЛМ дискретной ППЭ. Значения в скобках - результат полной оптимизации геометрии на том же уровне расчета (в случае, когда оптимизация сходится к той же структуре).

Рис. 2. Структура и геометрические параметры глобального минимума ППЭ комплекса H2O•O3, оптимизированного на уровне MP2/6-311++G(2d,2p).

240 300 360 2.5 3.0 3.5 4.0

Я? А

Рис. 3. Одномерные сечения ППЭ комплекса H2O•O3, проходящие через точку ГМ.

Рассчитанные полная энергия, энергия комплексо-образования, величины БЭвЕ и 7РЕ, а также корректированные величины энергии связи приведены в табл. 1.

Часть оптимизированных структур локальных минимумов дискретной ППЭ соответствовала истинным минимумам (отсутствовали мнимые колебательные частоты). Однако большая часть структур имела мнимые частоты, их величина обычно составляла 10-100 см-1. Такая ситуация является типичной при расчете слабосвязанных комплексов. Вследствие малой величины мнимых частот, нередко на основе результатов оптимизации, невозможно однозначно определить, является ли оптимизированная структура ложным минимумом (ошибочно найденным программой оптимизации вследствие того, что ППЭ системы почти плоская), либо мнимые частоты являются следствием численных ошибок при расчете гессиана в точке истинного, но слабовыра-женного минимума. Этот факт является еще одним аргументом в пользу использования предлагаемой здесь методики: наличие минимума на ППЭ подтверждает наличие локального минимума, даже если гармонические колебательные частоты являются мнимыми.

Рассчитанная ППЭ позволяет изучить профили ППЭ в окрестности точки глобального минимума. На рис.3 показаны одномерные профили ППЭ системы Н20 03, рассчитанные на уровне МР2/6-311++0(2<<,2р) в направлениях изменения всех координат. Как видно из приведенных данных, ППЭ системы является существенно ангармоничной, с сильным смешиванием колебательных и вращательных степеней свободы.

Для расчета термодинамических параметров данного комплекса использовалась интерполяция значений межмолекулярного потенциала линейной полиномиальной функцией по вышеописанному методу. При интегрировании по Е верхний предел интегрирования составлял 50 кДж моль-1, шаг интегрирования 0.5 кДж моль-1. Для каждого значения Е фазовый интеграл определялся методом Монте-Карло с простой случайной выборкой и подвижными пределами интегрирования по фазовым переменным, число точек интегрирования - 10000. Пределы интегрирования по Я выбирались на основе предварительного анализа ППЭ для каждого значения Е. Максимальное значение интервала интегрирования составляло [2.25-6.80А]. Как показывает рас-

Таблица 2

Молекулярные параметры для расчета ТД параметров реакции образования комплекса H2O■O3 (расчет MP2/6-311++G(2d,2p))

Молекулярный параметр I Н2О | О3 | Н2ОО3

Спиновая мультиплетность 1 1 1

Молекулярная масса, а.е.м. 18.0106 47.9847 65.9953

Главные моменты инерции, а.е. 6.3316 155.2853

4.1203 137.8867

2.2112 17.3985

Число вращательной симметрии 2 2 1

Масштабирующий фактор колеба-

тельных частот 0.9500 0.9500 0.9500

Колебательные частоты, см-1 1660.9 735.7 (74.8) *

3848.6 1154.0 (92.5)

3974.0 2305.5 (99.2) (112.5)

(163.6)

(245.6)

739.0

1151.8

1660.7

2298.7

3850.7

3969.0

* В скобках приведены частоты колебаний, которые соответствуют ровибрацион-ным степеням свободы в приближении КЛАРВ/ККГ.

Таблица 3

Стандартные ТД параметры комплекса H2O■O3, рассчитанные в различных приближениях

ТД параметр Приближение

ЖРГО КЛАРВ/КГК (без коррекции) Коррект.КЛАРВ/КГК ^о = 0.9075)

А£, кДж моль-1 -8.67 -8.67 -8.67

Лги°(298), кДж моль-1 -0.69 -5.10 -5.20

ЛГН°(298), кДж моль-1 -3.17 -7.58 -7.68

А^°(298), Дж К-1 моль-1 -67.05 -62.30 -63.20

ЛrG0(298), кДж моль-1 16.82 11.00 11.16

К°(298) 1.1310-3 1.18-10-2 1.11-10-2

чет, такой выбор параметров интегрирования приводит к попаданию в фазовый объем системы ~70% случайно выбираемых точек. Время расчета - от 10 до 40 мин (Се1егоп/1500 МГц) в зависимости от выбора других параметров расчета.

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

плекса, используемые для расчета ТД параметров, представлены в табл. 2. Результаты расчета термодинамических параметров при вариации параметров расчета показаны в табл. 3.

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

На рис.4 приведена зависимость стандартных констант комплексообразования Н2003 от температуры, рассчитанная в гармоническом приближении и на основе интерполированной ангармонической ППЭ при явном учете ровибрационных степеней свободы (приближение КЛАРВ/КГК). Как следует из приведенных данных, величина К0, вычисленная с учетом эффектов ангармоничости, значительно превышает результаты гармонического расчета, причем различие усиливается при

пони-жении температуры. Этот факт может иметь боль-

Рис. 4. Температурная зависимость стандартной константы равновесия комплекса H2O-O3, вычисленная в приближениях ЖРГО и КЛАРВ/КГК.

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

Выводы

Впервые квантовохимическим методом MP2/6-311++G(2d,2p) исследована полная межмолекулярная поверхность потенциальной энергии слабосвязанного бинарного молекулярного комплекса озона и воды H2O O3, включающая 315 900 точек, среди которых 14 016 уникальных по симметрии. Определена структура глобального минимума, вычислены его энергия связывания и колебательные частоты. С использованием приближения, учитывающего ангармонизм ППЭ, определены термодинамические функции данного комплекса в состоянии идеального газа, оценена его константа равновесия и ее температурная зависимость. Эти характеристики существенно отличаются от результатов расчетов, основанных на более простом приближении ЖРГО, что свидетельствует о важности учета ангармонизма при оценке термодинамических функций слабосвязанных комплексов и перспективности развиваемого метода.

Литература

1. Slanina Z, Uhlik F., Saito AT, Osawa E. Computing Molecular Complexes in Earth's and Other Atmosphères // Phys. Chem. Earth. (C). 2001. Vol. 26. No. 7. P. 505-511.

2. Headrick J.E., Vaida V. Significance of water complexes in the atmosphere // Phys. Chem. Earth. (C). 2001. Vol. 26. No.7. P. 479-486.

3. Horn A.B., Sodeau J.R. Interactions and Photochemistry of small molecules on ice surface: from atmospheric chemistry to astrophysics, in Water in Confining Geometries, Buch V., Delvin J.P., Editors. 2003. Springer Verlag: Berlin Heidelberg. P. 295-336.

4. Abbatt J.P.D. Interaction of atmospheric trace gases with ice surface: adsorption and reaction // Chem.Rev. 2003. Vol. 103. P. 4783-4800.

5. Implications for water monomer and dimer solar absorption from observations at Boulder, Colorado / J.S.Daniel, S.Solomon, R.W.Sanders, R.W.Portman, D.C.Miller, W.Madsen. // J. Geophys. Res. 1999. Vol. 104, No. 14. P. 16785-16791.

6. Vigasin A.A., Camy-Peyret C., eds. Weakly Interacting Molecular Pairs: Unconventional Absorbers of Radiation in the Atmosphere / NATO Science Series IV: Earth & Environmental Sciences. 2003. 308p.

7. Vaida V., Headrick J.E. Physicochemical Properties of Hydrated Complexes in the Earth's Atmosphere // J. Phys. Chem. A. 2000. Vol. 104. No. 23. P. 5401-5412.

8. Wittig C, Sharpe S., Beaudet RA. Photoinitiated reactions in weakly bonded complexes // Acc. Chem. Res. 1988. Vol. 21. No. 9. P. 341-347.

9. Gillies J.Z., Gillies C.W., Suenram R.D., Lovas F.J., Schmidt T, Cremer D. A microwave spectral and ab initio investigation of O3-H2O // J. Mol. Spectr. 1991. Vol. 146. No.2. P. 493-512.

10. Schriver L., Barreau C, Schriver A. Infrared spectroscopic and photochemical study of water-ozone complexes in solid argon // Chem. Phys. 1990. Vol. 140. No.3. Р. 429-438.

11. Молекулярный комплекс озона с водой. Неэмпирический расчет с учетом электронной корреляции / И.И.Захаров, О.И.Колбасина, Т.Н.Семенюк, Н.Ф.Тюпало, Г.М.Жидомиров // Журн. структ. химии, 1993. Вып. 34. №3. С. 26-30.

12. Tachikava H, Abe S. Ozone-water 1:1 complexes O3-H2O: An ab initio study // Inorg. Chem. 2003. Vol. 42. No.7. Р. 2188-2190.

13. Seinfeld J.H., Pandis S.N. Atmospheric Chemistry and Physics: From Air Pollution to Climate Change. 1998, N.Y.: John Wiley and Sons.

14. Wayne R.P. Chemistry of Atmospheres. Oxford: Clarendon Press, 1991. 645 p.

15. General atomic and molecular electronic structure system / M.W.Schmidt, K.K.Baldridge, J.A.Boatz, S.T.Elbert, M.S.Gordon, J.H.Jensen, S.Koseki, N.Matsunaga, K.A. N., S.J.Su, T.L.Windus, M.Dupuis, J.A.Montgomery // J. Comput. Chem. 1993. Vol. 14. No.11. Р. 1347-1363.

Статья поступила в редакцию 07.12.2010.