電子效應(yīng)調(diào)控丁烯二腈光分子開(kāi)關(guān)的CASSCF和MS-CASPT2理論研究

曹　丹, 李苑瑩, 蘇青青, 王　斌, 劉峰毅, 王文亮

(陜西省大分子科學(xué)重點(diǎn)實(shí)驗(yàn)室, 陜西師范大學(xué)化學(xué)化工學(xué)院, 西安 710062)

?

電子效應(yīng)調(diào)控丁烯二腈光分子開(kāi)關(guān)的CASSCF和MS-CASPT2理論研究

曹丹, 李苑瑩, 蘇青青, 王斌, 劉峰毅, 王文亮

(陜西省大分子科學(xué)重點(diǎn)實(shí)驗(yàn)室, 陜西師范大學(xué)化學(xué)化工學(xué)院, 西安 710062)

摘要采用多組態(tài)CASSCF方法和MS-CASPT2方法研究了丁烯二腈中性分子及陽(yáng)離子和陰離子的順-反異構(gòu)化機(jī)理. 結(jié)果表明, 中性分子和離子態(tài)的光順-反異構(gòu)化反應(yīng)經(jīng)歷不同的非絕熱躍遷途徑: 中性丁烯二腈受光激發(fā)至\%S\%1態(tài)后, 需克服一個(gè)不低于19.7 kJ/mol的能壘才有機(jī)會(huì)到達(dá)基態(tài)和激發(fā)態(tài)間的圓錐交叉(\%S\%1/\%S\%0-CI), 隨后經(jīng)非輻射躍遷回到基態(tài)， \%S\%1/\%S\%0-CI在結(jié)構(gòu)上偏離CC雙鍵旋轉(zhuǎn)路徑, 且能量較高, 因此會(huì)降低旋轉(zhuǎn)速度, 阻礙旋轉(zhuǎn)的單向性；丁烯二腈陽(yáng)離子和陰離子自由基的\%D\%0態(tài)和\%D\%1態(tài)旋轉(zhuǎn)勢(shì)能面在90°處相交, 優(yōu)化的\%D\%1/\%D\%0-CI與\%D\%1態(tài)中間體的結(jié)構(gòu)和能量均相近, 因此從\%D\%1態(tài)經(jīng)由\%D\%1/\%D\%0-CI無(wú)輻射躍遷到\%D\%0態(tài)的過(guò)程無(wú)勢(shì)壘, 在此過(guò)程中CC旋轉(zhuǎn)方向性得到最大限度的保持. 研究結(jié)果證實(shí)了電子誘導(dǎo)不僅能降低基態(tài)熱旋轉(zhuǎn)勢(shì)壘, 而且能夠調(diào)控光旋轉(zhuǎn)的非絕熱躍遷機(jī)理.

關(guān)鍵詞順-反異構(gòu); 圓錐交叉; 丁烯二腈; 分子開(kāi)關(guān)

1計(jì)算方法

在CASSCF計(jì)算中, 選擇分子的全部10個(gè)π型分子軌道作為活化軌道(圖S1, 見(jiàn)本文支持信息). 其中, 中性丁烯二腈分子的活化電子數(shù)為10個(gè), 活化空間CASSCF(10e, 10o); 而在其陽(yáng)離子和陰離子中分別減少和增加一個(gè)(活化)電子, 產(chǎn)生的活化空間為CASSCF(9e, 10o)和CASSCF(11e, 10o). 為了進(jìn)一步考慮動(dòng)態(tài)電子相關(guān)效應(yīng), 在CASSCF優(yōu)化的分子構(gòu)型和波函數(shù)基礎(chǔ)上采用多態(tài)二級(jí)微擾理論(MS-CASPT2)進(jìn)行單點(diǎn)能校正(部分結(jié)構(gòu)優(yōu)化), 以獲得更精確的能量(和結(jié)構(gòu))信息[22,23]. 所有計(jì)算均采用Molcas@UU 8.0程序[24,25]完成.

2結(jié)果與討論

