包含平面五配位和平面四配位碳化合物的結(jié)構(gòu)和穩(wěn)定性

吳林峰 李飛飛 張聰杰

(陜西師范大學(xué)化學(xué)與化工學(xué)院,陜西省大分子科學(xué)重點實驗室,西安710062)

包含平面五配位和平面四配位碳化合物的結(jié)構(gòu)和穩(wěn)定性

吳林峰 李飛飛 張聰杰*

(陜西師范大學(xué)化學(xué)與化工學(xué)院,陜西省大分子科學(xué)重點實驗室,西安710062)

在B3LYP/6-311+G＊＊水平上,研究了包含平面五配位碳(ppC)和平面四配位碳(ptC)的化合物CB5C2H2(C3B2)nC2H2CB5(n=1-5)的結(jié)構(gòu)、能隙、IR光譜、電子光譜、Wiberg鍵指數(shù)以及芳香性.計算結(jié)果表明,這五種化合物的能量最低結(jié)構(gòu)均位于勢能面的極小值點,HOMO-LUMO能隙在0.5到1.2 eV之間,第一電子的激發(fā)能在1780到2910 nm之間,且與分子尺寸呈現(xiàn)非單調(diào)變化趨勢.Wiberg鍵指數(shù)和鍵長表明這五種化合物的能量最低結(jié)構(gòu)均包含平面五配位和平面四配位碳.C3B2單元和右側(cè)CB5單元中三元環(huán)中心的核獨立化學(xué)位移(NICS(0))值均為負值,而左側(cè)CB5單元中的三元環(huán)中心的NICS(0)值中兩個為正值,另兩個為負值,NICS(1)值與NICS(0)值也一致,因此電子的局域離域?qū)Y(jié)構(gòu)的穩(wěn)定是很重要的.

密度泛函理論;平面五配位碳原子;平面四配位碳原子;芳香性

1 引言

1970年,Hoffmann等1通過雜化軌道分析了平面甲烷結(jié)構(gòu)的分子的化學(xué)鍵,認(rèn)為包含平面四配位碳原子(ptC)的化合物可能存在.40年來,在理論和實驗上,科研工作者對包含平面四配位化合物的結(jié)構(gòu)和穩(wěn)定性進行了廣泛的研究.2-28Keese13和Merino12等詳細闡述了ptC領(lǐng)域在實驗上和理論上的一些重要研究成果.在理論上得到的分子主要包含兩類,一類分子是ptC通過化學(xué)鍵被束縛在籠中.2-5另一類是ptC置于平面分子中,通過離域效應(yīng)得到穩(wěn)定構(gòu)型.6-28在實驗和理論相結(jié)合的基礎(chǔ)上,證實了存在包含ptC的Al4C-和MAl4(M=Si,Ge)團簇.6,7,22一些理論研究表明很多包含ptC的共軛體系滿足4n+2規(guī)則.23,24我們26發(fā)現(xiàn)在包含ptC的共軛分子中均存在共軛三元環(huán),所以我們系統(tǒng)研究了包含ptC共軛分子的穩(wěn)定化規(guī)則,而且得到了很好的驗證,并擴展到包含平面四配位氮原子體系.

在所有包含ptC的分子中,我們發(fā)現(xiàn)C3B2H4非常有趣,因為C3B2H4分子與實驗上得到了螺旋結(jié)構(gòu)的C5H4化合物的分子式非常接近,所以在理論上對C3B2H4分子的穩(wěn)定結(jié)構(gòu)和性質(zhì)進行了系統(tǒng)的研究,計算結(jié)果表明包含ptC的C3B2H4分子結(jié)構(gòu)比螺旋結(jié)構(gòu)穩(wěn)定,以C3B2H4分子為基礎(chǔ),我們研究了C3B2H4衍生物的結(jié)構(gòu)和性質(zhì),27以及C3B2H4的聚合行為,得到一系列包含ptC的化合物.28在包含ptC的研究基礎(chǔ)上,科研工作者還對包含平面五配位碳(ppC)和平面六配位碳(phC)化合物的結(jié)構(gòu)和穩(wěn)定性也進行了研究,29-31發(fā)現(xiàn)包含一個ppC的CB5H2分子以及包含兩個ppC的C2B10團簇是穩(wěn)定的.30,31

