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

25

на правах рукописи

БАБИН Юрий Владимирович

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

02.00.08 - химия элементоорганических соединений

02.00.04 - физическая химия

АВТОРЕФЕРАТ

диссертации на соискание учёной степени доктора химических наук

Владивосток -- 2007

Работа выполнена в ГОУ ВПО Тихоокеанский государственный экономический университет

Научный консультант	доктор химических наук, профессор Устынюк Юрий Александрович
Официальные оппоненты	Член-корреспондент РАН, профессор Нифантьев Эдуард Евгеньевич доктор химических наук, профессор Жидомиров Георгий Михайлович доктор химических наук, профессор Вовна Виталий Иванович
Ведущая организация	Российский химико-технологический университет им. Д.И. Менделеева

Защита состоится 05 ноября 2007 года в 1400 часов на заседании диссертационного совета Д 212.056.05 при Дальневосточном государственном университете по адресу: 690600, г.Владивосток, ул. Октябрьская, 27, ДВГУ, Химический факультет, тел./факс 8(4232) 457609, e-mail svistun@chem.dvgu.ru

С диссертацией можно ознакомиться в библиотеке Дальневосточного государственного университета.

Автореферат разослан « » ____________ 2007 года.

Учёный секретарь диссертационного совета, кандидат химических наук И.В. Свистунова

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность исследования.

Миграции атомов водорода - важнейший тип химических процессов. Подавляющее число кислотно-основных равновесий в растворах и твердых телах, многочисленные гомогенные и гетерогенные каталитические реакции и многие трансформации органических молекул включают перенос атомов или ионов водорода в качестве ключевых элементарных стадий. Протонные переносы в цепочках водородных связей, в том числе эстафетные множественные переносы, лежат в основе функционирования всех живых систем, определяют механизмы действия большинства ферментов, управляют процессами распознавания, связывания и транспорта биосубстратов в клетках. Хорошо известны и детально изучены как с экспериментальной, так и с теоретической точек зрения разнообразные типы прототропных таутомерных превращений, в которых атом водорода перемещается между атомами элементов-органогенов второго периода. Таковы, например, кетоенольная таутомерия карбонильных соединений (обратимая миграция СО), лактим-лактамная таутомерия (обратимая миграция ОN), прототропная таутомерия углеводородов циклопентадиенового и инденового рядов (обратимая миграция СС), движения протона в двуямном потенциале в енольных формах 1,3-дикарбонильных соединений (обратимая миграция ОО), изомеризации производных пиразола и имидазола (обратимые миграции NN) и многие другие процессы такого типа. Среди них следует особо отметить синхронные перемещения двух водородных атомов в потенциалах более сложной формы, такие как согласованные миграции протонных пар в димерах карбоновых кислот и в четырехямном потенциале в порфиринах. Атомы водорода в таких системах ведут себя как квантовые частицы, и поэтому такие миграции часто совершаются по туннельному механизму.

В отличие от переносов атомов водорода между атомами второго периода их миграции в гетероатомных системах, включающих атомы третьего и более высоких периодов исследованы значительно хуже. Среди таких миграций особенно мало изучены те из них, где атом водорода перемещается между тяжелыми атомами, существенно различающимися между собой по атомной массе и/или электроотрицательности. Например, в химии фосфорорганических соединений экспериментально обнаружено несколько интересных перегруппировок, протекающих с миграцией водорода между фосфором и кислородом, а также между фосфором и азотом. Однако до недавнего времени механизм большинства из них методами квантовой химии высокого уровня не изучался.

Создание и быстрое внедрение в практику методов металлокомплексного катализа несомненно стало одним из самых ярких достижений химической науки последних пяти десятилетий, оказавших революционизирующее влияние на лабораторный и промышленный органический синтез. Значительное число каталитических реакций, протекающих под влиянием комплексов поздних переходных металлов, в том числе реализованных в виде многотоннажных производств (гидрирование, гидросилилирование, гидрометаллирование алкенов и алкинов) также протекают с миграциями атомов водорода во внутренней координационной сфере тяжелого атома металла на ключевых стадиях, в том числе на стадиях окислительного присоединения и восстановительного элиминирования. Теоретическому изучению этих важнейших реакций современными квантовохимическими методами посвящены сотни работ, однако до сих пор остается дискуссионным вопрос о роли туннелирования в таких процессах. Между тем решение этой проблемы является важным шагом на пути создания новых типов каталитических процессов, в том числе для мягкой функционализации алканов, и вследствие этого представляет значительный интерес не только с фундаментальной, но и с прикладной точки зрения.

В последнее время химия гидрофосфорильных соединений (ГФС) получила новый мощный импульс в своем развитии благодаря тому, что вторичные фосфиноксиды и их таутомеры, фосфинистые кислоты, начали широко использоваться в качестве донорных лигандов в комплексах с поздними переходными металлами (Pd, Pt, Ni, Rh), которые выступают как эффективные катализаторы различных реакций кросс-сочетания, а также гидроформилирования и гидрофосфорилирования. В диссертации подробно исследованы механизмы гидроформилирования и гидрирования алкенов на комплексе платины с гидрофосфорильными лигандами 3.

3 (R,R' = H)

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

Представленная диссертационная работа обобщает результаты, полученные автором в цикле исследований, выполненном с 1985 по 2007 годы, посвященном процессам переноса атомов водорода в гетероатомных элементоорганических системах, включающих атомы, резко различающиеся по атомным массам и электроотрицательностям. В качестве объектов исследования выбраны две системы такого типа - таутомерия гидрофосфорильных соединений и каталитические реакции, протекающие с переносом водорода на металлоорганических комплексах с одним и двумя атомами переходных металлов.

Переносы атомов водорода между атомами кислорода и фосфора изучены на примере прототропной таутомерии большой серии гидрофосфорильных соединений. Эти соединения занимают особое положение в химии фосфорорганических соединений. Они обладают способностью претерпевать в растворах таутомерный переход между двумя формами - вторичным фосфиноксидом 1 и фосфинистой кислотой 2, которая, в свою очередь, может существовать в виде двух ротамеров - транс-2 (2а) и цис-2 (2б). Взаимопревращение 1 2 обеспечивает ГФС уникальное сочетание свойств и реакционной способности производных пятивалентного фосфора и производных трехвалентного фосфора.

(1)

1 2а 2б

R, R' = Алкил, арил, алкокси, арокси.

Переносы атомов водорода в системах, включающих один и два атома позднего переходного металла, изучены автором на модельных реакциях окислительного присоединения молекулярного водорода и алканов (метан, этан) к атому палладия, кластеру Pd2, металлоорганическим комплексам палладия.

Основные положения, выносимые на защиту:

ѕ модели образования сольватных комплексов гидрофосфорильных соединений с молекулами воды за счет водородных связей О-Н•••О и О-Н•••Р, которые предсказывают дополнительную стабилизацию фосфиноксидов в случае донорных заместителей на фосфоре и, напротив, дополнительную стабилизацию фосфинистых кислот в случае электроноакцепторных заместителей на фосфоре;

ѕ общие теоретические представления о димеризации фосфинистых кислот путем формирования одной водородной связи или О-Н•••Р или двух водородных связей О-Н••Р, а также о димеризации фосфиноксидов при образовании двух, трех и четырех водородных связей С-Н•••О в случае донорных заместителей на фосфоре, либо димеризации за счет диполь-дипольного взаимодействия в случае акцепторных заместителей;

ѕ теоретическое обоснование бимолекулярного механизма таутомерной перегруппировки гидрофосфорильных соединений путем синхронного переноса двух атомов водорода в цепочках водородных связей О-Н••Р в димерах;

ѕ оптимальный реакционный путь в каталитическом цикле гидроформилирования алкенов на модельном комплексе Pt(II) {[(Н2PO)2H]Pt(PH3)(Н)]} с бидентатным лигандом, образованным двумя молекулами фосфинистой кислоты, связанными водородной связью О-Н•••О;

