CC BY
55
УДК 570.1:532.7
СРЕДНЕЕ ВРЕМЯ ЖИЗНИ ВОДОРОДНЫХ СВЯЗЕЙ В ВОДЕ И МЕТАНОЛЕ ПО РЕЗУЛЬТАТАМ КОМПЬЮТЕРНОГО МОДЕЛИРОВАНИЯ
О.В. Ведь1, Н.А. Абакумова2, М.Л. Антипова1, В.Е. Петренко1
Институт химии растворов РАН, г. Иваново (1); кафедра «Химия», ГОУВПО «ТГТУ» (2)
Представлена членом редколлегии профессором Н.Ц. Гатаповой
Ключевые слова и фразы: водородные связи; времена жизни водородной связи; межмолекулярные взаимодействия.
Аннотация: Методом классической молекулярной динамики рассчитано среднее время жизни водородных связей в воде и метаноле с использованием различных подходов и разных модельных потенциалов. Показано, что времена жизни водородных связей характеризуют высокую корреляцию движений молекул ближайшего окружения.
Структура и динамические свойства ассоциированных жидкостей, таких как вода и спирты, простейшим из которых является метанол, определяются, главным образом, водородными связями (ВС). Также можно утверждать, что, в определенном смысле, сетка ВС и является структурой таких систем, поэтому изучение сетки ВС неотделимо от исследования жидкостей, характеризующихся ее наличием. В последнее время большая часть публикаций, посвященных компьютерному моделированию этих жидкостей, тем или иным образом связана с изучением сеток водородных связей.
В компьютерном моделировании пара молекул считается связанной, если она удовлетворяет заранее сформулированным требованиям, называемым критериями образования ВС. Известно несколько критериев, с помощью которых стараются отличить пары молекул, образующих водородную связь. Широко используемые определения ВС основаны, либо на геометрическом, либо на энергетическом критериях. В принципе считается, что ВС образуется, если энергия ВС (ЕВС) больше (алгебраически) некоторой выбранной величины, которая, как правило, изменяется в пределах -2...-5 ккал/моль. В случае геометрического критерия выбираются три предельных значения, которые должны удовлетворяться для того, чтобы ВС существовала (рис. 1). В качестве предельных обычно выбирают следующие значения, полученные из анализа дифракционных исследований (рентген и дифракция нейтронов) кристаллогидратов: Я00 = 3,6 А (или 3,3 А), Я0Н = 2,4 А и а = 30°. Очень часто в исследованиях ВС используется комбинированный критерий ВС, представляющий собой сочетание энергетического и геометрического.
Если топология характеризует взаимное расположение и ориентацию молекул - иными словами, саму форму сетки ВС, то динамический подход рассматривает длительность происходящих в жидкости процессов: движения и переориентации молекул, разрушения и восстановления связей между ними. В компьютерном эксперименте динамического типа, каким является классическая молекулярная динамика (МД), грань между этими двумя сторонами проявления феномена водородной связи практически стирается, поэтому данный метод позволяет параллельно применять оба подхода. С динамической точки зрения сетку водород-
Roh — 2,6 A
/ту*
— 3,3 a
ных связей начали изучать лишь около Донор
двадцати лет назад. Здесь главными характеристиками ВС являются вероятность ее образования и время ее жизни. В подавляющем большинстве случаев под време- Акцетгор
нем жизни водородной связи понимают пе- "■
риод, в течение которого в данной паре молекул соблюдается выбранный критерий существования ВС.
Различают два типа времени жизни Roo — ■
водородной связи [1]. Непрерывное
(continious) или короткое время iSp , кото- Рис 1 Схема> п°яшяющая
ВС геометрический критерий
рое длится до первого нарушения в данном образования водородной связи
димере выбранного критерия существования связи. Если данная связь разорвалась, а потом восстановилась, то она рассматривается как новая, вне зависимости от длительности периода разрыва. Мерцающее (intermittent) или продолжительное время жизни водородной связи тВс длится до тех пор, пока один из участников данной ВС не найдет себе партнера для образования новой связи. В этом случае промежуточные нарушения критерия существования связи не учитываются вне зависимости от их длительности.
На сегодняшний день существует несколько подходов к изучению динамики сетки водородных связей.
В первом и наиболее популярном методе строится автокорреляционная функция (АКФ) следующего вида
X (t0 )SIJ (t)
J
X Sfj (t0) 'Л tbc
С (t) = I,J ^ — = expl------— I, (1)
I, 3
где параметр Бц (/) = 1, если в момент / между молекулами I и 3 связь существует, в противном случае Бц (/) = 0 [1].
Автор работы [1] показал, что разброс текущих значений функции С(/), обусловленный выбором точки отсчета /о, составляет около 3 % в области С(/) = 0,5 и достигает 15 % при значениях С(/) = 0,1, где становится значительным вклад от статистического шума. Также им была изучена зависимость функции С(/) от значения энергетического параметра и от температуры. Количественных оценок в работе [1] не приводится, однако по графикам, содержащимся в ней, видно, что как при повышении температуры, так и при ужесточении критерия спадание функции С(/) становится более резким, а значит время жизни водородной связи уменьшается. С помощью данного метода были определены диапазоны возможных значений времен жизни ВС как для продолжительного, так и для короткого варианта. Первое из них тВс = 1.10 пс качественно согласуется со временем дипольной релаксации, а второе имеет величину на два порядка меньше твс = 0,05.0,3 пс. Это показывает, что молекула испытывает десятки флуктуаций взаимной ориентации относительно партнера по уже существующей связи, прежде чем образует новую. Аппроксимация АКФ существования ВС одной экспонентой является, по признанию самого автора [1], слишком грубым приближением, оправданным лишь на больших временах. Поэтому он предлагает оценивать «время полураспада ВС» т^, связанное с твс соотношением
Т1/2 = —. (2)
1/2 1п2
АКФ существования ВС были использованы при расчете времени жизни ВС и в работе [2], однако, вместо двоичного параметра Бц (/) был введен оператор заселенности водородной связи к(()
к(/) = 0(фс - ф, /)0(гс - г, /), (3)
где 0 - степ-функция Хевисайда. Как видно из формулы (3), для отбора ВС применен геометрический критерий. Критические значения выбраны из следующих
соображений: расстояние гс = 2,45 А удовлетворительно согласуется со значением гОГ (положение первого минимума функции радиального распределения кислород-водород), а угол фс = 30° выбран так, чтобы среднее число водородных связей в равновесном ансамбле из N молекул не зависело от значения фс. Затем авторы
[3] видоизменили формулу (1). Для времен, превышающих время перехода, она
имеет вид
С( )=Ц§=1 (4)
Как можно видеть, данная формула придает функции С(/) новый смысл - теперь это АКФ констант скоростей образования водородной связи к(/), которые рассчитываются по формуле
к(/) = “7^(ХХЧ (0)^(/)\. (5)
Начальное значение к(0) - это оценка константы скорости, полученная из теории переходного состояния, а отношение к (/ )/к (0) описывает отклонение реальной константы от данной оценки. Показано, что благодаря динамическим эффектам теория переходного состояния на порядок переоценивает значения констант скоростей образования связей, а значит, недооценивает времена их жизни.
В работах Г.Г. Маленкова [4, 5] среднее время жизни ВС определяется следующим образом. На основе выбранного критерия для каждой молекулы в ячейке строятся распределения водородных связей по времени их существования на протяжении компьютерного эксперимента. При анализе полученных распределений рассматриваются только те водородные связи, которые образовались после начала моделирования и разорвались до его окончания, что позволило авторам [4, 5] получить интервал значений времен жизни водородных связей Твс = 6 ± 2 пс (в расчетах использовался потенциал Маленкова-Полтева (МП)).
В работе [6] Б. 8сюг1то и 8.Ь. БогтИ рассматривали зависимость от времени отношения числа ВС в данный момент ^) к числу ВС в начальный момент N(0)
^ = N(/)/N(0). (6)
По истечении небольшого начального периода /Ъе§ ~ 3 пс зависимость 1пР(/) выходит на прямую линию, с помощью тангенса угла наклона которой можно определить время жизни водородной связи.
Фактически, аналогичный метод был использован авторами [7] при изучении динамики сетки ВС в метаноле. На каждом шаге моделирования подсчитывалось число связей Щ/„), которые разрушились по истечении ровно п временных шагов.
Если за время наблюдения за системой какая-либо водородная связь успевала разрушиться и восстановиться, то она рассматривалась как образованная заново,
вне зависимости от длительности периода разрыва. Таким образом, авторы [7]
были сторонниками концепции «непрерывного» времени жизни водородной связи. Затем на основе полученных значений N (/п) строилась функция вида
ТО / ТО
н (т ) = 2 N (п+1) 2 N (/п ), (7)
п=т / п=1
где Н (/т) представляет собой долю водородных связей, сохранившихся по истечении времени /т. Среднее время жизни водородной связи оценивали, интегрируя функцию Н (/т )
ТО
ТВС = ^ 2(Я (/п)+Н (/п+1)). (8)
0
Другой способ расчета времени жизни водородной связи основан на оценке вероятности ее существования (или сохранения) рв . Со временем жизни ВС эта величина связана следующими соотношениями [8, 9]:
РВ =-3+-; ^=-а^ (9)
тВС +т/ т^ 1 - Рв
Авторы [8, 9] предположили, что водородные связи каждой молекулы энергетически равноценны и образуются и разрываются по одной. В рамках модели пяти состояний была получена следующая система дифференциальных уравнений
= -4 Ар0 + Ярь
а/
= 4АР0 -(в + 3А)Рь
а/
% = 3Ар1 - 2(А + В)р1, (10)
а/
аР3 = 2(А + В)р2 -(А + 3В)р3 + 4Вр4; а/
^ = АР3 - 4BР4, а/
где р1 - вероятность перехода молекулы в состояние с I сохраненными связями; / -время. Вероятности перехода связаны со средним числом ВС на молекулу п следующим уравнением
п(/) = Р1(/) + 2Р2(/) + 3р3(/) + 4Р4(/) = 4р(/), (11)
где Р - средняя вероятность образования ВС. Скорость изменения этой величины можно записать, согласно системе (10), в следующем виде
0р = 1 аР1 + 1 ар2 , 3 Фэ + аР4 (12)
а/ 4 а/ 2 а/ 4 а/ а/
и, в то же время, как
Ф =_^+Ьр . (13)
dt ТВС т/
В результате подстановки скоростей из системы (10) с последующей перегруппировкой получим:
А = —; В = —^, (14)
т / тВС
где хвс - время жизни связи; т^ - время ее разрыва [10].
Таким образом, время жизни ВС можно получить, изучая динамику среднего числа ВС в расчете на молекулу
(л л \
dn
Bc
dt
11
- + -
TBC Tf
4
nBC +-------------------------------------. (15)
т f
Подобный подход использовали М. Ма1Бито1о и К.Е. виЪЫш [11] при изучении динамики ВС метанола.
Приведенные здесь методики являются прямыми, то есть они лишь констатируют факт образования или разрыва водородной связи, не объясняя причину этого явления. В этом кроются и их достоинства, и их недостатки: с одной стороны, им неважно, какой выбран критерий образования связи, но при этом, если критерий сменить, то результаты для той же системы при тех же условиях будут уже другими. Таким образом, время жизни также становится очень субъективной величиной, что усложняет сопоставление между литературными данными. Один из выходов - пересчитать все цитируемые результаты, используя один критерий образования связи.
Среднее время жизни водородной связи
При расчетах среднего времени жизни водородных связей в воде или метаноле, кроме выбора метода, оказывают весьма существенное влияние еще два момента: выбор потенциальной функции парного взаимодействия и выбор критерия, по которому оценивается существование ВС между молекулами. Это хорошо видно из табл. 1, в которой приведены времена жизни водородных связей в воде и метаноле, полученные в ряде работ с использованием различных методов расчета, критериев ВС и потенциалов взаимодействия.
Как можно видеть, разброс в результатах очень значителен. Он заранее предопределен не только многообразием подходов, способов расчета, но и выбором модельного потенциала и критерием ВС. Мы попытались выявить влияние каждого из этих факторов.
Зависимость времени жизни ВС от метода расчета проиллюстрирована в табл. 2 на примере модели Т1Р4Р.
Из табл. 2 видно, что время жизни ВС значительно зависит от метода расчета. Это вполне закономерно. Для получения оценок времен жизни ВС методом, описанным в работах [8, 9], данные, получаемые в ходе моделирования, необходимо численно дифференцировать по времени, а затем проводить их линейный регрессионный анализ, что проще и удобнее делать с помощью внешней программы. Проведенные нами исследования показали, что данный метод [8, 9] применим для оценки только короткого времени жизни ВС, поскольку при учете восстанавливающихся ВС, их общее число (и среднее на молекулу) колеблется около постоянного значения, и ни о какой линейной корреляции между его значениями и скоростями его изменения говорить не приходится. Если же разрушенную ВС исключать из анализа необратимо, то получаем спадающие функции, что и необ-
С I
Короткое (тВс) и продолжительное (твс) времена жизни водородных связей в воде и метаноле по литературным данным
Модель* Критерий ВС ТВС > пс ТВ С > пс Литература
Вода
NCC - Q,Q3...Q,Q7 - [12]
MCY ІЕВС| > 8,37 - 16,74 кДж/моль Q,Q5...Q,3Q 1...1Q [1]
SPC г00 < 3,4 А, Гон < 2,4 А Q,5Q 3,QQ [13, 14]
гоо < 3,6 А, Гон < 2,45 А 3,QQ [2]
SPC2 Гоо < 3,3 А, Гон < 2,6 А Q,68 - [15]
Пространственный критерий Q,76 -
МП Гоо < 3,3 А, Гон < 2,6 А б ± 2 [4]
Q,Q1... 16,QQ [16]
TIP4P гоо <3,5 А, і Евс|> 14,64 кДж/моль Q,65 - [8, 9]
ST2 Гоо < 3,5 А, I Евс |> 14,64 кДж/моль Q,Q2...Q,2Q - [17]
I Евс I >9,41 кДж/моль 1,QQ [18]
Метанол
н1 roo < 3,5 А, roн < 2,6 А 1,2Q -
н2 roo < 3,5 A, roн < 2,6 А 2,5Q - [7]
J1 roo < 3,5 А, ^н < 2,б А Q,7Q -
roo < 3,5 А, ^н < 2,б А 1,4Q -
J2 roo < 3,5 А, ^н < 2,б А 1,5Q 16,5Q [13, 14]
rc** < 3,8 А 1Q,7Q [11]
IЕВС I > 12,QQ кДж/моль 11,8Q
*
Параметры моделей для воды и метанола подробно приведены в [19 - 21].
** С
г - расстояние между центрами масс в молекуле.
Таблица 2
Времена жизни водородных связей в воде для модели ТІР4Р, рассчитанные различными способами. Критерий - геометрический: гоо < 3,3 А, г он < 2,6 А
Метод расчета по формуле ТВС > пс тВС, пс
(1)* Q,41 4,17
(б)** - 5,26
(9)*** Q,5Q9 -
*
Из АКФ существования ВС [1].
Через отношение числа ВС в данный момент к числу ВС в начальный момент [6]. Через вероятность существования ВС [8, 9].
ходимо в данном методе. В табл. 3 приведены оценки времен жизни для популярных модельных потенциалов для воды, полученные методом ВеАоИт й а1. [8, 9].
Из табл. 3 можно видеть, что полученные данным методом твс отличаются от результатов, полученных методом (1) [1] (результаты для соответствующих потенциалов приведены в табл. 4). Проанализировав полученные данные, мы остановились на том, что метод, основанный на построении АКФ существования ВС (1), более информативен, поскольку эти АКФ отражают как состояние сетки ВС в каждый момент моделирования, так и ее изменение со временем. Этот метод мы использовали в дальнейшем.
Мы провели расчет средних времен жизни ВС в воде и метаноле с использованием наиболее популярных модельных потенциалов. Был применен геометрический критерий ВС: для воды г00 < 3,3 А, г0н < 2,6 А, для метанола г00 < 3,5 А, г0н < 2,6 А. Использовался метод АКФ существования ВС по (1). Значения времен жизни рассчитаны интегрированием функций С(/) аналогично [1]. Полученные нами данные, иллюстрирующие влияние выбора потенциальной функции на результаты расчета, представлены в табл. 4.
Во всех случаях (для всех моделей) соотношения между временами жизни ВС сохраняются: продолжительное время на порядок больше короткого и для во-
С I
ды, и для спирта, при этом тВс для спирта - величина того же порядка, что и Твс
воды. На наш взгляд, это можно объяснить следующим образом. При построении АКФ для расчета короткого времени жизни допускается только один переход параметра от единицы к нулю, то есть разорванная связь больше не учитывается,
даже если в ходе моделирования она восстанавливается. В случае продолжительного времени жизни учитывается восстановление любой ВС, вне зависимости от
Таблица 3
Времена жизни водородной связи (9)
Модель SPC TIP4P SRWK2 TIPS2 BF TIP4P-HB*
Твс пс Q,952 Q,5Q9 1,Q1Q Q^5 Q,756 Q,66Q
*
Вид и параметры модели приведены [22, 23]
Таблица 4
Времена жизни водородных связей в воде и метаноле по результатам расчетов с использованием геометрического критерия г00 < 3,3 А, гон < 2,6 А - для воды и г0о < 3,5 А, г0н < 2,6 А - для метанола
Вода Метанол
Модель тВс , пс ТВс j пс Модель ТВс j пс ТВ с , пс
SPC Q,B9 4,84 H1 ЗДО 15,99
SPC2 Q^Q 4,2Q H2 5,2Q 24,91
TIP4P Q,41 4,17 J1 1,66 6,96
BF Q^B 4,1Q J2 З,69 17,Q2
SRWK2 Q,22 5,QB
TIPS2 Q,21 5,45
TIP4P-HB Q,46 4,QB
длительности периода ее отсутствия. Вероятность восстановления существовавшей связи очень высока, поэтому продолжительное время жизни ВС на порядок превышает короткое. При расчете продолжительного времени жизни существует вероятность (хотя и очень малая) того, что рассматриваемая связь разорвется в начале моделирования, а восстановится в конце, и тогда временем ее жизни будет считаться весь период наблюдения за системой.
При сопоставлении значений средних времен жизни ВС, приведенных в табл. 4, видно, что в отличие от продолжительных, короткие времена жизни ВС для разных потенциалов воды расходятся почти в два раза. Явный учет в потенциальной функции неэлектростатической составляющей энергии образования ВС (потенциал TIP4P-HB) приводит к некоторому снижению продолжительного времени жизни водородной связи и незначительному увеличению короткого времени жизни ВС. Аналогичные расхождения наблюдаются в метаноле. Более того, среди потенциалов для метанола модель H2 дает наибольшее время жизни ВС. К сожалению, нет принципиальной возможности сравнивать рассчитанные времена жизни водородной связи с экспериментальными, чтобы оценить корректность полученных результатов. Дело в том, что не существует определенного способа получения времени жизни ВС экспериментальным путем. Проблема заключается в том, что все экспериментальные методы очень чувствительны ко времени присутствия, а для такой системы как жидкая вода, частицы которой склонны к образованию одновременно нескольких водородных связей, то есть координационных оболочек, время присутствия намного превышает время жизни ВС. Тем не менее, существует ряд непрямых (indirect [24]) методов определения времени жизни, но они дают настолько расходящиеся между собой результаты, что использовать их в качестве опорных нецелесообразно. Например, в работе [25] получено значение 2,1 пс (использовался метод микроволновой спектроскопии), а в работе [26] - 0,6 пс (время жизни ВС оценивалось по спектрам деполяризованного молекулярного рассеяния света).
Таким образом, для дальнейших исследований динамики ВС жидкостей, очевидно надо принять решение относительно того, мерцающая, непрерывная или еще какая-то версия жизни ВС предпочтительнее для описания конкретных экспериментально определяемых величин. Кроме того, пороговое значение энергии в энергетическом и смешанном критерии должно быть заключено в достаточно узкий диапазон с помощью независимых методик.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований, проект № 07-03-00646, и Фонда содействия отечественной науке.
Список литературы
1. Rapaport, D.C. Hydrogen bonds in water. Network organization and lifetimes / D.C. Rapaport // Mol. Phys. - 1983. - Vol. 50, part 5. - P. 1151-1162.
2. Luzar, A. Resolving the hydrogen bond dynamics conundrum / A. Luzar // J. Chem. Phys. - 2000. - Vol. 113, part 23. - P. 10663-10675.
3. Luzar, A. Structure and hydrogen bond dynamics of water dimethylsulfoxide mixtures by computer simulations / A. Luzar, D. Chandler // J. Chem. Phys. - 1993. -Vol. 98, part 10. - P. 8160-8173.
4. Метод молекулярной динамики в физической химии (Method of molecular dynamics in physical chemistry) / под ред. Ю.К. Товбина. - М. : Наука, 1996. - 308 c.
5. Маленков, Г.Г. Динамический критерий водородной связи для анализа структуры водных кластеров / Г.Г. Маленков, Д.М. Тытик // Известия АН. - Сер. Физика. - 2000. - Т. 64, вып. 8. - С. 1469-1472.
6. Sciortino, F. Hydrogen bond cooperativity in simulated water: Time dependence analysis of pair interactions / F. Sciortino, S.L. Fornili // J. Chem. Phys. -1989. -Vol. 90, part 5. - P. 2786-2792.
7. Haughney, M. Molecular dynamics simulation of liquid methanol / M. Haugh-ney, M. Ferrario, I.R. McDonald // J. Chem. Phys. - 1987 - Vol. 91, part 19. -P. 4934-4940.
8. A “microscopic” model for the dynamics of water / D. Bertolini, M. Cassettari, M. Ferrario, G. Sallvetti // Chem. Phys. Lett. - 1985. - Vol. 62, part 6. - P. 548-553.
9. Bertolini, D. Diffusion effects of hydrogen bond fluctuations. The long time regime of the translation and rotational diffusion of water / D. Bertolini, M. Cassettari, M. Ferrario // J. Chem. Phys. - 1989. - Vol. 91, part 2. - P. 1179-1190.
10. Hydrogen bond statistics and dynamics in water: Self-diffusion and dielectric relaxation / D. Bertolini, M. Cassettari, M. Ferrario, G. Sallvetti // J. Chem. Phys. - 1984. - Vol. 81, part 12. - P. 6214-6223.
11. Matsumoto, M. Hydrogen bonding in liquid methanol / M. Matsumoto, K.E. Gubbins // J. Chem. Phys. - 1990. - Vol. 9, part 3. - P. 1981-1994.
12. Niesar, U. Preliminary observations on a new water-water potential / U. Niesar, G. Corongiu, M.J. Huang // Int. J. Quant. Chem.: Quant. Chem. Sym. -1989. - Vol. 23. - P. 421-443.
13. Padro, J.A. Hydrogen bonding in liquid alcohols: a computer simulation study / J.A. Padro, L. Saiz, E. Guardia // J. Mol. Struct. - 1997. - Vol. 416, part 416. -P. 243-248.
14. Dynamics in hydrogen bonded liquids: water and alcohols / E. Guardia, J. Marti, J.A. Padro and ets. // J. Mol. Liq. - 2002. - Vol. 96-97. - P. 3-17.
15. Bergman, D.L. Topological properties of the hydrogen - bond network in liquid water / D.L. Bergman // J. Chem. Phys. - 2000. - Vol. 253, part 2/3. -P. 267-282.
16. Malenkov, G.G. Hydrogen bonds in computer - simulated water / G.G. Malenkov, D.L. Tytik, E.A. Zheligovskaya // J. Mol. Liq. - 1999. - Vol. 82. -P. 27-38.
17. Structure and dynamics of the hydrogen-bond network in water by computer simulation / A. Geiger, P. Mausbach, J. Schnitker and ets. // J. Physique. - 1984. -Vol. 45, part 7. - P. 7 - 30.
18. Sciortino, F. Lifetime of the bond network and gel-like anomalies in supercooled water / F. Sciortino, P.H. Poole, H.E. Stanley // Phys. Rev. Lett. - 1990. -Vol. 64, part 14. - P. 1686-1689.
19. Кесслер, Ю.М. Вода: структура, состояние, сольватация. Достижения последних лет (Water: Structure, State, Solvation. Achievements of the Last Years) / Ю.М. Кесслер, В.Е. Петренко ; под ред. А.М. Кутепова (A.M. Kutepov). - М. : Наука, 2003. - 420 с.
20. Haughney, M. Molecular-Dynamics Simulation / M. Haughney, M. Ferrario, I.R. McDonald // J. Chem. Phys. - 1987. - Vol. 91, part 19. - P. 4934 - 4940.
21. Jorgensen, W.L. Transferable intermolecular potential functions for water, alcohols and ethers / W.L. Jorgensen // J. Amer. Chem. Soc. - 1981. - Vol. 103, part 2. - P. 345-340.
22. Petrenko, V.E. Sensivity of model water structure to hydrogen bond included / V.E. Petrenko, M.L. Antipova, A.V. Borovkov // J. Mol. Liq. - 2005. -Vol. 120. -P. 19-21.
23. Влияние на структуру модельной жидкости учета водородных связей в потенциале парного взаимодействия вода-вода / А.В. Боровков [и др.] // Журн. структ. химии. - 2005. - Том 45, вып. 3. - С. 678 -682.
24. Luzar, A. Water hydrogen-bond dynamics close to hydrophobic and hydrophilic groups / A. Luzar // Faraday Discuss. - 1996. - Vol. 103. - P. 29-40.
25. Matubayasi, N. Super and subcritical hydration of nonpolar solutes. Thermodynamics of hydration / N. Matubayasi, M. Nakahara // J. Chem. Phys. - 2000. -Vol. 112, part 18. - P. 8089-8109.
26. Conde, O. Hydrogen bond dynamics in water studied by depolarized scattering / O. Conde, J. Teixeira // J. Physique. - 1983. - Vol. 44. - P. 525-529.
Average Life Cycle of Hydrogen Links in Water and Methanol Based on Computer Simulation
O.V. Ved1, N.A. Abakumova2, M.L. Antipova1, V.E. Petrenko1
Institute of Chemistry of Solutions of RAS, Ivanovo (1);
Department “Chemistry ”, TSTU (2)
Key words and phrases: hydrogen links; intermolecular interactions; life cycle of hydrogen link.
Abstract: By the technique of classical molecular dynamics the average life cycle of hydrogen links in water and methanol is calculated through different approaches and different model potentials. It is shown that life cycle of hydrogen links characterizes high correlation of motion of nearby molecules.
Mittlere Lebensdauer der Wasserstoffverbindungen im Wasser und im Methanol nach den Ergebnissen der Computermodellierung
Zusammenfassung: Von der Methode der klassischen Molekulardynamik ist die mittlere Lebensdauer der Wasserstoffverbindungen im Wasser und im Methanol mit der Nutzung verschiedener Herangehen und verschiedener Modellpotentiale berechnet. Es ist aufgezeigt, dass die Lebensdauer der Wasserstoffverbindungen die hohe Korrela-tion der Bewegungen der Molekule der nachsten Umgebung charakterisiert.
Temps de vie moyen des liens d’hydrogene dans l’eau et le carbinol d’apres les resultats du modelage d’ordinateur
Resume: Par la methode de la dynamique moleculaire classique est calcule le temps de vie moyen des liens d’hydrogene dans l’eau et le carbinol avec l’utilisation de differentes approches et divers potentiels de modelage. Est montre que les temps de vie moyen des liens d’hydrogene caracterisent une haute correlation des mouvements des mollecules de l’environnement le plus proche.