由于共軛聚合物端位結(jié)構(gòu)的不同,將導(dǎo)致不同的共軛效應(yīng)和光譜性質(zhì).32所以本文將C3B2H4的聚合物兩端采用CB5H2結(jié)構(gòu)封閉,得到CB5C2H2(C3B2)nC2H2CB5(n=1-5)五種化合物,在理論上研究其穩(wěn)定結(jié)構(gòu)和電子性質(zhì).

2 模型設(shè)計和計算方法

首先設(shè)計了(C3B2)n兩端用CB5單元封閉,由于張力原因,在(C3B2)n與CB5之間采用用―CH―基連接,得到了包含平面五配位和平面四配位碳化合物CB5C2H2(C3B2)nC2H2CB5(n=1-5).在B3LYP/6-311+ G**水平上,我們優(yōu)化了CB5C2H2(C3B2)nC2H2CB5(n= 1-5)化合物結(jié)構(gòu),計算了它們的振動頻率.采用規(guī)范不變原子軌道(GIAO)方法計算了CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物能量最低結(jié)構(gòu)中三元環(huán)、六元環(huán)和八元環(huán)平面上兩個位置的核獨立化學(xué)位移(NICS),一個是平面中心(NICS(0)),一個是在中心上方0.1 nm處(NICS(1)).運用自然鍵軌道方法(NBO)預(yù)測了他們的Wiberg鍵指數(shù).33采用含時密度泛函理論TD-B3LYP和6-311+G**基組,我們計算CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物最低能量結(jié)構(gòu)的垂直激發(fā)能和振子強度.所有計算均采用Gaussian 03程序包34完成.

3 結(jié)果與討論

3.1 CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物的

結(jié)構(gòu)和穩(wěn)定性

在B3LYP/6-311+G**水平上優(yōu)化包含平面五配位和平面四配位碳化合物CB5C2H2(C3B2)nC2H2CB5(n=1-5)的結(jié)構(gòu),優(yōu)化的幾何結(jié)構(gòu)如圖1中所示.因為CB5C2H2(C3B2)2C2H2CB5中C3B2的取向不同,存在3個異構(gòu)體,如圖1中A2-a,A2-b和A2-c所示.同理, CB5C2H2(C3B2)3C2H2CB5存在4個異構(gòu)體(圖1中A3-a,A3-b,A3-c和A3-d).優(yōu)化這7個結(jié)構(gòu)及其頻率計算得到其電子態(tài)、能量、相對能、HOMO-LUMO能隙和最小振動頻率列于表1中.由表1可知,A2-b和A2-c的能量分別比A2-a高69.75和27.33 kJ· mol-1,A3-b、A3-c和A3-d的能量分別比A3-a高69.58、42.28和125.58 kJ·mol-1,表明C3B2單元在同方向的成鍵結(jié)構(gòu)是更有利的,與我們以前研究(C3B2)nH4的聚合物的結(jié)構(gòu)一致.27所以,當(dāng)n=4,5時,我們只研究了CB5C2H2(C3B2)nC2H2CB5化合物的C3B2單元的取向是相同方向的結(jié)構(gòu),且將其電子態(tài)、能量、能隙及最小振動頻率也列于表1中.由表1可知,A1、A2-a、A3-a、A4和A5的最小振動頻率都為正值,均是勢能面上的極小點.

圖1所示A1、A2-a、A3-a、A4和A5都是C2v點群.A1結(jié)構(gòu)中含有兩個ppCs和一個ptC.B―C、ptC—C和ptC—B的鍵長分別是0.1488、0.1399和0.1659 nm,,其中B—C和ptC—B鍵分別比C3B2H4中相應(yīng)的鍵增長0.0085和0.0044 nm,但ptC—C幾乎不變.28左邊CB5單元中ppC—B的鍵長介于0.1571和0.1665 nm之間,左邊CB5單元中ppC—B的鍵長相對右邊CB5單元中對應(yīng)鍵大約長0.002 nm,而且比C2B10團簇中相應(yīng)的鍵長也長.31但是右邊CB5單元中ppC—B鍵長與C2B10團簇中相應(yīng)的鍵長非常接近.這主要是因為—CH—C—導(dǎo)致左邊CB5上的p電子增加,使得左邊CB5單元中ppC—B鍵長增加,雖然右邊CB5單元也與—CH—相連,但由于—CH—左邊是B原子,所以左邊—CH—基對左邊CB5單元影響大.