ѕ основные факторы, ответственные за высокую каталитическую активность и селективность плоских квадратных 16е комплексов Pt этого типа в реакциях гидроформилирования алкенов;

ѕ общие закономерности окислительного присоединения и последующей активации С-Н связей простейших алканов на атомах, кластерах и комплексах палладия, включая оценки туннельных вкладов в константы скоростей в рамках подхода гамильтониана реакционного пути;

ѕ оптимальный реакционный путь в каталитическом цикле гидрирования этилена на кластере Pd2 и особенности протекания этой реакции.

Научная новизна.

Методом функционала плотности (РВЕ, ТZ2p базис) детально исследовано таутомерное равновесие (1) для представительного ряда гидрофосфорильных соединений RR'P(H)О (1) RR'POH (2) (R,R'=H, CH3, CF3, C6H5, C6F5, ОCH3, N(CH3)2, N(CF3)2 в различных сочетаниях) и выявлены эффекты электронодонорных и акцепторных заместителей, а также образования водородных связей с протонодонорными молекулами на его положение. Показано, что внутримолекулярный перенос протона в этом таутомерном равновесии запрещен, и для его описания необходимо рассматривать иные модели с участием протонодонорных молекул переносчиков. Установлено, что в инертных неполярных растворителях ГФС образуют прочные димерные ассоциаты различного строения за счет водородных связей типа О-Н···О, О-Н···Р и С-Н···О, в то время как связи Р-Н неполярны и в образовании водородных связей не участвуют. Синхронный перенос двух атомов водорода в димерах ГФС симметричного строения осуществляется через шестичленные плоские переходные состояния с низкими активационными барьерами, подобно тому, как это происходит в димерах карбоновых кислот. Таким образом в работе разрешен вопрос о механизме этого важного для химии ГФС процесса. В диссертации построена также динамическая модель внутреннего вращения вокруг связей Р-О в фосфинистых кислотах (2), приводящего к взаимопревращению их цис- и транс-конформеров.

В работе проведён анализ ППЭ реакций окислительного присоединения водорода, метана и этана по С-Н связи к атому Pd, кластеру Pd2 и дифосфиновому и этилендифосфиновому комплексам палладия в приближении ГРП с использованием метода ФП. Рассчитаны реакционные вероятности и термические полуклассические константы скоростей этих процессов. Оценены туннельные и активационные вклады в константы скоростей. Выявлены общие закономерности этих реакций: преобладание туннельного механизма при низких температурах для активации С-Н и Н-Н-связей, и бльшая активность и селективность кластера Pd2 по сравнению с атомарным палладием за счёт более интенсивного туннелирования, более низкой энергии активации и бльшего теплового эффекта реакции. Проведена оценка Н/D-кинетических изотопных эффектов в этих реакциях.

Теоретическое изучение таутомерии ГФС и механизма гидроформилирования алкенов на комплексах ГФС выполнено методом функционала плотности, с помощью которого были детально изучены стационарные точки ППЭ и построены связывающие их внутренние координаты реакции. При изучении этих систем были также построены динамические модели с использованием формализма гамильтониана реакционного пути (ГРП), что впервые позволило оценить вклады туннелирования в константы скорости.

Теоретическая и практическая значимость работы

Высокая каталитическая активность комплексов ГФС с переходными металлами побудила нас провести теоретическое моделирование механизмов гидроформилирования и гидрирования алкенов на комплексе Pt(II) с ГФС {[(Н2PO)2H]Pt(PH3)(Н)]} (3), в котором две молекулы фосфинистой кислоты образуют во внутренней координационной сфере металла необычный хелатирующий бидентатный моноанионный лиганд. При этом было показано, что протон в цепочке водородной связи -P-OH....О=Р- легко мигрирует от одного атома кислорода к другому и тем самым обеспечивает тонкую подстройку электронной плотности в каталитическом центре на каждом шаге реакции, выполняя таким образом роль молекулярного переключателя.

Полученные новые результаты теоретического анализа открывают путь к прогнозированию свойств ГФС в качестве лигандов в металлокомплексном катализе при создании новых эффективных каталитических систем.

Апробация работы. Результаты работы докладывались на Конференции «Современные проблемы АТР» (Владивосток, 1994), XVIII научной сессии Российского семинара по химии порфиринов и их аналогов. (Иваново, 1998), Всероссийских научно-технических конференциях ТОВМИ им. С.О. Макарова. (Владивосток, 1998-2006), Третьем международном симпозиуме “Chemistry and Chemical Education” (Владивосток , 2003), международной конференции «VII VOEVODSKY CONFERENCE “PHYSICS AND CHEMISTRY OF ELEMENTARY CHEMICAL PROCESSES” (Черноголовка, 2007); международной конференции “X INTERNATIONAL CONFERENCE ON THE PROBLEMS OF SOLVATION AND COMPLEX FORMATION IN SOLUTIONS” (Суздаль, 2007).

Публикации. По теме диссертации опубликовано 32 работы.

Структура и объём диссертации. Диссертация состоит из введения, четырех глав, выводов и списка цитируемой литературы. Работа изложена на 244 страницах, содержит 46 рисунков и 35 таблиц. Список цитируемой литературы включает 237 наименований.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении к диссертации обосновывается актуальность темы и формулируются цели и задачи, решаемые в диссертационном исследовании.

Глава 1. Методы теоретического исследования механизмов миграции атомов водорода

В этой главе дан краткий обзор квантовохимических методов, используемых для расчета молекулярного строения элементоорганических и металлоорганических соединений и процессов переноса атомов водорода в гетероатомных системах. Описан формализм гамильтониана реакционного пути и методы оценки туннельных вкладов в константы скорости химических реакций. На основе сравнительного анализа дано обоснование применения метода функционала плотности для выполнения диссертационного исследования.

