Н. Д. Зелинского к. Я. Бурштейн П. П. Шорыгин квантовохимические расчеты в органической химии и молекулярной спектроскопии

Таблица 1.9 Вычисленные энергии активации (кДж/моль) для реакций, приведенных на схемах I - V [56] Реакция МПДП КМПДП ХФ КХФ* 1 400 331 344 327 2 523 420 483 419 3 627 589 638 559 4 343 316 315 286 5 396 271 212 202 6 458 390 437 402 7 405 299 455 361 8 84 96 119 131 9 538 502 550 — 10 568 540 483 430 11 458 382 411 436 12 425 347 426 403 13 453 323 348 378 14 341 349 338 353 15 318 268 278 262 16 78 38 26 23 17 360 201 218 161 18 444 352 397 326 19 609 467 538 428 20 370 238 254 200 21 399 372 452 390 22 286 280 250 186 23 403 310 420 — 24 437 275 319 259 *При расчете энергии активации использована геометрия переходного состояния, предварительно вычисленная в приближении Хартри—Фока. Таким образом, мы видим, что для определения константы ско­рости реакции нет необходимости рассчитывать всю ППЭ. Достаточно найти на ней лишь стационарные точки. Геометрию исходных реа­гентов для любой реакции можно найти с помощью минимизации полной энергии. Эту задачу сравнительно легко решить. Гораздо сложнее найти геометрию переходного состояния. Переходные состоя­ния являются седловыми точками на ППЭ. Они похожи на минимумы, поскольку в них, как и в минимумах, градиент полной энергии по координатам равен нулю, но их нельзя найти с помощью минимизации полной энергии. В работе [55] было предложено искать переходные состояния с помощью минимизации квадрата модуля градиента, т.е. минимизировать сумму квадратов производных полной энергии по всем независимым координатам: Минимизация этого выражения действительно позволяет найти переходные состояния, однако сделать это достаточно сложно. Выражение для σ имеет очень много минимумов, так как не только стационарные точки дают минимумы для σ, но и любой изгиб ППЭ типа плеча приводит к появлению минимума для этой величины. Поэтому предложенная методика [55] пригодна лишь в непосред­ственной близости от седловой точки, и основная трудность при расчете геометрии переходных состояний этим методом заключается в том, что исходная (стартовая) геометрия, с которой начинается поиск, должна быть достаточно близка к конечной (истинной) гео­метрии переходного состояния. Определять стартовую геометрию приходится на основе интуитивных соображений, так как из экспе­римента получить хотя бы приближенную информацию о геометрии переходных состояний невозможно. Поэтому успешный поиск переходных состояний с помощью минимизации выражения для σ в значительной степени зависит от опыта исследователя. Наибольший интерес вызывает вопрос о точности, с которой можно рассчитать геометрию переходного состояния и энергию активации методами квантовой химии. К сожалению, мы не располагаем доста­точно полным набором данных эксперимента по энергиям активации для различных типов реакций с участием маленьких молекул, а экспе­риментальные значения геометрических параметров переходных состоя­ний получить невозможно. Поэтому в настоящем разделе нам пришлось ограничить рассмотрение вопроса о точности расчета переходных состояний методами квантовой химии сопоставлением результатов, которые были получены полуэмпирическими и неэмпирическими ме­тодами (см. [56, 57] и ссылки в них). В табл. 1.9 и 1.10 сопоставлены данные расчета параметров переходных состояний для реакций, изображенных на схемах I—V (см. стр. 40), методами МПДП и КМПДП (метод МПДП с учетом электрон­ной корреляции), неэмпирическим методом в приближении Хартри-Фока без учета и с учетом электронной корреляции. В таблицах неэмпирический расчет в приближении Хартри-Фока без учета электронной корреляции обозначен ХФ, с учетом электрон­ной корреляции - КХФ [56, 57]. Таблица 1.10 Статистический анализ результатов расчета геометрии переходных состояний и энергий активации методами МПДП и КМПДП [56, 57] Вычисленная величина Число расчетов Абсолютная ошибка*1 МПДП КМПДП Длина валентных связей*2, нм из них активных*3 пассивных*3 Валентный угол*2, град Торсионный угол*2, град Энергия активации, кДж/моль 112 70 42 58 20 24 0,0057 0,0078 0,0018 7,9 11,6 55 90*4 0,0056 0,0073 0,0025 6,2 7,9 45 36*4 *1 По сравнению с данными неэмпирического расчета в приближении Хартри-Фока. *2 Для переходного состояния. *3 Активные связи разрываются или образуются в ходе реакции, пассивные связи в ходе реакции формально остаются неизменными. *4 По сравнению с данными неэмпирических расчетов с учетом электронной корреляции. Для реакций, изображенных на схемах I - III, расчеты с оптимизацией геометрии в приближении Хартри-Фока проведены в базисе 6-31ГФ*, для реакций, изобра­женных на схемах IV и V, - в базисе 3-21ГФ или 4-31ГФ. Электронная корреляция учитывалась только при вычислении энергии активации. Эти расчеты проведены по теории возмущений Мёллера-Плезетта с точностью до четвертого порядка включительно или в приближении связанных электронных пар. Схема I Схема II Схемa III Схема IV Схемa V Из этих данных видно, что геометрические параметры пере­ходных состояний, вычисленные методами МПДП и КМПДП, находят­ся в хорошем согласии с данными неэмпирических расчетов без учета электронной корреляции. Однако для длин валентных связей, которые формально разрываются или образуются в ходе реакции (активные связи), расхождение результатов, полученных полуэмпирическими и неэмпирическими методами, заметно больше, чем для валентных связей, которые в ходе реакции остаются формально неизменными (пассивные связи). Сопоставить результаты расчета геометрии пере­ходных состояний методами МПДП и КМПДП с данными неэмпири­ческих расчетов с учетом электронной корреляции невозможно из-за почти полного отсутствия в литературе последних. При расчете энергий активации величины, полученные методами МПДП и КМПДП, отличаются на 15 - 20% от значений, вычисленных неэмпирическим методом без учета электронной корреляции. При сравнении с данными неэмпирических расчетов с учетом электронной корреляции оказывается, что метод КМПДП дает более близкие результаты. Это, по-видимому, связано с выбором реакций для тестирования. В отобранных реакциях учет электронной корреляции может заметно изменить результаты расчета, так как для сравнения авторам работ [56, 57] пришлось взять лишь те химические прев­ращения, для которых из литературы были известны результаты вычислений неэмпирическими методами с использованием достаточно большого базиса и с учетом электронной корреляции. Такие расчеты, естественно, были опубликованы лишь для тех реакций, в которых электронная корреляция играет существенную роль. Для других типов реакций литературные данные отсутствуют, и их получение связано с очень большим объемом вычислений. Таблица 1.11 Вычисленные тепловой эффект Q и энергия активации Н# реакции (25) (кДж/моль) Метод Q H# H#обр ХФ/3-21ГФ ХФ/6-31ГФ* МП2/6-31ГФ* МПЗ/6-31ГФ* МПЗ/6-311ГФ** МП4/6-311ГФ** С учетом нулевых колебаний 162 144 191 178 181 193 185 270 213 203 207 201 202 189 108 69 12 31 20 9 4 Примечание. Использованы следующие обозначения: ХФ — расчет в приближении Хартри-Фока; МП2, МПЗ и МП4 — расчеты с учетом электронной корреляции во втором, во втором и третьем и во втором, третьем и четвертом порядках теории возмущений Меллера-Плезетта. В методах МПЗ/6-31ГФ**, МПЗ/6-311ГФ** и МП4/6-311ГФ** использована геометрия, вычисленная предварительно методом МП2/6-31 ГФ*; Нобр - энергия активации обратной реакции. Некоторое представление о влиянии базиса и электронной кор­реляции на результаты неэмпирических расчетов высоты активационных барьеров можно получить из табл. 1.11 и 1.12, в которых приведены данные для реакций HC≡CH → H2C=C: (25), Для реакции (25) зависимость достаточно сильная [58-60], для (26) - слабая [61]. Существенное влияние электронной корреляции на результаты расчета высоты активационных барьеров можно ожидать для реакций, в ходе которых образуются бирадикалы, карбены и родственные им структуры. Кроме того, электронную корреляцию обязательно надо учитывать при расчете ППЭ для реакций гомолитического разрыва валентных связей. Таблица 1.12 Энергии активации реакций (26) (кДж/моль) [61] Метод (26а) (26б) (26в) Метод (26а) (26б) (26в) ХФ/3-21ГФ МП2/3-21ГФ ХФ/6-31ГФ МП2/6-31ГФ ХФ/6-31ГФ* 176 126 218 156 271 161 150 197 174 231 181 134 192 145 217 МП2/6-31ГФ* ХФ/6-31ГФ** МП2/6-31ГФ** МПЗ/6-31ГФ** МП4/6-31ГФ** 180 265 176 202 197 167 224 163 170 174 165 213 162 185 176 Примечание. Обозначения те же, что и в предыдущей таблице. Во всех расчетах использована геометрия молекул и переходных состояний, вычисленная предварительно методом ХФ/3-21ГФ. Вопрос о влиянии электронной корреляции на геометрию пере­ходных состояний был рассмотрен в работе [62]. В ней мето­дом МПДП без учета и с учетом электронной корреляции была рассчитана геометрия переходных состояний для реакций, изоб­раженных на схемах I - V, и показано, что учет электронной корреляции изменяет межатомные расстояния в среднем не более чем на 0,003 нм, при этом изменения длины активных валентных связей, т.е. тех, которые разрываются или образуются в ходе реакции, не превышают 0,004 нм, а пассивных, т.е. тех, кото­рые в ходе реакции остаются формально теми же, 0,001 нм, валентные углы меняются менее чем на 3°. Эти небольшие изме­нения геометрических параметров незначительно меняют структуру переходного состояния и практически не сказываются на резуль­татах расчета энергии активации. На основе этих данных в рабо­те [62] был сделан вывод о том, что корреляционные эффекты можно не учитывать при квантовохимических расчетах геометрии переходных состояний. Аналогичных исследований, выполненных не­эмпирическими методами, в настоящее время нет. 1.7. РАСЧЕТЫ ДЛЯ СИСТЕМ С ВОДОРОДНЫМИ СВЯЗЯМИ Системы с водородными связями играют важную роль в химии и биологии, поэтому для их изучения широко применяются как экспериментальные, так и теоретические методы. Квантовохимические расчеты для систем с водородными связями проводят полу­эмпирическими (ППДП/2, ППДП/БУ, МЧПДП/3Н и МПДП/Н) и неэмпирическими методами. Наибольшее число работ, опубликован­ных за последние 20 лет, выполнено методами ППДП/2 и ОСТ-3 ГФ [63]. Однако эти методы имеют ряд существенных недостатков. Для большинства систем они в 1,5 - 2 раза переоценивают энергию образования водородных связей и на 0,02—0,03 нм занижают их длину. Кроме того, при решении прикладных задач, наряду с вы­числением энергии и длины водородной связи обычно приходится рассчитывать термохимические параметры соединений и комплек­сов с водородными связями. Методы ППДП/2 и ОСТ-3ГФ не позволя­ют этого сделать, в результате область их применения при изу­чении реакционной способности систем с водородными связями ограничена. Ими пользовались в течение длительного времени только потому, что более точных и столь же экономичных мето­дов не было. Таблица 1.13 Вычисленные энергии водородных связей (кДж/моль) [64] Комплекс ППДП/БУ ППДП/2 ОСТ-ЗГФ 4-31ГФ Экспсримен Н2О...ОН 184 302 294 172 142 Н2О...О2Н3 87 155 155 126 97 Н3О+...ОН2 113 252 250 185 136 - 155 В работе [64] была сделана попытка использовать для изуче­ния реакционной способности систем с водородными связями ме­тод ППДП/БУ. Оказалось, что по сравнению с методами ППДП/2 и ОСТ-ЗГФ он дает заметно более точные результаты при рас­чете энергий водородных связей. Это особенно хорошо видно на примере прочности водородных связей в заряженных системах (табл. 1.13), точность получается примерно такая же, как в не­эмпирических расчетах с использованием базиса 4-31ГФ [65]. Однако метод ППДП/БУ существенно занижает длину водородных связей (немного сильнее, чем методы ППДП/2 и ОСТ-ЗГФ) и дает боль­шие ошибки при расчете термохимических параметров органи­ческих соединений. Поэтому метод ППДП/БУ, так же как и методы ППДП/2 и ОСТ-ЗГФ, вряд ли целесообразно использовать при количественном решении прикладных задач. Для изучения реакционной способности органических и био­органических систем с водородными связями необходим метод, который позволял бы одновременно вычислять как термохимичес­кие параметры молекул и переходных состояний, так и парамет­ры систем с водородными связями. В настоящее время этим тре­бованиям удовлетворяют методы МЧПДП/3Н и МПДП/Н [66 -71], которые являются модифицированными вариантами широко извест­ных методов МЧПДП/3 и МПДП. Основной недостаток методов МЧПДП/3 и МПДП заключается в том, что они не могут быть исполь­зованы для изучения систем с водородными связями, хотя послед­ние очень часто встречаются при решении прикладных задач, осо­бенно когда мы имеем дело с биоорганическими соединениями. В методах МЧПДП/3Н и МПДП/Н этот недостаток устранен за счет не­большой модификации некоторых эмпирических формул, которые используются для оценки матричных элементов. Кроме этих двух методов, для расчета параметров систем с водородными связями можно использовать метод AM1, но по точности он заметно уступает методу МПДП/Н. Неэмпирические расчеты для систем с водородными связями в настоящее время чаще всего проводят с использованием валентно-расщепленных базисов (обычно 3-21ГФ или 4-31ГФ). Их выбор связан с необходимостью идти на компромисс между точ­ностью расчета и затратами машинного времени. Расчеты в этих базисах немного завышают энергию образования водородных свя­зей, занижают их длину (расстояние между гетероатомами), а так­же существенно переоценивают стабильность циклических структур с водородными связями. Так, при использовании базисов 3-21ГФ и 4-31ГФ циклический димер воды получается стабильнее линейного [65]. Циклическая структура комплекса аниона с водой стабильнее линейной [72] и т.д., хотя из эксперимента известно, что линейные структуры являются более стабильными. Включение в базис поляризационных орбиталей устраняет эти недостатки [73]. Таким образом, мы видим, что при использовании неэмпирических методов расчеты для систем с водородными связями желательно проводить в базисах с поляризационными орбиталями, т.е. исполь­зовать базис 6-31ГФ* или 6-31ГФ**, однако сделать это практически очень трудно, так как системы с водородными связями, которые представляют интерес для органической химии и биохимии, обычно имеют большие размеры и их свойства не удается моделировать с помощью простых систем, состоящих из небольших молекул. В табл. 1.14 - 1.16 сопоставлены результаты квантовохимических расчетов параметров систем с водородными связями методами ППДП/2, МПДП/Н и 4-31ГФ. Из этих данных видно, что мето­ды МПДП/Н и 4-31ГФ дают приблизительно одинаковые резуль­таты. Следует, однако, отметить, что метод МПДП/Н в отличие от 4-31ГФ правильно передает относительную стабильность ли­нейных и циклических систем с водородными связями. Метод ППДП/2 существенно уступает по точности методам МПДП/Н и 4-31ГФ. Таблица 1.14 Энергии водородных связей (кДж/моль) [70] У—Н...Х ППДП/2 МПДП/ВС 4-31ГФ Экспери­мент У—Н X Н2О Н2О* 37 18 34 21 Н2О ОН - 302 134 172 142 Н2О ОН- • Н2О 155 88 126 97 Н2O ОН- • 2Н2О 105 63 97 76 Н3О+ Н2О 252 130 185 137 - 155 Н3О+ • Н2О Н2О 147 97 130 92 - 97 Н3О+ • 2Н2О Н2О 105 80 109 71 NН3 Н2О 19 8 15 — Н2О NH3 46 12 33 — NH3 NH3 19 4 14 19 NH4+ Н2О — 82 118 71 NH4+• Н2О Н2О — 71 — 63 NН4+ • 2Н2О Н2О — 63 — 55 NH4+ • 3Н2О Н2O — 47 — 50 NH4+ NH3 — 155 130 105 NH4+ • NH3 NH3 — 97 — 76 NH4+ • 2NH3 NH3 — 71 — 59 NH4+ • 3NH3 NH3 — 59 — 50 FH Н2О 63 38 59 30 Н2О FH 13 21 21 — FH NH3 84 50 63 — NH3 FH 6 9 13 — FH FH* 39 37 33 21 - 29 FH F - 428 197 — 168 - 252 Н2О F -* 334 130 — 97 Н2О F - • Н2О* — 97 — 71 Н2О F - • 2Н2О* — 76 — 59 Н2О F -• ЗН2О* — 50 — 55 *Данные для наиболее стабильной линейной структуры. При решении прикладных задач в случае систем с водород­ными связями часто возникает необходимость вычислять сечения ППЭ для реакций переноса протона между гетероатомами, так как эти процессы сопровождают многие органические и биоорганиче­ские реакции. Расчет сечения ППЭ переноса протона провести относительно легко, но при этом необходимо учитывать сильную зависимость результата от длины водородной связи. Так, в си­стемах H2O...OH и Н2O...Н3О+ при расстояниях О...О меньше 0,25 нм на сечении ППЭ получается только один минимум, который соот­ветствует положению протона посредине между атомами кисло­рода. При расстоянии О...О больше 0,25 нм на сечении ППЭ будут уже два минимума, соответствующие положению протона либо у одного атома кислорода, либо у другого. Потенциальный барь­ер между этими минимумами очень быстро возрастает с увеличением расстояния между атомами кислорода [64, 74, 75]. В ра­боте [64] показано, что в ходе переноса протона расстояние O...O может немного уменьшиться и при этом активационный барьер переноса протона будет снижаться. Рассмотрим в качестве примера результаты расчетов методами МПДП/Н и 4-31ГФ сечения ППЭ для синхронного переноса двух протонов в циклическом димере му­равьиной кислоты. Таблица 1.15 Вычисленные расстояния (нм) между гетероатомамн (X—Y) (длины водородных связей) [70] У—Н...Х ППДП/2 МПДП/Н 4-31ГФ У—Н X Н2O Н2O 0,254 0,262 0,283 Н2O OH - 0,234 0,242 0,245 Н3О+ H2O 0,234 0,244 0,236 NН3 H2O 0,269 0,290 0,319 Н2O NН3 0,258 0,300 0,298 NH3 NН3 0,274 0,340 0,336 NH4+ Н2O — 0,273 0,265 NH4+ NH3 — 0,254 0,275 HF Н2О 0,242 0,254 0,260 Н2O HF 0,255 0,268 0,292 HF NH3 0,246 0,253 0,272 NH3 HF 0,275 0,286 0,319 HF PA 0,244 0,256 0,270 Таблица 1.16 Перенос электронной плотности с X на Y—Н при образовании водородных связей У—Н...Х (единицы заряда электрона) [70] У—Н...Х ППДП/2 МПДП/Н 4-31 ГФ У—Н Х Н2O Н2O 0,036 0,007 0,028 Н2О NН3 0,059 0,007 0,037 H2O HF — 0,014 0,028 HF Н2O 0,064 0,036 0,036 HF NН3 0,097 0,064 0,057 H2O OH - 0,365 0,207 0,210 H2O ОН - • Н2O — 0,132 0,130 Н2O ОН - • 2Н2O — 0,109 0,100 Н3О+ Н2O 0,289 0,139 0,140 Н3О+ • Н2O Н2O — 0,086 0,080 Н3О+• 2Н2O Н2O — 0,073 0,070 Из эксперимента известно, что высота активационного барьера этой реакции мала, однако в ряде опубликованных работ квантовохимические расчеты привели к противоположному результату [76 - 83] Это связано с тем, что не проводилась полная опти­мизация геометрии и не учитывалось изменение длины водородных связей в переходном состоянии. Высота активационного барьера рассматриваемой реакции, вычисленная с полной оптимизацией гео­метрии методами МПДП/Н и 4-31ГФ, составила менее 40 кДж/моль [69, 72]. Полученные геометрические параметры приведены в табл. 1.17. Из этих данных видно, что в ходе переноса протонов в переходном состоянии происходит сжатие димера. Длина водородной связи (рас­стояние между атомами кислорода) уменьшается на 0,03 - 0,04 нм. За счет этого снижается высота активационного барьера. Если зафик­сировать расстояние между атомами кислорода в водородных связях, то высота активационного барьера увеличится до 280 кДж/моль. Таблица 1.17 Межатомные расстояния (нм) димера муравьиной кислоты и переходного состояния реакции синхронного переноса двух протонов [69] Связь Димер Переходное состояние МПДП/Н 4-31ГФ эксперимент МПДП/Н 4-31ГФ С=O 0,122 0,122 0,122 0,128 0,126 С—О 0,133 0,131 0,132 0,128 0,126 O...O 0,261 0,270 0,270 0,237 0,239 С—Н 0,110 0,107 0,108 0,110 0,107 Из приведенного в настоящей главе материала читатель смог убедиться, что квантовохимические расчеты позволяют получать много полезной информации о реакционной способности органи­ческих соединений, которая представляет интерес для специали­стов в области тонкого органического синтеза и физической орга­нической химии. Это прежде всего индексы реакционной способ­ности. В некоторых случаях они дают возможность определить направление реакции, а иногда и выход конечных продуктов. Бо­лее детальную информацию о механизме реакции удается получить на основе расчетов тепловых эффектов и энергий активации. Расчеты индексов реакционной способности обычно проводят по­луэмпирическими методами (чаще всего пользуются методом ППДП/2), при этом либо используют геометрию, полученную эксперимен­тально, либо полагают, что все межатомные расстояния и валентные углы равны своим стандартным значениям. При вычислении термодинамических параметров необходимо находить оптимальную геометрию, соответствующую минимуму полной энергии. Расчеты с оптимизацией геометрии целесообразно проводить методами МПДП или 3-21 ГФ. Эти методы позволяют вычислять геометрию как стабиль­ных соединений, так и короткоживущих интермедиатов и переход­ных состоянии. Для уточнения значений термодинамических пара­метров иногда достаточно провести расчет с фиксированной гео­метрией в базисе 6-31ГФ. Для получения более надежных результа­тов необходимо включить в базис поляризационные орбитали на всех атомах и учесть энергию корреляции. Однако практически сделать это можно лишь для систем, содержащих небольшое число атомов. Для соединений, которые представляют интерес для органической и биоорганической химии, приходится ограничиться рас­четами методом МПДП. В этом случае необходимо предварительно на небольших модельных системах, для которых доступны данные хороших неэмпирических расчетов, убедиться, что метод МПДП дает для данного класса реакций правильные результаты.

