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

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

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

Ранее неоднократно изучались строение, конформационные переходы таутомеров 1,3 - дикетонов и их производных различными экспериментальными и теоретическими квантовохимическими полуэмпирическими и неэмпирическими методами [1], [2], [3]. Между тем, наличие такого атома третьего периода, как сера, может оказывать существенное влияние как на сами свойства таутомеров 1,3 - тиокарбонильного соединения, так и на процесс таутомеризации в целом. В частности, помимо енольных форм (V, VI, VII), возможна также енолизация исходного тиодикетона I по атому серы, с образованием тиольных таутомеров II, III, IV (рис. 1).

Рис. 1. Конформации и таутомерные формы тиоацетилацетона

тиоацетилацетон конформер таутомер квантохимический

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

В настоящей работе ставилась задача исследования неэмпирическими методами квантовой химии геометрического строения конформеров тиоацетилацетона, как простейшего представителя 1,3 - тиокарбонильного соединения, связь конформационных форм его таутомеров с энергией водородной связи, а также изучалось влияние расчетного базиса и полученной в ходе различных приближений энергии электронной корреляции на расчетную величину энергии молекулы и ее пространственное строение. Все расчеты проводились в корреляционно-согласованном базисе cc-pVDZ, и улучшенном aug-cc-pVDZ. Для учета корреляционных поправок использовался метод возмущений Меллера-Плессета второго порядка (MP2) и метод функционала электронной плотности (DFT) с применением гибридного функционала B3LYP. Для расчетов использовалась программа MPQC [4] гибридного функционала B3LYP. Для расчетов использовалась программа MPQC [4]

Рис. 2. Структура исходного тиокетона I

Рис. 3. Структура енольного таутомера II

Структуры таутомерных молекул и их конформеров оптимизированные в расчетах b3lyp/cc-pVDZ представлены на рис. 2-8. Для вывода графических моделей использовалась программа Gabedit [5].

Следует отметить, что структурные параметры молекул, полученные в различных методах различаются незначительно, однако видно (табл. 1), что расширение базиса cc-pVDZ до aug-cc-pVDZ приводит к некоторому завышению межатомных расстояний, тогда как расчеты с учетом электронной корреляции напротив приводят к их заметному сокращению.

Рис. 4. Структура енольного таутомера III

Рис. 5. Структура енольного таутомера IV

Рис. 6. Структура енольного таутомера V

Рис. 7. Структура енольного таутомера VI

Таблица 1. Равновесные межатомные расстояния S-O, Е

Рис. 8. Структура енольного таутомера VII

Длина водородной связи в тех таутомерах, где она может образовываться (II и V), также обнаруживает те же тенденции при расширении базиса и учете электронной корреляции (табл. 2). Кроме того, вполне ожидаемо ее увеличение для атома серы, однако оно крайне незначительно, и не превышает 0,1 - 0,15 Е.

Таблица 2. Ддина водородной связи, Е

Из таблицы 3, в которой приведены относительные энергии рассчитанных соединений, видно, что учёт корреляционных эффектов оказывает весьма существенное влияние на результат расчета энергетических характеристик молекул. Так, если расширение числа базисных функций со 142 (базис cc-pVDZ) до 237 (базис aug-cc-pVDZ) приводит к выигрышу порядка 40 - 45 кДж/моль, то энергия корреляции, полученная методом возмущений MP2, составляет более трех с половиной, а метод функционала электронной плотности показывает корреляционную энергию более шести тысяч кДж/моль. Очевидно, что величины такого масштаба не могут не оказывать влияния на качественный результат проводимых расчетов, что и наблюдается при сравнении относительной энергии таутомеров (табл. 4). В частности, одноэлектронные расчеты в базисе cc-pVDZ указывают, что енолизация по по атому серы (структура II) энергетически выгодней енолизации по кислороду (структура V). Улучшение базиса до aug-cc-pVDZ хотя качественно меняет этот результат, однако количественный энергетический выигрыш в этом случае невелик. В то же время, учет энергии корреляции приводит к явному предпочтению енолизации тиокетона по кислороду, а таутомеры VI и VII являются энергетически невыгодными. Также можно видеть, что енолизация по атому серы без образования водородной связи (структуры III и IV) также не приводит к сколько-нибудь значительному энергетическому выигрышу по сравнению с неенолизованой формой (I). К настоящему времени разработано и продолжает разрабатываться [см. напр. 7] множество методов учета электронной корреляции в рамках различных приближений, однако предпочтительность того или иного метода пока может оцениваться только на основе сравнительного анализа с экспериментальными данными.

Таблица 3. Энергия электронной корреляции таутомерных и енольных форм тиоацетилацетона, кДж/моль

Несмотря на то, что в виду таких особенностей ab initio расчетов, как учет многоцентровых взаимодействий не существует достаточно строгого метода, позволяющего выделить отдельно двухатомную энергию и энергию связи в частности, тем не менее, расчет в данном случае позволяет сделать приблизительную оценку энергии водородной связи, сравнивая энергию соответствующих таутомеров с водородной связью и без нее. Так, по данным DFT-расчетов конформеры II и III различаются приблизительно на 13 кДж/моль, а V и VI на более, чем 70 кДж/моль (табл. 4). Однако такой неожиданный результат может объясняться не только более высокой прочностью водородной связи S-H, но и относительно невыгодным взаимным расположением других атомов, в частности, близким расположением атомов кислорода и серы и, вследствие этого, их взаимным отталкиванием.

Таблица 4. Относительная молекулярная энергия таутомеров тиоацетилацетона, кДж/моль

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

Список литературы

1. Шокова Э.А. 1,3 - дикетоны. Синтез и свойства / Э.А. Шокова, Дж.К. Ким, В.В. Ковалев // Журнал органической химии. - 2015. - Т. 51. - Вып. 6 - С. 773-847.

2. Иванов Ю.В. Исследование таутомерии иминопроизводных оксиацетофенона неэмпирическими методами / Ю.В. Иванов // Современные проблемы химической науки и образования, Чебоксары, 2012. - 19-20 апреля. - С. 84-85.

3. Кострюкова Т.С. Получение фторированных в-дикетонов / Т.С. Кострюкова, Н.В. Васильев // Вестник Московского государственного областного университета серия «Химия». - - №1. - С. 150-156.

4. The Massively Parallel Quantum Chemistry Program (MPQC), Version 2.3.1 / Curtis L. Janssen, Ida B. Nielsen, Matt L. Leininger, Edward F. Valeev, Joseph P. Kenny, Edward T. Seidl, Sandia National Laboratories, Livermore, CA, USA, 2008.

5. Allouche, A.-R. Gabedit. A graphical user interface for computational chemistry softwares // Journal of Computational Chemistry. - 2011 - №32 - 174-182.

6. Аникин Н.А. Новая аппроксимация DFT для быстрых расчетов био- и наноструктур / Н.А. Аникин, М.Б. Кузьминский // Информационно-вычислительные технологии в науке - 2011, Сб. Тезисов. - М.: 2011 - С. 29-30.

Размещено на Allbest.ru