Квантовохимические расчеты выполнены методом функционала плотности (ФП) с помощью комплекса программ ПРИРОДА при использовании градиент-корректированного функционала РВЕ и полноэлектронного трехэкспонентного набора базисных функций. Энергии стационарных точек ППЭ были рассчитаны с учетом поправки на нулевые колебания (ZPVE), определение типа стационарных точек проводилось путем анализа гессианов. Для проверки сопряжения обнаруженных переходных состояний, исходных соединений и продуктов использовался метод построения внутренней координаты реакции (IRC). Пространственное строение всех исследованных молекул и димеров полностью оптимизировалось. Для расчета термодинамических функций использовались статистические формулы жесткого ротатора и гармонического осциллятора. Вид граничных орбиталей ГФС рассчитан с помощью программы MOLDEN. Заряды на атомах рассчитывались по Хиршфельдту.

Для расчётов, кроме метода ФП, использовался спин-ограниченный метод ХФР с учётом корреляционных эффектов по теории возмущений Меллера-Плессета II порядка (ОХФ-МП2). Такие расчёты произведены по программе GAMESS.

Глава 2. Исследование молекулярного строения гидрофосфорильных соединений, их гидратов и димеров, и механизма переноса протонов в них методом функционала плотности

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

2.1 Влияние заместителей у атома фосфора на положение таутомерного равновесия

Положение таутомерного равновесия (1) зависит от электронных свойств заместителей на атоме фосфора: согласно данным эксперимента донорные заместители смещают равновесие в сторону формы оксида 1, а сильные электроноакцепторные (R=R'=CF3) - в сторону формы фосфинистой кислоты 2. Имеются примеры экспериментального наблюдения равновесия обеих таутомерных форм (R=R'=C6F5). В диссертации приведены результаты расчетов молекулярного строения большой группы ГФС, а также переходных состояний для их перегруппировок по мономолекулярному и бимолекулярному механизму. Всего было исследовано более 40 пар таутомеров ГФС вида:

RR'P(H)OH RR'P(O)H2 (2)

(RR')N-P(H)OH (RR')N-P(O)H2 (3)

(CH3)(RR'N)POH (CH3)(RR'N)PHO (4)

(RR?N)(R?R??N)POH (RR?N)(R?R??N)PHO (5)

где R, R', R?, R??- метильные группы, в которых атомы водорода последовательно заменялись на атомы фтора, от CH3 до CF3. При этом исследовались тенденции изменения геометрических, электронных и энергетических параметров молекулярного строения.

В таблице 1 приводятся рассчитанные значения энергий таутомеризации для 8 пар таутомеров ГФС в сравнении с известными данными эксперимента и расчетов. Из данных таблицы видно, что наш расчет методом ФП с функционалом PBE хорошо воспроизводит экспериментальные тенденции при согласии с результатами расчетов методами ФП и ab initio. Можно видеть, что индуктивный эффект и эффект сопряжения дают приблизительно равные вклады.

Наблюдается хорошее согласие рассчитанных геометрических параметров с экспериментом. Расхождение длин связей не превышает 0,04 A. Дипольные моменты молекул оксидов во всех случаях существенно выше, чем у таутомерных фосфинистых кислот, что обусловлено высокой полярностью связи Р>О. Очевидно, что полярные растворители будут стабилизировать оксиды и смещать в их сторону положение таутомерного равновесия. А влияние протонодонорных растворителей на это равновесие необходимо рассматривать отдельно, поскольку основным фактором при этом может стать образование водородных связей.

Таблица 1.Относительные энергииа таутомеров ГФС (ккал/моль) по данным квантово-химических расчетов (условия газовой фазы)

Заместители	RR'P-OH (цис)	RR'P-OH (транс)	Примечание
R	R'	ДEo	ДG298	ДEo	ДG298
H	H	-3,64	-3,82	-3,85	-4,17
H	H	-2,47		-2,25		Расчет МР4/6-311G** (Kwiatkowski, Leszczyn?ski 1992)
H	H		0,23		-1,14	Расчет В3LYP (Chesnut 2000)
H	H	1,15		1,00		Расчет cc-pVTZ/CCSD(T) (Wesolovski 2002)
CH3	CH3	4,72	4,68	3,72	3,53
CF3	CF3	-8,18	-7,77	-6,94	-6,64
CF3	CF3	-3,35				Расчет B3PW91/6-311G(3d,p) (Hoge 2005)
C6H5	C6H5	3,12	3,50	2,36	2,01
C6F5	C6F5	-3,19	-2,84	-1,24	-1,30
C6F5	C6F5	- 0,41				Расчет B3LYP/6-311G(2d,p) (Hoge 2005)
CH3	N(CH3)2	3,81	3,74	2,41	2,32
CH3	N(CF3)2	-1,83	-1,74	-0,89	-0,98
ОCH3	ОCH3	3,71	3,58	4,14	3,96
ОCH3	ОCH3	6,5				Эксперимент (Pietro, Hehre 1982)

За нуль отсчета принята энергия соответствующего фосфиноксида, ДЕ0 - относительная энергия ротамеров кислот с учетом поправок на энергию нулевых колебаний, ДG298 -свободная энергия Гиббса при Т = 298 К

2.2 Внутримолекулярный перенос протона в гидрофосфорильных соединениях

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

Мы рассчитали высоту активационных барьеров для мономолекулярного механизма более 20 таутомерных перегруппировок гидрофосфорильных соединений в таутомерных парах (5) и (6). Активационные энергии находятся в интервале 43,3 - 52,2 ккал/моль. Переходное состояние для мономолекулярного механизма переноса атома водорода содержит напряженный структурный треугольный фрагмент Р-О-Н:

Для модельной перегруппировки простейшего фосфиноксида H3POH2POH был найден еще более высокий активационный барьер: 60 ккал/моль (ab initio, Wesolowski 2002).

Это обстоятельство заставляет рассматривать альтернативные межмолекулярные механизмы такой перегруппировки с участием молекул-переносчиков протона. В качестве таковых могут выступать молекулы протонодонорного растворителя, например воды или спиртов. В отсутствие протонодонорных молекул в неполярных апротонных растворителях перегруппировка должна осуществляться как бимолекулярная реакция с участием не менее двух молекул ГФС. Вследствие этого изучение комплексов ГФС с протонодонорными молекулами, а также их димеров, образующихся за счет водородных связей, представляется весьма актуальным.

2.3 Исследование строения модельных моногидратов гидрофосфорильных соединений