圖1 包含ptCs和ppCs的CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物的結(jié)構(gòu)Fig.1 Optimized geometries of CB5C2H2(C3B2)nC2H2CB5(n=1-5)compounds with ptCs and ppCs bond length in nm

在A2-a、A3-a、A4和A5結(jié)構(gòu)中,ppC—B鍵長在0.1562-0.1674 nm之間,ptC—C和ptC—B鍵長分別在0.1422-0.1448 nm和0.1607-0.1668 nm之間. ptC—C和ptC—B鍵長較A1結(jié)構(gòu)中的ptC—C和ptC—B變化較大,這是因為隨著n的增大,中間部分將形成非常強的局部離域效應(yīng).27由于A2-a、A3-a、A4和A5結(jié)構(gòu)的CB5與C3B2間存在—CH—基,所以這四個結(jié)構(gòu)中兩端CB5部分的離域效應(yīng)主要由(CH)2CB5形成,兩邊對應(yīng)的ppC—B鍵長很接近,比如,A2-a和A3-a中左邊ppC—B鍵長最大差值為0.008 nm,右邊對應(yīng)的ppC—B鍵長最大差值為0.007 nm.

在B3LYP/6-311+G**水平上,計算了能量最低結(jié)構(gòu)A1,A2-a,A3-a,A4和A5的Wiberg鍵指數(shù)(WBI),列于表2中.從表2中可以看出,這5個結(jié)構(gòu)中ppC—B鍵的WBI介于0.6579-0.9793之間, ptC—B鍵的WBI位于0.5440-0.7114之間,表明這5種結(jié)構(gòu)中均存在ppC和ptC原子.而且ppC和ptC的總WBI接近于4,滿足八隅規(guī)則.

表1 含ptCs和ppCs的CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物的能量最低結(jié)構(gòu)的對稱性、電子態(tài)、總能量(E)、相對能(Er)、HOMO和LUMO能量(EHOMO和ELUMO),HOMO-LUMO能隙(Egap)和最小振動頻率(Freq)Table 1 Symmetries,electronic states,total energies(E),relative energies(Er),HOMO and LUMO energies(EHOMOand ELUMO),HOMO-LUMO energy gaps(Egap),and the minimum vibrational frequencies(Freq)of the lowest energy structures of CB5C2H2(C3B2)nC2H2CB5(n=1-5)compounds with ptCs and ppCs

3.2 CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物能量最低結(jié)構(gòu)的芳香性

圖2給出CB5C2H2(C3B2)nC2H2CB5(n=1-5)團簇的穩(wěn)定構(gòu)型的所有三元環(huán)、六元環(huán)和八元環(huán)環(huán)中心NICS(0)以及中心上方0.1 nm處位置的NICS(1)值.從圖2中可以看出,結(jié)構(gòu)中八元環(huán)平面中心的NICS(0)和NICS(1)值都是正值,表明都是反芳香性.六元環(huán)中心NICS(0)的值均為正值,除化合物A2-a為正值外,NICS(1)值都略為負值,表明A3-a、A4和A5的六元環(huán)具有比較弱的芳香性.左側(cè)CB5單元中的三元環(huán)中心的NICS(0)值兩個為正值(上下兩三元環(huán)),另兩個為負值(左邊兩三元環(huán)),而右側(cè)CB5單元和C3B2單元中三元環(huán)中心的NICS(0)值均為負值.兩端CB5單元中的三元環(huán)中的上下兩三元環(huán),隨著離環(huán)中心距離的增加而減小,而左邊兩三元環(huán),隨著離環(huán)中心距離的增加而增大.說明在CB5單元中上下兩三元環(huán)的芳香性越來越強,而左邊兩三元環(huán)的芳香性卻越來越弱.但對于整個結(jié)構(gòu)來說,隨著n的增加,C3B2單元中的三元環(huán)相同位置的NICS(0)值均逐漸減小,因此C3B2單元中三元環(huán)的芳香性逐漸增強,而左右側(cè)CB5單元中三元環(huán)的NICS(0)值隨著n增加,變化均很小,所以,所有三元環(huán)芳香性的變化與鍵長的變化趨勢是吻合的.兩端CB5三元環(huán)中心的NICS值的差異性主要是由于與這兩個單元相連的—CH—C—和—B—CH—單元的離域效應(yīng)的差異性導(dǎo)致的.三元環(huán)中NICS(1)值均大于NICS(0),除A1外,其它三元環(huán)中心NICS(1)值都為較小的負值.所以,三元環(huán)的芳香性表明局部離域?qū)ζ矫娼Y(jié)構(gòu)的穩(wěn)定是非常重要的,這與我們27以前的結(jié)論是一致的.