2.1丁烯二腈順-反異構(gòu)體的基態(tài)與激發(fā)態(tài)電子結(jié)構(gòu)

圖3還示出了中性分子、 陽(yáng)離子和陰離子態(tài)丁烯二腈的前線(xiàn)軌道能量隨基態(tài)順-反異構(gòu)化進(jìn)程的變化曲線(xiàn). 在整個(gè)旋轉(zhuǎn)路徑上, 軌道能的順序?yàn)殛庪x子>中性分子>陽(yáng)離子, 說(shuō)明電子占據(jù)會(huì)導(dǎo)致成鍵π和反鍵π*軌道能量的整體上升. 同時(shí), 在3個(gè)體系中, MO 20(π)和MO 21(π*)的能量差均在θ=90°處有極小值, 一方面反映出在該區(qū)域附近體系波函數(shù)具有明顯的多態(tài)特征; 另一方面, 也預(yù)示著基態(tài)和激發(fā)態(tài)在該處可能存在交叉或避免交叉.

2.2丁烯二腈中性分子的光致順-反異構(gòu)化反應(yīng)與非絕熱躍遷機(jī)理

圖4(A)～(C)分別為采用CASSCF，MS-CASPT2//CASSCF和MS-CASPT2方法計(jì)算得到的中性丁烯二腈的\%S\%0態(tài)和\%S\%1態(tài)的絕熱旋轉(zhuǎn)勢(shì)能曲線(xiàn)(曲線(xiàn)a和c). 以反式結(jié)構(gòu)的能量為零點(diǎn), 采用MS-CASPT2//CASSCF和MS-CASPT2方法計(jì)算的順式結(jié)構(gòu)的相對(duì)能量分別為3.4和3.3 kJ/mol, 具有典型的雙穩(wěn)態(tài)特征; 同時(shí)，從勢(shì)能曲線(xiàn)上獲取的從反式結(jié)構(gòu)出發(fā)的能壘高達(dá)232.2和225.2 kJ/mol,較全優(yōu)化得到的過(guò)渡態(tài)能壘略高(二者在MS-CASPT2//CASSCF水平的能量差為6.2 kJ/mol), 說(shuō)明中性丁烯二腈的熱順-反異構(gòu)化較難發(fā)生. 分析圖4的激發(fā)態(tài)勢(shì)能曲線(xiàn)發(fā)現(xiàn), CASSCF和MS-CASPT2優(yōu)化的\%S\%1-MEP(曲線(xiàn)c)均呈“U”形, 因此在θ=90°處可能存在\%S\%1態(tài)中間體(\%S\%1-min); 但由于在該處\%S\%1-MEP與其在\%S\%0態(tài)的投影(曲線(xiàn)b)在能量上很接近, 無(wú)法優(yōu)化\%S\%1-min. 此外, 比較MS-CASPT2//CASSCF和MS-CASPT2的計(jì)算結(jié)果還可發(fā)現(xiàn), 二者在定量上相吻合, 因此在后續(xù)計(jì)算中, 主要采用效率較高的MS-CASPT2//CASSCF方法.

2.3丁烯二腈陽(yáng)離子和陰離子的光異構(gòu)化反應(yīng)與非絕熱躍遷機(jī)理

圖6給出了丁烯二腈陽(yáng)離子和陰離子的基態(tài)和激發(fā)態(tài)絕熱旋轉(zhuǎn)路徑, 其中CASSCF水平優(yōu)化得到的陽(yáng)離子和陰離子體系的\%S\%0-MEP和\%S\%1-MEP分別見(jiàn)圖6(A)和(B), 對(duì)應(yīng)的MS-CASPT2校正后的能量曲線(xiàn)分別見(jiàn)圖6(C)和(D).

從以上分析可見(jiàn), 在丁烯二腈的成鍵π軌道上電離一個(gè)電子, 或在反鍵π*軌道加入一個(gè)電子, 形成的陽(yáng)離子和陰離子自由基不僅熱異構(gòu)化勢(shì)壘大大降低, 而且更為重要的是, 光異構(gòu)化非絕熱反應(yīng)途徑由有勢(shì)壘的、 多反應(yīng)模式參與的旋轉(zhuǎn)變?yōu)闊o(wú)勢(shì)壘、 單坐標(biāo)的旋轉(zhuǎn)過(guò)程.

