[DMAPH]Br催化CO2和环氧氯丙烷反应的机理

第40卷第5期肇庆学院学报V ol.40,No.5 2019年9月JOURNAL OF ZHAOQING UNIVERSITY Sep.2019
[DMAPH]Br催化CO2和环氧氯丙烷反应的机理
宋振艳1，冯华杰2，李志伟1
（1.肇庆学院环境与化学工程学院，广东肇庆526061；2.海南师范大学化学与化工学院，海南海口571158）
摘要：用密度泛函理论和前线分子轨道理论研究了4-二甲胺基吡啶氢溴酸盐([DMAPH]Br)催化CO2和环氧氯丙烷反应生成环状碳酸酯的反应机理.研究了2种可能的反应途径，对不同途径各个步骤的中间体和过渡态进行了理论计算，得出了最优反应途径.
关键词：离子液体；CO2；环氧氯丙烷；密度泛函理论
中图分类号：O641-O649文献标志码：A文章编号：1009-8445（2019）05-0025-06
近年来CO2加剧排放引发的全球变暖和环境恶化问题受到越来越多的关注，将大气中CO2固定转化为高附加值精细化学品，无论从资源回收利用还是从环境保护方面，都具有巨大的现实意义和应用前景.与其他常见碳源相比，CO2无毒价廉、取之不尽，被广泛应用于基础研究和应用研究领域[1-2].但相对于其巨大的储量，CO2的利用非常有限，原因在于它很高的热力学稳定性和动力学惰性，需要高效的催化剂才能将其活化.
在CO2固定活化催化剂中，离子液体(IL)有着举足轻重的地位[3-6].离子液体催化具有活性高、工艺简单、收率高、能循环利用、环境友好的优点.环碳酸酯是CO2固定转化中最流行的目标产物之一[7-9]，为更深入地理解此类反应，以便开发高效、价廉、绿色、环保的离子液体催化体系，有必要对催化反应机理、反应途径，离子液体在反应中的作用进行剖析[10-11].最近，Zhang[12]等报道了离子液体[DMAPH]Br催化CO2与环氧氯丙烷生成环碳酸酯的反应，产物收率和选择性都很高，为明确其反应机理，通过理论计算从微观角度进行探讨.
1计算方法
首先在B3LYP/6-31+G(d,p)计算水平下，对反应过程中不同步骤的反应物、过渡态、中间体和产物进行结构优化，然后对结构进行频率计算，确保反应物、中间体、产物是稳定的，而过渡态有且只有一个合理的虚频.为确保所有过渡态的准确性，即可通过势能面连接相应的反应物和产物，用内禀反应坐标(IRC)方法对反应路径进行验证.此外，根据前线分子轨道理论分析反应物间的成键情况.所有的计算都在Gaussian09[13]程序包中完成，能量均采用校正后的吉布斯自由能.
2结果与讨论
2.1反应机理
计算体系如图1所示，离子液体催化剂为4-二甲胺基吡啶氢溴酸盐.首先提出2个反应途径A和B，反应物、中间体、过渡态及产物的结构分别如图2和图3所示，而能量如表1和图4所示，相对能量的零点为相同
收稿日期：2019-01-26
基金项目：肇庆学院大学生创新创业训练计划项目（201710580109）；广东省自然科学基金项目（2015B030313705）
作者简介：宋振艳（1995-），女，陕西榆林人，肇庆学院环境与化学工程学院2015级本科生.
通信作者：李志伟（1977-），男，湖南汝城人，肇庆学院环境与化学工程学院教授，博士.
肇庆学院学报第40卷
计算水平下所得的环氧氯丙烷、CO 2以及催化剂[DMAPH]Br 的能量之和
.
图1[DMAPH]Br 催化CO 2和环氧化合物的环加成反应
路径A 中，环氧化物通过其氧原子与[DMAPH]Br 的质子之间形成氢键而活化，导致C-O 键极化，[DMAPH]Br 中的Br 攻击环氧化物中与氯甲基相连的碳原子(C1)，经过相对能量为115.12kJ/mol 的过渡态TS1A ，使其开环并形成中间体A1，这一步的活化能为70.64kJ/mol ，放热91.47kJ/mol.在TS1A 中，将要断键的O1-C1被拉长至0.1714nm ，将要成键的O1-H1，Br-C1分别缩短至0.1394和0.2762nm ，TS1A 有唯一虚频(-439.36icm –1)，其振动方式为O1和C1伸长、H1和O1缩短、C1和Br 缩短，正是断键和成键的方向，可见这一步骤的过渡态是正确的.
随后A1的O 原子与CO 2的C 原子结合，越过75.52kJ/mol 的能垒，经过过渡态TS2A ，形成中间体A2，这一步骤吸热67.39kJ/mol.在TS2A 中，O1与C3的距离缩短至0.1782nm ，逐渐靠近成键，H1与O1的相互作用减弱，键长拉长至0.1334nm ，该过渡态有且只有一个虚频(-317.07icm –1)，振动方式为O1和C1靠近成键、H1和O1远离断键的方向，最终生成中间体A2，成功捕获CO 2分子，同时Br-C1键长增加，作用减弱.
中间体A2继续反应，经过过渡态TS3A ，能垒为101.68kJ/mol ，闭环形成环状碳酸酯，放热45.92kJ/mol.在TS3A 中，Br-C1键长增加至0.2574nm ，O2和C1的距离缩短至0.1935nm ，且有唯一虚频(-418.47icm –1)，振动方式为Br 与C1伸长断键、C1和O2靠近成键的方向，最终闭环生成产物，同时离子液体催化剂[DMAPH]Br 复原.
表1
相关中间体的相对吉布斯自由能(ΔG)结构A0TS1A A1A1+CO 2TS2A A2TS3A PC
虚频/icm –1
-439.36-317.07-418.74ΔG/(kJ·mol -1)
44.48115.12-46.99-35.9439.5831.45133.13-14.47
结构B0TS1B B1B1+CO 2TS2B B3TS3B PC
虚频/icm –1
-378.98-246.11-434.83ΔG/(kJ·mol -1)
41.10115.54-54.36-44.7234.3538.82109.89-14.47
26
第5期宋振艳等：[DMAPH]Br 催化CO 2和环氧氯丙烷反应的机理
+
O1
C1
C2H1O1-H1
0.2021O1-C10.1453Br-C10.3679N-H10.1020
H1
C1
C2O1
O1-H10.1394O1-C10.1714Br-C10.2762N-H10.1131
Epichlorohydrin [DMAPH]Br
TS1A
H1O1C1
C2O1-H10.0990O1-C10.2445Br-C10.1997N-H10.1824
2
H1O1
C1C2C3
O1-H10.1039O1-C30.2004N-H10.1484?C3160.61°C3
H1
O1C1
C2O1-H10.1334O1-C30.1782N.-H10.1169?C3147.00°TS2A
A1
O1-H10.1621O1-C30.1534O2-C10.3091N-H10.1060?C3137.14°H1O1
C2C1
C3
O2
C1
O1
C3O2
O1-C30.1534O2-C10.1935
O1
C3O2C1
O1-C30.1360O2-C10.1443
PC
TS3A A0
Br
N
N
Br N
Br
N
Br
N
Br
N
A2
A1+CO 2
图2
反应路径A 的反应物、中间体、过渡态和产物的结构和参数(nm 或
°)
图3反应路径B 的反应物、中间体、过渡态和产物的结构和参数(nm 或°)
27
肇庆学院学报第40卷28
2.2路径B的反应机理
与A不同，反应路径B中的Br攻击了环氧化物中与氧原子相连的另一个碳原子(C2).首先，由[DMAPH]Br 的H原子使环氧氯丙烷中的氧原子活化，而Br离子亲核攻击C2原子，在2种因素的协同作用之下，经过能垒为74.44kJ/mol的过渡态TS1B生成B1，放热95.46kJ/mol.在TS1B中，H1和O1相互靠近，距离缩短至0.1428nm，Br与C2的距离缩短至0.2805nm，均有成键的趋势，而O1-C2键被拉长至0.1655nm，有断键趋势.过渡态TS1B有唯一虚频(–378.98icm–1)，且振动沿着H1与O1、Br与C2靠近成键，O1-C2键伸长断开的方
向进行，最终生成B1结构，完成环氧氯丙烷的开环
.
图4路径A和B的相对吉布斯自由能图
然后，B1中插入CO2分子，经过能垒为79.07kJ/mol的过渡态TS2B得到B2.TS2B中，被H1活化的O1逐渐远离H1，O1-H1伸长至0.1371nm，呈断键趋势.同时，O1与CO2分子的C3逐渐靠近，距离缩短至0.1805nm，CO2的空间结构也发生一定的变化，变为夹角为147.46°的非直线型，TS2B有且只有一个虚频(-246.11icm–1)，振动沿着O1-H1键拉长断键、O1和C3距离缩短成键的方向进行，分别指向TS2B前后的反应物和产物.经过过渡态TS2B结构生成B2，其相对能量为38.82kJ/mol，完成CO2的捕捉，同时Br-C2键长增加，作用减弱.需要指出的是，如果不考虑零点能，则过渡态TS2B比中间体B2的能量高，考虑零点能后，B2的能量比TS2B 略高，这是个较常见的问题[14-15].
第三步，中间体B2中的C2与CO2分子中的O2相互作用，经过过渡态TS3B，越过71.07kJ/mol的能垒，闭环形成环状碳酸酯，放热53.29kJ/mol.TS3B中，Br-C2键伸长至0.2522nm，呈断键趋势，而C2与O2距离靠近缩短至0.1896nm，呈成键趋势.过渡态TS3B有唯一虚频(-434.83icm–1)，振动沿着旧键Br-C2断裂、新键C2-O2生成的方式进行，最终闭环完成反应，同时，离子液体催化剂[DMAPH]Br复原.
由上述分析和图4可知，路径A中，所有过渡态中TS3A相对能量最高，为133.13kJ/mol，即路径A的开环、CO2插入、闭环3个步骤中，最后一步闭环生成产物所要跨越的能垒最高，是反应的决速步，而总的能垒为180.12kJ/mol.路径B中，过渡态TS1B相对能量最高，为115.54kJ/mol，即路径B开环、CO2插入、闭环3个反应步骤中，第一步开环生成B1所要跨越的能垒最高，是路径B的决速步，反应的总能垒为169.90kJ/mol.由总能垒的大小可知，路径A的总反应能垒要比路径B总反应能垒高，因此反应路径B为更优反应路径.
第5期宋振艳等：[DMAPH]Br催化CO2和环氧氯丙烷反应的机理
29 2.3前线分子轨道理论分析
图5给出了环氧氯丙烷、CO2、环氧氯丙烷和催化剂配合物A0和B0的前线轨道HOMO和LUMO以及相应的能量.从图5中可看出，电子从环氧氯丙烷的HOMO轨道跃迁到CO2的LUMO轨道所要跨越的能垒为818.29kJ/mol，而电子从CO2的HOMO轨道跃迁到LUMO轨道所要跨越的能垒相对较高，为980.68kJ/mol，
.
但两者都很大，所以很难自发反应，需要催化体系来活化Array
图5前线分子轨道的能量(kJ/mol)
在催化剂的作用下，催化剂和环氧氯丙烷先形成配合物A0和B0，两者的LUMO和HOMO的能量相差不大，但与CO2相比则相差很大.根据前线分子轨道理论，A0和B0的LUMO与CO2的HOMO能量相差较大，分别为843.62和842.36kJ/mol，而A0和B0的HOMO与CO2的LUMO能量相差较小，分别为527.95和527.05kJ/mol，可见反应由A0和B0提供电子，CO2提供空轨道，且B0的HOMO与CO2的LUMO能量差与A0相比略低，从另一角度证实了B路径更优.
3结论
在密度泛函理论水平B3LYP/6-31+G(d，p)下，计算了离子液体[DMAPH]Br催化CO2和环氧氯丙烷反应生成环碳酸酯的反应机理.设计了2个反应路径A和B，对不同路径所涉及的反应物、中间体、过渡态以及产物的相关结构参数和能量进行了理论计算.根据计算的结果，B路径的总能垒较低，所以路径B为更优反应路径.前线轨道理论分析表明：反应发生时环氧氯丙烷的HOMO轨道与CO2的LUMO轨道发生对称性重叠，电子从环氧氯丙烷的HOMO轨道流向CO2的LUMO轨道从而发生化学键的生成和断裂.
参考文献：
[1]潘星.锌系催化剂催化二氧化碳/环氧化合物共聚机理研究[D].上海:华东理工大学,2013.
[2]王勤.若干金属配合物活化二氧化碳机理的理论研究[D].山西:山西师范大学,2015.
[3]ZHANG X,SU D,XIAO L F,et al.Immobilized protic ionic liquids:Efficient catalysts for CO2fixation with epoxides[J].J.
CO2Utilization，2017,17:37-42.
[4]JUTZ F,ANDANSON J M,BAIKER A.Ionic Liquids and Dense Carbon Dioxide:A Beneficial Biphasic System for Catalysis
肇庆学院学报第40卷30
[J].Chem.Rev.2011,111:322-353.
[5]YANG Z Z,ZHAO Y N,HE L N.CO2chemistry:task-specific ionic liquids for CO2capture/activation and subsequent conver-
sion[J].RSC Advances,2011,1:545-567.
[6]YUE S,HAO X J,WANG P P,et al.Amino acid-based ionic liquids for CO2conversion to form cyclic carbonate under sol-
vent-free conditions[J].Molecular Catalysis，2017,433:420-429.
[7]LIU M S,LI X,LIANG L,et al.Protonated triethanolamine as multi-hydrogen bond donors catalyst for efficient cycloaddi-
tion of CO2to epoxides under mild and cocatalyst-free conditions[J].J.CO2Utilization，2016,16:384-390.
[8]KIM H,BABU R,ROSHAN R,et al.Catalytic performance of metal azolate frameworks in the solventless synthesis of cy-
clic carbonates from CO2and epoxides[J].Applied Catalysis A:General，2017,538:59-65.
[9]LIU D,LI G,LIU H O.Functionalized MIL-101with imidazolium-based ionic liquids for the cycloaddition of CO2and epox-
ides under mild condition[J].Applied Surface Science，2018,428:218-225.
[10]WANG L,HUANG T F,CHEN C,et al.Mechanism of hexaalkylguanidinium salt/zinc bromide binary catalysts for the fixa-
tion of CO2with epoxide:A DFT investigation[J].J.CO2Utilization，2016,14:61-66.
[11]CHEN C,MA Y,ZHENG D N,et al.Insight into the role of weak interaction played in the fixation of CO2catalyzed by the
amino-functionalized imidazolium-based ionic liquids[J].J.CO2Utilization，2017,18:156-163.
[12]ZHANG Z G,FAN F J,XING H B,et al.Efficient Synthesis of Cyclic Carbonates from Atmospheric CO2Using a Positive
Charge Delocalized Ionic Liquid Catalyst[J].ACS Sustainable Chem.Eng.2017,5:2841-2846.
[13]FRISCH M J,TRUCKS G W,SCHLEGEL H B,et al.Gaussian09,Revision E.01,Gaussian Inc.,Wallingford CT,2013.
[14]FEYEL S,DÖBLER J,HOECKENDORF R,et al.Activation of methane by oligomeric[(Al2O3)x+(x=3,4,5)]:The role of Ox-
ygen-Centered radicals in thermal Hydrogen-Atom abstraction[J].Angew.Chem.Int.Ed.2008,47:1946-1950.
[15]ZHAO Y X,WU X N,MA,J B,et al.Experimental and theoretical study of the reactions between Vanadium-Silicon hetero-
nuclear oxide cluster anions with N-Butane[J].J.Phys.Chem.C，2010,114:12271-12279.
Mechanism for the Cycloaddition of the CO2with Epichlorohydrin
Catalyzed by[DMAPH]Br
SONG Zhenyan1，FENG Huajie2，LI Zhiwei1
(1,School of Environmental and Chemical Engineering,Zhaoqing University,Zhaoqing,Guangdong526061；
2,School of Chemistry and Chemical Engineering,Hainan Normal University,Haikou,Hainan571158,China) Abstract:The cycloaddition mechanism between carbon dioxide and epichlorohydrin catalyzed by4-dimeth-ylaminopyridine hydrobromide([DMAPH]Br)was studied by density functional theory and frontier molecular orbital theory.Two possible reaction pathways were designed and calculated,and the optimal reaction pathway was obtained.
Keywords:ionic liquid;carbon dioxide;epichlorohydrin;density functional theory
（责任编辑：张宝杰）。