В этом разделе приводятся результаты теоретического анализа структур модельных комплексов, которые образуют диметилфосфинистая кислота (4) и диметилфосфиноксид (5), а также их фторированные аналоги - бис(трифториметил)фосфинистая кислота (6) и бис(трифториметил)-фосфиноксид (7) с молекулой воды. Мы исследовали влияние такого комплексообразования на положение таутомерных равновесий (6) и (7), а также изучили структуры димеров этих ГФС и механизмов переноса атомов водорода в них. Экспериментально надежно установлено, что равновесие (6) полностью смещено в сторону диметилфосфиноксида 4, а равновесие (7) - в сторону бис(трифторметил)фосфинистой кислоты 7.

 (CH3)2POH (CH3)2P(O)H (6)

4 5

(CF3)2POH (CF3)2P(O)H (7)

6 7

При образовании водородных связей энергии и узловые свойства граничных орбиталей молекул донора и акцептора имеют первостепенное значение. На рис.1 представлен вид верхних занятых молекулярных орбиталей (ВЗМО) 4a, 4b и 5. Для 6a, 6b и 7 эти орбитали имеют сходный вид (рис.2). ВЗМО оксидов локализована на фосфорильном атоме кислорода, она вырождена и содержит вклады двух неподеленных электронных пар. В фосфинистых кислотах ВЗМО локализована на атоме фосфора и представлена его неподеленной парой. Основные вклады в нее дают 3p- и 3s-АО фосфора, но имеется также небольшой вклад р-орбиталей кислорода, а в случае бис(трифторметил)-фосфиноксида небольшой вклад дают также 2р-АО фтора. Такая структура ВЗМО способствует образованию водородных связей типа Р>О···Н-Х для оксида и Х-Н···Р либо Х-Н··О для фосфинистой кислоты с протонодорными молекулами Н-Х.

Из данных табл. 2, в которой приведены рассчитанные значения энергий граничных орбиталей еВЗМО и еНВМО для изучаемых соединений, видно, что замещение донорных метильных групп СF3-группами при фосфоре приводит к понижению энергии ВЗМО как оксидов, так и кислот на 0.5 - 1.5 эВ.

цис-(CН3)2РОН (4a)	транс-(CН3)2РОН (4b)	(CН3)2Р(О)Н (5)

Рис. 1. Структура верхних занятых молекулярных орбиталей 4a, 4b и 5



цис-(CF3)2РОН (6a)	транс-(CF3)2РОН (6b)	(CF3)2Р(О)Н (7)

Рис. 2. Структура верхних занятых молекулярных орбиталей 6a, 6b и 7

Разность энергий ВЗМО и НВМО достаточно велика (4 ? 5 эВ в кислотах и 4 ? 6 эВ в оксидах), поэтому учет электронных возбуждений для анализа таутомерных превращений ГФС и других их реакций, протекающих в термических условиях, не требуется.

Таблица 2. Энергии граничных орбиталей еВЗМО и еНВМО и их разность Де (эВ) в таутомерах гидрофосфорильных соединений

Таутомер ГФС	еВЗМО	еНВМО	Де
4a	-4.93	-0.58	4.35
4b	-4.62	-0.50	4.12
5	-6.10	-0.42	5.69
6a	-7.30	-2.17	5.13
6b	-6.89	-2.22	4.66
7	-7.73	-1.75	5.98

Большинство экспериментальных данных о таутомерном равновесии (1) получено для жидкостей или растворов. Ранее вопрос о зависимости положения этого таутомерного равновесия (1) от растворителя исследовался на теоретическом уровне лишь на примере простейшего фосфиноксида H3PO и его транс-таутомера Н2POH. Модель поляризованного континуума предсказывает относительную стабилизацию фосфиноксида по сравнению с кислотой на 3,5 ккал/моль, при учете образования водородной связи фосфиноксида с одной молекулой воды, выступающей в качестве донора протона, относительная стабилизация оксида составила 4.0 ккал/моль (Chesnut, 2000). Сольватация H3PO и кислоты Н2POH молекулами воды исследовалась методом Монте-Карло. Относительная стабилизация фосфиноксида составила около 9 ккал/моль (Wesolowski, 2002).

Мы исследовали относительную стабильность комплексов, образуемых диметилфосфиноксидом (CH3)2P(O)H (5) и диметилфосфинистой кислотой (CH3)2POH (4), а также их фторированными аналогами (CF3)2P(O)H (7) и (CF3)2POH (6), с молекулой воды с целью выявить вклад водородных связей различного типа в относительную устойчивость таутомеров этих ГФС.

Структуры наиболее устойчивых моногидратов (CH3)2P(O)H и (CH3)2POH, в которых вода выступает как донор водородной связи (8 и 9), и данные по энергиям их образования ДЕ представлены на рис. 3.

8	9a	9b
R(O···H) = 1.830 A	R(O···H) = 1.940 A	R(O···H) = 1.896 A
ДЕ = -9.58 ккал/моль	ДЕ = -6.09 ккал/моль	ДЕ = -6.26 ккал/моль

Рис. 3. Структуры моногидратов диметилфосфиноксида и таутомерных цис- и транс-форм диметилфосфинистой кислоты, в которых молекула воды играет роль донора водородной связи

Координация воды с оксидом (CH3)2P(O)H проходит по кислороду фосфорильной группы. Несмотря на то, что ВЗМО фосфинистых кислот локализована на атоме фосфора, координация молекулы воды с ними происходит за счет неподеленных электронных пар кислорода группы Р-ОН: жесткий протон координируется по более жесткому кислородному центру. Детальное исследование ППЭ для моногидрата цис-ротамера этой кислоты в окрестности, где возможна координация молекулы воды по атому фосфора, позволило обнаружить локальный минимум, соответствующий структуре 9с (рис. 4).

9с
R(H···P) = 2.435 A
ДE = -4.95 ккал/моль

Рис. 4. Структура моногидрата цис-ротамера диметилфосфинистой кислоты, в котором молекула воды играет роль донора водородной связи при координации по атому фосфора

Большое расстояние R(H···P) = 2.435 A и относительно небольшая энергия гидратации ДE = -4.95 ккал/моль указывают на образование слабой водородной связи О-H···P.Вода как донор водородной связи стабилизирует оксид сильнее, чем кислоту, на 3.32 ккал/моль за счет образования более сильной водородной связи O···H-О по фосфорильной группе, где атом кислорода несет более высокий отрицательный заряд ( qO= -0.35e), чем атом кислорода кислоты (qO= -0.22e).

Структуры наиболее устойчивых моногидратов 10 - 11, в которых вода выступает по отношению к (CH3)2P(O)H и (CH3)2POH как акцептор водородной связи, донируя для этого свою электронную пару, представлены на рис. 5.