3結(jié)論

支持信息見(jiàn)http://www.cjcu.jlu.edu.cn/CN/10.7503/cjcu20150966.

參考文獻(xiàn)

[1]Feringa B. L.,J.Org.Chem., 2007, 72(18), 6635—6652

[2]Kazaryan A., Kistemaker J. C. M., Schafer L. V., Browne W. R., Feringa B. L., Filatov M.,J.Phys.Chem.A, 2010, 114(15), 5058—5067

[3]Cnossen A., Kistemaker J. C. M., Kojima T., Feringa B. L.,J.Org.Chem., 2014, 79(3), 927—935

[4]Lu X. Y., Zhu W. H.,Chin.J.Org.Chem., 2007, 27(11), 1352—1357(陸欣宇, 朱為宏. 有機(jī)化學(xué), 2007, 27(11), 1352—1357)

[5]Fang C., Oruganti B., Durbeej B.,RSCAdvances, 2014, 4(20), 10240—10251

[6]Filatov M.,Chem.Phys.Chem., 2011, 12(17), 3348—3353

[7]Guo X., Zhou J. W., Siegler M. A., Bragg A. E., Katz H. E.,Angew.Chem.Int.Ed., 2015, 54(16), 4782—4786

[8]Freund L., Klessinger M.,Int.J.QuantumChem., 1998, 70, 1023—1028

[9]Barbatti M., Paier J., Lischka H.,J.Chem.Phys., 2004, 121(23), 11614—11624

[10]Barbatti M., Granucci G., Persico M., Lischka H.,Chem.Phys.Lett., 2005, 401(S1—S3), 276—281

[11]Barbatti M., Ruckenbauer M., Lischka H.,J.Chem.Phys., 2005, 122(17), 174307

[12]Ben-Nun M., Martínez T. J.,Chem.Phys.Lett., 1998, 298(1), 57—65

[13]Michael F., Massimo O.,J.Org.Chem., 2014, 79(8), 3587—3600

[14]Klaiman S., Cederbaum L. S.,Angew.Chem.Int.Ed., 2015, 127(36), 10616—10619

[15]Finley J., Malmqvist P. ?., Roos B. O., Serrano-Andrés L.,Chem.Phys.Lett., 1998, 288(S2—4), 299—306

[16]Andersson K., Malmqvist P. ?., Roos B. O.,J.Chem.Phys., 1992, 96(2), 1218—1226

[17]Malmqvist P. ?., Roos B. O., Schimmelpfennig B.,Chem.Phys.Lett., 2002, 357(3), 230—240

[18]Andersson K.,Theor.Chim.Acta, 1995, 91(1), 31—46

[19]Levine B. G., Ko C., Quenneville J., Martínez T. J.,Mol.Phys., 2006, 104(5), 1039—1051

[20]Levine B. G., Coe J. D., Martínez T. J.,J.Phys.Chem.B, 2008, 112(2), 405—413

[21]Maeda S., Ohno K., Morokuma K.,J.Chem.TheoryComput., 2010, 6(5), 1538—1545

[22]Qu Z. X., Gao J. L.,Chem.J.ChineseUniversities, 2015, 36(11), 2236—2240(曲澤星, 高加力. 高等學(xué)校化學(xué)學(xué)報(bào), 2015, 36(11), 2236—2240)

[23]Song Y. L., Liu Y. J.,Chem.J.ChineseUniversities, 2015, 36(11), 2163—2170(宋艷麗, 劉亞軍. 高等學(xué)?；瘜W(xué)學(xué)報(bào), 2015, 36(11), 2163—2170)

[24]Aquilante F., de Vico L., Ferré N., Ghigo G., Malmqvist P. ?., Neogrády P., Pedersen T. B., Pitoňák M., Reiher M., Roos B. O., Serrano-Andrés L., Urban M., Veryazov V., Lindh R.,J.Comput.Chem., 2010, 31(1), 224—247

