Государственное образовательное учреждение
высшего профессионального образования
КАЗАНСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
им В. И. Ульянова-Ленина
На правах рукописи
ЮЛЬМЕТОВ Айдар Рафаилевич
СТРУКТУРА И МАГНИТНОРЕЗОНАНСНЫЕ ПАРАМЕТРЫ
МОЛЕКУЛЯРНЫХ СИСТЕМ НА ОСНОВЕ МЕТОДОВ
МОЛЕКУЛЯРНОЙ МЕХАНИКИ, КВАНТОВОЙ ХИМИИ
И СПЕКТРОСКОПИИ ЯМР
01.04.07 — физика конденсированного состоянияАВТОРЕФЕРАТ
диссертации на соискание учёной степени кандидата физико-математических наук
Казань —
Работа выполнена на кафедре общей физики и в лаборатории исследования структуры органических соединений Государственного образовательного учреждения высшего профессионального образования «Казанский государственной университет им. В. И. Ульянова-Ленина»
Научный руководитель: доктор химических наук, профессор Аминова Роза Мухаметовна
Научный консультант: доктор химических наук, профессор Аганов Альберт Вартанович
Официальные оппоненты: доктор физико-математических наук, профессор Жидомиров Георгий Михайлович доктор физико-математических наук, профессор Овчинников Игорь Васильевич
Ведущая организация: Ростовский государственный университет
Защита состоится «2» ноября 2006 г. в «1430 » на заседании Диссертационного совета Д 212.081.15 при Казанском государственном университете им. В. И. Ульянова-Ленина по адресу: 420008, г. Казань, ул. Кремлёвская, 18, КГУ.
С диссертацией можно ознакомится в научной библиотеке им. Н. И. Лобачевского Казанского государственного университета.
Автореферат разослан « » сентября 2006 г.
Учёный секретарь диссертационного совета, доктор физико-математических наук, профессор М.В. Ерёмин
Общая характеристика работы
Актуальность темы исследования. Актуальными задачами современной химической и молекулярной физики и молекулярной биологии являются установление строения и изучение свойств природных и синтетических супрамолекулярных систем. Комплексное использование данных спектроскопии ЯМР высокого разрешения для растворов и твёрдых тел (метод кросс-поляризации магнитно разбавленных спинов в сочетании с вращением образца под магическим углом, ЯМР КП/ВМУ), а также результатов теоретических расчётов констант ядерного магнитного экранирования (ЯМЭ) позволяет значительно расширить круг решаемых задач. В частности, становится возможным сравнительный анализ структуры молекулы в растворе и в твёрдом теле, что позволяет выявить характер влияния ближнего и дальнего окружения молекулы на её строение и свойства и установить определённые закономерности такого влияния.
Интерпретация экспериментально наблюдаемых магнитнорезонансных параметров супрамолекулярных систем представляет собой весьма сложную задачу, поскольку включает в себя необходимость выявления и учёта слабых нековалентных как внутри-, так и межмолекулярных взаимодействий (вандерваальсовы взаимодействия, нетрадиционные водородные связи и т.д.), которые, как стало очевидно к настоящему времени, оказывают существенное влияние на формирование новых физико-химических свойств молекул. Константы ЯМЭ, особенно для тяжёлых ядер, в супрамолекулярных системах не могут быть корректно описаны аддитивным учетом магнитно-анизотропных и электрических вкладов удаленных групп, как это было принято ранее согласно работе Сайки и Слихтера (Saika A, Slichter C. P., J. Chem. Phys., 1985,V. 22, p. 26). В таких супрамолекулах необходимо принимать во внимание вариации электронного распределения, обусловленные всеми типами нековалентных взаимодействий во всей молекулярной системе в целом. Большое количество опубликованных в научной литературе работ посвящено выявлению влияния нековалентных взаимодействий на энергетику образования молекулярных кластеров, однако гораздо меньше внимания уделено анализу влияния таких взаимодействий на магнитнорезонансные параметры (в частности, на константы ЯМЭ). Существующий термин «эффекты упаковки» понимается разными исследователями неоднозначно и, в целом, означает некий феноменологический подход к описанию отличий между экспериментальными результатами для разных фаз вещества без объяснения природы таких отличий.
Целью работы является выработка комплексного подхода к анализу экспериментально наблюдаемых констант ядерного магнитного экранирования (ЯМЭ) 13 C и 1 H для изучаемого класса соединений в растворе и кристаллической фазе. Выполнение работы потребовало разработки метода, позволяющего моделировать наиболее оптимальные (с энергетической точки зрения) возможные пространственные структуры молекулярного кристалла и сольватной оболочки молекулы и проведения на основе полученных геометрических параметров кластеров квантовохимического расчёта констант ЯМЭ.
Объекты исследования. Для решения поставленной задачи выбран ряд соединений, имеющих в своём составе фенильные фрагменты, вандерваальсовы взаимодействия с участием которых в значительной мере стабилизируют супрамолекулярные системы.
Научная новизна. C целью адаптации метода молекулярной механики к задачам данного исследования предложены потенциалы для учёта вкладов дальнодействующих взаимодействий, эффектов диполь-дипольных взаимодействий, специфических взаимодействий типа слабых водородных связей, а также учёта симметрии пространственной ячейки и периодичности кристаллической решётки. Предложены новые алгоритмы и разработана программа для расчётов геометрической структуры молекулярных кристаллов. На основе анализа различных алгоритмов оптимизации предложены методы для улучшения сходимости результатов и ускорения вычислений.
Впервые предложен способ моделирования элементарной ячейки молекулярного кристалла на основе анализа возможных структур, соответствующих наиболее распространённым пространственным группам симметрии.
Впервые предложен способ и определены граничные условия для сшивки результатов расчёта поверхностей потенциальной энергии при одновременном учёте эффектов связанных с образованием водородных связей и вандерваальсовых взаимодействий (как внутри-, так и межмолекулярных).
Впервые применён комплексный подход для расчёта параметров спектров ЯМР 13 C и 1 H в кристаллических образцах и в растворах, сочетающий результаты предложенного модифицированного метода молекулярной механики с методами квантовой химии.
Проведены систематические вычисления констант ЯМЭ молекулярных систем в газовой фазе и молекулярном кристалле, а также молекулярных кластеров, включающих молекулы растворителя. Проведены квантовохимические расчёты констант ЯМЭ с учётом влияния различных растворителей в рамках поляризованной континуальной модели.
Впервые для вандерваальсовых кластеров выявлены закономерности, характеризующие изменение констант ЯМЭ и электронного распределения в зависимости от межмолекулярного расстояния. Установлены закономерности в изменении констант ЯМЭ от вариации распределения зарядов.
Практическая новизна и значимость. Предложенный комбинированный подход, позволяющий на основе методов молекулярной механики и квантовомеханических методов моделировать структуру молекулярных кристаллов и вычислять теоретические значения констант ЯМЭ, может оказаться полезным при анализе спектров ЯМР КП/ВМУ. Метод построения молекулярных кластеров конечного размера, моделирующих структуру сольватных оболочек и других сложных молекулярных систем может быть востребован при анализе строения супра- и наномолекулярных систем.
Обоснованность и достоверность полученных результатов обеспечивалась комплексным характером выполненных исследований, использованием современных экспериментальных и теоретических методов и концепций с постоянным сравнением полученных данных с имеющимися в литературе. Результаты исследований неоднократно представлялись на отечественных и международных научных конференциях, где они анализировались на предмет достоверности.
На защиту выносятся положения, сформулированные в выводах.
Апробация работы. Основные результаты докладывались и обсуждались на VIII, IX, X, XI Всероссийских конференциях «Структура и динамика молекулярных систем» (Яльчик, 1999, 2000, 2001, 2002), Молодёжной научной школе «Актуальные проблемы магнитного резонанса и его приложений» (Казань, 2002), III Всероссийской конференции «Новые достижения ЯМР в структурных исследованиях» (Казань, 2000), I Научной конференции молодых ученых, аспирантов и студентов «Материалы и технологии XXI века», (Казань, 2000), международном семинаре «Современное развитие магнитно-резонансной томографии спектроскопии. Физические основы и применение в медицине и биологии» (Казань, 2001), I и III международных конференциях по новым технологиям и приложениям физико-химических методов для изучения окружающей среды (Ростов-наДону, 2001, 2005), VI и VII международных семинарах по магнитному резонансу (спектроскопия, томография и экология) (Ростов-на-Дону, 2002, 2004), Юбилейной научной конференции физического факультета КГУ (Казань, 2004), Международной конференции «Nuclear Magnetic Resonance in Condensed Matter» (Санкт-Петербург, 2005).
Публикации результатов исследований По материалам диссертации опубликованы 18 работ, из них 4 статьи в центральной печати, 3 — в сборниках статей, 11 — тезисы докладов.
Структура диссертации Диссертационная работа состоит из введения, четырёх глав, выводов, приложения и списка литературы (150 наименований).
Во введении обоснована актуальность работы, рассмотрены цели и задачи исследования, изложены структура и содержание работы.
В первой главе дано описание метода ЯМР кросс-поляризации разбавленных спинов и вращения образца под магическим углом (КП/ВМУ).
В первой главе также приведён обзор научных публикаций, посвящённых применению метода ЯМР КП/ВМУ для изучения структуры молекулярных кристаллов органических соединений. Анализ литературных данных показывает, что в настоящее время методом ЯМР высокого разрешения возможно изучение структуры широкого класса органических молекул в кристаллической фазе. Тем не менее, анализ полученных данных ЯМР для твердотельных образцов часто затруднён из-за усложнённой, по сравнению со спектрами жидкостей, структуры спектров. Также можно сделать вывод, что для корректной интерпретации особенностей кристаллической структуры обязательным является учёт межмолекулярных взаимодействий, как специфических (водородные связи), так и неспецифических (вандерваальсовы взаимодействия).
Во второй главе рассматриваются предложенные изменения в методе молекулярной механики (ММ), позволившие провести моделирование не только структуры изолированных молекул, но и структуры молекулярных кристаллов, а также сольватных оболочек.
Метод ММ был нами дополнен и усовершенствован для учёта специфических межмолекулярных взаимодействий, в частности, для одновременного учёта вандерваальсовых взаимодействий и водородных связей, а также для учёта несвязанных взаимодействий в кристаллической решётке.
С целью адекватного описания межмолекулярных вандерваальсовых взаимодействий нами была использована потенциальная функция Бэкингема, т.к. при рассмотрении плотноупакованных молекулярных структур, в частности, молекулярных кристаллов, становится важным более точное описание энергии взаимодействия на малых расстояниях:
Для учёта водородных связей нами был использован потенциал вида который широко применяется для моделирования структуры различных биологических объектов, в частности, протеинов.
Потенциал вандерваальсовых взаимодействий (1) на расстояниях, меньших 2, 5 A описывает сильное отталкивание между атомами, в то время как именно в этой области могут образовываться водородные связи с локальным минимумом энергии взаимодействия (2) на расстояниях порядка 1, 5 2, 0 A. Поэтому, при разработке предложенного модифицированного метода MM, специально вводились граничные условия, позволяющие одновременно учитывать эти эффекты.
«Сшивка» потенциалов в нашей программе производилась путём использования дополнительной взвешивающей функции s(r), такой, что:
где rHbond и rV dW — расстояния между атомами, соответствующие минимуму потенциальной энергии для потенциалов (2) и (1), соответственно.
В качестве функции, удовлетворяющей указанному условию, был использован кубической полином, коэффициенты которого находили исходя из граничных условий (3). Полученный в результате «сшивки» потенциал для несвязанных межатомных взаимодействий можно записать в виде:
Enonbonded = s(r) · EHbond (r) + (1 s(r)) · EV dW (r). (4) Молекулярный кристалл нами рассматривался в виде совокупности кристаллических ячеек, внутри каждой из которых помещались одна, две или несколько изучаемых молекул в зависимости от выбранной пространственной группы симметрии ячейки. Положение атомов внутри каждой такой элементарной ячейки задавалось декартовыми координатами, отсчитываемыми относительно её геометрического центра. Проводился учёт взаимодействия между молекулами, находящимися в разных ячейках, координаты атомов которых получались трансляционными преобразованиями из координат атомов исходной ячейки. Также была добавлена возможность моделирования кристаллических структур для заданных пространственных групп симметрии ячейки. Поскольку все относительные внутримолекулярные расстояния для молекул, связанных преобразованиями симметрии, идентичны, оказалось возможным лишь однократное вычисление внутримолекулярных (всех связанных и части несвязанных) взаимодействий для каждой такой группы молекул, что позволяет значительно ускорить ход вычислений.
В второй главе также приводится краткий обзор методов расчёта констант ЯМЭ и обосновывается выбор метода градиентно-инвариантных атомных орбиталей (ГИАО) с использованием волновых функций, рассчитанных в рамках метода функционала плотности (DFT) на уровне UB3LYP/6– 31G(d,p).
Таблица 1. Размеры элементарных ячеек, соответствующих разным пространственным группам симметрии, рассчитанных для молекулы нафталина (линейные размеры — A, углы —, объём ячейки — A3 ), а также данные РСА.
В третьей главе приведены результаты апробации предложенного в гл. 2 метода на примере ряда соединений (I, II) для которых имелись литературные данные о кристаллической структуре, полученные методом РСА, а также результаты применения разработанного подхода для анализа спектров ЯМР поликристаллических образцов соединений (III), (IV) и (V).
несколько возможных структур элементарной ячейки кристалла с различными пространственными группами симметрии и разным количеством независимых молекул в ячейке. Геометрические характеристики полученных кристаллических структур приведены в табл. 1, из которой I видно, что наибольшей плотностью (наименьшим объёмом элементарной ячейки, приходящейся на одну молекулу V /Z) обладают структуры, содержащие по две молекулы в элементарной ячейке и с симметрией ячейки P 1 (V /Z 179, 3 A3 ) и P 21 (V /Z 177, 0 A3 ).
Такие результаты хорошо согласуются с литературными данными РСА — из эксперимента следует, что кристалл нафталина обладает симметрией P 21. Рассчитанные размеры ячейки для группы P 21 также близки к экспериментально измеренным (табл. 1).
Энергетические параметры рассчитанных структур приведены в табл. 2.
Из анализа полученных значений можно сделать вывод, что наиболее стабильной из всех рассмотренных структур также является кристаллическая ячейка c симметрией P 21.
С целью апробации, предложенного подхода нами были проведены вычисления геометрической структуры кристаллических ячеек молекулы 5-(pметилбензоил)-1,3-диоксана (II), для различных пространственных групп Таблица 2. Результаты вычислений энергетических характеристик кристаллических структур (энергия взаимодействия атомов внутри элементарной ячейки кристалла Eintra, энергия взаимодействия атомов одной ячейки с атомами соседних ячеек Eintra, энтальпия сублимации Hs, энергия напряжения изолированной молекулы Eg, удельные значения, ккал/моль) для молекулы нафталина симметрии. Наилучшее согласие геометрических параметров (размеры ячейки, длины связей, валентные и двугранные углы) с имеющимися данными РСА наблюдается для структуры с пространственной группой симметрии P 21 /c.
Для проверки предложенного комбинированного подхода к вычислению единение (II), было нами изучено с по- мощью спектроскопии ЯМР 13 C высоO кого разрешения в поликристалле. В таблице 3 приведены экспериментальные значения констант ЯМЭ ядер 13 C ми функциями рассчитанными в рам- пользовались геометрические структуры молекулярных кристаллов, полученных методом, описанным ранее в гл. 2. Наблюдается достаточно удовлетворительное согласие между экспериментальными значениями констант ЯМЭ и рассчитанными теоретически для элементарной ячейки кристалла с симметрией P 21 /c, определённой ранее как предпочтительной.
Таблица 3. Результаты вычислений констант ЯМЭ ядер 13 C для нескольких рассчитанных кристаллических структур для соединения II. Приведены также экспериментальные значения, полученные методом ЯМР КП/ВМУ Конформационная структура и внутримолекулярная подвижность 2,9,10триметил-1,3-дитиа-5,6-бензоциклогептена (III) в растворе были хорошо изучены ранее с помощью метода динамического ЯМР. Из низкотемпературных (188-193 К) спектров ЯМР 1 H удалось выяснить, что в растворах соединение (III) существует в виде смеси двух конформеров — кресло-е и ванна-е:
Для установления пространственной структуры молекул данного соединения в кристалле соединение (III) было также изучено нами методом ЯМР 13 C КП/ВМУ, и на основе сравнения величин химических сдвигов ядер углерода было показано, что в поликристалле существует только конформация кресло-е. Для более детального анализа возможных вариантов Таблица 4. Результаты вычисления констант ядерного магнитного экранирования ядер 13 C для модельных структур кристаллической ячейки с различной симметрией для соединения III и экспериментальные значения, полученные методом ЯМР 13 C КП/ВМУ упаковки молекул данного соединения в поликристалле нами были промоделированы структуры молекулярных кристаллов обоих конформеров соединения (III) для нескольких наиболее распространённых пространственных групп симметрии. Для полученных пространственных структур также были рассчитаны константы ЯМЭ (табл. 4).
Сравнение рассчитанных значений констант ЯМЭ с экспериментально измеренными позволяет сделать вывод, что наиболее вероятной схемой упаковки молекул в кристалле для соединения (III) является упаковка конформера кресло-е в ячейке с пространственной симметрией P 21. Следует отметить, что этой структуре соответствует также максимальная плотность упаковки (минимальный объем, приходящийся на одну молекулу).
Спектры ЯМР 13 C 2-фенилэтинилдиметил-1,3-диоксана(IV), записанные ранее для раствора и по- 14 ли рассчитаны возможные структуры кристаллических ячеек, содержащих одну и, особенно, две кристаллографически неэквивалентные молекулы (т.е. не связанных ни одним из возможных для данной пространственной группы преобразований симметрии) соединения (IV). Размеры рассчитанных кристаллических ячеек приведены Рис. 1. Спектры ЯМР 13 C 2-фенилэтинил-5,5-диметил-1,3-диоксана в растворе (CD3(CD3 )2 CO и )2 CO (а) CКП/ВМУ поликристалла (б).
в табл. 5. Так же, как и для предыдущих соединений, для полученных в результате геометрических параметров молекул были рассчитаны значения констант ЯМЭ.
Полученные результаты позволяют сделать вывод о возможном способе упаковки молекул соединения (IV) в кристалле, как «голова к хвосту»
(рис 2), и позволяют объяснить особенности спектра ЯМР поликристалла данного соединения: элементарная ячейка молекулярного кристалла содержит две молекулы с различной геометрической структурой и различным магнитным окружением. Это отличие, которое наиболее велико для метильных групп в положении 5, приводит к разным значениям констант ЯМЭ для этих групп в каждой их этих молекул. Поэтому в спектре ЯМР Таблица 5. Геометрические параметры (линейные размеры —, A, углы —, объём ячейки — 3 ) элементарных ячеек молекулярных кристаллов, рассчитанных для молекулы соединения (IV) Ванна-е P 21 /c Кресло-е P 21 /c Рис. 2. Модельная кристаллическая структура соединения IV. Симметрия ячейки P 1. Две независимых молекулы в ячейке.
Энергия связывания E, ккал/моль Рис. 3. Зависимости энергии связывания, электронной плотности и констант ЯМЭ для кластера из молекул метана и бензола.
должны наблюдаться 2 группы из 2 линий (соответствующих аксиальному и экваториальному положению этих метильных групп) каждая. При этом значения теоретически полученных констант ЯМЭ для двух из этих четырёх линий весьма близки (33,8 м.д. и 33.7 м.д.) что при обычной для спектров КП/ВМУ ширине линии ЯМР должно приводить к наблюдению в спектре трёх линий, что полностью соответствует полученным ранее экспериментальным данным. Для анализа механизмов влияния окружения молекулы (IV) на её магнитнорезонансные параметры нами были рассчитаны зависимости констант ЯМЭ и электронных плотностей q для кластеров, состоящих из одной, трех и девяти кристаллических ячеек.
Для анализа роли слабых вандерваальсовых взаимодействий нами были проведены расчёты для некоторых систем (кластеры молекулы бензола с молекулами метана и ацетилена), которые моделируют фрагменты молекулярных кластеров, рассмотренных выше. Были получены и проанализированы зависимости энергии кластера, электронной плотности на ядрах C и 1 H и констант ЯМЭ от расстояния между молекулами. Некоторые из полученных зависимостей приведены на рис. 3. Анализ полученных зависимостей позволяет сделать вывод, что в данном случае (для вандерваальсовых кластеров) применяемое обычно разделение вкладов в константы ЯМЭ на диамагнитный, парамагнитный и т.п. не применимо, поскольку используемые в таком разложении параметры не трансферабельны, и их значения сами по себе зависят от магнитного окружения.
Спектры ЯМР КП/ВМУ третбутилкаликс[4]арена (V) и комплекса третбутилкаликс[4]арена с молекулой толуола приведены на рис. 5. Cигнал от толуо- CH ла оказывается сдвинутым относительно своего по- ложения в растворе в сторону сильных магнитных полей, причём сигнал от метильной группы (CH3 ) смещён сильнее, чем сигнал от фенильного кольца комплексов показывают, что атомы молекулы то- OH луола находятся от атомов молекулы каликсарена на расстояниях, характерных для вандервальсовых комплексов. Расчёты констант ЯМЭ для изолированной молекулы толуола и полученных структур комплекса «гость»–«хозяин»
позволяют сделать вывод, что наблюдаемое экспериментально смещение сигналов ЯМР молекулы тоулола можно объяснить влиянием фенильных колец молекулы каликсарена на ядра 13 C толуола-«гостя», помещённого в полость каликсарена-«хозяина».
каликсарена(V) и толуола. Атомы водорода условно не показаны.
Рис. 5. Спектры ЯМР 13 C КП/ВМУ трет-бутилкаликс[4]арена(V) (сверху) и его комплекса с молекулой тоулола (внизу).
Таблица 6. Экспериментально полученные константы ядерного магнитного экранирования 13 C для молекулы толуола в растворе, поликристаллическом комплексе с молекулой каликсарена(V) и рассчитанные теоретически значения для изолированной молекулы и комплекса Четвёртая глава посвящена учёту влияния растворителя на константы ЯМЭ. Учёт влияния растворителя проводился нами двумя методами: в рамках супермолекулярного подхода и в рамках поляризованной континуальной модели (PCM). Расчёты электронной и пространственной структуры молекулярных систем в рамках PCM проводились нами с использованием метода UB3LYP/6–31G(d,p) по программе GAMESS. Расчёты констант ЯМЭ с использованием полученной таким образом геометрии показали, что применение континуальной модели не позволяет достаточно успешно описать влияние растворителя на параметры ЯМР растворенной молекулы.
Результаты отражают только тенденцию изменения химических сдвигов в большинстве случаев.
Для расчётов в приближении «супермолекулы» сольватная оболочка для различного числа молекул растворителя конструировалась нами в соответствии с описанным в главе 2 алгоритмом. Из рассмотрения полученных структур (рис. 6), можно видеть, что расстояния от растворенной молекулы до ближайших атомов имеют порядок 2, 6 3 A, что указывает на хорошее согласие как с экспериментальными данными для димеров двухатомных молекул N2 и O2, так и с результатами неэмпирических квантовомеханических расчётов высокого уровня точности для димеров молекул средних размеров.
Пять молекул ацетона ориентированы вокруг молекулы метана по поверхности, близкой к сферической (рис. 6, a). При этом молекулы ацетона расположены друг от друга на расстояниях порядка 2, 012, 0 A, что позволяет говорить о наличии слабых водородных связей типа C H· · · O между молекулами ацетона.
Для кластеров, образованных молекулами бензола (рис. 6, б) можно видеть, что ближайшие расстояния от атомов молекулы CH4 до атомов Рис. 6. a) Молекула метана в окружении пяти молекул ацетона; б) молекула метана в окружении трёх молекул бензола.
турными данными о межмолекулярных расстояниях в кластерах со слабым нековалентным взаимодействием.
Анализ полученных результатов показывает, что наблюдаемое изменение химических сдвигов в таком супермолекулярном подходе обусловлено, в основном, влиянием сольватной оболочки, окружающей растворенную молекулу. С изменением числа молекул растворителя длина связи C H растворённой молекулы меняется незначительно (не более чем на 0.003 A для всех рассмотренных кластеров). Изменения валентных углов также недостаточно велики, чтобы объяснить изменение констант ЯМЭ только за счёт изменения геометрической структуры молекулы.
Проведённые нами расчёты геометрической структуры (табл. 7) кластеров молекулы соединения (III) с молекулами растворителя (сероуглерод), для обоих экспериментально наблюдаемых конформеров показали, что под влиянием сольватной оболочки происходит незначительная деформация молекулы в сторону «твист»-конформации, причём для конформации «ванна» эта деформация проявляется сильнее, чем для конформации «кресло», т.е. конформация «кресло» является более стабильной, что хорошо согласуется с накопленными экспериментальными данными (доля «кресло-е»
конформера в сероуглероде 0,85).
Для полученных пространственных структур были также рассчитаны константы ЯМЭ для ядер 13 C (табл. 8). Рассчитанные значения находятся в хорошем согласии с экспериментально наблюдаемыми, причём увеличение числа молекул растворителя позволяет получить рассчитанные значения, более близкие к экспериментальным.
Таблица 7. Молекула соединения (III) в сероуглероде. Геометрические параметры (валентные и торсионные углы, ).
Таблица 8. Молекула (III) в сероуглероде. Константы ядерного магнитного экранирования.
1. На основании комплексного использования модифицированного метода молекулярной механики и квантовохимических расчётов констант ядерного магнитного экранирования изучены спектры ЯМР ряда соединений в кристаллической фазе и в растворах.
2. Предложены следующие усовершенствования метода ММ: потенциалы для учёта вкладов дальнодействующих взаимодействий, эффектов диполь-дипольных взаимодействий, специфических взаимодействий типа водородных связей. Реализована возможность учёта взаимодействия молекул кристаллической ячейки с ближайшим окружением, а так же учёта пространственных групп симметрии при моделировании кристаллической ячейки. На основе анализа различных алгоритмов оптимизации предложены методы для улучшения сходимости результатов и ускорения вычислений. Предложен метод, позволяющий производить построение молекулярных кластеров конечного размера, позволяющий моделировать структуру сольватных оболочек и других сложных молекулярных систем.
3. Апробация предложенного метода, проведённая для кристаллических структур нафталина и 5-(p-метилбензоил)-1,3-диоксана, и сравнение с литературными данными РСА позволяет сделать вывод, что предложенный метод даёт возможность получать структуры, адекватные экспериментально наблюдаемым. Предложенный алгоритм поиска строения кластеров позволяет выбрать наиболее подходящие структуры как на основе сравнения энергии полученных кластеров, так и на основе принципа плотной упаковки (сравнения плотности полученных модельных кристаллов). Были получены спектры ЯМР 13 C КП/ВМУ для (p-метилбензоил)-1,3-диоксана и 2,9,10-триметил-1,3-дитиа-5,6бензоциклогептена, сравнение экспериментальных и рассчитанных значений констант ЯМЭ показало, что наилучшее их согласие наблюдается для структур с минимальной энергией и максимальной плотностью.
4. Для соединения 2-фенилэтинил-5,5-диметил-1,3-диоксана с помощью предложенного метода удалось объяснить сложную структуру спектров ЯМР КП/ВМУ. Наиболее вероятной (как с энергетической так и с точки зрения максимальной упаковки) оказалась структура с двумя кристаллографически неэквивалентными молекулами в кристаллической ячейке. Вычисления, проведённые с учётом нескольких ячеек улучшает согласие с экспериментом и свидетельствует о влиянии межмолекулярных взаимодействий на на магнитные свойства молекул в кристаллах.
5. Проведены расчёты пространственной структуры и констант ЯМЭ молекулы трет-бутил каликс[4]арена и комплекса этой молекулы с молекулой толуола. На основании сравнения этих данных с полученными спектрами 13 C ЯМР КП/ВМУ можно сделать вывод, что молекула толуола в комплексе ориентирована метильной группой внутрь молекулярной полости каликсарена.
6. Рассчитаны пространственные структуры молекулярных кластеров молекул метана и 2,9,1-триметил-1,3-дитиа-бензоциклогептена (для конформаций «кресло»-е и «ванна»-е) в окружении разного количества молекул растворителя. Для полученных структур проведены расчёты констант ядерного магнитного экранирования с учётом электронной корреляции квантовохимическими методами в рамках метода функционала плотности уровня B3LYP/6-31(d,p) и использования GIAO.
7. На основе полученных результатов моделирования растворителя на магнитнорезонансные параметры можно сделать вывод, что применение континуальной модели не даёт результатов, адекватных экспериментальным, тогда как применение предложенного подхода построения больших молекулярных кластеров при увеличении количества молекул растворителя приводит к лучшему согласию рассчитанных и экспериментальных значений.
Публикации по теме диссертации 1. Julmetov, A. R. Structure of molecules in crystalline lattice obtained by a modified method of molecular mechanics: calculation of 13 C chemical shifts / A. R. Julmetov, R. M. Aminova, A. V. Aganov // Applied magnetic resonance. — 2005. — Vol. 26. — Pp. 317–326.
2. Julmetov, A. R. Calculations of intermolecular interactions and its influences on structure and 13 C, 1 H magnetic shielding constants of solute molecules / A. R. Julmetov, R. M. Aminova, A. V. Aganov // Applied magnetic resonance. — 2005. — Vol. 29. — Pp. 277–285.
3. Клочков, В. В. Конформации 1,5-дигидро-3,7,8-триметил-2,4-бензодитиепина в различных растворителях и твердой фазе. Данные ЯМР. / В. В. Клочков, А. И. Рахматуллин, А. Р. Юльметов, А. В. Аганов // Ж. общ.
химии. — 1998. — Т. 68. — С. 837–841.
4. Каратаева, Ф. Х. Структура и внутримолекулярная подвижность nтио)фосфорил(тио)амидов. XIII. Структура N–фенил–N’(диизопропокситиофосфорил)тиомочевины. / Ф. Х. Каратаева, А. Р. Юльметов, А. В. Забиров и др. // Ж. общ. химии. — 2005. — Т. 75, № 6. — С. 908–911.
5. Юльметов, А. Р. Расчет геометрических параметров и констант ядерного магнитного экранирования 2-фенилэтинил-5,5-диметил-1,3-диоксана в кристаллической решетке / А. Р. Юльметов, Р. М. Аминова, А. В. Аганов // Новости ЯМР в письмах. — 2002. — Т. 1-2. — С. 1263–1268.
6. Юльметов, А. Р. Моделирование влияния среды на геометрические и магнитные свойства молекул в рамках континуальной модели с учетом специфической сольватации / А. Р. Юльметов, Р. М. Аминова, А. В. Аганов // Структура и динамика молекулярных систем. Cборник статей. Выпуск XI. — Казань, 2004. — Т. 2. — С. 94–96.
7. Юльметов, А. Р. Расчет геометрических параметров и констант ядерного магнитного экранирования 2-фенилэтинил-5,5-диметил-1,3-диоксана в кристаллической решетке / А. Р. Юльметов, Р. М. Аминова, А. В. Аганов // Новые аспекты применения магнитного резонанса. Сб. трудов молодёжной научной школы «Актуальные проблемы магнитного резонанса и его приложений». — Казань, 2002. — С. 170–174.
8. Брусько, В. В. Комплексы тиофосфорилированных тиомочевин. Данные ЯМР / В. В. Брусько, А. И. Рахматуллин, А. Р. Юльметов и др. // Новые достижения ЯМР в структурных исследованиях. Тезисы докладов. — Казань, 2000. — С. 105.
9. Каратаева, Ф. Х. Экспериментальное и теоретическое изучение строения N-фенил,N’(диизопропокситиофосфорил-тио)мочевины / Ф. Х. Каратаева,, А. Р. Юльметов и др. // III Всероссийская конференция «Новые достижения ЯМР в структурных исследованиях». Тезисы докладов. — Казань, 2000. — С. 108.
10. Юльметов, А. Р. Моделирование структуры молекулярных кластеров методом молекулярной механики / А. Р. Юльметов, А. И. Рахматуллин, А. В. Аганов // Международная конференция по новым технологиям и приложениям современных физико-химических методов (ядерный магнитный резонанс, хроматография/масс-спектрометрия, ИКФурье спектроскопия и их комбинации) для изучения окружающей среды, включая секции молодых ученых Научно-образовательных центров России. Материалы конференции. Ростов-на-Дону, 25-28 сентября 2001 г. — Ростов-на-Дону, 2001. — С. 134–135.
11. Юльметов, А. Р. Программа для вычисления структуры молекулярных кристаллов в рамках метода молекулярной механики / А. Р. Юльметов, А. И. Рахматуллин, А. В. Аганов // VIII Всероссийская конференция «Структура и динамика молекулярных систем». — Казань– Москва–Йошкар-Ола–Уфа, 2001. — С. 82.
12. Karataeva, F. K. NMR experimental and theoretical studing of the intramolecular mobility of N–phenil–N’(diizopropoxythiophosphorilthio)urea / F. K. Karataeva, A. R. Julmetov, A. V. Aganov // Международный семинар «Современное развитие магнитно-резонансной томогорафиии спектроскопии. Физические основы и применение в медицине и биологии», 11-13 июня, Казань». — Казань, 2001. — С. 63.
13. Юльметов, А. Р. Спектры КП/ВМУ и расчеты геометрических параметров и констант ядерного магнитного экранирования для ряда молекулярных кристаллов / А. Р. Юльметов, Р. М. Аминова, А. В. Аганов // X Всероссийская конференция «Структура и динамика молекулярных систем», cборник тезисов докладов и сообщений, 30 июня – 4 июля 2003 г. — Казань–Москва–Йошкар-Ола–Уфа, 2003. — С. 82.
14. Юльметов, А. Р. Моделирование влияния среды на геометрические и магнитные свойства молекул в рамках континуальной модели с учетом специфической сольватации / А. Р. Юльметов, Р. М. Аминова, А. В. Аганов // X Всероссийская конференция «Структура и динамика молекулярных систем», сборник тезисов докладов и сообщений, июня – 4 июля 2003 г. — Казань–Москва–Йошкар-Ола–Уфа, 2003. — 15. Юльметов, А. Р. Межмолекулярные взаимодействия и их влияние на структуру и константы ядерного магнитного экранирования 13 C и H / А. Р. Юльметов, Р. М. Аминова, А. В. Аганов // Юбилейная научная конференция физического факультета. Программа. Сборник тезисов. Казань, 10 ноября 2004 г. — Казань, 2004. — С. 88.
16. Юльметов, А. Р. Моделирование влияния среды на структуру и константы ядерного магнитного экранирования / А. Р. Юльметов, Р. М. Аминова, А. В. Аганов // IV Всероссийская конференция «Новые достижения ЯМР в структурных исследованиях».Сборник тезисов. Казань, 4–7 апреля 2005 года. — Казань, 2005. — С. 103.
17. Аганов, А. В. Теоретическое изучение влияния слабых невалентных взаимодействий на структуру и параметры ЯМР ассоциатов / А. В. Аганов, Р. М. Аминова, А. А. Нафикова, А. Р. Юльметов // III Международная конференция по новым технологиям и приложениям современных физико-химических методов (ядерный магнитный резонанс, хроматография/масс-спектрометрия, ИК-Фурье спектроскопия и их комбинации) для изучения окружающей среды, включая секции молодых ученых Научно-образовательных центров России. Материалы конференции. Ростов-на-Дону, 21-25 марта 2005 г. — Ростов-наДону, 2005. — С. 57–58.
18. Julmetov, A. R. Van der Waals and H-bond interactions and their influences on the structure and chemical shifts of a molecule in solution / A. R. Julmetov, R. M. Aminova, A. V. Aganov // International Symposium and Summer School in “Nuclear Magnetic Resonance in Condensed Matter”, 20nd Meeting “NMR in Life Sciences”, 11-15 July 2005, SaintPetersburg. — Saint-Petersburg, 2005. — P. 79.