• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    第一性原理研究富勒烯衍生物PCBM的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜

    2015-09-03 07:46:31王廣偉孫紹濤宋秀能
    物理化學(xué)學(xué)報(bào) 2015年8期
    關(guān)鍵詞:富勒烯譜峰碳原子

    馬 勇 王廣偉 孫紹濤 宋秀能,*

    (1山東師范大學(xué)物理與電子科學(xué)學(xué)院,濟(jì)南 250014; 2山東信息通信技術(shù)研究院管理中心,濟(jì)南 250101)

    第一性原理研究富勒烯衍生物PCBM的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜

    馬 勇1王廣偉1孫紹濤2宋秀能1,*

    (1山東師范大學(xué)物理與電子科學(xué)學(xué)院,濟(jì)南 250014;2山東信息通信技術(shù)研究院管理中心,濟(jì)南 250101)

    富勒烯衍生物[6,6]-苯基-C61-丁酸甲酯(PCBM)在有機(jī)聚合物太陽(yáng)能電池的電子輸運(yùn)方面扮演著非常重要的角色.利用密度泛函理論計(jì)算了PCBM的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜及未占據(jù)分子軌道.通過(guò)對(duì)比計(jì)算得到的PCBM分子中不等價(jià)碳原子的譜線(xiàn)總和,將該分子的主要共振吸收峰進(jìn)行了標(biāo)定.我們分析了第一個(gè)π*共振吸收峰高能區(qū)右肩吸收峰的來(lái)源,并確定了該吸收峰主要來(lái)自于C60籠子中碳原子能量較高的未占據(jù)分子軌道的躍遷.

    密度泛函理論; X射線(xiàn); 精細(xì)結(jié)構(gòu); 富勒烯衍生物; 完全芯態(tài)空穴

    1 引 言

    球形富勒烯C60及其衍生物[6,6]-苯基-C61-丁酸甲酯(PCBM)通常被用來(lái)作為有機(jī)聚合物光伏器件的活性電子受體材料.1,2在小分子光伏器件中,太陽(yáng)能電池是通過(guò)在真空中沉積一層C60和電子受體材料而形成.在聚合物光伏器件中,活化層是由一層電子受體材料PCBM和電子給體型共軛聚合物的混合溶液覆蓋而成的薄膜.這些共軛分子的電子結(jié)構(gòu)在激子的解離和分離過(guò)程中具有重要作用.C60和PCBM分子的最低未占據(jù)分子軌道(LUMO)在電子轉(zhuǎn)移以及將電子從"給體-受體"交界面?zhèn)鬏數(shù)骄哂休^低功函數(shù)電極的過(guò)程中扮演著至關(guān)重要的角色.

    近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜是表面科學(xué)領(lǐng)域一種非常重要的技術(shù).該技術(shù)可以探測(cè)特定元素中電子從1s軌道到未占據(jù)分子軌道的躍遷,從而對(duì)有機(jī)分子表面材料的各種性質(zhì)進(jìn)行研究.3–5近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜目前已經(jīng)被廣泛應(yīng)用于光伏器件材料的研究,6–9以測(cè)定異質(zhì)結(jié)太陽(yáng)能電池活性層的化學(xué)組成.10,11共振吸收后,近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜對(duì)應(yīng)的末態(tài)要經(jīng)歷電子的弛豫過(guò)程,其態(tài)密度與該吸收對(duì)應(yīng)的未占據(jù)分子軌道在體系處于基態(tài)(GS)時(shí)的態(tài)密度是不一樣的.碳元素K殼層的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜在研究共軛多聚物12–14以及一些由碳元素形成的材料如富勒烯、富勒烯衍生物、15納米碳管16和石墨烯17等方面具有非常重要的意義.在C60和C70分子特性以及其分子間的弱相互作用被發(fā)現(xiàn)之后,二者截然不同的共振吸收就被通過(guò)近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜進(jìn)行了研究.Terminello等18利用Hückel理論較符合地模擬了它們的譜線(xiàn).近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜曾經(jīng)被廣泛地用以研究PCBM分子的電子結(jié)構(gòu)和幾何結(jié)構(gòu),6,12,19–21Germack等6的工作中因其較低的分辨率而無(wú)法解釋所有的π*共振,而近期的一些工作則研究了較為豐富的π*軌道的精細(xì)結(jié)構(gòu).12,19–21目前,在PCBM分子X(jué)射線(xiàn)共振吸收的研究方面有很多報(bào)道,Bazylewski等19對(duì)理論計(jì)算的PCBM譜線(xiàn)與實(shí)驗(yàn)測(cè)得PCBM薄膜的譜線(xiàn)進(jìn)行了對(duì)比,并將PCBM側(cè)鏈中的苯環(huán)以及富勒烯碳籠子對(duì)譜線(xiàn)的貢獻(xiàn)做了區(qū)分.Richter等20詳細(xì)地分析了PCBM的X射線(xiàn)吸收譜中的共振.現(xiàn)有PCBM的X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜的研究中,關(guān)注度最高的是在實(shí)驗(yàn)譜線(xiàn)中測(cè)得的第一個(gè)π*共振吸收峰的右肩峰(能量比主峰高大約0.5 eV).12,19–21這一特征峰在C60的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜中是不存在的,Bazylewski等19認(rèn)為PCBM側(cè)鏈的存在會(huì)導(dǎo)致C60最低未占據(jù)分子軌道的強(qiáng)離域,進(jìn)而使得軌道的對(duì)稱(chēng)性破缺是產(chǎn)生這一特征峰的原因.Patnaik等22研究金表面上的C60功能化自組裝單層膜的X射線(xiàn)吸收譜時(shí)也發(fā)現(xiàn)了類(lèi)似的高能肩峰,通過(guò)對(duì)比石墨譜線(xiàn)中相似的譜峰,這個(gè)肩峰的出現(xiàn)則被認(rèn)為是C60出現(xiàn)的聚集態(tài)所導(dǎo)致的.到目前為止,關(guān)于這個(gè)特征譜峰的來(lái)源還沒(méi)有確切的結(jié)論,因此我們需要深入的理論研究來(lái)解釋這個(gè)譜峰.

    本工作中,我們利用密度泛函理論(DFT)對(duì) PCBM的碳原子1s軌道的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜進(jìn)行了研究.計(jì)算中我們采用孤立的PCBM分子,同時(shí)對(duì)分子中不同組成部分的貢獻(xiàn)分別進(jìn)行了研究.詳細(xì)地研究了譜線(xiàn)中的第一個(gè)π*共振的精細(xì)結(jié)構(gòu).結(jié)果顯示第一個(gè)π*共振吸收是由富勒烯和苯環(huán)中碳原子的1s軌道到LUMO的躍遷組成的,而另外兩個(gè)主要的共振則主要來(lái)源于碳籠子中碳原子的貢獻(xiàn).第一個(gè)π*共振高能區(qū)的肩峰由三個(gè)小的吸收峰組成,這三個(gè)吸收峰主要來(lái)自碳籠子中部分特定碳原子的貢獻(xiàn).

    2 計(jì)算方法

    PCBM分子的幾何結(jié)構(gòu)(圖1所示)是在Gaussian 09程序包23中通過(guò)完全放開(kāi)所有原子進(jìn)行優(yōu)化的.在優(yōu)化分子結(jié)構(gòu)的計(jì)算過(guò)程中采用雜化的B3LYP泛函,24對(duì)于所有的原子,基函數(shù)都選用6-31G(d,p).25在優(yōu)化得到PCBM分子幾何結(jié)構(gòu)的基礎(chǔ)上,利用StoBe程序包26計(jì)算了碳元素1s 軌道的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜.近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜對(duì)原子的局域化學(xué)環(huán)境非常敏感,而非局域化學(xué)環(huán)境對(duì)其影響不大,因此分子間的弱相互作用對(duì)該譜線(xiàn)的作用可以忽略.Luo及其合作者27–30利用孤立分子模型計(jì)算的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜與實(shí)驗(yàn)測(cè)量固相薄膜樣品的譜線(xiàn)符合較好.有機(jī)分子對(duì)X射線(xiàn)的吸收過(guò)程對(duì)應(yīng)于兩個(gè)電子態(tài):基態(tài)和核電子激發(fā)態(tài).根據(jù)末態(tài)定則,31,32X射線(xiàn)吸收譜可以利用核電子激發(fā)態(tài)的波函數(shù)進(jìn)行計(jì)算得到.對(duì)于核電子被激發(fā)后的電子態(tài),本工作中采用完全芯態(tài)空穴(FCH)近似來(lái)描述,29該近似方法已經(jīng)廣泛地應(yīng)用于X射線(xiàn)譜的理論計(jì)算,其計(jì)算結(jié)果與實(shí)驗(yàn)譜線(xiàn)符合較好.29,33–37X射線(xiàn)吸收的振子強(qiáng)度由下式計(jì)算:

    這里,Iif是從躍遷振子強(qiáng)度,m是電子質(zhì)量,?是約化普朗克常數(shù),Ψif是躍遷過(guò)程中涉及到的兩個(gè)分子軌道,ωif是兩個(gè)軌道的能級(jí)間隔.x、y、z三個(gè)方向分量取平均值則是由于分子取向的隨機(jī)分布導(dǎo)致的統(tǒng)計(jì)效果.得到上述的能級(jí)間隔和吸收振子強(qiáng)度之后,需要利用第一激發(fā)能,也就是從1s軌道到LUMO的激發(fā)能去標(biāo)定整個(gè)譜線(xiàn).我們利用ΔKohn-Sham(ΔKS)方法38完全考慮第一激發(fā)過(guò)程的電子弛豫,得到第一激發(fā)能:

    這里N指系統(tǒng)里包含的電子數(shù).然后,將標(biāo)定之后的激發(fā)能和吸收振子強(qiáng)度根據(jù)不等價(jià)碳原子的權(quán)重利用高斯分布函數(shù)卷積,即可得到近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜.在利用高斯函數(shù)卷積時(shí),采用分段的半高全寬(FWHM)技術(shù).高于電離能(IP)的能量范圍屬于連續(xù)態(tài),我們采用Stieltjes方法38,39處理,低于電離能的能量范圍,則根據(jù)實(shí)驗(yàn)譜線(xiàn)分辨率不同分別采用0.4和0.1 eV作為FWHM.分子中碳元素1s軌道的電離能是利用ΔKS方法38進(jìn)行計(jì)算的,該電離能在數(shù)值上等于核電離態(tài),即FCH與GS之間的能量差I(lǐng)P=N–1EFCH–NEGS.這里N是基態(tài)的電子數(shù)目,N–1則對(duì)應(yīng)1s軌道的一個(gè)電子被電離之后的核電子電離態(tài)的電子數(shù)目.計(jì)算得到的富勒烯碳籠子中碳原子1s軌道的電離能范圍在289.8–290.6 eV之間,而苯環(huán)中碳原子的電離能在289.6–289.9 eV之間,側(cè)鏈上其他碳原子的電離能在289.9–293.0 eV之間.

    X射線(xiàn)譜的計(jì)算同樣也是利用密度泛函理論進(jìn)行的,我們采用的是包含梯度校正的Becke 交換泛函40和Perdew關(guān)聯(lián)泛函.41為了描述核電子被激發(fā)的碳原子,計(jì)算中我們采用了Kutzelnigg,Fleischer和Schindler發(fā)展的IGLO-III基函數(shù).42對(duì)于其余的碳原子,在StoBe程序中采用了描述4個(gè)電子的有效核勢(shì).對(duì)于其他原子,在計(jì)算中采用的是TZVP基函數(shù).43

    圖1 PCBM分子的幾何結(jié)構(gòu)Fig.1 Geometry structure of PCBM molecule

    3 結(jié)果與討論

    PCBM分子優(yōu)化后的結(jié)構(gòu)如圖1所示.結(jié)果顯示該分子屬于Cs對(duì)稱(chēng)性,其對(duì)稱(chēng)面為C2和C3兩原子的中垂面,因此PCBM分子具有42個(gè)不等價(jià)碳原子(碳籠子中32個(gè),側(cè)鏈中10個(gè)).對(duì)比富勒烯C60分子,PCBM的碳籠子部分出現(xiàn)較大的形變,在B3LYP/6-31G(d,p)計(jì)算水平下的形變能約為119.32 kJ·mol–1.從分子結(jié)構(gòu)上看,碳籠子中遠(yuǎn)離側(cè)鏈部分的碳碳鍵的鍵長(zhǎng)變化不大,而靠近側(cè)鏈部分的碳碳鍵長(zhǎng)的變化達(dá)到0.01–0.02 nm.該形變已經(jīng)打破了球形富勒烯C60的Ih高對(duì)稱(chēng)性.計(jì)算結(jié)果顯示PCBM分子中碳籠子部分的密里根電荷為–0.1e,表明當(dāng)C60分子與側(cè)鏈相連后存在從側(cè)鏈到碳籠子的電荷轉(zhuǎn)移.對(duì)稱(chēng)性的破壞以及上述的電荷轉(zhuǎn)移都會(huì)對(duì)近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜造成影響.

    圖2 理論計(jì)算的PCBM分子及其組成部分的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜與實(shí)驗(yàn)譜線(xiàn)12的對(duì)比Fig.2 Comparison of the calculated near-dege X-ray absorption fine structure spectra of PCBM molecule and contributions from selected moieties with the experimental spectrum12

    我們計(jì)算的PCBM分子碳原子1s軌道的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜(FWHM=0.4 eV)和Watts等12測(cè)量的實(shí)驗(yàn)譜線(xiàn)如圖2所示.根據(jù)PCBM分子的對(duì)稱(chēng)性以及碳原子化學(xué)環(huán)境的不同,我們將碳原子分為六類(lèi),每類(lèi)碳原子的譜線(xiàn)也列在圖2中用以研究實(shí)驗(yàn)譜線(xiàn)中各共振吸收峰的來(lái)源.我們將實(shí)驗(yàn)譜線(xiàn)中的譜峰標(biāo)定了數(shù)字(1–7),各譜峰的能量分別為284.5 eV(1),285.8 eV(2),286.2 eV(3),288.4 eV(4),290.8 eV(5),292.7 eV(6)和295.3 eV(7).1號(hào)譜峰高能區(qū)肩峰中心能量大約為285.0 eV,被標(biāo)定為S.結(jié)果顯示,我們理論計(jì)算的PCBM分子譜線(xiàn)與實(shí)驗(yàn)譜線(xiàn)吻合較好.而與C60富勒烯的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜18相比,PCBM分子譜線(xiàn)中的特征譜峰S主要來(lái)源于碳籠子中碳原子的貢獻(xiàn).為了研究實(shí)驗(yàn)譜線(xiàn)中各共振吸收峰的來(lái)源,我們將計(jì)算得到的譜線(xiàn)與其進(jìn)行了詳細(xì)的對(duì)比.1號(hào)譜峰與實(shí)驗(yàn)值符合較好,該譜峰包含碳籠子中碳原子(圖1中所示1–60號(hào)碳原子)以及苯環(huán)中碳原子(圖1中所示62–67號(hào)碳原子)的貢獻(xiàn).2號(hào)和3號(hào)譜峰包含碳籠子中碳原子的貢獻(xiàn).而側(cè)鏈部分碳原子的貢獻(xiàn)集中在從4號(hào)譜峰開(kāi)始的高能量范圍,以及4號(hào)譜峰低能區(qū)大約在287.4 eV的肩峰.理論計(jì)算得到的在287.4 eV處的肩峰以及4號(hào)譜峰與5號(hào)譜峰之間的部分譜峰在實(shí)驗(yàn)譜線(xiàn)中并未發(fā)現(xiàn).另外計(jì)算得到的2號(hào)譜峰與3號(hào)譜峰之間的強(qiáng)度明顯比實(shí)驗(yàn)譜線(xiàn)中二者的強(qiáng)度比要大.

    圖3 (a)理論計(jì)算和(b)實(shí)驗(yàn)測(cè)量12的π*共振區(qū)域的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜Fig.3 (a) Calculated and(b) experimental12near-edge X-ray absorption fine structure spectra of PCBM in the π*-resonance region

    圖4 理論計(jì)算的(a) PCBM;(b) 富勒烯;(c) C8;(d) C29以及(e) C62原子的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜Fig.4 Calculated near-edge X-ray absorption fine structure spectra of(a) PCBM;(b) the fullerene;(c) C8,(d) C29,and(e) C62 atoms

    本工作中,我們比較感興趣的是實(shí)驗(yàn)譜線(xiàn)中1號(hào)譜峰高能區(qū)的肩峰S,該肩峰在實(shí)驗(yàn)測(cè)量中可以較清晰地分辨出來(lái),然而在理論計(jì)算結(jié)果中卻不明顯.從計(jì)算的譜線(xiàn)中可以看出第一個(gè)π*共振峰的高能區(qū)部分是不對(duì)稱(chēng)的.實(shí)驗(yàn)譜線(xiàn)中該肩峰的強(qiáng)度高于我們計(jì)算得到的1號(hào)譜峰高能區(qū)不對(duì)稱(chēng)部分.為了準(zhǔn)確地分析計(jì)算得到的該肩峰的來(lái)源,我們將該π*共振吸收區(qū)域放大并應(yīng)用了更小的半高全寬(FWHM=0.1 eV),以便更清晰地顯示該譜峰的不對(duì)稱(chēng)性,如圖3所示.我們計(jì)算的1號(hào)譜峰非對(duì)稱(chēng)部分主要由介于284.6和285.2 eV之間三個(gè)獨(dú)立的譜峰導(dǎo)致.在圖4中,將這三個(gè)獨(dú)立的譜峰標(biāo)定為S1、S2和S3.為了研究該三個(gè)譜峰的來(lái)源,首先研究了所有碳原子的1s軌道在該能量區(qū)域的譜線(xiàn),來(lái)確定分子中對(duì)該譜峰提供主要貢獻(xiàn)的部分碳原子,然后我們根據(jù)計(jì)算的這些共振吸收對(duì)應(yīng)的躍遷能再進(jìn)行推斷.圖4中給出了PCBM分子總的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜和富勒烯碳籠子的貢獻(xiàn),以及62號(hào)碳原子和碳籠子中另外兩個(gè)碳原子C8和C29的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜.從圖2中可見(jiàn)除了C2和C3以外,碳籠子中所有的碳原子都對(duì)1號(hào)譜峰高能區(qū)部分有貢獻(xiàn).碳籠子中大部分碳原子對(duì)S1和S2譜峰有貢獻(xiàn),其中具有主要貢獻(xiàn)的碳原子(C18,C19,C23,C25,C27,C28,C29,C39,C40,C43,C45,C46,C47,C52,C53,C54)來(lái)自分子的左右兩個(gè)半球,如圖1所示.上述碳原子都屬于碳籠子,它們的化學(xué)環(huán)境比較類(lèi)似,其譜線(xiàn)也基本相同,因此將計(jì)算的C29譜線(xiàn)作為該系列碳原子的例子在圖4中給出.對(duì)譜峰S3貢獻(xiàn)較多的原子(C4,C5,C6,C7,C8,C9,C10,C11,C12,C13)主要集中于靠近側(cè)鏈連接點(diǎn)的部分,如圖1所示.圖4中給出了C8的譜線(xiàn)作為該部分碳原子的示例.單個(gè)碳原子π*共振能量偏移的不同主要是由于其1s基態(tài),也就是化學(xué)環(huán)境的差異.C61是側(cè)鏈與碳籠子之間的紐帶碳原子,側(cè)鏈中的碳原子除了C62,也就是苯環(huán)中跟C61相連的碳原子之外,其他碳原子對(duì)S1、S2和S3三個(gè)譜峰都沒(méi)有貢獻(xiàn).C62對(duì)S1譜峰有較強(qiáng)的貢獻(xiàn),而譜峰S1也正是來(lái)自該碳原子最低能量的共振,也就是從C62的1s軌道到LUMO之間的躍遷.由于其與C61相連,相比于苯環(huán)中其他的碳原子,C62與三個(gè)碳原子成鍵,因此C62的1s軌道能量被降低了,導(dǎo)致從1s軌道到LUMO的躍遷往高能方向偏移.苯環(huán)中其他碳原子最低的躍遷則對(duì)1號(hào)譜峰的主峰有主要貢獻(xiàn).表1中分別列出了碳籠子和苯環(huán)中碳原子1s軌道到LUMO,LUMO+1或 LUMO+2軌道躍遷的躍遷能及其對(duì)應(yīng)譜峰的標(biāo)定.從表1中可以看到,1號(hào)譜峰是由碳籠子和苯環(huán)中碳原子的1s軌道到LUMO的躍遷導(dǎo)致的,而S1譜峰是由富勒烯碳籠子中1s軌道到LUMO+1的躍遷以及C62的1s軌道到LUMO的躍遷貢獻(xiàn)的.S2譜峰主要由富勒烯中碳原子1s軌道到LUMO+2和一部分1s軌道到LUMO+ 1的躍遷組成.而肩峰的第三部分S3譜峰,僅由富勒烯碳籠子中一部分碳原子的1s軌道到LUMO+2的電子躍遷導(dǎo)致,因此該譜峰的強(qiáng)度比前兩個(gè)譜峰要低.

    表1 理論計(jì)算的1號(hào)譜峰及其肩峰(S1、S2、S3)對(duì)應(yīng)的富勒烯和苯環(huán)中碳原子的平均電離能(IPs)、譜峰的平均躍遷能(Eexc)以及對(duì)應(yīng)的軌道躍遷Table 1 Calculated average ionization potentials(IPs) and excitation energies(Eexc) of the fullerene and the phenyl components of peak 1 and the shoulder peaks(S1,S2,S3),with the assignment to the corresponding transition

    我們的計(jì)算結(jié)果顯示富勒烯碳籠子與側(cè)鏈之間的碳原子C61對(duì)1號(hào)譜峰、2號(hào)譜峰以及譜峰S都沒(méi)有貢獻(xiàn).實(shí)際上,C61的1s軌道到未占據(jù)分子軌道的第一個(gè)躍遷處在高能區(qū),并對(duì)3號(hào)譜峰有貢獻(xiàn),如圖2所示.C61原子譜線(xiàn)中第二個(gè)譜峰對(duì)4號(hào)譜峰的左肩峰有較小的貢獻(xiàn),而C61譜線(xiàn)中的第三個(gè)譜峰僅對(duì)4號(hào)譜峰有貢獻(xiàn).同時(shí)發(fā)現(xiàn)C2和C3碳原子對(duì)該肩峰S沒(méi)有貢獻(xiàn).另外,計(jì)算結(jié)果顯示C2和C3這兩個(gè)與C61化學(xué)成鍵的碳原子譜線(xiàn)與C61的譜線(xiàn)非常接近.相比C61的譜線(xiàn),C2和C3的譜線(xiàn)只是整體地往低能方向平移了0.3 eV.計(jì)算得到的C61、C2和C3的譜線(xiàn)與富勒烯碳籠子中其他碳原子之間譜線(xiàn)的區(qū)別主要?dú)w結(jié)于其軌道雜化的不同.C61、C2和C3有四個(gè)共價(jià)鍵,因此它們是sp3雜化的,而碳籠子中的其他碳原子都是sp2雜化的.44

    我們理論計(jì)算的譜線(xiàn)和他人實(shí)驗(yàn)測(cè)量的譜線(xiàn)均具有一個(gè)非常有意義的特點(diǎn),就是第一個(gè)譜峰的強(qiáng)度非常大.該譜峰的高強(qiáng)度主要是由于在C60和C70等物質(zhì)中發(fā)現(xiàn)的激子效應(yīng).29激子是一種由電子和空穴之間的庫(kù)侖相互作用導(dǎo)致的準(zhǔn)粒子.電子空穴相互作用會(huì)使得未占據(jù)分子軌道的能量降低,同時(shí)使得原子中產(chǎn)生芯態(tài)空穴附近的分子軌道更加局域化,從而使1s軌道到LUMO的躍遷強(qiáng)度增大.苯環(huán)中碳原子對(duì)1號(hào)譜峰有較強(qiáng)的貢獻(xiàn),主要是因?yàn)楸江h(huán)中也有類(lèi)似的激子效應(yīng).苯環(huán)中碳原子上產(chǎn)生芯態(tài)空穴后導(dǎo)致的LUMO在苯環(huán)附近的局域性,甚至比LUMO在碳籠子中的局域性更大.當(dāng)芯態(tài)空穴出現(xiàn)在酯基中時(shí),類(lèi)似的激子效應(yīng)也可以在計(jì)算的譜線(xiàn)中看到,該效應(yīng)使得4號(hào)譜峰的肩峰強(qiáng)度增強(qiáng).

    4 結(jié) 論

    利用密度泛函理論,結(jié)合完全芯態(tài)空穴近似方法對(duì)富勒烯衍生物PCBM的近邊X射線(xiàn)吸收精細(xì)結(jié)構(gòu)譜進(jìn)行了理論研究.理論計(jì)算的譜線(xiàn)與實(shí)驗(yàn)譜線(xiàn)符合較好.通過(guò)深入地分析PCBM分子中各組成部分單獨(dú)的譜線(xiàn),我們對(duì)實(shí)驗(yàn)譜線(xiàn)中的各個(gè)譜峰都進(jìn)行了理論標(biāo)定.實(shí)驗(yàn)譜線(xiàn)中的第一個(gè)π*共振譜峰對(duì)應(yīng)于富勒烯和苯環(huán)中碳原子的1s軌道到LUMO軌道的躍遷.第二個(gè)和第三個(gè)譜峰則源自富勒烯碳籠子中碳原子1s軌道到高能未占據(jù)分子軌道的躍遷,而側(cè)鏈中的碳原子1s軌道到未占據(jù)分子軌道的躍遷則對(duì)實(shí)驗(yàn)測(cè)量的第四個(gè)譜峰有較大的貢獻(xiàn).實(shí)驗(yàn)譜線(xiàn)中第一個(gè)譜峰對(duì)應(yīng)的1s軌道到LUMO的躍遷強(qiáng)度具有最大值是由激子效應(yīng)導(dǎo)致的.實(shí)驗(yàn)譜線(xiàn)中觀測(cè)到的在第一個(gè)譜峰右側(cè)能量大約為285.0 eV的肩峰主要是來(lái)自碳籠子中碳原子的貢獻(xiàn).

    (1)Zhuo,Z.L.;Zhang,F.J.;Xu,X.W.;Wang,J.;Lu,L.F.;Xu,Z.Acta Phys.-Chim.Sin.2011,27(4),875.[卓祖亮,張???許曉偉,王 健,盧麗芳,徐 征.物理化學(xué)學(xué)報(bào),2011,27(4),875.] doi:10.3866/PKU.WHXB20110414

    (2)Li,D.;Liang,R.;Yue,H.;Wang,P.;Fu,L.M.;Zhang,J.P.;Ai,X.C.Acta Phys.-Chim.Sin.2012,28(6),1373.[李 丹,梁然,岳 鶴,王 鵬,付立民,張建平,艾希成.物理化學(xué)學(xué)報(bào),2012,28(6),1373.] doi:10.3866/PKU.WHXB201204061

    (3)St?hr,J.NEXAFS Spectroscopy;Springer Verlag:Berlin,1996;pp 1–3.

    (4)Zhu,M.Q.;Pan,G.;Liu,T.;Li,X.L.;Yang,Y.H.;Li,W.;Li,J.;Hu,T.D.;Wu,Z.Y.;Xie,Y.N.Acta Phys.-Chim.Sin.2005,21(12),1378.[朱孟強(qiáng),潘 綱,劉 濤,李賢良,楊玉環(huán),李薇,李 晉,胡天斗,吳自玉,謝亞寧.物理化學(xué)學(xué)報(bào),2005,21(12),1378.] doi:10.3866/PKU.WHXB20051210

    (5)Guo,H.L.;Wang,J.Y.;Wu,Z.H.;Jiang,S.C.Acta Polymerica Sinica 2014,No.2,179.[郭慧龍,王佳怡,吳忠華,蔣世春.高分子學(xué)報(bào),2014,No.2,179.]

    (6)Germack,D.S.;Chan,C.K.;Hamadani,B.H.;Richter,L.J.;Fischer,D.A.;Gundlach,D.J.;DeLongchamp,D.M.Appl.Phys.Lett.2009,94,233303.doi:10.1063/1.3149706

    (7)Germack,D.S.;Chan,C.K.;Kline,R.J.;Fischer,D.A.;Gundlach,D.J.;Toney,M.F.;Richter,L.J.;DeLongchamp,D.M.Macromolecules 2010,43(8),3828.doi:10.1021/ma100027b

    (8)Xue,B.;Vaughan,B.;Poh,C.H.;Burke,K.B.;Thomsen,L.;Stapleton,A.;Zhou,X.;Bryant,G.W.;Belcher,W.;Dastoor,P.C.J.Phys.Chem.C 2010,114(37),15797.doi:10.1021/jp104695j

    (9)Tillack,A.F.;Noone,K.M.;MacLeod,B.A.;Nordlund,D.;Nagle,K.P.;Bradley,J.A.;Hau,S.K.;Yip,H.L.;Jen,A.K.Y.;Seidler,G.T.;Ginger,D.S.ACS Appl.Mater.Interfaces 2011,3(3),726.doi:10.1021/am101055r

    (10)Anselmo,A.S.;Dzwilewski,A.;Svensson,K.;Moons,E.J.Polm.Sci.Part B:Polym.Phys.2013,51(3),76.

    (11)Anselmo,A.S.;Lindgren,L.;Rysz,J.;Bernasik,A.;Budkowski,A.;Andersson,M.;Svensson,K.;van Stam,J.;Moons,E.Chem.Mater.2011,23(9),2295.doi:10.1021/cm1021596

    (12)Watts,B.;Swaraj,S.;Nordlund,D.;Lüning,J.;Ade,H.J.Chem.Phys.2011,134,024702.doi:10.1063/1.3506636

    (13)DeLongchamp,D.M.;Lin,E.K.;Fischer,D.A.Proc.SPIE 2005,5940,59400A.

    (14)McNeill,C.R.;Ade,H.J.Mater.Chem.C 2013,1(2),187.doi:10.1039/C2TC00001F

    (15)Mikoushkin,V.M.;Shnitov,V.V.;Bryzgalov,V.V.;Gordeev,Y.S.;Boltalina,O.V.;Gol'dt,I.V.;Molodtsov,S.L.;Vyalykh,D.V.Nanotubes Carbon Nanostruct.2008,16(5–6),588.doi:10.1080/15363830802286574

    (16)Tang,Y.H.;Sham,T.K.;Hu,Y.F.;Lee,C.S.;Lee,S.T.Chem.Phys.Lett.2002,366(5–6),636.doi:10.1016/S0009-2614(02)01620-2

    (17)Pacilé,D.;Papagno,M.;Rodríguez,A.F.;Grioni,M.;Papagno,L.;Girit,C.?.;Meyer,J.C.;Begtrup,G.E.;Zettl,A.Phys.Rev.Lett.2008,101,066806.doi:10.1103/PhysRevLett.101.066806

    (18)Terminello,L.J.;Shuh,D.K.;Himpsel,F.J.;Lapiano-Smith,D.A.;St?hr,J.;Bethune,D.;Meijer,S.G.Chem.Phys.Lett.1991,182(5),491.doi:10.1016/0009-2614(91)90113-N

    (19)Bazylewski,P.F.;Kim,K.H.;Forrest,J.L.;Tada,H.;Choi,D.H.Chem.Phys.Lett.2011,508(1–3),90.doi:10.1016/j.cplett.2011.04.017

    (20)Richter,M.H.;Friedrich,D.;Schmeiber,D.BioNanoSci 2012,2(1),59.doi:10.1007/s12668-011-0034-1

    (21)Friedrich,D.;Henkel,K.;Richter,M.;Schmeiber,D.BioNanoSci 2011,1(4),218.doi:10.1007/s12668-011-0025-2

    (22)Patnaik,A.;Okudaira,K.K.;Kera,S.;Setoyama,H.;Mase,K.;Ueno,N.J.Chem.Phys.2005,122(15),154703.doi:10.1063/1.1880952

    (23)Frisch,M.J.;Trucks,G.W.;Schlegel,H.B.;et al.Gaussian 09,Revision A.01;Gaussian,Inc.:Wallingford,CT,2009.

    (24)Becke,A.D.J.Chem.Phys.1993,98(7),5648.doi:10.1063/1.464913

    (25)Rassolov,V.;Pople,J.A.;Ratner,M.;Windus,T.L.J.Chem.Phys.1998,109(4),1223.doi:10.1063/1.476673

    (26)Hermann,K.;Pettersson,L.;Casida,M.;et al.StoBe Version 3.0;StoBe Software:Stockholm,Sweden,2007.

    (27)Luo,Y.;?gren,H.;Keil,M.;Friedlein,R.;Salaneck,W.R.Chem.Phys.Lett.2001,337,176.doi:10.1016/S0009-2614(01)00181-6

    (28)Hellgren,N.;Guo,J.;S?the,C.;Agui,A.;Nordgren,J.;Luo,Y.;?gren,H.;Sundgren,J.E.Appl.Phys.Lett.2001,79(26),4348.doi:10.1063/1.1428108

    (29)Nyberg,M.;Luo,Y.;Triguero,L.;Pettersson,L.G.M.;?gren,H.Phys.Rev.B 1999,60,7956.doi:10.1103/PhysRevB.60.7956

    (30)Carlegrim,E.;Gao,B.;Kanciurzewska,A.;de Jong,M.P.;Wu,Z.;Luo,Y.;Fahlman,M.Phys.Rev.B 2008,77,054420.doi:10.1103/PhysRevB.77.054420

    (31)von Barth,U.;Grossman,G.Solid State Commun.1979,32(8),645.doi:10.1016/0038-1098(79)90719-1

    (32)von Barth,U.;Grossman,G.Phys.Rev.B 1982,25,5150.doi:10.1103/PhysRevB.25.5150

    (33)Luo,Y.;?gren,H.;Gelmukhanov,F.;Guo,J.;Skytt,P.;Wassdahl,N.;Nordgren,J.Phys.Rev.B 1995,52,14479.

    (34)Gao,B.;Liu,L.;Wang,C.;Wu,Z.;Luo,Y.J.Chem.Phys.2007,127(16),164314.doi:10.1063/1.2800028

    (35)Zhao,T.;Gao,B.;Liu,L.;Ye,Q.;Chu,W.S.;Wu,Z.Y.Chin.Phys.C 2009,33(11),954.doi:10.1088/1674-1137/33/11/005

    (36)Qi,J.;Hua,W.;Gao,B.Chem.Phys.Lett.2012,539–540,222.

    (37)Song,X.;Ma Y.;Wang,C.;Dietrich,P.D.;Unger,W.E.S.;Luo,Y.J.Phys.Chem.C 2012,116(23),12649.doi:10.1021/jp302716w

    (38)Triguero,L.;Pettersson,L.G.M.;?gren,H.Phys.Rev.B 1998,58,8097.doi:10.1103/PhysRevB.58.8097

    (39)Triguero,L.;Plashkevych,O.;Pettersson,L.G.M.;?gren,H.J.Electron Spectrosc.Relat.Phenom.1999,104(1–3),195.doi:10.1016/S0368-2048(99)00008-0

    (40)Becke,A.D.Phys.Rev.A 1988,38,3098.doi:10.1103/PhysRevA.38.3098

    (41)Perdew,J.P.Phys.Rev.B 1986,33,8822.doi:10.1103/PhysRevB.33.8822

    (42)Kutzelnigg,W.;Fleischer,U.;Schindler,M.NMR:Basic Principles and Progress;Springer Verlag:Berlin Heidelberg,1990;Vol.213.

    (43)Sch?fer,A.;Huber,C.;Ahlrichs,R.J.Chem.Phys.1994,100(8),5829.doi:10.1063/1.467146

    (44)Dresselhaus,M.S.;Dresselhaus,G.;Eklund,P.C.Science of Fullerenes and Carbon Nanotubes;Academic Press:London,1996.

    First-Principles Study on the Near-Edge X-ray Absorption Fine Structure Spectroscopy of the Fullerene-Derivative PCBM

    MA Yong1WANG Guang-Wei1SUN Shao-Tao2SONG Xiu-Neng1,*
    (1School of Physics and Electronics,Shandong Normal University,Jinan 250014,P.R.China;
    2Administration Center,Shandong Academy of Information and Communication Technology,Jinan 250101,P.R.China)

    Fullerene-derivative [6,6]-phenyl-C61-butyric acid methyl ester(PCBM) plays an important role in terms of electron transport in polymer solar cells.The electronic structure of PCBM is of much importance to investigate.In this study,the near-edge X-ray absorption fine structure spectroscopy and unoccupied orbitals of PCBM were researched with density functional theory.By comparing the calculated sum spectra of nonequivalent carbon atoms,we assigned the main resonances of PCBM.The origin of the shoulder in the right side of the first π* resonance was analyzed,and the results showed that this absorption peak was mainly contributed by the transitions to higher unoccupied orbitals of the unmodified carbons in the C60cage.

    Density functional theory; X-ray; Fine structure; Fullerene-derivative; Full core hole

    January26,2015;Revised:May 22,2015;Published on Web:May 25,2015.

    O641

    icle]

    10.3866/PKU.WHXB201505251 www.whxb.pku.edu.cn

    *Corresponding author.Email:xiuneng@sdnu.edu.cn;Tel:+86-13953138878.

    The project was supported by the National Natural Science Foundation of China(21303096,11374195),Promotive Research Fund for Young and Middle-Aged Scientists of Shandong Province,China(BS2013CL016),China Postdoctoral Science Foundation(2013M541951),and Scientific Research Foundation for the Returned Overseas Chinese Scholars,Ministry of Education,China.

    國(guó)家自然科學(xué)基金(21303096,11374195),山東省優(yōu)秀中青年科學(xué)家科研獎(jiǎng)勵(lì)基金(BS2013CL016),中國(guó)博士后科學(xué)基金(2013M541951)和教育部留學(xué)回國(guó)人員科研啟動(dòng)基金資助項(xiàng)目

    ? Editorial office of Acta Physico-Chimica Sinica

    猜你喜歡
    富勒烯譜峰碳原子
    建構(gòu)模型認(rèn)知 突破有機(jī)物分子中原子共平面問(wèn)題
    連續(xù)波體制引信多譜峰特性目標(biāo)檢測(cè)方法
    X射線(xiàn)光電子能譜復(fù)雜譜圖的非線(xiàn)性最小二乘法分析案例
    基于無(wú)基底扣除的數(shù)據(jù)趨勢(shì)累積譜峰檢測(cè)算法
    色譜(2021年6期)2021-05-06 02:18:56
    巖性密度測(cè)井儀工作原理與典型故障分析
    科技資訊(2020年12期)2020-06-03 04:44:20
    鉆石級(jí)抗衰 諾貝爾獎(jiǎng)光環(huán)揭開(kāi)“富勒烯”的神秘面紗
    碳鏈異構(gòu)有關(guān)的同分異構(gòu)體書(shū)寫(xiě)補(bǔ)遺
    富勒烯基聚合物制備研究的進(jìn)展
    有機(jī)化合物命名易錯(cuò)題直擊
    新型富勒烯的合成
    a级毛片在线看网站| 国产色婷婷99| 精品亚洲乱码少妇综合久久| 久热这里只有精品99| 免费大片黄手机在线观看| 永久网站在线| 少妇人妻 视频| 多毛熟女@视频| 91精品伊人久久大香线蕉| 简卡轻食公司| tube8黄色片| 欧美老熟妇乱子伦牲交| 欧美xxⅹ黑人| 美女主播在线视频| 一级毛片电影观看| 亚洲国产色片| 国产成人精品一,二区| 亚洲性久久影院| 欧美xxⅹ黑人| av免费在线看不卡| 91精品伊人久久大香线蕉| 久久99热6这里只有精品| 久久久久久久亚洲中文字幕| 成人亚洲欧美一区二区av| 亚洲精品亚洲一区二区| 亚洲国产精品专区欧美| 国产高清有码在线观看视频| 国产无遮挡羞羞视频在线观看| 亚洲三级黄色毛片| 国产免费一级a男人的天堂| 又大又黄又爽视频免费| 成人18禁高潮啪啪吃奶动态图 | 亚洲三级黄色毛片| av.在线天堂| 免费观看a级毛片全部| 一级毛片电影观看| 午夜激情福利司机影院| 80岁老熟妇乱子伦牲交| av网站免费在线观看视频| 少妇人妻久久综合中文| 久久 成人 亚洲| 久久人人爽人人片av| tube8黄色片| 大香蕉久久网| 欧美人与善性xxx| 精品国产一区二区三区久久久樱花| 五月伊人婷婷丁香| 国产日韩一区二区三区精品不卡 | 大香蕉久久网| 草草在线视频免费看| 97在线视频观看| 亚洲国产最新在线播放| 国产在线一区二区三区精| 十八禁网站网址无遮挡| 亚洲熟女精品中文字幕| 综合色丁香网| 亚洲欧美日韩另类电影网站| 精品人妻在线不人妻| 高清午夜精品一区二区三区| 日日撸夜夜添| 亚洲欧美清纯卡通| av在线老鸭窝| 亚洲人成77777在线视频| 国产一区二区三区av在线| 欧美少妇被猛烈插入视频| 春色校园在线视频观看| 涩涩av久久男人的天堂| 国产av码专区亚洲av| 少妇丰满av| 亚洲精品久久午夜乱码| 日韩亚洲欧美综合| 如何舔出高潮| 亚洲精品一二三| 高清视频免费观看一区二区| 超色免费av| 精品午夜福利在线看| 国模一区二区三区四区视频| 中文字幕人妻熟人妻熟丝袜美| 国产探花极品一区二区| 国产午夜精品久久久久久一区二区三区| 午夜激情久久久久久久| 亚洲精品国产av成人精品| 亚洲国产成人一精品久久久| 国产免费现黄频在线看| 色吧在线观看| 国内精品宾馆在线| 久久婷婷青草| 这个男人来自地球电影免费观看 | 日韩在线高清观看一区二区三区| 亚洲av国产av综合av卡| 亚州av有码| 国产免费一区二区三区四区乱码| 91成人精品电影| 波野结衣二区三区在线| 日韩大片免费观看网站| 曰老女人黄片| 热99久久久久精品小说推荐| 免费人妻精品一区二区三区视频| 99九九线精品视频在线观看视频| 亚洲精品中文字幕在线视频| 国产在线免费精品| 婷婷色综合大香蕉| 亚洲成人手机| 国产精品麻豆人妻色哟哟久久| 欧美日韩综合久久久久久| 亚洲伊人久久精品综合| 能在线免费看毛片的网站| 亚洲不卡免费看| 国产精品99久久久久久久久| 欧美人与性动交α欧美精品济南到 | 欧美精品一区二区免费开放| 免费观看性生交大片5| 另类精品久久| 九草在线视频观看| 国产成人91sexporn| 少妇猛男粗大的猛烈进出视频| 如何舔出高潮| 国产男女超爽视频在线观看| 爱豆传媒免费全集在线观看| 国产精品嫩草影院av在线观看| 亚洲图色成人| 在线天堂最新版资源| 五月伊人婷婷丁香| 亚洲情色 制服丝袜| 国产片特级美女逼逼视频| 一个人免费看片子| 国产精品久久久久久久电影| 大码成人一级视频| 午夜91福利影院| 国产乱人偷精品视频| 成人亚洲精品一区在线观看| 爱豆传媒免费全集在线观看| 国产精品熟女久久久久浪| 黑人巨大精品欧美一区二区蜜桃 | 久久久久精品久久久久真实原创| 能在线免费看毛片的网站| 在线观看人妻少妇| av.在线天堂| 自拍欧美九色日韩亚洲蝌蚪91| 在现免费观看毛片| 亚洲精品久久成人aⅴ小说 | 国产精品久久久久久精品电影小说| 久久久久久久亚洲中文字幕| 蜜臀久久99精品久久宅男| 欧美性感艳星| 亚洲av不卡在线观看| 久久精品国产亚洲av涩爱| 91aial.com中文字幕在线观看| 黄片播放在线免费| 国产精品久久久久久精品电影小说| 乱人伦中国视频| 国产精品人妻久久久久久| 91精品三级在线观看| 国产成人一区二区在线| 国产成人午夜福利电影在线观看| 91成人精品电影| av在线老鸭窝| tube8黄色片| 黄色配什么色好看| 久久亚洲国产成人精品v| 亚洲欧美成人精品一区二区| 欧美激情国产日韩精品一区| 人人妻人人爽人人添夜夜欢视频| 国产高清三级在线| 亚洲四区av| 简卡轻食公司| 香蕉精品网在线| 亚洲欧美色中文字幕在线| 日本-黄色视频高清免费观看| 亚洲国产av新网站| 老司机影院成人| 精品一区二区三卡| 免费观看a级毛片全部| 黄色毛片三级朝国网站| av在线老鸭窝| 赤兔流量卡办理| 大香蕉久久成人网| av卡一久久| 久热久热在线精品观看| 99国产综合亚洲精品| 精品久久蜜臀av无| 国产乱来视频区| 免费大片黄手机在线观看| 人人妻人人澡人人看| 日韩免费高清中文字幕av| 精品久久久噜噜| 母亲3免费完整高清在线观看 | 肉色欧美久久久久久久蜜桃| 精品久久国产蜜桃| 亚洲精品日本国产第一区| 午夜视频国产福利| 一区在线观看完整版| 午夜福利视频精品| 考比视频在线观看| 亚洲伊人久久精品综合| 国产精品国产三级国产专区5o| 欧美日本中文国产一区发布| 一区在线观看完整版| 亚洲国产最新在线播放| 国产精品人妻久久久久久| 欧美日韩视频精品一区| 大香蕉97超碰在线| 午夜免费男女啪啪视频观看| 人人妻人人澡人人看| 国产精品熟女久久久久浪| 精品一品国产午夜福利视频| 99国产综合亚洲精品| 高清黄色对白视频在线免费看| 亚洲国产av新网站| 人人澡人人妻人| 一区二区三区乱码不卡18| 好男人视频免费观看在线| 91精品伊人久久大香线蕉| 中文字幕精品免费在线观看视频 | 狂野欧美白嫩少妇大欣赏| 韩国高清视频一区二区三区| 亚洲av综合色区一区| 欧美一级a爱片免费观看看| 精品久久久噜噜| 久久久欧美国产精品| 久久久精品免费免费高清| 国产在线免费精品| 久久久国产一区二区| 天堂俺去俺来也www色官网| 伊人久久精品亚洲午夜| 女人精品久久久久毛片| 亚洲av综合色区一区| 欧美日韩综合久久久久久| 午夜久久久在线观看| 在线观看www视频免费| 亚洲成人av在线免费| 97在线视频观看| 中文字幕制服av| 三级国产精品欧美在线观看| 久久97久久精品| 午夜影院在线不卡| 午夜av观看不卡| 欧美日韩国产mv在线观看视频| 国产淫语在线视频| 国产国语露脸激情在线看| 男人爽女人下面视频在线观看| 日韩在线高清观看一区二区三区| 久久久久视频综合| 丝袜在线中文字幕| 边亲边吃奶的免费视频| 国产成人av激情在线播放 | 亚洲精品日韩av片在线观看| 女性被躁到高潮视频| 亚洲人成网站在线观看播放| 亚洲一级一片aⅴ在线观看| 韩国av在线不卡| 国产高清不卡午夜福利| 丰满迷人的少妇在线观看| 十八禁高潮呻吟视频| 久久热精品热| 亚洲婷婷狠狠爱综合网| 亚洲无线观看免费| 国产精品国产三级国产av玫瑰| 亚洲欧美日韩另类电影网站| 中文字幕人妻丝袜制服| 久久 成人 亚洲| 一级毛片aaaaaa免费看小| 精品久久久久久久久亚洲| 大香蕉97超碰在线| av有码第一页| 久久ye,这里只有精品| 亚洲不卡免费看| 亚洲激情五月婷婷啪啪| 亚洲国产最新在线播放| 亚洲精品美女久久av网站| 久久久久视频综合| 欧美亚洲日本最大视频资源| 亚洲激情五月婷婷啪啪| 国内精品宾馆在线| 亚洲经典国产精华液单| 18禁在线播放成人免费| 亚洲丝袜综合中文字幕| 极品人妻少妇av视频| 精品人妻熟女av久视频| 精品国产一区二区三区久久久樱花| 搡老乐熟女国产| 亚洲精品国产av蜜桃| 我的老师免费观看完整版| 国产综合精华液| 天天影视国产精品| www.色视频.com| 美女国产高潮福利片在线看| 黄色一级大片看看| 成人影院久久| 久久精品熟女亚洲av麻豆精品| av福利片在线| 久久精品国产亚洲网站| 欧美日韩综合久久久久久| 欧美性感艳星| 一区二区日韩欧美中文字幕 | 精品久久国产蜜桃| 国产日韩欧美亚洲二区| 国产欧美日韩综合在线一区二区| 一边摸一边做爽爽视频免费| 日本黄大片高清| 久久精品久久精品一区二区三区| av视频免费观看在线观看| 亚洲国产av新网站| 欧美成人精品欧美一级黄| 午夜福利网站1000一区二区三区| 好男人视频免费观看在线| 亚洲欧美成人综合另类久久久| 水蜜桃什么品种好| 在线天堂最新版资源| 久久99热这里只频精品6学生| 黄色欧美视频在线观看| 日本wwww免费看| 青春草国产在线视频| 人体艺术视频欧美日本| 亚洲色图综合在线观看| 亚洲国产精品999| 各种免费的搞黄视频| 国产一级毛片在线| freevideosex欧美| 精品人妻熟女av久视频| 男女免费视频国产| 国模一区二区三区四区视频| 亚洲中文av在线| 青春草亚洲视频在线观看| 精品亚洲成a人片在线观看| 亚洲精品日韩av片在线观看| 久久99热6这里只有精品| 久久99一区二区三区| 91久久精品国产一区二区成人| 久久毛片免费看一区二区三区| 一本色道久久久久久精品综合| 国产精品久久久久久久久免| 日产精品乱码卡一卡2卡三| 夜夜骑夜夜射夜夜干| 国产一区二区三区综合在线观看 | 最近2019中文字幕mv第一页| 亚洲欧美精品自产自拍| 午夜福利影视在线免费观看| 亚洲av二区三区四区| 中国国产av一级| 男人添女人高潮全过程视频| 国产成人免费观看mmmm| 中文字幕免费在线视频6| av网站免费在线观看视频| 一区二区三区乱码不卡18| 纯流量卡能插随身wifi吗| 国产精品久久久久久久电影| 91精品国产国语对白视频| 亚洲图色成人| 欧美日韩成人在线一区二区| 99热6这里只有精品| 高清黄色对白视频在线免费看| 王馨瑶露胸无遮挡在线观看| 久久久久国产网址| 人人澡人人妻人| 午夜视频国产福利| 色哟哟·www| 91成人精品电影| 欧美日韩视频精品一区| 另类精品久久| 国产午夜精品久久久久久一区二区三区| 水蜜桃什么品种好| 久久久久久久久久久久大奶| 国产精品久久久久成人av| 女性生殖器流出的白浆| 啦啦啦视频在线资源免费观看| 国产不卡av网站在线观看| av又黄又爽大尺度在线免费看| 久久精品国产鲁丝片午夜精品| 美女主播在线视频| 免费看av在线观看网站| 最新的欧美精品一区二区| 国产精品嫩草影院av在线观看| 欧美日韩综合久久久久久| 黄色毛片三级朝国网站| 精品酒店卫生间| 日韩成人av中文字幕在线观看| av视频免费观看在线观看| 五月玫瑰六月丁香| 欧美一级a爱片免费观看看| 免费播放大片免费观看视频在线观看| 免费日韩欧美在线观看| 日日爽夜夜爽网站| 久久久久精品性色| 秋霞伦理黄片| 亚洲精品亚洲一区二区| 一边亲一边摸免费视频| 我要看黄色一级片免费的| 成人综合一区亚洲| 18在线观看网站| 精品一区在线观看国产| 色婷婷久久久亚洲欧美| 中文字幕久久专区| 男的添女的下面高潮视频| 亚洲精品乱久久久久久| 免费观看性生交大片5| 丝瓜视频免费看黄片| 尾随美女入室| 日本与韩国留学比较| 久久影院123| 久久久久久久大尺度免费视频| 看十八女毛片水多多多| 国产精品一区二区在线不卡| 国产精品不卡视频一区二区| 成人漫画全彩无遮挡| 97超碰精品成人国产| 有码 亚洲区| 亚洲第一av免费看| 午夜激情久久久久久久| 秋霞伦理黄片| 伦理电影免费视频| 国产精品人妻久久久久久| a级毛片免费高清观看在线播放| 午夜免费观看性视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧美色中文字幕在线| 桃花免费在线播放| 亚洲欧美色中文字幕在线| 亚洲精品乱久久久久久| 久久国产精品男人的天堂亚洲 | 夫妻午夜视频| 国产又色又爽无遮挡免| 激情五月婷婷亚洲| 午夜免费男女啪啪视频观看| 在线观看国产h片| 少妇的逼水好多| av一本久久久久| 日韩欧美一区视频在线观看| 精品久久国产蜜桃| 国产午夜精品一二区理论片| 久久久久人妻精品一区果冻| 亚洲欧美成人精品一区二区| 欧美性感艳星| 国产一区二区三区综合在线观看 | 亚洲av不卡在线观看| 精品久久久久久久久亚洲| 婷婷成人精品国产| 最黄视频免费看| 蜜桃国产av成人99| 国产欧美日韩综合在线一区二区| 久久久久久久久久久丰满| 亚洲熟女精品中文字幕| 中文字幕免费在线视频6| 亚洲天堂av无毛| 天堂俺去俺来也www色官网| 街头女战士在线观看网站| 妹子高潮喷水视频| 久久精品国产自在天天线| 欧美成人精品欧美一级黄| 日本色播在线视频| 国产精品久久久久久精品古装| 汤姆久久久久久久影院中文字幕| 亚洲精品成人av观看孕妇| 超色免费av| 建设人人有责人人尽责人人享有的| 日韩av不卡免费在线播放| 丝袜脚勾引网站| 中文天堂在线官网| 国产精品国产三级国产专区5o| 街头女战士在线观看网站| 美女视频免费永久观看网站| 少妇人妻 视频| 国产日韩欧美亚洲二区| 国产男人的电影天堂91| 大又大粗又爽又黄少妇毛片口| 亚洲欧洲日产国产| kizo精华| 校园人妻丝袜中文字幕| 色94色欧美一区二区| 丰满少妇做爰视频| 久久人妻熟女aⅴ| 国产精品熟女久久久久浪| 三级国产精品欧美在线观看| 极品少妇高潮喷水抽搐| 久久99蜜桃精品久久| 丝袜在线中文字幕| 成人国产麻豆网| 精品亚洲成a人片在线观看| 日韩中文字幕视频在线看片| 少妇精品久久久久久久| 日韩中文字幕视频在线看片| tube8黄色片| 精品久久久久久久久av| 久久婷婷青草| 搡老乐熟女国产| 国产在线视频一区二区| 久久韩国三级中文字幕| 精品国产露脸久久av麻豆| 丰满少妇做爰视频| 欧美三级亚洲精品| 日韩亚洲欧美综合| 一级,二级,三级黄色视频| 久久精品国产亚洲网站| 免费观看a级毛片全部| 亚洲不卡免费看| 各种免费的搞黄视频| h视频一区二区三区| 国产片内射在线| 一个人看视频在线观看www免费| 91精品伊人久久大香线蕉| 国产av码专区亚洲av| 亚洲国产色片| 久久99精品国语久久久| 欧美 日韩 精品 国产| 日韩不卡一区二区三区视频在线| 日本91视频免费播放| 有码 亚洲区| 人妻一区二区av| 久久韩国三级中文字幕| 免费观看av网站的网址| 最近2019中文字幕mv第一页| 欧美日韩在线观看h| 国产深夜福利视频在线观看| 色吧在线观看| 久久这里有精品视频免费| 18禁在线播放成人免费| 女人久久www免费人成看片| 能在线免费看毛片的网站| 成人毛片60女人毛片免费| 午夜激情福利司机影院| 69精品国产乱码久久久| 精品亚洲成a人片在线观看| 日本猛色少妇xxxxx猛交久久| 久久久欧美国产精品| 91在线精品国自产拍蜜月| 免费av不卡在线播放| 五月天丁香电影| 欧美亚洲 丝袜 人妻 在线| 三上悠亚av全集在线观看| 亚洲精品乱久久久久久| 午夜视频国产福利| 啦啦啦中文免费视频观看日本| 99国产综合亚洲精品| 精品人妻熟女av久视频| 午夜av观看不卡| 极品少妇高潮喷水抽搐| 国产亚洲午夜精品一区二区久久| 日日摸夜夜添夜夜添av毛片| 免费观看无遮挡的男女| 欧美xxⅹ黑人| 国产视频内射| 九九久久精品国产亚洲av麻豆| 久久精品国产a三级三级三级| 国产亚洲欧美精品永久| 久热这里只有精品99| 国产精品女同一区二区软件| 国产欧美亚洲国产| 99九九在线精品视频| 两个人的视频大全免费| 超碰97精品在线观看| 亚洲国产av新网站| 亚洲精品久久成人aⅴ小说 | 久久毛片免费看一区二区三区| 一区二区三区四区激情视频| 国产精品.久久久| 欧美 亚洲 国产 日韩一| 蜜臀久久99精品久久宅男| 亚洲不卡免费看| 九色成人免费人妻av| 亚洲性久久影院| 国模一区二区三区四区视频| 日韩一本色道免费dvd| 纯流量卡能插随身wifi吗| 亚洲国产欧美在线一区| 亚洲精品一二三| 精品卡一卡二卡四卡免费| 精品人妻偷拍中文字幕| 蜜桃国产av成人99| 亚洲精品乱久久久久久| 久久av网站| av电影中文网址| 蜜臀久久99精品久久宅男| 99久久人妻综合| 2018国产大陆天天弄谢| 久久久久久久久久成人| 春色校园在线视频观看| 只有这里有精品99| 涩涩av久久男人的天堂| 成人影院久久| 大陆偷拍与自拍| 丰满少妇做爰视频| 免费观看av网站的网址| 插阴视频在线观看视频| 成人免费观看视频高清| 丰满少妇做爰视频| 午夜福利在线观看免费完整高清在| 中文字幕av电影在线播放| 黄色配什么色好看| 美女福利国产在线| 性高湖久久久久久久久免费观看| 久久99热这里只频精品6学生| 熟女av电影| 国产国拍精品亚洲av在线观看| 免费少妇av软件| xxx大片免费视频| 久久久a久久爽久久v久久| 免费不卡的大黄色大毛片视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产av影院在线观看| 精品亚洲成a人片在线观看| 九九久久精品国产亚洲av麻豆| 18禁观看日本| 人人妻人人澡人人看| 高清在线视频一区二区三区| 色吧在线观看| 中文字幕久久专区| 97在线视频观看| 日日摸夜夜添夜夜添av毛片| 校园人妻丝袜中文字幕| 精品酒店卫生间| 王馨瑶露胸无遮挡在线观看| 色网站视频免费| 一本一本综合久久| 我的老师免费观看完整版| 欧美xxxx性猛交bbbb|