УДК 548.737:541.632
ДИМЕРНЫЙ ПОДХОД В ОПИСАНИИ МЕЖМОЛЕКУЛЯРНЫХ ВЗАИМОДЕЙСТВИЙ В КРИСТАЛЛАХ ПОЛИФТОРИРОВАННЫХ СОЕДИНЕНИЙ
Е.В. Барташевич, ЕЛ. Шманина
Выполнено теоретическое исследование межмолекулярных взаимодействий с участием атомов фтора для полифторированных диазафлуоренов, енаминоиминоинданов и енаминоиминоинденов, образующих молекулярные кристаллы. Произведена оценка энергий межмолекулярных взаимодействий димеров, представляющих собой пары соседних в кристалле молекул. Выполнено сопоставление энергий взаимодействий в димере с квантово-топологическими характеристиками электронной плотности (ЭП) в критических точках связей, отвечающих межмолекулярным взаимодействиям, преимущественно для F...F, F...H и О...Н контактов.
Ключевые слова: полифторированные ароматические соединения, молекулярные кристаллы, межмолекулярные взаимодействия, квантово-топологичес-кий анализ электронной плотности, теория AIM.
Введение
Соединения фтора всегда представляли интерес, как для фундаментальной науки, так и для реализации ее прикладных аспектов. В области фторорганических соединений интенсивное развитие получила химия фтораренов. Актуальность исследований этих соединений базируется на уникальной возможности применения полифтораренов в качестве ключевых «строительных блоков» в фармацевтике - антибиотики, противораковые препараты, при создании жидких кристаллов и функционально обогащенных полимеров. Перспективными они остаются в дизайне и синтезе молекулярных проводников и магнетиков, инженерии кристаллов с заданными свойствами. В связи с разработкой новых органических материалов и высокоэффективных лекарств, стали популярными работы, связанные с заменой атомов Н на атомы F [1]. Атомы фтора могут имитировать атомы водорода в живом организме. Вследствие этого многие частично фторированные органические соединения обладают физиологической активностью. Так, гексафторбензол исследован как ингаляционный анестетик в ветеринарии; фторбензол и дифторбензол, полицикличе-ские фторароматические соединения - полупродукты в синтезе фармацевтических и пестицид-ных препаратов, 5-фторурацил используется в хемотерапии рака [2].
Кроме того, интерес к перфторированным ароматическим соединениям обусловлен нетривиальными особенностями их электронной структуры, богатством структурных типов. Расположение молекул в кристаллах определяется межмолекулярными взаимодействиями. Наиболее распространенными являются кристаллические мотивы, образованные за счет водородных связей. Как известно, введение атома фтора часто резко изменяет свойства органических соединений. Это обусловлено высокой электроотрицательностью, малым размером атома фтора. Связь C-F характеризуется малым межатомным расстоянием, энергия связи превышает энергию связей С-Н и С-С1 [3].
Исследования строения фторированных и перфторированных органических соединений продемонстрировали возможность образования новых мотивов, связанных с формированием F...F взаимодействий в кристаллах. В ряде работ отмечалось, что для веществ, имеющих два и более атома фтора в молекуле, оказывается возможным формирование даже F-агрегатов [4]. Подобные данные являются результатом подробного изучения строения и принципов упаковки полифторированных молекул в кристаллах. Они находятся в противовесе с устоявшимся мнением, что фтор не образует коротких межмолекулярных расстояний F...F. Обсуждение характеристик контактов в кристаллах фторарильных соединений, например, в работе [5] не дают однозначного ответа на вопрос, как должны соотноситься кристаллохимические радиусы F и Н, исходя из характеристик
контактов С Н. ..Ни CH.. .F. Таким образом, проблема интерпретации данных о сильных и слабых взаимодействиях с участием фтора на сегодняшний день остается актуальной.
Целью нашего теоретического исследования явилось описание межмолекулярных взаимодействий с участием фтора на примере представленного ряда соединений с помощью димерного подхода. Димерный подход подразумевает рассмотрение всех пар соседних молекул в кристалле. Это удобная структурная единица, позволяющая перечислить полный набор атом-атомных взаимодействий, если последовательно восстанавливается всё окружение молекулы, которая изначально выбирается за опорную. Для количественного описания взаимодействий, кроме энергетического критерия, который в таком приближении дает лишь некоторую относительную оценку, нами использованы электронные критерии. Для оценки характеристик межмолекулярного связывания удобно использовать квантово-топологический анализ электронной плотности, базирующийся на теории Бейдера AIM [6].
Следует отметить особую роль квантово-топологических характеристик электронной плотности в построении эмпирических зависимостей, пригодных для теоретической оценки физикохимических свойств соединений и их реакционной способности. Так, например, показаны возможности привлечения квантово-топологических характеристик для вычисления индексов порядков связей, позволяющих производить оценку направлений химических реакций [7]. Наиболее популярным привлечение расчетных квантово-топологических характеристик остается для прогноза энергетических характеристик взаимодействий в атомно-молекулярных системах. Например, для слабых водородных связей локальная плотность потенциальной энергии У(гь) может быть использована как дескриптор для определения энергии контакта согласно эмпирическому уравнению Есот - -■ 1 /2 V(rh), предложенному в работе [8]. Данная эмпирическая зависимость выведена с привлечением экспериментальных данных по энергиям диссоциации слабых водородных связей. Однако применение данных зависимостей, вычисленных для типичных и слабых водородных связей, с нашей точки зрения, не является обоснованным, если мы рассматриваем взаимодействия другого сорта, например, контакты с участием атома фтора. Зависимость параметров таких эмпирических уравнений от сорта атомов, участвующих в межмолекулярных взаимодействиях, очевидна и вид количественных описаний таких зависимостей должен быть, по крайней мере, скорректирован на диапазоны значений электронной плотности для новых рассматриваемых объектов.
Кроме того, подобная проблема корректировки параметров не является единственной. Следует отметить, что необходимым условием применения квантово-топологической теории является использование равновесной геометрии системы. Между тем, молекулярные структуры, полученные из данных рентгеноструктурного анализа, не обязательно отвечают равновесию. Как правило, это структуры, в которых положения атомов усреднены по времени, а легкие атомы водорода чаще всего локализованы механически. Кроме того, геометрия трехмерной упорядоченной кристаллической структуры отлична от равновесной геометрии фрагмента этого же кристалла, представленного димером - парой соседних молекул. Выделенная из кристаллического континуума пара молекул может не сохранять прежнюю взаимную ориентацию после локализации равновесной геометрии, т. е. в отсутствии взаимодействий с другими молекулами. Могут также наблюдаться существенные конформационные изменения и, как правило, удлиняться контакты межмолекулярных взаимодействий. В связи с этим, в круг поставленных задач вошло рассмотрение применимости результатов оценки межмолекулярных взаимодействий в димерном приближении для интерпретации свойств, связанных с межмолекулярным взаимодействием в кристаллическом состоянии.
Таким образом, оценка возможностей димерного подхода в описании межмолекулярных взаимодействий с участием фтора в полифторированных молекулярных кристаллах включает следующие задачи:
-произвести локализацию равновесной геометрии димеров, сопоставить равновесную геометрию димера и геометрию соответствующего фрагмента кристаллической структуры, в том числе, посредством вычисления среднеквадратичных смещений координат атомов - RMSD;
- использовать волновые функции B3LYP/6-311G для вычисления локальных свойств электронной плотности, представленных квантово-топологическими характеристиками в критических точках связей; установить количественные закономерности изменения этих величин в зави-
симости от межъядерных расстояний, скорректировать вычисленные значения с учетом изменений структуры в результате оптимизации геометрии димеров;
- сравнить энергии взаимодействия в димерах и величины характеристик ЭП в критических точках связей в межмолекулярном пространстве;
- количественно оценить межмолекулярные взаимодействия с участием фтора и сравнить характеристики связей для фрагментов С-F.. .X, где X = Н, С, N, О между собой и с водородными связями N-H.. .N N H.. .О, реализующимися в данных системах.
Экспериментальная часть
Объектами нашего теоретического исследования выступили шесть кристаллических структур полифторированных енаминоиминоинданов, енаминоиминоинденов и диазафлуорена, структура которых подтверждена методами PC А [9-11]. Интерес к данным структурам обусловлен особенностями формирования разнообразных специфических взаимодействий в кристаллах с участием атомов фтора.
Поскольку молекулярные структуры, полученные из данных рентгеноструктурного анализа, содержали элементы разупорядоченности, в особенности для атомов водорода, предварительно была произведена оптимизация геометрии димеров. Общее количество сформированных димеров составило 49. В 30 из них в межмолекулярных взаимодействиях принимают участие атомы водорода. Локализация равновесной геометрии димеров и изолированных молекул произведена посредством квантово-химических расчетов методом Кона-Шэма в приближении B3LYP/6-311G с помощью пакета программ PC GAMESS / Firefly version 7.1.G [12] на суперкомпьютерном кластере СКИФ-Урал.
Энергию взаимодействий В димере ВЫЧИСЛЯЛИ как AEinier = Edimer - 2Етопотег, где энергия димера и мономера представляет собой полную энергию, получаемую в результате локализации равновесной в димере геометрии.
Количественная оценка геометрических изменений, произошедшая в димерах при локализации равновесной геометрии произведена с помощью удобного критерия - среднеквадратичного отклонения расстояний, на которые произошло смещение атомов в результате оптимизации геометрии. Величина RMSD, позволяет сравнивать изменения позиций атомов при равенстве всех их весов, w = 1, RMSDjnt - изменения позиций только тех атомов, которые участвуют в межмолекулярных взаимодействиях. При этом веса остальных атомов при вычислении RMSDint принимались равными нулю: w = 0.
В том же приближении, что и оптимизация геометрии, для всех рассматриваемых структур димеров выполнены расчеты волновых функций, которые использовались для вычисления квантово-топологических характеристик ЭП, согласно теории Бейдера AIM [13]. Данный подход включает топологический анализ функции распределения электронной плотности р(г), полученной с помощью квантово-химических расчетов. Найдены экстремумы электронной плотности, являющиеся критическими точками связей. Критические точки ЭП, характеризующие связь, обнаруживаются в структурах димеров в области межмолекулярных взаимодействий и определяются как взаимодействия по типу закрытых электронных оболочек. Локальная потенциальная энергия V(r) определяет среднее эффективное поле, действующее на электрон в многоэлектронной системе. Лапласиан электронной плотности V2p(r) связан с плотностями кинетической и потенциальной энергии в каждой точке системы: -1/4 V2p(г) = 2G(r) + V(r). Для исследуемых пар димеров рассчитаны значения ЭП рО/,), лапласиана ЭП V2p(r/,), собственные числа матрицы Гессе ЭП Я.ь Х2, а также плотности полной H(rh) = V(rh) + G(rb), потенциальной V(rb) и кинетической G{rh) энергии электронов в экстремумах ЭП - критических точках, характеризующих связь. Все вычисления проведены с помощью пакета программ AIMPAC [6]. Для каждого димера найдены максимальные и средние значения квантово-топологических характеристик ЭП. Вместе со значениями энергий взаимодействий АЕШег и среднеквадратичными смещениями координат атомов в результате оптимизации RMSDi эти данные представлены в табл.1.
Таблица 1
Энергетические и квантово-топологические характеристики димеров пролифторированных соединений
№ Соединение Типы межмолекулярных контактов в димерах ЯМБО, ^ Ду'/)та\ У(Гь) шах АЕшег
1 2 3 4 5 6 7 8
1 3-(1-амино-2,2,2-трифторэтилиден)-2-метилимино-1,1,4,5,6,7-гексафториндан Р...Р,Р...Н 0,335 0,007 0,007 -0,003 -2,68
2 Б...Б, Б...С 0,320 0,005 0,006 -0,003 -2,05
3 Р../И, Р...С, Р...Н 0,552 0,006 0,005 -0,001 -4,61
4 Р...Р, Р...С 0,313 0,005 0,005 -0,003 -1,34
5 Р...Р,Р...С 0,320 0,005 0,006 -0,003 -2,05
6 Р...Р,К..Н 0,310 0,021 0,015 -0,001 -6,05
7 3-метил амино-2-метилимино-1,1,4,5,6,7-гексафториндан Р...Р,Р...Н 0,321 0,006 0,007 -0,002 -1,53
8 Р...Р,Р...Н 0,278 0,005 0,005 -0,003 -1,33
9 Р...Р,Р...С, Р...Н, С...С 0,297 0,008 0,007 -0,001 -4,03
10 Р...Р,Р...С, Р...Н,С...С 0,298 0,008 0,007 -0,001 -4,02
11 Р...Р,Р...:И, Р...С,Р...Н, С...с 0,386 0,005 0,006 -0,002 -2,75
12 Р...Р 0,308 0,006 0,007 -0,001 -2,16
13 Р...Н 0,324 0,006 0,006 -0,003 -2,86
14 2-амино-З -гексафто- ринденилтрифторме- тилкетон О...Н,Р...О, Р...Н 1,152 0,026 0,023 -0,003 10,44
15 О...Н,Р...О, Р...Н 0,722 0,027 0,023 -0,004 10,12
16 Р...Р,Р...С, Р...Н 1,225 0,013 0,012 -0,003 -5,69
17 Р...Р, Р...С,Р...Н, о...с,с...с 0,787 0,007 0,007 -0,002 -3,02
18 Р...Р,Р...С 0,384 0,005 0,005 -0,001 -3,02
19 Р...Р 0,089 0,005 0,006 -0,002 -0,73
20 Р...Р,Р...С 0,474 0,006 0,006 -0,002 -1,66
21 Р.-Т 0,549 0,005 0,006 -0,003 -0,73
22 Р...Н 0,770 0,010 0,009 -0,006 -3,00
23 Р...Р,Р...С, Р...Н 1,592 0,012 0,011 -0,003 -4,06
24 Р...Р,Р...С 1,520 0,005 0,005 -0,003 -1,52
25 Р...Н,С...О 1,351 0,011 0,010 -0,004 -5,00
26 Р...Р,Р...С 0,748 0,007 0,008 -0,002 -2,37
1 2 3 4 5 6 7 8
27 3-(1-амино- 2,2,2-трифторэтили- ден)-2-имино- 1,1,4,5,6,7- гексафториндан Р...Р, Р...с, с../ы 0,610 0,004 0,004 -0,002 -2,56
28 Р...Р,Р...Н 0,647 0,012 0,010 -0,003 -4,36
29 Р...Р,Р...Н 0,649 0,005 0,006 -0,001 -1,86
30 Р...Р,Р...С 1,472 0,006 0,006 -0,003 -2,18
31 Р...Н 1,679 0,010 0,008 -0,007 -2,60
32 Р...Р,Р...С, К.-Н 0,655 0,008 0,008 -0,002 -4,02
33 Р...Р,Р...Н 0,655 0,012 0,010 -0,003 ^1,36
34 Р...Р,Р...Н 1,121 0,012 0,010 -0,003 -3,44
Окончание табл. 1
1 2 3 4 5 6 7 8
35 3-(1-амино-2,2,2-трифторэтилиден)-1,1,4,5,6,7-гексафто-риндан-2-он О...Н, Р...Н 0,210 0,021 0,017 -0,003 -6,69
36 О...Н,Р...Н 0,920 0,021 0,017 -0,003 -6,76
37 Р...Р,Р...Н 2,138 0,013 0,012 -0,003 -4,34
38 Р...Р 0,983 0,006 0,007 -0,003 -2,86
39 Р...Р 0,613 0,006 0,007 -0,003 -2,85
40 Р...Р,Р...Н 2,144 0,013 0,012 -0,003 —4,34
41 Р...Р,Р...С 4,030 0,005 0,005 -0,002 -3,53
42 Р...Р, Р..ДР...С, о...н,с...с 0,612 0,006 0,005 -0,002 -5,02
43 (2,4-бис(трифторметал)-5,6,7,8,9,9-гексафтор-1,3-диазафлуорен к13221 р...р,р...с 0,828 0,006 0,006 -0,001 -4,94
44 Р...Р, Р...1*,Р...С 0,891 0,005 0,006 -0,002 -3,90
45 Р..Л7,Р...с 0,609 0,005 0,005 -0,003 -1,90
46 Р...Р,Р...С 1,621 0,006 0,005 -0,002 -*,62
47 Р...Р,Р...С 0,868 0,006 0,006 -0,002 -3,12
48 Р...Р, Р...К,Р...С 1,921 0,005 0,005 -0,001 -4,02
49 Р...Р,Р...С 1,579 0,006 0,006 -0,002 -3,82
Обсуждение результатов
Как известно, галогены в молекулярных кристаллах образуют агрегаты, в которых каждый атом имеет в своем ближайшем окружении один или несколько таких же атомов из соседних молекул. В исследованиях авторов [14] описаны следующие структурные особенности галоген-агрегации. Они могут быть конечными -/^агрегаты (двухатомные/(2), трехатомные ДЗ) и т. д.) и бесконечными в одном измерении - /'-агрегаты (цепочечные, ленточные, стержневые), в двух измерених - /2-агрегаты (слои разнообразного строения) или трех измерениях - г3-агрегаты (каркасные). При этом контактирующими считаются атомы, расстояние между которыми превышает удвоенное значение ван-дер-ваальсова радиуса галогена более чем на 0,5 А. Следует отметить, что при описании особенностей агрегации в кристаллах используют следующие характеристики: координационное число атома галогена по галогену и отклонение среднего расстояния до ближайших соседей от удвоенного ван-дер-ваальсова радиуса [15].
В самом простом приближении межмолекулярные взаимодействия в димерах исследуемых нами полифторированных структур можно классифицировать как внутрислоевые, внутристопоч-ные и межстопочные (рис. 1). В первом случае, это взаимодействия, при которых обе молекулы принадлежат одному слою; во втором - одной стопке, в том числе, когда реализуются 7Г-7Г-взаимо-действия. Межстопочные взаимодействия в нашем случае подразумевают пару молекул, принадлежащих соседним стопкам, либо расположенных друг к другу под углом, значимо отличным от нуля. Как демонстрирует рис. 1, а, внутрислоевые взаимодействий могут реализовываться либо за счет водородных связей, что характерно для производных инданона, либо за счет Р...Р-взаи-модействий, присущих любой паре рассматриваемых молекул.
Другой способ классификации подразумевает выделение взаимодействий для атомов фтора, ковалентно связанных с углеродом ароматических колец САг--Р...X, либо входящих в состав метиленовой группы С Р2.. .X, либо трифторметильной группы СР3.. .X, где X = Н, С, N. О — атомы соседней молекулы.
Обсуждая проблему стабильности структуры агломерата или, в нашем случае, димера, «вынутого» из контекста кристаллической упаковки, а также принципы формирования взаимного расположения молекул, нельзя не коснуться вопросов соотношения таких величин, как энергия взаимодействия и «жесткость» контактов в структуре. По мнению авторов [16], «жесткость» контакта определяется не столько энергией межмолекулярного взаимодействия, сколько наличием достаточно высоких потенциальных барьеров, препятствующих существенному изменению их взаимного расположения. Оценку «жесткости» контактов авторы предлагают производить на основе степени геометрических изменений взаимного расположения контактирующих молекул, возникающих при минимизации энергии их взаимодействия.
6)
Рис. 1. Примеры внутрислоевых (а), внутристопочных и межстопочных (б) взаимодействий в димерах
Следует отметить, что величины АВ = Д,. х - /)р...х&У, где Ог хС'у - межъядерные расстояния в кристаллических структурах, установленные РСА, для большинства взаимодействий в нашей выборке являются отрицательными. Это связано с тем, что при оптимизации геометрии в ряде димеров наблюдается существенная перестройка взаимного расположения молекул и контактов, отвечающих таким взаимодействиям, становится больше. Однако, как показывает гистограмма на рис. 2, такая перестройка затрагивает контакты Р...X и в существенно меньшей степени характерна для контактов Б... Б. Изменения длин межмолекулярных контактов при оптимизации для исследуемых нами структур димеров полифторароматических соединений нельзя назвать критическими, поскольку для 70 % контактов выборки значение АО по абсолютной величине не превышает 0,5 А.
-5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 ДО
Рис. 2. Распределение изменений межъядерных расстояний (А) при оптимизации геометрии в сравнении для контактов Р...Р и Р...Х
Мы воспользовались возможностью произвести количественную оценку геометрических изменений, произошедшую в димерах при локализации равновесной геометрии. В среднем по выборке, величины таких отклонений составляют ЯМБО] = 0,84 и Н.М8Г),т = 0,74 А, соответственно. Качественно сопоставление изменений геометрии структур для предельного (4,00 А) и среднего значений ЬШБОы демонстрирует рис. 3. Выявлено, что не обнаруживается предпочтений и тенденций в большей или меньшей степени сохранять жесткость ориентации при сравнении внутри-слоевых, внутристопочных или межстопочных взаимодействий в димерах.
Рис. 3. Примеры изменений геометрии структуры димеров в случае малого значения [ЗМвОм = 0,26 (а) и для максимального наблюдаемого значения = 4,00 (б). Темным оттенком отмечена структура,
наблюдаемая в кристалле, светлым - локализованная равновесная геометрия в димере
Все более популярным в последнее время для описания и оценки межмолекулярных взаимодействий, в том числе, с участием атомов фтора [17], становится использование квантовотопологических характеристик ЭП [18]. Вычисленные значения и характеристики ЭП в критических точках в межмолекулярном пространстве, являются зависимыми от сорта атомов, участвующих в этом связывании. В современных исследованиях [19] приведены количественные выражения характеристик ЭП в зависимости от межъядерных расстояний для разных рядов водородных связей, с меняющимися атомами-акцепторами протона. Для межмолекулярных взаимодействий авторами приводятся эмпирические зависимости, выражающие экспоненциальный характер спадания ЭП в критических точках связей с ростом межъядерного расстояния.
В ряду изучаемых соединений нами исследовано поведение величин электронной плотности р(г/,), плотностей потенциальной и кинетической энергий У(гь), (}(гь) в критических точках, отвечающих взаимодействиям Р...Б, Б...О, Б..,1М, Б...С, Б...Н. На рис. 4 показано индивидуальное для каждого сорта рассматриваемых контактов поведение величин р(г/,), V(гь), 0(гл) от межъядерных расстояний. Наибольший диапазон расстояний от 2,12 до 3,48 А, в котором изменяются наблюдаемые величины, обнаруживается для Р...Р и Р... Н взаимодействий. Важно отметить, что изменение квантово-топологических характеристик ЭП для этих взаимодействий мы наблюдаем в довольно широком общем для них диапазоне межъядерных расстояний: от 2,65 до 3,19 А. Установленные значения коэффициентов в соответствующих уравнениях вида (1)—(3), а также статистические характеристики линейных зависимостей наблюдаемых и вычисленных величин приведены в табл. 2.
а)
б)
р(гь) = а ехр ( Ь Г)Р .х), в(гь)= а ехр ( -/) В1ХХ У(гь)= а ехр ( Ь Ог х)-
(1)
(2)
(3)
0.0.16 о
0.014 %
0.012 * в 1 , °
0.010 » » '
0.00S «?* о ♦
0.006 # §э
0.004 Ф %
0.002 ф о® » о о
0.000
2.0 22 2.4 2.6 2.8 3,0 3 2 3.4 3.6 3.8 4.0 4.2 4.4 4.6 D F...X
а)
* F...H
■ F...C
■ F...N
* F...O
о F...F
б)
в)
Рис. 4. Изменение электронной плотности ЭП (а), плотности потенциальной энергии (б), плотности кинетической энергии электронов (в) в критических точках связей, отвечающих межмолекулярным взаимодействиям С-Р...Х
Значения параметров и статистические характеристики взаимосвязи квантово-топологических характеристик ЭП и межъядерных расстояний в уравнениях
Таблица 2
Контакт, (число контактов) Уравнение, отклик Диапазон расстояний ¥...Х, А a b Коэффициент линейной корреляции І? Диапазон отклонений Д, ат. ед.
Б...Б (125) (1) рЫ 2,65- -3,48 8,991 -2,515 0,958 О о <о o' 1 -0,0007
Р...Н (47) (1) Pin) 2,12- -3,19 0,749 -1,896 0,985 о о о о -0,0024
Б...Б (125) (2) G(rb) 2,65- -3,48 4,639 -2,252 0,972 -0,0009- -0,0007
Р...Н (47) (2) G(rb) 2,12- -3,19 0,576 -1,840 0,980 -0,0011- -0,0018
Б...Б (125) (3) V(n,) 2,65- -3,48 -16,303 -2,746 0,972 1 1 о О 1 о: О: 0\ ; ; 0,0009
Р...Н (47) (3) V(rb) 2,12- -3,19 -0,700 -2,001 0,973 -0,0021- -0,0011
Скорость роста электронной плотности р(гь) и плотности кинетической энергии (Лп) в результате сокращения межъядерного расстояния Б...Б ощутимо выше, чем в случае Р...Н контактов. Наибольшего значения величина р(г/,)р...р = 0,012 ат. ед. достигает при межъядерном расстоянии 2,65 А, на этом же расстоянии величина р(г;,)Р н = 0,005 ат. ед., т. е. почти вдвое ниже. Высокие коэффициенты линейных корреляций исходных и прогнозируемых по уравнениям (1)-(3) значений квантово-топологических характеристик ЭП позволяют обоснованно использовать
предложенные уравнения для экстраполяции данных на новые диапазоны межъядерных расстояний. Диапазоны отклонений Д, представляющие собой разницу найденных и вычисленных по данным уравнениям значений, низки и лежат на границе точности экспериментального определения данных величин.
Предлагаемая методология представляется полезной, поскольку наблюдаемые в кристалле расстояния не всегда сохраняются в процессе локализации равновесной геометрии. Таким образом, полученные величины квантово-топологических характеристик ЭП могут быть скорректированы с учетом межмолекулярных расстояний в кристалле. Полученные значения величин, в свою очередь, могут быть полезны для более точного описания характеристик взаимодействий и восстановления взаимосвязей с полезными экспериментальными данными.
Значения вычисленных энергий взаимодействий в димерах &ЕШег лежат в диапазоне от -0,73 до -10,44 ккал/моль. Данные расчеты велись без коррекции ошибки суперпозиции базисного набора ВБЗЕ [20]. Наиболее низкие отрицательные значения энергий взаимодействий принадлежат димерам с типичными водородными связями О...Н, ]Ч...Н. В среднем, для них величина межъядерных расстояний Д-) Н в нашей выборке составляет 1,96±0,08 А. В случае самых энергетически выгодных взаимодействий с участием атомов фтора Р...Р, Р...Н, энергия АЕт1сг достигает -5,69 ккал/моль. Наикратчайший наблюдаемый в данных структурах контакт Р...Р составляет 2,84 А и принадлежит взаимодействиям Сл,—Р.. .Р—САг двух соседних фторарильных колец. Наикротчайшее Р.. Т расстояние не принадлежит димеру с наименьшей энергией взаимодействия.
Полученные величины хорошо соотносятся с описанными в литературе данными о расчетных и экспериментально полученных энергиях взаимодействий с участием фтора. Сравнение энергий межмолекулярных взаимодействий с помощью димерного подхода для арилнезамещен-ных и перфторированных структур дисульфиниламиноаренов, выполненное нами ранее, дает тот же диапазон значений [21]. В работе [22] приводятся расчетные значения энергии взаимодействия для заряженных систем (СбРб)г и (СбРб)г+, в том числе при различной взаимной ориентации префторбензольных колец. Так, например, для стопочных димеров с симметрией Ду„ стопочных со смещением (С2и) и Т-образных димеров (С2„) вычисленные значения энергий изменяются от -8,9 до -13,6 ккал/моль. Эти данные сопоставляются с экспериментальными значениями энтальпий ассоциаций, измеренных [23] и лежащих для этих же систем в диапазоне от -7,2 до -10,4 ккал/моль. В этой же работе для незаряженных систем приводятся менее высокие значения расчетных энергий взаимодействия, изменяющиеся в зависимости от взаимной ориентации перфторированных колец бензола или нафталина от -1,60 до -5,82 ккал/моль. Малые значения полных энергий для систем со слабыми взаимодействиями Н3С-Р...Н-СНэ обстоятельно обсуждаются в работе [24]. Авторами приводится значение -0,46 ккал/моль, вычисленное с поправкой на В88Е в приближении МР2/а1щ-сс-рУТ2 при существенном для данного сорта связи межъядер-ном расстоянии Р...Н 2,56 А. Сравнительно низкие значения энергии стабилизации в комплексах со множественными взаимодействиями Р...Н приводятся в работе [25]. Вычисленные в приближении МР2/аг^-сс-р’УТ2, энергии стабилизации в комплексах КРЬР СНР),, и ЫНРг-ШР'),,, где п = 2 или 3, а фтор - донор протона, составляют от —4,1 до -9,4 ккал/моль. Этот результат наилучшим образом согласуется с диапазоном полученных нами значений энергий, при этом также для систем, включающих от трех и более межмолекулярных контактов Р.. .Н и Р.. .Р.
Отметим, что вычисленная величина энергии взаимодействия для рассматриваемых структур никак не связана с типом димера, определяющим взаимное расположение молекул в паре. Например, межстопочные взаимодействия могут включать водородную связь и характеризоваться более низкими значениями АЕт1ег, а взаимодействия, соответствующие внутрислоевым - более высокими, чем внутристопочные. Таким образом, вероятно, энергия взаимодействия для пары молекул в большей степени определяется набором атом-атомных сокращенных контактов, по сравнению с типом мотива и укладки в кристалле, из контекста которого была выделена рассматриваемая пара молекул. Также установлено, что величины ЯМБО, и КМ5Е),пЬ призванные косвенно характеризовать жесткость контактов, количественно не связаны с вычисленными энергиями взаимодействий в димерах АЕШег.
Для выявления наличия количественной взаимосвязи локальных характеристик электронной плотности для взаимодействий Р...Х и Х...Н, где X = Р, О, N с величинами энергий АЕт1ег в димерах проанализированы суммы, средние и максимальные значения величин р(гД У(гь), 0{гь).
Суммирование или выбор максимального значения велись по количеству критических точек ЭП в межмолекулярном пространстве каждого димера, т. е. среди всех сортов наблюдаемых контактов. В результате регрессионного анализа установлено, что наиболее высокий коэффициент корреляции наблюдается для взаимосвязи Ет,ег и максимального из наблюдаемых в каждом димере значения плотности кинетической энергии электронов С{гь)тйХ в критических точках связей, отвечающих межмолекулярным взаимодействиям. Величина коэффициента корреляции составила 0,83, графическое выражение зависимости приведено на рис. 5. Среди локальных плотностей кинетической энергии 0(гь) естественно выделяются значения, соответствующие типичным водородным связям О...Н, N...11. Такие связи характеризуются низкими значениями АЕт1ег и сравнительно высокими значениями С(гь). Водородные связи с участием фтора Р...Н имеют значения С(гь) в диапазоне от 0,006 до 0,012 ат. ед. Плотность локальной кинетической энергии в критических точках взаимодействий Р.-.Б характеризуется еще более низким диапазоном: от 0,004 до
0,008 ат. ед. В это же время, диапазон энергий взаимодействия АЕШег в димерах для Р..Т и Б.-.Н контактов остается общим.
^(Гь)тах
Рис. 5. Взаимосвязь величин энергий взаимодействий в димерах (ккап/моль) и максимальных значений плотностей кинетической энергии электронов (ат. ед.) в критических точках в межмолекулярном пространстве
Следует признать, что количественная взаимосвязь топологических характеристик ЭП с энергией взаимодействия АЕ,тег прослеживается основным образом для димеров с водородной связью. Энергия АЕт1сг в димерах, включающих только Р...Р контакты, становится зависимой от других факторов: например, от количества таких контактов между соседними молекулами. Этот тезис можно проиллюстрировать на примере сравнения взаимодействий, включающих до десятка Р...Р контактов в димере З-метиламино-2-метилиминогексафториндана, где АЕт,сг= -2,16 ккал/моль и трех контактов Р...Р в димере 2-амино-З-гексафтор-инденил-трифторметилкетона с АЕЫег = -0,73 ккал/моль. В последнем случае мы наблюдаем слабые взаимодействия между триф-торметильной и дифторметиленовой группами. Тем не менее, здесь сравнительно длинные контакты Р.. Т сохраняются в ходе оптимизации и характеризуется этот димер, как имеющий самую высокую энергию в выборке. Таким образом, количественное описание энергии взаимодействия в димерах, основанное на использовании локальных свойств ЭП, должно уточняться для каждого сорта рассматриваемого взаимодействия и, вероятно, не только на уровне параметров линейных зависимостей.
Заключение
На примере димерного подхода для ряда полифторированных молекулярных кристаллов показано, что проблема корректного учета электронной плотности в критических точках межмолекулярных взаимодействий Р...Р и Р...Н, связанная с изменением геометрии при оптимизации, может быть решена путем использования количественных зависимостей топологических характеристик электронной плотности от межъядерных расстояний. Параметры в уравнениях количественных зависимостей могут быть установлены для каждого вида химической связи или межмо-
лекулярного взаимодействия. Предложены виды таких уравнений для F...F и F...H контактов в молекулярных кристаллах полифторированных енаминоиминоинданов, енаминоиминоинденов и диазафлуоренов.
Выявлено, что вычисленные величины энергии взаимодействия в димерах не связаны с «жесткостью» контакта и типом димера - расположением пары молекул в одной плоскости или принадлежащих одной стопке в кристалле.
Продемонстрирована взаимосвязь между энергией взаимодействия в димере, которая изменяется в рассмотренном ряду от -0,73 до -10,44 ккал/моль и максимальными значениями энергетических топологических характеристик экстремумов электронной плотности, отвечающих меж-молекулярным взаимодействиям F.. .F, F.. .Н и О.. .Н, N.. .Н.
Таким образом, димерный подход, при котором для оценки межмолекулярных взаимодействий рассматривается вся совокупность пар соседних в кристалле молекул, позволяет получить представления об энергетических характеристиках взаимодействий в молекулярном кристалле и не является препятствием для корректного учета топологических характеристик электронной плотности при использовании процедуры оптимизации геометрии димера.
Литература
1.Leroux, F. a-Fluorinated Ethers, Thioethers, and Amines: Anomerically Biased Species / F. Leroux, P. Jeschke, M. Schlosser // Chemical Reviews. - 2005. - Vol. 105 (3). - P. 827-856.
2. Reichenbacher, K. Fluorine in crystal engineering - “the little atom that could” / K. Reichen-bacher, H.I. Suss, J. Hulliger // Chem. Soc. Rev. - 2005. - Vol. 34. - P. 22-30.
3. Зоркий, П.М. Критический взгляд на основные понятия химии. Химическое и физическое взаимодействие атомов. Специфические молекулярные контакты / П.М. Зоркий // Рос. хим. журн. -1996.-Т. 40.-С. 5-25.
4. Roy Choudhury, A. How Realistic Are Interactions Involving Organic Fluorine in Crystal Engineering? Insights from Packing Features in Substituted Isoquinolines / A. Roy Choudhury, T.N. Guru Row // Crystal Growth & Design. - 2004. - Vol. 4. - № 1. - C. 47-52.
5. Dunitz, Jack. D. Molecular Pair Analysis: C-H...F Interactions in the Crystal Structure of Fluoro-benzene and Related Matters / Jack. D. Dunitz, W. Bemd Schweizer // Chem. Eur. J. - 2006. - Vol. 12. -P. 6804-6815.
6. Bader, F.W. AIM2000 Program Package Ver. 2.0 / F.W. Bader, F. Biegler-Konig, J. Schonbohm. Bielefield, Germany: 2002. - University of Applied Sciences.
7. Determination of covalent bond orders and atomic valence indices using topological features of the experimental electron density / V.G. Tsirelson, E.V. Bartashevich, A.I. Stash, V.A. Potemkin // Acta Cryst. - 2007. - Vol. 63. - P. 142-150.
8. Hydrogen bond strengths revealed by topological analyses of experimentally observed electron densities / E. Munoz, A.M. Benito, L.C. Estepa et al. // Chem. Phys. Lett. - 1998. - Vol. 285, № 3. -P170-173.
9. Синтез, молекулярная и кристаллическая структура полифтор-1,3-диазафлуоренов / В.М. Карпов, В.Е. Платонов, И.П. Чуйков и др. // Журн. орган, химии. - 2004. - Т. 40, № 3. - С. 448-452.
10. Рентгеноструктурное исследование монокристаллов комплексов 3-( 1-амино-2,2,2-три-фторэтилиден)-2-имино-1,1,4,5,6,7-гексафториндана с диоксаном, пиразином, пиридином / Т.В. Рыбалова, В.М. Карпов, В.Е. Платонов, Ю.В. Гатилов // Журн. структур, химии. - 2007. -Т. 48, №2.-С. 318-324.
11. Синтез, молекулярная и кристаллическая структура 3-(1-амино-2,2,2-триторэтилиден)-2-имино-1,1,4,5,6,7-гексафториндана и его N-метил- и N,N’-диметилпроизводных / В.М. Карпов,
B.Е. Платонов, Т.В. Рыбалова, Ю.В. Гатилов // Журн. структур, химии. - 2006. - Т. 47, № 3. -
C. 532-539.
12. Granovsky, А. -2003. -http://classic.chem.msu.su/gran/gamess/index.html
13. In Atoms in Molecules: A Quantum Theory; The International Series of Monographs of Chemistry / R.F.W. Bader, J. Halpen, M.L.H. Green et al. - Clarendon Press: Oxford, UK, 1990.
14. Гринева, O.B. Изоструктурность и неизоструктурность в рядах галогенорганических кристаллических веществ. Строение галоген-агрегатов / О.В. Гринева, П.М. Зоркий // Журн. структур. химии. - 2001. - Т. 42, № 1. - С. 21-30.
15. Гринева, О.В. Энергия межмолекулярного взаимодействия в кристаллах хлорсодержащих органических соединений / О.В. Гринева, П.М. Зоркий // Журн. физ. химии. - 1998. - Т. 72, № 4. -С.714-720.
16. Исследование устойчивости молекулярных агломератов в молекулярных кристаллах / А.В. Малеев, Б.Б. Седов, И.К. Житков, В.Г. Рау // Журн. структур, химии. - 2007. - Т. 48, № 1. —
С.148-152.
17. About the evaluation of the local kinetic, potential and total energy densities in closed-shell interactions / E. Espinosa, I. Alkorta, I. Rozas et al. // Chem. Phys. Lett. - 2001. - Vol. 336. - P. 457-461.
18. Checinska, L. F-H...F-C hydrogen bonds - The influence of hybridization of carbon atom connected with F-acceptor on their properties / L. Checinska, S.J. Grabowski // Chemical Physics. - 2006. -Vol. 327. - P. 202-208.
19. Universal Features of the Electron Density Distribution in Hydrogen-Bonding Regions: A Comprehensive Study Involving H...X (X = H, C, N, O, F, S, Cl, n) Interactions / I. Mata, I. Alkorta, E. Molins, E. Espinosa // Chem. Eur. J. - 2010. - Vol. 16. - P. 2442-2452.
20. Gianinetti, E. Modification of Roothan Equations to Exclude BSSE from Molecular Interaction Calculations / E. Gianinetti, M. Raimondi, E. Tomaghi // J. Quantum Chem. - 1996. - Vol. 60. - P. 157.
21. Теоретическое исследование межмолекулярных взаимодействий в кристаллах замещенных дисульфиниламиноаренов / Е.В. Барташевич, Д.К. Никулов, Е.А. Шманина, А.И. Локтев // Вестник ЮУрГУ. Серия «Химия». - 2009. - Вып. 1. - № 12. - С. 10-17.
22. Lorenzo, S. Supramolecular potentials and embraces for fluorous aromatic molecules /
S. Lorenzo, G.R. Lewis, la. Dance // New J. Chem. - 2000. - Vol. 24. - P. 295-304.
23. Hiraoka, K. Stabilities and structures of hexafluorobenzenes (C6,F6XC6F6) and С6[-64(С6Р6)) / К. Hiraoka, S. Mizuse and S. Yamabe // J. Phys. Chem. - 1990. - Vol. 94. - P. 3689.
24. D’Oria, E. On the hydrogen bond nature of the C-H...F interactions in molecular crystals. An exhaustive investigation combining a crystallographic database search and ab initio theoretical calculations / E. D’Oria and JJ. Novoa // CrystEngComm. - 2008. - Vol. 10. - P. 423-436.
25. Karpfen, A. The interaction of fluoramines, fluorophosphines and fluoroarsines with hydrogen fluoride clusters (HF)n: Model studies on blue-shifted hydrogen bonds / A. Karpfen // Journal of Molecular Structure: THEOCHEM. - 2005. - Vol. 757. - P. 203-215.
Поступила в редакцию 6 сентября 2010 г.
DIMERIC APPROACH IN DESCRIPTION OF INTERMOLECULAR INTERACTIONS IN CRYSTALS OF POLYFLUORINATED STRUCTURES
Theoretical research of intermolecular interactions with participation of atoms of fluorine for poly-fluorinated diazafluoren, enaminoindens and enaminoindans which forming molecular crystals is fulfilled. Energy of intermolecular interactions in pairs of neighbor molecules in a crystal is estimated. Correlation of interactions energy in dimers with quantum-topological characteristics of electronic density in bond critical points for F.. .F, F.. .H and О.. .H contacts is executed.
Keywords: polyfluorinated structures, molecular crystals, intermolecular interactions, quantum-topological analyzes of electron density, theory AIM.
Bartashevich Ekaterina Vladimirovna - PhD (Chemistry), Associate professor, Organic Chemistry Subdepartment, South Ural State University. 76, Lenin avenue, Chelyabinsk, 454080.
Барташевич Екатерина Владимировна - кандидат химических наук, доцент, кафедра органической химии, ГОУ ВПО «Южно-Уральский государственный университет». 454080, г. Челябинск, проспект им. В.И. Ленина, 76.
E-mail: kbartash@yandex.ru
Shmanina Elena Alexandrovna - Posrgraduate Student, Physical Chemistry Subdepartment, South Ural State University. 76, Lenin avenue, Chelyabinsk, Russia, 454080.
Шманина Елена Александровна - аспирант, кафедра физической химии, ГОУ ВПО «Юж-но-Уральский государственный университет». 454080, Челябинск, проспект им. В. И. Ленина, 76.
E-mail: shmelenok_1808@mail.ru