3.3 CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物能量最低結(jié)構(gòu)的IR光譜和電子光譜

表2 含ptCs和ppCs的CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物的能量最低結(jié)構(gòu)中的ptC―B和ppC―B鍵的Wiberg鍵指數(shù)(WBIC―B),ptCs和ppCs總的Wiberg鍵指數(shù)(WBIC)Table 2 WBIs of ptC―B and ppC―B bonds(WBIC―B)and total WBIs(WBIC)of the ptCs and ppCs of the lowest energy structures of CB5C2H2(C3B2)nC2H2CB5(n=1-5)with compounds ptCs and ppCs

圖2 包含ptCs和ppCs的CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物能量最低結(jié)構(gòu)中環(huán)中心的NICS(0)和NICS(1)值Fig.2 NICS(0)and NICS(1)values in the center of the rings of the lowest energy structures of CB5C2H2(C3B2)nC2H2CB5(n=1-5)compounds with ptCs and ppCsNICS(0):in plain,NICS(1):in italic

圖3給出了穩(wěn)定構(gòu)型的IR光譜圖,從圖3中可以看出,這五種構(gòu)型的強IR活性的振動大概位于800、1000以及1600 cm-1處,對應(yīng)于ppC和ptC與其周圍B和C的對稱和不對稱伸縮振動,其中C—ptC的對稱和不對稱伸縮振動具有較強的IR活性.最強的IR峰大概在1600 cm-1左右.

在TD-B3LYP/6-311+G**水平上,我們得到了CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物最低能量結(jié)構(gòu)的垂直激發(fā)能和振子強度,結(jié)果列于表3.從表3可以看出5個化合物的第一激發(fā)均來自X1A1→A1B2躍遷,激發(fā)能分別為:2907.87、1788.92、1782.69、2032.25和2317.37 nm,與分子尺寸呈現(xiàn)非單調(diào)遞增的變化趨勢.為了進一步探討這種現(xiàn)象,圖4給出了表3中所示的躍遷軌道的圖形.由圖4可知,A1和 A2-a分子中電子從基態(tài)到第一激發(fā)態(tài)都是HOMO到LUMO,但是A3-a、A4和A5從基態(tài)到第一激發(fā)態(tài)涉及的電子躍遷都來自HOMO-1到LUMO,而且A3-a、A4和A5的HOMO-1和LUMO圖形非常相似,所以A3-a、A4和A5的第一激發(fā)能隨著分子尺寸的增大而逐漸減小.

圖3 包含ptCs和ppCs的CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物能量最低結(jié)構(gòu)的紅外光譜Fig.3 IR spectra of the lowest energy structures of CB5C2H2(C3B2)nC2H2CB5(n=1-5)compounds with ptCs and ppCs

表3 含ptCs和ppCs的CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物的能量最低結(jié)構(gòu)的第一激發(fā)能(λ)和振子強度(f)Table 3 Vertical electron transition energies(λ)and oscillator strengths(f)of the lowest energy structures ofCB5C2H2(C3B2)nC2H2CB5(n=1-5)compounds with ptCs and ppCs

圖4 A1、A2-a、A3-a、A4和A5結(jié)構(gòu)的占據(jù)π分子軌道和LUMO圖Fig.4 Occupied π molecular orbitals and LUMO of the structuresA1,A2-a,A3-a,A4,andA5

4 結(jié)論