Похожие:

	Тематическое планирование по химии 11 класс Предмет органической химии. Краткий исторический очерк: виталистическая теория. Связь органической и неорганической химии. Органические...		Программа по формированию навыков безопасного поведения на дорогах... Повторение курса органической химии. Строение и номенклатура органических веществ. Вывод молекулярной формулы вещества
	Программа по формированию навыков безопасного поведения на дорогах... Знать : определение органической химии, что изучает данная наука; различия между органическими и неорганическими веществами, особенности...		Тема Форма проведения Свойство атомов углерода образовывать прямые, разветвленные и замкнутые цепи, одинарные и краткие связи. Гомология, изомерия. Значение...
	Ключевые слова Ключевые слова: структура молекул, водородная связь, ик спектроскопия, спектроскопия ямр, квантовохимические расчеты, молекулярные...		Учебно-методический комплекс учебной дисциплины методика изучения... Программа учебной дисциплины обсуждена и утверждена на заседании кафедры органической химии Протокол №9 от 23. 05. 2012 г
	Основные результаты, полученные в исследованиях по оптике и спектроскопии,... Методические указания разработаны кандидатом физико-математических наук, доцентом Нойкиным Ю. М		Рабочая программа по дисциплине ен. Ф. 04. 2 Органичская химия Целью дисциплины является освоение системы знаний о фундаментальных законах, теориях, фактах органической химии необходимых для понимания...
	Характеристика молекулярной патологии в увеальных меланомах Работа выполнена в лаборатории молекулярной генетики человека нии молекулярной медицины Московской медицинской академии им. И. М....		Программа дисциплины дпп. Ф. 10 Органический синтез цели и задачи... Курс «Органический синтез» проводится после изучения систематического курса органической химии и выполнения практических работ малого...
	Магнитной радиоспектроскопии в биофизических и медико-биологических исследованиях Диамагнетизм и парамагнетизм атомов и молекул. Ядерный магнетизм. Основы эпр-спектроскопии и ямр-спектроскопии		Программа по формированию навыков безопасного поведения на дорогах... Элективный курс предназначен для учащихся 10-х классов, изучающих химию на базовом уровне. Курс рассчитан на 34 часа. Введение данного...
	Программа по формированию навыков безопасного поведения на дорогах... Элективный курс предназначен для учащихся 10-х классов, изучающих химию на базовом уровне. Курс рассчитан на 34 часа. Введение данного...		Рабочая программа дисциплины Теоретические основы органической химии... Целями освоения дисциплины Теоретические основы органической химии биологически активных добавок являются
	Реакции гетероциклизации с участием цианаминопиранов ...		Урок по химии по разделу «Химические реакции в органической химии» Раннев П. П. Философия религии : Учебно-метод комплекс. Ростов н/Д.: Изд-во рсэи, 2011. 63 с