Связь Р-Н в 5 неполярна, поэтому нам не удалось обнаружить устойчивой структуры ассоциата молекул (СН3)2P(O)H и воды, в котором атом кислорода воды координируется по атому водорода этой связи. Единственный локальный минимум ППЭ в этом случае соответствует ассоциату 10, в котором атом кислорода молекулы воды сближен с атомами водорода двух метильных групп на расстояние R(O···H) = 2.55 A, что можно рассматривать как образование слабых водородных связей С-Н···O за счет дефицита электронной плотности на атомах водорода метильных групп (qн=0.06e). Ротамеры фосфинистой кислоты как доноры водородной связи образуют с водой значительно более прочные ассоциаты.

10	11a	11b
R(O···H) = 2.55 A	R(O···H) = 1.870 A	R(O···H) = 1.928 A
ДЕ = -6.27 ккал/моль	ДЕ = -6.16 ккал/моль	ДЕ = -8.16 ккал/моль

Рис. 5. Структуры моногидратов диметилфосфиноксида, а также цис- и транс-ротамеров диметилфосфинистой кислоты, в которых молекула воды играет роль акцептора водородной связи

Ассоциат транс-ротамера 11b на 2.0 ккал/моль более устойчив, чем ассоциат цис-ротамера 11a.

Из полученных результатов следует, что из всех ассоциатов, образуемых таутомерами (CH3)2P(O)H и (CH3)2POH с молекулой воды, наиболее устойчивым является моногидрат оксида 8. Вследствие этого относительная устойчивость оксида (CH3)2P(O)H 5 по отношению к фосфинистым кислотам (4a, b) возрастет. Этот результат находится в качественном согласии с результатами неэмпирических расчетов, полученными для пары фосфиноксид H3PO - транс-фосфинистая кислота Н2POH (Chesnut, 2000). Модели ассоциатов с молекулой воды были построены также для 6 и 7 (рис. 6 и 7).

Из данных по энергии гидратации следует, что способность этих соединений образовывать ассоциаты с молекулой воды в качестве донора водородной связи сильно понижается по сравнению с таковой для их нефторированных аналогов, что выражается как в понижении энергий связывания, так и в увеличении расстояний О···Н.

12	13a
R(O···H) = 2.139 A	R(O···H) = 2.133 A
ДЕ = -6.70 ккал/моль	ДЕ = -2.84 ккал/моль

Рис. 6. Структуры моногидратов бис-(трифторметил)фосфиноксида и цис- ротамера бис(трифторметил)фосфинистой кислоты, в которых молекула воды играет роль донора водородной связи

Для транс-ротамера 6b не удалось обнаружить минимума, соответствующего ассоциату этого типа, а цис-ротамер 6a образует слабый ассоциат 13a с энергией связывания 2.84 ккал/моль. Наиболее прочный ассоциат 12 образует молекула оксида, но энергия связывания в нем (-6.70 ккал/моль) существенно понижена по сравнению с таковой в 8.

Напротив, замена атомов водорода электроноакцепторными атомами фтора увеличивает способность ротамеров HOP(CF3)2 выступать по отношению к молекуле воды в качестве доноров водородной связи. Из данных по энергиям связывания ДЕ (рис. 7) видно, что наиболее прочный ассоциат 15a образует цис-ротамер кислоты, но короткие (1.664 A и 1.705 A) и прочные водородные связи этого типа образуют оба ротамера. Оксид OPH(CF3)2, в отличие от 5, способен образовыывать непрочный комплекс 14 за счет слабой водородной связи Р-H···О, поскольку введение двух мощных электроноакцепторных групп CF3 приводит к заметной поляризации связи Р-Н и появлению на атоме водорода небольшого положительного заряда.



14	15a	15b
R(O···H) = 2.136 A	R(O···H) = 1.664 A	R(O···H) = 1.705 A
ДЕ = -4.47 ккал/моль	ДЕ = -10.79 ккал/моль	ДЕ = -10.50 ккал/моль

Рис. 7. Структуры моногидратов 6 и 7, в которых молекула воды выступает в качестве акцептора водородной связи

Детальное изучение ППЭ для моногидратов этих кислот при координации атома водорода молекулы воды по атому фосфора показало, что локальный минимум, сходный по структуре с 9c (рис.4), в этом случае отсутствует. Очевидно, что в данном случае при наличии мощных электроноакцепторных заместителей на атоме фосфора электронная плотность на нем понижена настолько, что становится невыгодным образование даже слабой водородной связи О-Н···Р.

Из представленных данных следует, что из всех таутомеров OPH(CF3)2 наиболее прочный сольватный комплекс с молекулой воды образует цис-ротамер бис(трифторметил)фосфинистой кислоты (10.79 ккал/моль). Эта форма и будет в наибольшей степени стабилизироваться за счет образования водородных связей с полярными донорными растворителями. Этот результат согласуется с экспериментальными данными по влиянию растворителей на положение таутомерного равновесия бис-(пентафторфенил)фосфинистой кислоты и соответствующего оксида (Hoge, 2005).

2.4 Строение димеров гидрофосфорильных соединений и механизмы переноса атомов водорода в них

По данным наших теоретических расчетов, ГФС образуют несколько типов димеров c одной и двумя водородными связями, которые описаны в диссертационной работе. Здесь мы представляем результаты детального изучения структуры димеров симметричного строения и механизма синхронного переноса протонов в них.

Две молекулы 4b образуют достаточно прочный димер 16 (ДЕb = 9.07 ккал/моль) с двумя водородными связями О-Н···Р, в котором оба фрагмента Р-О-Н лежат в одной плоскости (рис. 8).

	R(Р-O) = 1.679 A R(P···H) = 2.346 A R(O-H) = 0.995A
16	ДЕb = -9.07 ккал/моль

Рис.8. Строение и некоторые параметры димера 16

Тщательный анализ ППЭ показал, диметилфосфиноксид V не образует димеров симметричного строения с водородными связями О···Н-Р, поскольку, как уже отмечалось выше, связь Р-Н в 5 неполярна и не способна участвовать в образовании таких водородных связей. Однако на ППЭ удается зафиксировать несколько локальных минимумов, соответствующих димерным ассоциатам, в которых фрагменты Р(О)Н не копланарны, и вместо водородных связей Р-Н···О образуются межмолекулярные водородные связи С-Н···О с атомами водорода метильных групп. Структуры и некоторые параметры димеров этого типа 17 - 19 представлены на рис. 9.

В димере 17 имеются две водородные связи С-Н···О. Энергия димеризации ДЕ = 6.64 ккал/моль, расстояния R(Н···О) = 2.38 A. Известно (Steiner, 2002), что среднее расстояние R(Н···О) при образовании водородных связей С-Н···О составляет около 2.36 A, а энергия таких связей может достигать 4 ккал/моль.

В образование димера 17 вносит некоторый вклад также диполь-дипольное взаимодействие антипараллельно расположенных поляризованных связей Р>О, однако по силе оно уступает взаимодействиям С-Н···О.