在B3LYP/6-311+G**水平上,得到了包含平面五配位和平面四配位碳化合物CB5C2H2(C3B2)nC2H2CB5(n=1-5)的能量最低結(jié)構(gòu)及其成鍵規(guī)律,結(jié)合鍵長和Wiberg鍵指數(shù)證實了這五種化合物的能量最低結(jié)構(gòu)均包含平面五配位和平面四配位碳.計算CB5C2H2(C3B2)nC2H2CB5(n=1-5)化合物的能量最低構(gòu)型的所有三元環(huán)、六元環(huán)和八元環(huán)中心位置的NICS(0)和NICS(1)值，表明負NICS值主要在三元環(huán)中心,因此,三元環(huán)的芳香性對平面結(jié)構(gòu)的穩(wěn)定是非常重要的.這五種化合物的能隙在0.5到1.2 eV之間,這幾個化合物的第一電子躍遷所產(chǎn)生的吸收光譜呈現(xiàn)非單調(diào)變化,第一激發(fā)能位于1780到2910 nm之間.

(1) Hoffmann,R.;Alder,R.W.;Wilcox,C.F.J.Am.Chem.Soc. 1970,92,4992.

(2) McGrath,M.P.;Radom,L.J.Am.Chem.Soc.1993,115,3320.

(3) Wang,Z.X.;Schleyer,P.v.R.J.Am.Chem.Soc.2001,123, 994.

(4) Wang,Z.X.;Schleyer,P.V.R.J.Am.Chem.Soc.2002,124, 11979.

(5) Rasmussen,D.R.;Radom,L.Angew.Chem.Int.Edit.1999,38, 2876.

(6) Li,X.;Wang,L.S.;Boldyrev,A.I.;Simons,J.J.Am.Chem. Soc.1999,121,6033.

(7) Wang,L.S.;Boldyrev,A.I.;Li,X.;Simons,J.J.Am.Chem. Soc.2000,122,7681.

(8) Liang,J.X.;Jia,W.H.;Zhang,C.J.;Cao,Z.X.Acta Phys.-Chim.Sin.2009,25,1847. [梁錦霞,賈文紅,張聰杰,曹澤星.物理化學(xué)學(xué)報,2009,25,1847.]

(9) Sun,W.X.;Zhang,C.J.Acta Phys.-Chim.Sin.2008,24,32. [孫文秀,張聰杰.物理化學(xué)學(xué)報,2008,24,32.]

(10) Merino,G.;Mendez-Rojas,M.A.;Vela,A.J.Am.Chem.Soc. 2003,125,6026.

(11) Merino,G.;Mendez-Rojas,M.A.;Beltran,H.I.;Corminboeuf, C.;Heine,T.;Vela,A.J.Am.Chem.Soc.2004,126,16160.

(12) Merino,G.;Mendez-Rojas,M.A.;Vela,A.;Heine,T. J.Comput.Chem.2007,28,362.

(13) Keese,R.Chem.Rev.2006,106,4748.

(14) Li,S.D.;Guo,J.C.;Miao,C.Q.;Ren,G.M.J.Phys.Chem.A 2005,109,4133.

(15)Li,S.D.;Guo,Q.L.;Miao,C.Q.;Ren,G.M.Acta Phys.-Chim.Sin.2007,23,743.[李思殿,郭巧凌,苗常青,任光明.物理化學(xué)學(xué)報,2007,23,743.]

(16) Roy,D.;Corminboeuf,C.;Wannere,C.S.;King,R.B.; Schleyer,P.v.R.Inorg.Chem.2006,45,8902.

(17) Zhang,C.J.;Sun,W.X.;Cao,Z.X.J.Am.Chem.Soc.2008, 130,5638.

(18) Wang,Z.X.;Zhang,C.G.;Chen,Z.;Schleyer,P.v.R.Inorg. Chem.2008,47,1332.

(19) Collins,J.B.;Dill,J.D.;Jemmis,E.D.;Apeloig,Y.;Schleyer, P.v.R.;Seeger,R.;Pople,J.A.J.Am.Chem.Soc.1976,98, 5419.

(20) Boldyrev,A.I.;Wang,L.S.J.Phys.Chem.A 2001,105,10759.

(21) Priyakumar,U.D.;Reddy,A.S.;Sastry,G.N.Tetrahedron Lett. 2004,45,2495.