[25]Aquilante F., Autschbach J., Carlson R. K., Chibotaru L. F., Delcey M. G., de Vico L., Galván I. F., Ferré N., Frutos L. M., Gagliardi L., Garavelli M., Giussani A., Hoyer C. E., Manni G. L., Lischka H., Ma D., Malmavist P. ?., Müller T., Nenov A., Olivucci M., Pedersen T. B., Peng D., Plasser F., Pritchard B., Reiher M., Rivalta I., Schapiro I., Segarra-Martí J., Stenrup M., Truhlar D. G., Ungur L., Valentini A., Vancoillie S., Veryazov V., Vysotskiy V. P., Weingart O., Zapata F., Lindh R.,J.Comput.Chem., 2016, 37, 506—541

[26]Liu F. Y., Morokuma K.,J.Am.Chem.Soc., 2012, 134(10), 4864—4876

[27]Ohmine L.,J.Chem.Phys., 1985, 83(5), 2348—2362

[28]Ben-Nun M., Martínez T. J.,Chem.Phys., 2000, 259(2), 237—248

[29]Krawczyk R. P., Viel A., Manthe U., Domcke W.,J.Chem.Phys., 2003, 119(3), 1397—1411

[30]Yamazaki S., Kato S.,J.Chem.Phys., 2005, 123(11), 114510

(Ed.: Y, Z, S)

? Supported by the National Natural Science Foundation of China(Nos.21473107, 21473108), the Natural Science Basic Research Plan in Shaanxi Province of China(No.2015JM2056) and the Fundamental Research Funds for the Central Universities, China(No.GK201502002).

CASSCF and MS-CASPT2 Studies on an Electron-tunable,1,2-Dicyanoethylene-based Optical Molecular Switch?

CAO Dan, LI Yuanying, SU Qingqing, WANG Bin, LIU Fengyi*, WANG Wenliang

(KeyLaboratoryforMacromolecularScienceofShaanxiProvince,SchoolofChemistry&ChemicalEngineering,ShaanxiNormalUniversity,Xi’an710062,China)

KeywordsCis-transisomerization; Conical intersection; 1,2-Dicyanoethylene; Molecular switch

AbstractCASSCF and MS-CASPT2 mechanistic studies were carried on thecis-transisomerization processes of 1,2-dicyanoethylene in its neutral, cationic and anionic forms. The results confirmed the importance of electron-induction in reducing the reaction barrier, and more importantly, revealed the different nonadiabatic channels tuned by electron attachment/detachment. The \%S\%1→\%S\%0decay of neutral dicyanoethylene need overcome a mild barrier(≥19.7 kJ/mol) to reach a H-migration-type, namely the \%S\%1/\%S\%0-CI, which is away from the CC torsional coordinates and may slow down the speed of CC rotation and hurt its directionality; while in cationic and anionic isomerization processes, the \%D\%1and \%D\%0PESs intersect along the rotary path, therefore, the nonadiabatic \%D\%1→\%D\%0decay is barrierless, as result, the directionality of CC rotation is maintained. The study revealed the role of electron induction in tuning thecis-transphotoisomerization, and shed light on the design of light-driven molecular rotary motors.

收稿日期:2015-12-20. 網(wǎng)絡(luò)出版日期: 2016-05-26.

基金項(xiàng)目:國(guó)家自然科學(xué)基金(批準(zhǔn)號(hào): 21473107， 21473108)、 陜西省自然科學(xué)基礎(chǔ)研究計(jì)劃項(xiàng)目(批準(zhǔn)號(hào): 2015JM2056)和中央高?；究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金(批準(zhǔn)號(hào): GK201502002)資助.

中圖分類(lèi)號(hào)O641

文獻(xiàn)標(biāo)志碼A

聯(lián)系人簡(jiǎn)介: 劉峰毅, 男, 博士, 教授, 博士生導(dǎo)師, 主要從事理論與計(jì)算化學(xué)研究. E-mail: FengyiLiu@snnu.edu.cn