17	18	19
R(О···H) = 2.38 A R(О···Р) = 3.30 A	R(О···H1) = 2.34 A R(О···H2) = 2.45 A R(О···Р) = 3.45 A	R(О···H) = 2.45 A R(О···Р) = 3.50 A
ДЕ = -6.64 ккал/моль	ДЕ = -7.46 ккал/моль	ДЕ= -9.36 ккал/моль

Рис. 9. Структуры димеров оксида (СН3)2Р(О)Н с образованием нескольких водородных связей С-Н···О

Мы обнаружили также димеры с тремя (18) и четырьмя (19) водородными связями С-Н···О. В 18 расстояния R(О···Н1) = 2.34 A для связи с одной метильной группой и R(О···Н2) = 2.45 A для связей с двумя метильными группами. В 19 все четыре расстояния R(Н···О) = 2.45 A, и энергия димеризации ДЕ составляет 9.36 ккал/моль. Всем димерам молекул диметилфосфиноксида соответствуют локальные минимумы на ППЭ, что подтверждено расчетом градиентов и частот нормальных колебаний.

Рис. 10. Расположение связевых критических точек типа (3,-1) по Бейдеру в структуре димера диметилфосфиноксида 18. Видно наличие трех межмолекулярных водородных связей С-Н···О.

Вопрос о том, в какой мере относительно слабые (1 - 4 ккал/моль) межмолекулярные взаимодействия С-Н....О могут рассматриваться как нормальные водородные связи, был предметом оживленных дискуссий. Обстоятельный теоретический анализ с использованием неэмпирических расчетов показал, что такое описание вполне обосновано. Мы провели теоретический AIM-анализ структуры димера 18 по Бейдеру и установили (рис. 10) наличие в ней трех связевых критических точек типа (3,-1) (f=0.0097, h=(-0.0088, -0.0078, 0.0468)), соответствующих образованию трех водородных связей С-Н....О.

Транс-ротамер бис(трифторметил)фосфинистой кислоты 6b образует симметричный димер 20 с двумя водородными связями О-Н···Р, в котором оба фрагмента Р-О-Н лежат в одной плоскости. В образование симметричного димера бис(трифторметил)фосфиноксида 21 помимо двух относительно слабых взаимодействий О-Н···Р основной вклад дает диполь-дипольное взаимодействие антипараллельно ориентированных полярных групп Р>О (рис. 11).

20	21
R(Р-O) = 1.655 A R(P···H) = 2.501 A R(O-H) = 0.986 A	R(О···H) = 2.523 A R(PH) = 1.423 A R(O=Р) = 1.505 A
ДЕ = -6.36 ккал/моль	ДЕ = -4.04 ккал/моль

Рис. 11. Структуры и некоторые параметры симметричных димеров 20 и 21

Координата реакции для перегруппировки симметричного димера диметилфосфинистой кислоты 16, протекающей с синхронным переносом обоих протонов в цепочках водородных связей Р-О-Н···Р (уравнение 8), представлена на рис. 12.

(СН3)2РОН + (СН3)2РОН > (СН3)2Р(О)Н + (СН3)2Р(О)Н (8)

Она отчетливо распадается на две ветви. Сначала симметричное переходное состояние ПС1 достигается при плавном возрастании энергии за счет увеличения межатомных расстояний R(O-H) от 0.995 A в 16 до 1.344 A в ПС1. Отметим, что на этом участке реакционного пути и в переходном состоянии

	ПС1
	R(PO) = 1.593 A R(P···H) = 1.637 A R(O···H) = 1.344 A R(PC) = 1.842 A

Рис.12. Координата реакции синхронного переноса двух атомов водорода при перегруппировке димера диметилфосфинистой кислоты 16 в димер диметилфосфиноксида 17 и некоторые характеристики ПС1 сохраняется копланарность фрагментов Р-О-Н обеих молекул

Траектория спуска из ПС1 имеет интересную особенность. На начальном участке спуска копланарность фрагментов Р-О-Н сохраняется, и энергия плавно уменьшается вплоть до того момента, когда межатомные расстояния Р......Н уменьшаются до 1.431 A. Межатомные расстояния О---Н при этом увеличиваются до 2.453 A. В этой точке в гессиане системы имеется мнимая частота, соответствующая изменению двугранного угла Н-Р...Р-Н. Дальнейший спуск происходит за счет уменьшения этого угла от 1800 до 710 при практически неизменном расстоянии Р···Н. Молекулы оксида разворачиваются относительно друг друга вокруг связей Р>О, возникают 2 водородные связи С-Н···О между фосфорильными атомами кислорода и атомами водорода метильных групп, и образуется димер 17, уже описанный выше. Высота активационного барьера на пути такой бимолекулярной перегруппировки кислоты в оксид составляет ДЕ0 = 4.84 ккал/моль, что на порядок меньше высоты барьера для мономолекулярного переноса атома водорода. Таким образом, расчет однозначно показывает, что синхронный перенос двух протонов в димере диметилфосфинистой кислоты 16 осуществляется очень легко, что и объясняет протекание таутомерного перехода (2) в мягких условиях в отсутствие молекул-переносчиков.

Перегруппировка димера бис(трифторметил)фосфинистой кислоты 20 в димер бис(трифторметил)фосфиноксида 21 (уравнение 9) с синхронным переносом обоих протонов в цепочках водородных связей Р-О-Н···Р, протекает через переходное состояние ПС2.

(СF3)2РОН + (СF3)2РОН > (СF3)2Р(О)Н + (СF3)2Р(О)Н (9)

Координата реакции для этой перегруппировки представлена на рис. 13. Она адекватно описывается изменением одного параметра - межатомных расстояний в цепочках Р-О-Н···Р. Плоскость симметрии, в которой лежат фрагменты Р-О-Н···Р, сохраняется вдоль всего реакционного пути, и в этом случае не наблюдается особенности, отмеченной для перегруппировки (8), связанной с разворотом молекул, поскольку симметричному димеру 21 соответствует отчетливый минимум ППЭ. Достаточно низкий активационный барьер процесса (5.41 ккал/моль) показывает, что и в этом случае таутомерные переходы между 6 и 7 могут осуществляться в мягких условиях в отсутствие молекул-переносчиков, что соответствует экспериментальным данным.



	ПС2
	R(PO) = 1.562 A R(P···H) = 1.592 A R(O···H) = 1.408 A R(PC) = 1.923 A

Рис.13. Координата реакции синхронного переноса двух протонов в димере бис(трифторметил)фосфинистой кислоты 20 и некоторые характеристики переходного состояния ПС2

2.5 Полуклассическая оценка константы равновесия для перегруппировки фосфинистой кислоты в фосфиноксид

С целью оценки точности вычисления констант равновесия прототропной перегруппировки по уравнению изотермы реакции мы построили динамическую модель в приближении гамильтониана реакционного пути (ГРП) для модельной реакции перегруппировки фосфинистой кислоты в фосфиноксид H3POH2POH.