(22) Pancharatna,P.D.;Mendez-Rojas,M.A.;Merino,G.;Vela,A.; Hoffmann,R.J.Am.Chem.Soc.2004,126,15309.

(23) Perez,N.;Heine,T.;Barthel,R.;Seifert,G.;Vela,A.; Mendez-Rojas,M.A.;Merino,G.Org.Lett.2005,7,1509.

(24) Esteves,P.M.;Ferreira,N.B.P.;Correa,R.J.J.Am.Chem. Soc.2005,127,8680.

(25) Perez-Peralta,N.;Sanchez,M.;Martin-Polo,J.;Islas,R.;Vela, A.;Merino,G.J.Org.Chem.2008,73,7037.

(26) Zhang,C.J.;Jia,W.H.;Cao,Z.X.J.Phys.Chem.A 2010,114, 7960.

(27) Zhang,C.J.;Wang,P.;Liang,J.X.;Jia,W.H.;Cao,Z.X. J.Mol.Struct.-Theochem 2010,941,41.

(28) Sun,W.X.;Zhang,C.J.;Cao,Z.X.J.Phys.Chem.C 2008, 112,351.

(29) Exner,K.;Schleyer,P.v.R.Science 2000,290,1937.

(30) Wang,Z.X.;Schleyer,P.v.R.Science 2001,292,2465.

(31) Liang,J.X.;Zhang,C.J.Acta Chim.Sin.2010,68,7. [梁錦霞,張聰杰.化學(xué)學(xué)報,2010,68,7.]

(32) Jamie,K.;Robert,M.;Michael,J.F.;Rik,R.T.Org.Lett.2008, 10,2163.

(33) Wiberg,K.B.Tetrahedron 1968,24,1083.

(34) Frisch,M.J.;Trucks,G.A.;Schlegel,H.B.;et al.Gaussian 03, Revision C.02;Gaussian Inc.:Wallingford,CT,2004.

December 21,2011;Revised:February 21,2012;Published on Web:March 7,2012.

Structure and Stability of Compounds with Planar Pentacoordinate Carbons and Planar Tetracoordinate Carbons

WU Lin-Feng LI Fei-Fei ZHANG Cong-Jie*
(Key Laboratory of Macromolecular Science of Shaanxi Province,School of Chemistry and Chemical Engineering, Shannxi Normal University,Xi?an 710062,P.R.China)

We have investigated the structures,gaps,IR spectra,electronic spectra,Wiberg bond indices (WBIs),and aromaticity of CB5C2H2(C3B2)nC2H2CB5(n=1-5)with planar pentacoordinate carbons(ppC)and planar tetracoordinate carbons(ptC)at the B3LYP/6-311+G＊＊level.Calculations indicate that the five compounds with the lowest energies are located at the minima of the potential energy surfaces.The energy gaps between the highest occupied molecular orbital(HOMO)and the lowest unoccupied molecular orbital (LUMO)vary between 0.5 and 1.2 eV;the first electronic transition wavelengths are between 1780 and 2910 nm and depend non-monotonically on the size of the compounds.WBIs of the five compounds show that they contain both ppC and ptC.The nucleus-independent chemical shift(NICS(0))values of the centers of the three-membered rings of the CB5sections on the right side of these compounds,as well as the C3B2sections,are negative,while the NICS(0)values of only two centers of the three-membered rings of the CB5sections on the left side are negative.In addition,since the NICS(0)of the centers of the three-membered rings are consistent with those of NICS(1),then local delocalization of the π electrons must play an important role in stabilizing these compounds.

Density functional theory;Planar pentacoordinate carbon atom;Planar tetracoordinate carbon atom;Aromaticity

10.3866/PKU.WHXB201203071

?Corresponding author.Email:zcjwh@snnu.edu.cn.Tel:+86-29-81530726.

The project was supported by the Fundamental Research Funds for the Central Universities,China(GK201002013)and State Key Laboratory for Physical Chemistry of Solid Surface(Xiamen University),China(2010).

2010年度中央高?；究蒲袠I(yè)務(wù)費專項資金(GK201002013)和廈門大學(xué)固體表面物理化學(xué)國家重點實驗室開放課題(2010)資助

O641