Согласно нашим результатам, разница в энергиях Гиббса структур фосфинистой кислоты и фосфиноксида равна -3.1 ккал/моль, разница энтальпий -2.5 ккал/моль.

На рис. 14 приведена потенциальная функция V0(s) (пунктир) и функции частот поперечных колебаний в зависимости от естественной координаты реакции s, выраженной в единицах A? (mH -- масса протона, равная 1847 а.е.м.).

Левый минимум реакционного пути соответствует двум ассоциированным молекулам фосфиноксида, правый -- двум молекулам фосфинистой кислоты. Димеризация фосфинистой кислоты сопровождается понижением энтальпии образования на 3.2 ккал/моль, фосфиноксида -- на 2.3 ккал/моль с учётом энергии нулевых колебаний.

На всём протяжении реакционного пути молекулярная система сохраняет симметрию C2h. Оба атома фосфора, оба атома кислорода и мигрирующие атомы водорода лежат в одной плоскости.

Рис. 14. Потенциальная функция перегруппировки V0(s) (пунктир) и функции частот поперечных колебаний в зависимости от естественной координаты реакции s

Из рис. 14 видно, что лишь немногие функции щ(s) претерпевают существенные изменения в ходе реакции. В частности, это частоты щ2,3(s), которые в левой части РП представляют собой симметричные и антисимметричные колебания q+(P-H) и q-(P-H), а в правой -- q+(O-H) и q-(O-H), соответственно симметрии A и B. Высота барьера прямой реакции без учёта нормальных колебаний составляет 5.6 ккал/моль, с учётом нормальных колебаний -- 4.1 ккал/моль.

Оценка констант скорости прямой kI и обратной kII реакции показывает, что для прямой реакции значения >2 при T<165K, для обратной -- при T<190K.

Значения находятся в пределах от 1,4 до 1,3 в температурном интервале от 300 до 400K. Это означает, что при использовании, например, уравнения изотермы химической реакции происходит систематическая недооценка значения константы равновесия (см. рис. 15).

Таким образом, при исследовании таутомерного равновесия между фосфинистой кислотой H2POH и фосфиноксидом H3PO при комнатной и повышенной температуре необходимо учитывать туннельный эффект, так как при использовании уравнения изотермы реакции константа равновесия оказывается недооцененной примерно на ?.

Рис. 15. Температурные функции логарифмов констант равновесия таутомерной перегруппировки фосфиноксида в фосфинистую кислоту, рассчитанных как отношения полных констант скоростей kCUM, надбарьерных вкладов с учётом надбарьерного отражения kARR и вычисленных по уравнению изотермы реакции

1 -- , 2 -- , 3 --

Глава 3. Построение и исследование моделей каталитических циклов на металлоорганических комплексах, связанных с процессами переноса атома водорода в реакциях гидроформилирования и гидрирования

3.1 Теоретические исследования механизмов переноса атома водорода в процессе гидроформилирования алкенов

Приводятся основные известные теоретические и экспериментальные данные о каталитической реакции гидроформилирования на металлорганических комплексах кобальта, родия и платины. Подобно третичным фосфинам, фосфинистые кислоты RR'POH образуют с поздними переходными металлами прочные комплексы (Roundhill, 1975), которые выступают как активные катализаторы в реакциях кросс-сочетания (Li, 2002), гидроформилирования алкенов (Van Leeuven, 1986) и гидрофосфинилирования алкинов (Tanaka, 1996). Однако, в отличие от комплексов третичных фосфинов, эти комплексы стабильны на воздухе вследствие высокой оксидативной устойчивости ГФС.

3.2 Исследование альтернативных каталитических циклов гидроформилирования этилена на комплексах платины с гидрофосфорильными лигандами

В последнее десятилетие химия ГФС получила новый мощный импульс в своем развитии. Ван Лёвен с сотр. показали, что гидроформилирование гептена-1 на комплексе [(Ph2PO)2H]Pt(PPh3)(Н)] (22) протекает с региоселективностью более 90% в очень мягких условиях (уравнение 10). Мы провели теоретическое моделирование механизма этой реакции для молекулы этилена на комплексе [(H2PO)2H]Pt(PH3)(Н)] (3), в котором для сокращения вычислительных затрат Ph-радикалы при атомах фосфора заменены атомами водорода.

Для расчетов использованы метод функционала плотности, градиент-корректированный функционал PBE, трехэкспоненциальный набор базисных функций гауссовского типа TZ2P для описания валентных электронов и псевдопотенциалы SBK-JC для остовных электронов. Энергии стационарных точек рассчитаны с учетом поправки на нулевые колебания ZPVE.

Табл.3.Некоторые геометрические параметры комплексов 3 и 23 по данным расчета и РСА (Tanaka, 1996)

	Параметр	23 X=Ph, Y= Et	23 (РСА) X= Ph, Y= Et	3 X=Y=H
	r(Pt1-P1) A	2,341	2,280	2,315
	r(Pt1-P2) A	2,366	2,313	2,353
	r(Pt1-P3) A	2,344	2,301	2,323
	r(O1-O2) A	2,455	2,317	2,470
	r(P1-O1) A	1,573	1,532	1,563
	r(P2-O2) A	1,614	1,569	1,607

С5Н11CH=CH2 + CO + H2 > С5Н11CH2CH2CHO (10)

Проверка сопряжения обнаруженных переходных состояний с исходными соединениями и продуктами осуществлялась построением координаты реакции IRC. Все расчеты выполнены с использованием программы PRIRODA.

Сравнение геометрических параметров, рассчитанных для 3 и 23 с данными РСА для 23 показывает (табл.3), что расчет удовлетворительно воспроизводит геометрию координационного узла. Прочность водородной связи P(2)H2O(2)H....О(1)Р(1) в 3 составляет 13.1 ккал/моль. Эта оценка была получена расчетом структуры при повороте группы P(2)H2OH в 3 относительно связи Pt-P(2) на 1800, которой соответствует локальный минимум на ППЭ.

Присоединение молекулы этилена к 3 относительно легко протекает как энтропийно контролируемая реакция (ДH298 = -5.15 ккал/моль, ДS298 = -41.58 э.е., ДG298 = 7.25 ккал/моль) с образованием 24, в котором этилен связан с атомом платины по з2-типу. Весьма примечательно, что при координации этилена атом водорода в цепочке P(2)-O(2)-H.....O(1)-P(1) перемещается к кислороду O(1).

{[(Н2PO)2H]Pt(PH3)(Н)]} + CH2=CH2 > {[(Н2PO)2H]Pt(PH3)(Н)з2(CH2=CH2)]} (11)

3 24

Перегруппировка 24 в з1-этильный комплекс 25 (уравнение 12) протекает через раннее переходное состояние ПС3 с энергией активации 14.15 ккал/моль.

{[(Н2PO)2H]Pt(PH3)(Н)з2(CH2=CH2)]} > {[(Н2PO)2H]Pt(PH3)CH2CH3]} (12)

24 25

24	ПС3	25

Процесс экзотермичен (ДH298 = -12.26 ккал/моль, ДG298= -13.89 ккал/моль). Внедрение этилена по связи Pt-H сопровождается обратным переносом атома водорода в цепочке водородной связи от О(1) к О(2).

Присоединение оксида углерода к 25 (уравнение 13) протекает безбарьерно (ДH298 =-10.18 ккал/моль, ДG298 =-0.41 ккал/моль) с образованием комплекса 26. Оно снова сопровождается переносом протона в цепочке водородной связи.

{[(Н2PO)2H]Pt(PH3)CH2CH3]} + СО > {[(Н2PO)2H]Pt(PH3)(CO)CH2CH3]} (13)

25 26

Вслед за этим осуществляется экзотермическая изомеризация (ДH298 =
-8.98 ккал/моль, ДG298 = -8.86 ккал/моль) комплекса 26 в плоский квадратный з1-ацильный комплекс 27 (уравнение 14), что связано с преодолением самого высокого барьера на этом реакционном пути, соответствующего переходному состоянию ПС4 (энергия активации ДEa0= 17.60 ккал/моль). Процесс также сопряжен с переносом протона в цепочке водородной связи.

26	ПС4	27

{[(Н2PO)2H]Pt(PH3)(CO)CH2CH3]}>{[(Н2PO)2H]Pt(PH3)COCH2CH3]} (14)

26 27

Окислительное присоединение молекулы водорода к 27 (уравнение 15) приводит к октаэдричесому цис-дигидриду 28. Это энтропийно контролируемый процесс (ДH298 = -2.01 ккал/моль, ДG298 = 9.69 ккал/моль, ДEa0= 12.41 ккал/моль), идущий через переходное состояние ПС5.

{[(Н2PO)2H]Pt(PH3)COCH2CH3]} + H2 > {[(Н2PO)2H]Pt(H)2(PH3)COCH2CH3]} (15)

27 28

Замыкает каталитический цикл экзотермическая (ДH298 = -0.86 ккал/моль, ДG298 = -14.25 ккал/моль) стадия восстановительного элиминирования пропаналя (уравнение 16), которая проходит через переходное состояние ПС6 с энергией активации ДEa0= 15.51 ккал/моль.

{[(Н2PO)2H]Pt(H)2(PH3)COCH2CH3]} > {[(Н2PO)2H]Pt(H)(PH3)} + C2H5COH (16)

3

ПС5	28	ПС6

Детальное изучение ППЭ и возможных альтернативных путей осуществления каталитического процесса показало, что описанный выше каталитический цикл (Рис. 16) является оптимальным по кинетическим и термодинамическим характеристикам. Он согласуется с экспериментальными данными (Van Leeuven, 1986).

Известно, что часто в качестве побочного процесса при гидроформилировании алкенов выступает их гидрирование. Оно может осуществлятся за счет окислительного присоединения молекулы водорода к 25 c последующим элиминированием молекулы этана. Расчет показал, что присоединение Н2 к 25 идет c энергией активации 14.15 ккал/моль, в результате чего образуется дигидрид 29. Процесс контролируется энтропийным фактором и эндотермичен (ДH298 = 1.44 ккал/моль, ДS298 = -32.55 э.е., ДG298 = 11.14 ккал/моль).

{[(Н2PO)2H]Pt(PH3)CH2CH3]} + Н2 > {[(Н2PO)2H]Pt(PH3)( Н)2CH2CH3]} (17)

25 29

Отщепление этана от 29 (уравнение 18) далее осуществляется с энергией активации 13.93 ккал/моль как сильно экзотермическая реакция (ДH298 = -18.61 ккал/моль, ДG298 = -30.25 ккал/моль).

Рис.16. Каталитический цикл гидроформилирования этилена на комплексе 3. Стадии, сопровождающиеся переносом протона в цепочке водородной связи P(2)H2O(2)H....О(1)Р(1) отмечены звездочкой *

{[(Н2PO)2H]Pt(PH3)( Н)2CH2CH3]} > {[(Н2PO)2H]Pt(PH3) Н]} + С2Н6 (18)

Поскольку оксид углерода присоединяется к 25 безбарьерно и экзотермично (уравнение 13), в то время как присоединение молекулярного водорода (уравнение 17) связано с преодолением заметного барьера (14.15 ккал/моль) и эндотермично, вклад гидрирования как конкурирующего процесса в мягких условиях не должен быть значительным, если в реакционной среде отсутствует избыток водорода. Это заключение вполне подтверждается экспериментальными данными (Van Leeuven, 1986).

Таким образом результаты теоретического моделирования показывают, что высокая эффективность комплексов типа 22 и 3 как катализаторов гомогенного гидроформилирования алкенов обусловлена следующими факторами:

Наличие свободного координационного места в этих плоских квадратных 16е комплексах платины обеспечивает возможность координации алкена на первой стадии без энергетически невыгодной предварительной диссоциации одной из связей металл - лиганд.

Высокая прочность водородной связи P(2)H2O(2)H....О(1)Р(1) приводит к формированию в координационной сфере металла бидентатного лиганда, который достаточно жестко фиксирует геометрию каталитического центра, что, как известно, увеличивает региоселективность процесса.

Протон в цепочке P(2)H2O(2)H....О(1)Р(1), обладая высокой подвижностью, легко перемещается в ней и тем самым обеспечивает тонкую подстройку электронной плотности в каталитическом центре на каждом шаге реакции, выполняя таким образом роль молекулярного переключателя.

Эти закономерности могут быть использованы для оптимизации структуры каталитических комплексов типа 22 и 3.

Глава 4. Окислительное присоединение водорода, метана, этана и этилена к каталитическим системам на основе палладия

В этой главе кратко изложены основные экспериментальные факты, показывающие наличие квантового туннелирования в активации С-Н и Н-Н - связей на атомах, кластерах и кристаллических поверхностях переходных металлов. Обобщены имеющиеся литературные данные по квантовохимическим расчётам активации углеводородов на переходных металлах. Описаны результаты исследований разработанных автором моделей каталитической активации связей Н-Н и С-Н.

4.1 Окислительное присоединение водорода, метана и этана к атому Pd

В рамках приближения гамильтониана реакционного пути исследованы реакции присоединения молекул водорода, метана и этана к атому палладия:

Pd + H2 H-Pd-H (19)

Pd + CH4 H-Pd-CH3 (20)

Pd + C2H6 H-Pd-C2H5 (21);

Следует упомянуть о том, что присоединение этана к палладию используется для моделирования разрыва С-С-связи, поэтому реакция (21) на неэмпирическом уровне никем до сих пор не исследовалась (разрыв С-Н связи исследуется обычно только на примере метана).