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

    分子篩限域孔道中吡啶的吸附結(jié)構(gòu)和能量

    2012-12-21 06:32:50褚月英鄭安民
    物理化學(xué)學(xué)報(bào) 2012年2期
    關(guān)鍵詞:鍵長(zhǎng)色散孔道

    韓 冰 褚月英 鄭安民 鄧 風(fēng)

    (1中國(guó)科學(xué)院武漢物理與數(shù)學(xué)研究所,波譜與原子分子物理國(guó)家重點(diǎn)實(shí)驗(yàn)室,武漢核磁共振中心,武漢430071; 2中國(guó)科學(xué)院研究生院,北京100049)

    分子篩限域孔道中吡啶的吸附結(jié)構(gòu)和能量

    韓 冰1,2褚月英1,2鄭安民1,*鄧 風(fēng)1,*

    (1中國(guó)科學(xué)院武漢物理與數(shù)學(xué)研究所,波譜與原子分子物理國(guó)家重點(diǎn)實(shí)驗(yàn)室,武漢核磁共振中心,武漢430071;2中國(guó)科學(xué)院研究生院,北京100049)

    以吸附于ZSM-5孔道中的吡啶分子為例,利用量子化學(xué)理論方法考察了計(jì)算模型和密度泛函方法的選擇對(duì)吡啶吸附結(jié)構(gòu)和吸附能的影響,從而為準(zhǔn)確計(jì)算分子篩限域孔道中客體分子吸附態(tài)結(jié)構(gòu)和能量參數(shù)提供了依據(jù).計(jì)算結(jié)果表明,吡啶吸附能隨著所選用的分子篩的計(jì)算模型(從8T到128T)增大而增大,當(dāng)選用的孔道結(jié)構(gòu)能夠?qū)⒄麄€(gè)分子篩的孔道結(jié)構(gòu)完全包括進(jìn)來(lái)的時(shí)候(72T)達(dá)到收斂.與常規(guī)的密度泛函方法(B3LYP和M06-2X)相比較,考慮到色散作用校正的B97D泛函方法能夠很好地處理分子篩體系中主客體間的長(zhǎng)程相互作用和弱相互作用,計(jì)算得到的能量數(shù)據(jù)與實(shí)驗(yàn)結(jié)果符合得很好.

    量子化學(xué)計(jì)算;分子篩;限域效應(yīng);吸附能

    1 引言

    分子篩催化劑因?yàn)榄h(huán)境友好而在石油化工中被廣泛應(yīng)用.作為固體酸催化劑,它具有不腐蝕設(shè)備、容易與液相反應(yīng)體系分離、環(huán)境污染小、選擇性高等特點(diǎn).1由于分子篩具有獨(dú)特的孔道特征、較好的熱穩(wěn)定性和優(yōu)越的離子交換能力,2因此在裂解、質(zhì)子化、異構(gòu)化和聚合等反應(yīng)中有廣泛的應(yīng)用.3分子篩的酸性特征直接與催化活性相關(guān)聯(lián),研究分子篩的表面酸性對(duì)于深入理解其催化活性及催化機(jī)理具有重要意義.4,5吡啶作為一種堿性探針?lè)肿颖粡V泛地用來(lái)表征固體酸催化劑的酸強(qiáng)度,例如在程序升溫吸脫附(TPD)方法中,通過(guò)比較吡啶脫附溫度的高低來(lái)確定酸強(qiáng)度的大小;而在核磁共振(NMR)探針?lè)肿臃椒ㄖ?通過(guò)比較吡啶吸附前后的1H化學(xué)位移變化來(lái)定量測(cè)量酸強(qiáng)度.6通過(guò)吡啶的TPD和NMR測(cè)定方法,人們可以很直觀地比較催化劑的酸強(qiáng)度大小,而這些譜學(xué)參數(shù)與吡啶在酸位表面的吸附構(gòu)型和吸附熱參數(shù)是密切相關(guān)的.7由此看來(lái),研究吡啶與分子篩酸位的主客體相互作用的結(jié)構(gòu)和能量(吸附熱)參數(shù)對(duì)于在原子分子水平下解釋實(shí)驗(yàn)譜學(xué)參數(shù)顯得特別重要,然而這些重要信息并不能通過(guò)實(shí)驗(yàn)方法來(lái)直接獲取.

    理論計(jì)算的方法作為實(shí)驗(yàn)方法有力的補(bǔ)充,已經(jīng)被廣泛用來(lái)分析固體催化劑活性中心的結(jié)構(gòu)和催化反應(yīng)機(jī)理.彭少逸等8研究了B、Al、Ga等同晶取代進(jìn)入絲光沸石骨架后可能存在的位置,考察了雜原子進(jìn)入骨架對(duì)絲光沸石Br?nsted酸性的影響,孫淮等9研究了磷在P-ZSM-5沸石中存在的形態(tài), Yang等10研究了La離子在分子篩骨架中的存在形態(tài)和Lewis酸性特征,Zhou等11研究了鉬摻雜的ZSM-5分子篩上甲烷脫氫芳構(gòu)化反應(yīng)路徑.在眾多的量子化學(xué)理論方法中,密度泛函理論(DFT)是目前應(yīng)用最廣泛的理論計(jì)算方法之一,該方法不需要很大的計(jì)算資源卻得到不錯(cuò)的精確度,因此被廣泛應(yīng)用到有機(jī)體系熱力學(xué)和動(dòng)力學(xué)的預(yù)測(cè).12,13

    分子篩催化劑由于具有納米量級(jí)的三維孔道結(jié)構(gòu)而在石油化工領(lǐng)域有著廣泛的應(yīng)用,其獨(dú)特產(chǎn)物的選擇性源于客體分子(反應(yīng)物、過(guò)渡態(tài)、產(chǎn)物或反應(yīng)中間體)與分子篩限域孔道中的吸附或排斥作用.孔道限域效應(yīng)主要表現(xiàn)為吸附分子與分子篩孔道之間的長(zhǎng)程主客體相互作用(如范德華作用).眾所周知,經(jīng)典的DFT泛函并不能用來(lái)精確描述范德華作用力以及氫鍵這類(lèi)弱相互作用,理論預(yù)測(cè)的吸附能往往較實(shí)驗(yàn)值偏低.14因此,建立新的理論方法來(lái)準(zhǔn)確描述氫鍵和范德華作用力是當(dāng)前理論研究的熱點(diǎn),其中Becke小組15在1998年開(kāi)發(fā)出混合梯度近似泛函如B97系列,這種泛函對(duì)于研究短程作用為主的體系有較高的計(jì)算精度.后來(lái),Truhlar等16于2006年研發(fā)了M06系列泛函(M06、M06-L、M06-2X),研究結(jié)果表明這一系列泛函能夠更好地描述非鍵相互作用.這是因?yàn)樯⒘ψ鳛榉兜氯A作用力的主要來(lái)源,在研究分子間的相互作用時(shí)非常重要而不能被忽視.可以設(shè)想,如果在理論計(jì)算中能夠?qū)⑸⑾嗷プ饔每紤]進(jìn)來(lái),那么計(jì)算得到的吸附結(jié)構(gòu)和吸附能參數(shù)將更接近于實(shí)驗(yàn)值.Grimme17于2004年提出了DFT-D算法,利用半經(jīng)驗(yàn)的-C6/R6形式的色散項(xiàng)來(lái)表示范德華作用力.在DFT-D算法中,系統(tǒng)的總能量(EDFT-D)表示為DFT計(jì)算得到的能量(EDFT)和長(zhǎng)程相互作用的色散能貢獻(xiàn)(Edisp)之和,即:EDFT-D=EDFT+Edisp.其中長(zhǎng)程相互作用的色散能貢獻(xiàn)Edisp對(duì)于短程作用為零,則回歸為DFT泛函的原始形式,不影響泛函本身在短程相關(guān)效應(yīng)描述較好的原有優(yōu)點(diǎn).DFT-D解決了DFT泛函不能很好地描述弱相互作用及長(zhǎng)程作用中的范德華作用的不足,因而引起了人們?cè)絹?lái)越多的興趣.基于B3LYP而改進(jìn)的B3LYP-D算法,是利用在B3LYP方法優(yōu)化的基礎(chǔ)上,在最后加上色散校正項(xiàng)而得到最終的單點(diǎn)能參數(shù);同樣B97-D方法與B3LYP-D方法類(lèi)似,在單點(diǎn)能部分加上色散作用的校正.兩種方法的優(yōu)勢(shì)還在于占用的計(jì)算資源不大,計(jì)算精度也比較高.而B(niǎo)97D方法則與B97-D方法有所不同,它是在結(jié)構(gòu)優(yōu)化的每一步都將色散相互作用考慮進(jìn)來(lái).18

    準(zhǔn)確描述有機(jī)物在分子篩孔道中的吸附行為的理論研究對(duì)于解釋實(shí)驗(yàn)現(xiàn)象、揭示反應(yīng)機(jī)理是非常重要的.目前B3LYP、M06-2X、B97D等方法廣泛用于研究分子篩的酸性特征,預(yù)測(cè)分子篩催化機(jī)理和過(guò)渡態(tài)特征.19-22然而,哪種泛函方法得到的結(jié)構(gòu)參數(shù)和吸附能與實(shí)驗(yàn)更為接近,文獻(xiàn)一直沒(méi)有相關(guān)的報(bào)導(dǎo).所以有必要通過(guò)系統(tǒng)考察B3LYP、B971、M06-2X以及考慮色散校正作用的B3LYP-D、B971-D和B97D等6種理論方法的計(jì)算結(jié)果,與實(shí)驗(yàn)結(jié)果對(duì)比來(lái)確定能夠準(zhǔn)確預(yù)測(cè)有機(jī)物種在分子篩孔道中吸附結(jié)構(gòu)和能量的泛函方法.本文擬從分子篩計(jì)算模型和理論方法的選擇兩個(gè)方面入手,探索出能夠準(zhǔn)確預(yù)測(cè)有機(jī)分子在分子篩孔道中結(jié)構(gòu)參數(shù)和吸附能量的準(zhǔn)確方法.

    2 計(jì)算方法

    2.1 模型的選取

    為了考察分子篩孔道結(jié)構(gòu)對(duì)Br?nsted酸位和吡啶吸附的影響,分別選取8T、30T、46T、60T、72T和128T團(tuán)簇模型來(lái)表示分子篩ZSM-5的孔道結(jié)構(gòu)(見(jiàn)圖1).需要指出的是,8T模型僅將Br?nsted酸活性中心的局域結(jié)構(gòu)考慮進(jìn)來(lái),隨著模型的逐漸增大,分子篩完整的孔道結(jié)構(gòu)將被重現(xiàn)出來(lái).基于客體分子的易接近性以及酸位的穩(wěn)定性,選取分子篩ZSM-5中Si12-O24(H)-Al12位點(diǎn)來(lái)表示其酸性中心.23在搭建模型中,末端Si原子用H原子來(lái)飽和,H原子的取向與原晶體中被取代的氧原子取向一致,并且所有末端Si-H鍵長(zhǎng)設(shè)置為0.147 nm.

    2.2 DFT計(jì)算方法

    所有的理論計(jì)算都采用Gaussian 09軟件包24完成.隨著模型選擇的增大,體系原子增多(例如128T模型,結(jié)構(gòu)式為Si127Al1O219H75,共422個(gè)原子),對(duì)于如此龐大的體系如果選擇單純的DFT方法因計(jì)算量太大而在現(xiàn)有的計(jì)算水平下顯然無(wú)法實(shí)現(xiàn).為了提高計(jì)算效率,在這里采用分子力學(xué)和量子力學(xué)結(jié)合的ONIOM(our own N-layer integrated MO and MM method)方法,25,26該方法已成功應(yīng)用于研究分子篩的催化反應(yīng)機(jī)理.27在計(jì)算中,將體系分為兩層,其中活性中心采用比較精確的量子力學(xué)方法(比如B3LYP、M06-2X和B97),體系中其余的部分則用較低水平的半經(jīng)驗(yàn)(MNDO)方法加以處理.這樣可以保證較好的計(jì)算精度,同時(shí)也能減小計(jì)算量,提高計(jì)算效率.在30T-128T簇模型體系中,高層為8T區(qū)域(SiO)3Si-OH-Al(OSi)3與吸附的吡啶分子,優(yōu)化分別采用B3LYP、B971、B97D和M062X方法,基組為6-31G(d,p);剩余的部分為低層,采用MNDO半經(jīng)驗(yàn)方法.在結(jié)構(gòu)優(yōu)化過(guò)程中,只將分子篩的活性中心部位(O3Si-OH-AlO3)與吡啶分子進(jìn)行全優(yōu)化,剩下的原子都被固定在原來(lái)的晶格位置,這樣可以保持分子篩模型在結(jié)構(gòu)優(yōu)化的過(guò)程中不失去其特有的晶體結(jié)構(gòu).對(duì)于小的8T ZSM-5結(jié)構(gòu)優(yōu)化采用同樣的密度泛函方法以及基組進(jìn)行計(jì)算,在結(jié)構(gòu)優(yōu)化過(guò)程中對(duì)O3Si-OH-AlO3酸中心和吡啶分子進(jìn)行優(yōu)化,而末端-SiH3基團(tuán)都固定在原來(lái)的晶體位置.單點(diǎn)能計(jì)算分別選用B3LYP、B971、B97D和M06-2X方法,基組為6-31G(d,p).而B(niǎo)971-D和B3LYP-D的單點(diǎn)能則是在B971和B3LYP方法得到的單點(diǎn)能基礎(chǔ)上加上相應(yīng)的Edisp校正.

    圖1 B3LYP/6-31G(d,p)水平下8T-128T ZSM-5模型優(yōu)化得到的Br?nsted酸位結(jié)構(gòu)Fig.1 Optimized structures of Br?nsted acid sites with 8T-128T ZSM-5 models at B3LYP/6-31G(d,p)level (a)8T,(b)30T,(c)46T,(d)60T,(e)72T(f)128T;distances in nm

    2.3 色散相互作用(Edisp)的計(jì)算

    Grimme17和Sumpter28提出長(zhǎng)程相互作用的色散能貢獻(xiàn)可以用公式(1)表示:在這里Nat為體系中的原子數(shù),S6為校正系數(shù),通常不同的泛函對(duì)應(yīng)不同的校正系數(shù),比如PBE的校正系數(shù)為0.75,B3LYP為1.05,B97D為1.25;17C6為原子的色散系數(shù),其大小與原子的電離電勢(shì)和靜態(tài)偶極化率相關(guān);fdamp(Rij)是與原子間距離Rij相關(guān)的阻尼函數(shù),其表達(dá)式可以用公式(2)表示:其中R0為原子的半徑,i和j分別為體系中i和j兩個(gè)原子的原子序號(hào),Rij為i和j兩原子間的距離.表1給出了H、C、O、N、Al和Si原子的色散系數(shù)C6和原子半徑R0等有關(guān)參數(shù).25基于公式(1)和公式(2),色散相互作用(Edisp)的計(jì)算由本課題組編寫(xiě)的C程序完成.

    3 結(jié)果與討論

    表1 色散能計(jì)算中H、C、N、O、Al和Si原子的色散系數(shù)(C6)和原子半徑(R0)25Table 1 Dispersion coefficient(C6)and radii(R0)for elements H,C,N,O,Al,and Si in the calculations of dispersion energy25

    3.1 分子篩酸性特征的理論研究

    分子篩的酸性特征將直接影響到吸附的有機(jī)物種的結(jié)構(gòu)和催化反應(yīng)活性.van Santen29和Blowers30等研究了固體酸酸強(qiáng)度對(duì)烷烴活化反應(yīng)(如氫交換、裂解、脫氫等)活化能壘的影響,發(fā)現(xiàn)催化活性隨著酸強(qiáng)度的增大而顯著增強(qiáng).在前期的研究工作中,我們6,31,32發(fā)現(xiàn)固體酸的酸強(qiáng)度也將影響到探針?lè)肿拥奈胶妥V學(xué)參數(shù),固體酸的酸性越強(qiáng),酸性質(zhì)子向吡啶、三甲基膦氧等堿性探針?lè)肿愚D(zhuǎn)移的程度就越顯著,同時(shí)探針?lè)肿拥幕瘜W(xué)位移也表現(xiàn)出規(guī)律性的變化.為了討論固體酸模型大小對(duì)酸性中心的結(jié)構(gòu)和酸強(qiáng)度的影響,基于8T、30T、46T、60T、72T和128T模型的計(jì)算得到的Br?nsted酸中心的結(jié)構(gòu)參數(shù)和質(zhì)子的Mulliken電荷(QH)列于表2中.可以發(fā)現(xiàn)不論采用哪種泛函方法,隨著模型的增大結(jié)構(gòu)參數(shù)呈現(xiàn)出相同的變化規(guī)律.以B3LYP方法為例,隨著模型的增大,30T、46T模型的Si-O鍵長(zhǎng)分別從8T模型的0.1662 nm減小為0.1658和0.1654 nm,而Al-O鍵長(zhǎng)則在0.1833 nm的基礎(chǔ)上分別減小了0.0006和0.0013 nm,說(shuō)明模型的增大(從8T到46T)對(duì)酸位的局域骨架結(jié)構(gòu)有一定的影響;而46T增大到128T模型的結(jié)構(gòu)參數(shù)變化則不明顯,說(shuō)明對(duì)于46T以上模型,其結(jié)構(gòu)參數(shù)已經(jīng)達(dá)到收斂,模型的進(jìn)一步增大對(duì)骨架結(jié)構(gòu)的影響很小.同樣,對(duì)于羥基O-H鍵長(zhǎng),8T模型為0.0969 nm,30T和46T分別為0.0971和0.0973 nm,46T以上模型O-H鍵長(zhǎng)維持在0.0973 nm不變;Q(H)即H的電荷從0.388|e|逐漸增加到0.401|e|后保持不變,同樣在46T呈現(xiàn)收斂的趨勢(shì).

    表2 不同計(jì)算方法下8T-128T ZSM-5模型去質(zhì)子化能(DPE)、H的電荷(QH)和主要結(jié)構(gòu)參數(shù)Table 2 Deprotonation energy(DPE),Mulliken charge for H(QH),and main geometry parameters of 8T-128T ZSM-5 models calculated with different methods

    分子篩的催化反應(yīng)活性是與其酸性特征密切相關(guān)的,通過(guò)rO-H和Br?nsted質(zhì)子QH的變化可以定性地觀察酸強(qiáng)度的變化趨勢(shì).從表2給出的rO-H和QH的數(shù)據(jù)來(lái)看,分子篩的酸強(qiáng)度隨著選用模型的增加而緩慢增強(qiáng).除此之外,去質(zhì)子化能(DPE)是一種有效表征酸強(qiáng)度大小的定量方法.33DPE定義為酸位脫去一個(gè)質(zhì)子所要吸收的能量(AH→H++A-),在數(shù)值上表現(xiàn)為酸位(AH)的能量與失去質(zhì)子后形成共軛堿(A-)的能量差值:DPE=E(A-)-E(AH).較小的DPE代表較強(qiáng)的酸性,在酸催化反應(yīng)中質(zhì)子轉(zhuǎn)移就較容易.從表2中發(fā)現(xiàn)8T、30T、46T、60T、72T和128T的DPE計(jì)算值(由B3LYP方法計(jì)算得到)分別是1257.5、1256.2、1235.7、1235.6、1235.6和1235.6 kJ· mol-1,呈先降低后收斂的趨勢(shì);不同的固體酸模型的DPE值與ZSM-5分子篩的實(shí)驗(yàn)值范圍1217.8-1255.5 kJ·mol-134接近.如表2所示,這種變化趨勢(shì)在使用其他方法計(jì)算的結(jié)果中也有類(lèi)似的體現(xiàn).這說(shuō)明分子篩的酸強(qiáng)度主要取決于酸性中心的局域結(jié)構(gòu),而孔道結(jié)構(gòu)的影響不大.

    3.2 不同計(jì)算模型研究吡啶在分子篩孔道中的吸附

    圖2 B3LYP/6-31G(d,p)方法下8T-128T ZSM-5模型吸附吡啶后優(yōu)化得到的結(jié)構(gòu)Fig.2 B3LYP/6-31G(d,p)optimized structures of pyridine adsorption complexes on 8T-128T ZSM-5 models (a)8T,(b)30T,(c)46T,(d)60T,(e)72T,(f)128T;distance in nm

    分子篩孔道和被吸附分子之間的主客體相互作用將直接影響到吸附分子的構(gòu)型以及其物理化學(xué)性質(zhì).被吸附分子與分子篩孔道之間主要是靜電吸引和范德華作用,當(dāng)被吸附分子逐漸被孔道包圍時(shí),這種穩(wěn)定作用也逐漸增強(qiáng),即“孔道限域效應(yīng)”.在研究了不同大小的分子篩模型對(duì)其固有酸強(qiáng)度的影響之后,我們將進(jìn)一步研究孔道效應(yīng)對(duì)吡啶結(jié)構(gòu)的影響.圖2給出了B3LYP計(jì)算方法下ZSM-5不同酸模型的吡啶吸附結(jié)構(gòu)圖.比較圖1和圖2中分子篩和吡啶的構(gòu)型變化后可以發(fā)現(xiàn),在吡啶吸附前后rO-H鍵長(zhǎng)的變化是非常明顯的.以72T模型為例,吡啶吸附前,O-H鍵長(zhǎng)為0.0973 nm,吸附吡啶后, O-H鍵長(zhǎng)增大為0.1559 nm,而N-H鍵長(zhǎng)為0.1045 nm,這表明吡啶被質(zhì)子化形成了吡啶離子.相應(yīng)的,吡啶吸附前后Br?nsted酸位上的O原子與Si和Al的距離也呈現(xiàn)一定的趨勢(shì):對(duì)于8T模型, O-Al和O-Si鍵長(zhǎng)分別為0.1833和0.1662 nm;在吡啶探針?lè)肿幼饔孟?O-Al和O-Si鍵長(zhǎng)分別減小為0.1749和0.1594 nm;對(duì)于72T模型,吡啶吸附前后O-Al和O-Si鍵長(zhǎng)分別從0.1820和0.1653 nm減小為0.1744和0.1589 nm.表3給出了在不同模型和方法下計(jì)算得到質(zhì)子的Mulliken電荷(QH)和結(jié)構(gòu)參數(shù).研究結(jié)果表明不論采用哪種泛函方法,隨著模型的增大,Eads、QH以及結(jié)構(gòu)參數(shù)都呈現(xiàn)出相同的變化規(guī)律.以B3LYP方法為例,隨著酸模型的增大,N-H鍵長(zhǎng)逐漸從0.1075 nm減小到0.1045 nm,而后者與氣態(tài)吡啶離子的N-H鍵長(zhǎng)0.1017 nm是比較接近的.另外通過(guò)比較不同模型間的參數(shù)變化發(fā)現(xiàn),8T和30T模型間的參數(shù)變化為ΔrN-H=0.0026 nm,ΔrO-H=0.0005 nm;30T和46T之間變化很小ΔrN-H=0.0003 nm,ΔrO-H=0.0002 nm,而60T到128T之間則基本沒(méi)有變化.值得注意的是,在8T模型中, QH為0.420|e|,而在72T模型中QH減小為0.400|e|.這說(shuō)明隨著酸模型的增大,越來(lái)越多的正電荷被分散到吡啶分子的共軛環(huán)中,增大了吡啶離子的穩(wěn)定性.這也進(jìn)一步說(shuō)明了孔道結(jié)構(gòu)對(duì)吡啶吸附起到了穩(wěn)定的作用.不同方法計(jì)算得到的吡啶吸附參數(shù)有一定的差別,說(shuō)明不同的計(jì)算方法對(duì)吡啶在分子篩表面吸附的構(gòu)型是有影響的;但從8T到128T模型參數(shù)的變化趨勢(shì)上看,無(wú)論哪一種泛函方法,其整體變化趨勢(shì)都是一致的.

    表3 不同計(jì)算方法下8T-128T ZSM-5模型吸附吡啶后H的電荷和主要結(jié)構(gòu)參數(shù)Table 3 Mulliken charge for H and main geometry parameters of pyridine adsorption complexes on 8T-128T ZSM-5 models calculated with different methods

    表4給出了吡啶在不同的酸模型表面上的吸附能.吸附能計(jì)算公式如下:

    ΔEads=(EZOH+Epyridine)-Epyridine-ZOH(3)其中Epyridine-ZOH為復(fù)合物的能量,EZOH和Epyridine分別為分子篩模型和吡啶分子的單點(diǎn)能.在這里需要提出的是,計(jì)算得到的吸附能是體系在絕對(duì)零度的相互作用能,其值與實(shí)驗(yàn)得到的吸附熱接近.35,36表4給出的吡啶吸附能結(jié)果分析發(fā)現(xiàn),隨著模型的增大,吸附能呈增大趨勢(shì)并最后達(dá)到收斂.對(duì)于B3LYP方法,8T模型的吸附能為45.4 kJ·mol-1,它包含了酸位的局域結(jié)構(gòu),但沒(méi)考慮孔道;30T和46T模型吸附能分別為85.4和88.9 kJ·mol-1,與8T模型相比吸附能分別升高了40和43.5 kJ·mol-1,這是由于30T和46T模型既包含了酸位的局域結(jié)構(gòu),又包含了部分孔道;60T、72T和128T模型的吸附能分別為90.0、90.4和90.6 kJ·mol-1,模型之間能量只有0.6 kJ· mol-1的變化,說(shuō)明孔道結(jié)構(gòu)已經(jīng)逐漸完整,對(duì)于吡啶的吸附影響變小,能量達(dá)到收斂.這說(shuō)明對(duì)于復(fù)雜的ZSM-5分子篩體系,72T模型可以完全將孔道對(duì)吸附分子之間的相互作用都包括進(jìn)來(lái),可以用來(lái)反映ZSM-5分子篩真實(shí)的孔道結(jié)構(gòu).對(duì)于其他方法計(jì)算得到的吡啶吸附能(表4)也同樣表現(xiàn)出隨著計(jì)算模型的增大而增大,在72T模型時(shí)能量達(dá)到收斂.

    3.3 不同泛函方法研究吡啶在分子篩孔道中的吸附

    基于優(yōu)化得到的吡啶分子吸附結(jié)構(gòu),我們?cè)诒?中給出了不同模型和計(jì)算方法下吡啶在ZSM-5分子篩表面的吸附能.在前面的研究中我們已經(jīng)確定了72T模型可以用來(lái)反映分子篩ZSM-5的孔道結(jié)構(gòu),所以在這里以72T模型為例比較不同計(jì)算方法得到的吸附能.在B3LYP/6-31G(d,p)水平下計(jì)算得到的吸附能為90.4 kJ·mol-1,而實(shí)驗(yàn)值為200 kJ· mol-1.37眾所周知對(duì)于質(zhì)子化的吡啶離子與分子篩孔道之間的相互作用,除了電子相關(guān)效應(yīng)外,靜電吸引和范德華作用也起到穩(wěn)定吡啶離子的作用. B3LYP方法僅僅能對(duì)電子相關(guān)效應(yīng)作出較好的描述,而不能夠很好地描述范德華作用力,因此計(jì)算得到的吸附能與實(shí)驗(yàn)值誤差很大.B971/6-31G(d,p)水平下計(jì)算得到的吸附能為110.4 kJ·mol-1,與B3LYP方法相比計(jì)算精度稍有提高,這是因?yàn)锽97方法在電子相關(guān)方面進(jìn)行了一定的修正,對(duì)于研究短程作用為主的體系有較高的計(jì)算精度,但由于對(duì)范德華作用力的描述依舊不是很全面,所以計(jì)算得到的吸附能與實(shí)驗(yàn)值仍有很大的誤差.M06-2X方法考慮了弱相互作用力,在計(jì)算范德華作用方面性能比B3LYP有所改善,計(jì)算得到的吸附能量為144.3 kJ·mol-1,與B3LYP相比計(jì)算精確明顯提高.19但對(duì)于主客體間普遍存在的長(zhǎng)程色散作用,M06-2X未能體現(xiàn)出色散作用對(duì)吸附能量的影響,對(duì)于分子間的色散作用描述不是很充分,計(jì)算結(jié)果與實(shí)驗(yàn)值仍有約30%的誤差.B3LYP-D和B971-D分別是在B3LYP和B971優(yōu)化后對(duì)色散作用進(jìn)行校正,因此得到的能量比B3LYP和B971方法計(jì)算的結(jié)果進(jìn)一步提高,分別為122.4和122.5 kJ·mol-1,但這兩種校正方法的最大問(wèn)題在于僅僅在單點(diǎn)能計(jì)算部分增加了一項(xiàng)色散校正,并沒(méi)有將弱相互作用的影響引入結(jié)構(gòu)優(yōu)化過(guò)程中,雖然相對(duì)于純的密度泛函計(jì)算得到的能量有所改進(jìn),但是與實(shí)驗(yàn)值相差仍然約80 kJ·mol-1.B97D方法是在B97的基礎(chǔ)上,每一步結(jié)構(gòu)優(yōu)化的過(guò)程中都進(jìn)行色散校正,因此可以預(yù)期該種方法得到的結(jié)構(gòu)和能量精度將會(huì)有很大的提高.從表3給出的不同方法下吡啶吸附的結(jié)構(gòu)參數(shù)可以發(fā)現(xiàn),B97D方法計(jì)算得到的N-H和O-H鍵長(zhǎng)分別為0.1047和0.1563 nm,比B971方法計(jì)算得到的N-H鍵長(zhǎng)增加了0.0002 nm,O-H鍵長(zhǎng)減小了0.0004 nm;基于B97D水平下得到的結(jié)構(gòu)參數(shù),我們進(jìn)一步計(jì)算得到其吸附能為204.7 kJ·mol-1,與實(shí)驗(yàn)值200 kJ·mol-1吻合得較好,這也進(jìn)一步驗(yàn)證了B97D方法在計(jì)算色散作用力方面有較大的優(yōu)勢(shì).38圖3給出了不同團(tuán)簇模型和不同方法下吡啶的吸附能以及實(shí)驗(yàn)值,從圖中可以更直觀地看出B97D方法與實(shí)驗(yàn)值吻合得最好.

    表4 不同計(jì)算方法下得到的8T-128T ZSM-5模型吡啶吸附能(Eads)Table 4 Pyridine adsorption energies(Eads)on 8T-128T ZSM-5 models calculated with different methods

    3.4 分子篩體系理論計(jì)算模型和方法的選擇

    分子篩催化劑的優(yōu)越性主要體現(xiàn)在其具有一定的酸性和獨(dú)特的孔道特征,因此在選擇計(jì)算模型的時(shí)候需要充分考慮這兩個(gè)因素.通過(guò)系統(tǒng)地研究發(fā)現(xiàn)分子篩的酸強(qiáng)度主要取決于酸性中心的局域結(jié)構(gòu),孔道結(jié)構(gòu)對(duì)分子篩的固有酸強(qiáng)度影響不大.從圖2給出的吡啶吸附結(jié)構(gòu)圖和表4給出的吡啶吸附能結(jié)果可以發(fā)現(xiàn),隨著模型選擇的增大,孔道結(jié)構(gòu)逐漸完整,“孔道限域效應(yīng)”越來(lái)越明顯,吸附能呈增大趨勢(shì)并在60T-128T達(dá)到收斂.另外需要注意的是,隨著模型的增大,體系原子增多,計(jì)算量也相應(yīng)增大.基于以上因素考慮,我們選用72T模型來(lái)進(jìn)行分子篩體系的理論計(jì)算,因?yàn)?2T模型既能表現(xiàn)出ZSM-5的酸性特征,又具備完整的孔道結(jié)構(gòu),可以充分地反映ZSM-5的催化性能;另外72T模型原子數(shù)量適中,計(jì)算量比128T小,計(jì)算精度比小模型高,因此選用72T模型來(lái)代表ZSM-5分子篩的孔道結(jié)構(gòu)應(yīng)用于一些理論計(jì)算可以得到很好的結(jié)果.

    圖3 六種計(jì)算方法下8T-128T ZSM-5模型吡啶吸附能Fig.3 Adsorption energies of pyridine on 8T-128T ZSM-5 models calculated with six different methods The experimental value(200 kJ·mol-1)37is shown by the dotted line.

    分子篩體系中的主客體間相互作用非常復(fù)雜,除了電子相關(guān)效應(yīng)外,靜電吸引、范德華作用以及色散作用等長(zhǎng)程相互作用對(duì)于有機(jī)物種在分子篩孔道中的吸附都有不同程度的影響.因此選擇合適的計(jì)算方法用于準(zhǔn)確描述有機(jī)物種在分子篩孔道中的吸附結(jié)構(gòu)和能量,對(duì)于解釋實(shí)驗(yàn)現(xiàn)象、揭示反應(yīng)機(jī)理是非常重要的.通過(guò)研究不同泛函方法計(jì)算吡啶在分子篩孔道中吸附的結(jié)果發(fā)現(xiàn),B97D方法可以很好地處理范德華作用等弱相互作用,在計(jì)算色散作用力方面也有較大的優(yōu)勢(shì),是一種能夠準(zhǔn)確預(yù)測(cè)有機(jī)分子在分子篩孔道中結(jié)果參數(shù)和吸附能量的準(zhǔn)確方法.

    4 結(jié)論

    通過(guò)量子化學(xué)計(jì)算的方法系統(tǒng)考察了計(jì)算中分子篩模型和泛函方法對(duì)分子篩固有酸強(qiáng)度與吡啶吸附結(jié)構(gòu)和能量的影響.結(jié)果表明,分子篩的酸強(qiáng)度主要取決于酸性中心的局域結(jié)構(gòu).而在研究反應(yīng)物在分子篩限域孔道中吸附和轉(zhuǎn)化過(guò)程中,為了考察吸附分子和分子篩孔道之間的主客體相互作用,有必要采用大的模型將分子篩完整的孔道結(jié)構(gòu)包括進(jìn)來(lái).由于經(jīng)典的B3LYP泛函方法未能將色散作用力考慮進(jìn)來(lái),采用該方法計(jì)算的吸附能量明顯低于實(shí)驗(yàn)值.為了獲得精確的結(jié)構(gòu)和能量數(shù)據(jù), B97D泛函方法在結(jié)構(gòu)優(yōu)化的每一步都作出色散作用力的校正,計(jì)算結(jié)果與實(shí)驗(yàn)值能夠很好的吻合.References

    (1) Weitkamp,J.;Traa,Y.Catal.Today 1999,49,193.

    (2) Slagtern,A.;Dahl,I.M.;Jens,K.J.;Myrstad,T.App.Catal.AGen.2010,375,213.

    (3)Luzgin,M.V.;Parmon,V.N.Angew.Chem.Int.Edit.2008,47, 4559.

    (4)Yu,Z.W.;Zheng,A.M.;Wang,Q.;Huang,S.J.;Deng,F.;Liu, S.B.Chin.J.Magn.Reson.2010,27,485.[喻志武,鄭安民,王 強(qiáng),黃信炅,鄧 風(fēng),劉尚斌.波譜學(xué)雜志,2010,27,485.]

    (5) Zheng,A.M.;Huang,S.J.;Deng,F.;Liu,S.B.Phys.Chem. Chem.Phys.2011,13,14889.

    (6)Zheng,A.M.;Zhang,H.L.;Chen,L.;Yue,Y.;Ye,C.H.;Deng, F.J.Phys.Chem.B 2007,111,3085.

    (7) Coma,A.Chem.Rev.1995,95,559.

    (8)Yuan,S.P.;Wang,J.G.;Li,Y.W.;Peng,S.Y.Acta Physico-Chimica Sinica 2001,17,811.[袁淑萍,王建國(guó),李永旺,彭少逸.物理化學(xué)學(xué)報(bào),2001,17,811.]

    (9)Yang,J.;Sun,Y.X.;Zhao,L.F.;Sun,H.Acta Physico-Chimica Sinica 2011,27,1823. [楊 靜,孫迎新,趙立峰,孫 淮.物理化學(xué)學(xué)報(bào),2011,27,1823.]

    (10)Yang,G.;Wang,Y.;Zhou,D.H.;Zhuang,J.Q.;Liu,X.C.;Han, X.W.;Bao,X.H.J.Chem.Phys.2003,119,9765.

    (11)Zhou,D.H.;Ma,D.;Liu,X.C.;Bao,X.H.J.Chem.Phys. 2001,114,9125.

    (12) Hohenberg,P.;Kohn,W.Phys.Rev.1964,136,864.

    (13)Kohn,W.;Becke,A.D.;Parr,R.G.J.Phys.Chem.1996,100, 12974.

    (14) Kristyan,S.;Pulay,P.Chem.Phys.Lett.1994,229,175.

    (15)Schmider,H.L.;Becke,A.D.J.Chem.Phys.1998,108,9624.

    (16) Zhao,Y.;Schultz,N.E.;Truhlar,D.G.J.Chem.Theory Comput.2006,2,364.

    (17)Grimme,S.J.Comput.Chem.2004,25,1463.

    (18)Grimme,S.J.Comput.Chem.2006,27,1787.

    (19) Zhao,Y.;Truhlar,D.G.J.Phys.Chem.C 2008,112,6860.

    (20) Pidko,E.A.;Hensen,E.J.M.;van Santen,R.A.J.Phys.Chem. C 2008,112,19604.

    (21) Boronat,M.;Martinez,C.;Corma,A.Phys.Chem.Chem.Phys. 2011,13,2603.

    (22) Boekfa,B.;Choomwattana,S.;Khongpracha,P.;Limtrakul,J. Langmuir 2009,22,12990.

    (23) Vankoningsveld,H.;Van Bekkum,H.;Jansen,J.C.Acta Crystallogr.B 1987,43,127.

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

    (25) Maseras,F.;Morokuma,K.J.Comput.Chem.1995,16,1170.

    (26) Humbel,S.;Siebe,R.S.;Morokuma,K.J.Chem.Phys.1996, 105,1959.

    (27) Lesthaeghe,D.;Speybroeck,V.V.;Marin,G.B.;Waroquier,M. Chem.Phys.Lett.2006,417,309.

    (28) Sumpter,B.G.J.Chem.Theory Comput.2010,6,727.

    (29) Rigby,A.M.;Kramer,G.J.;van Santen,R.A.J.Catal.1997, 170,1.

    (30) Zheng,X.;Blowers,P.J.Phys.Chem.A 2006,110,2455.

    (31) Zheng,A.M.;Zhang,H.L.;Lu,X.;Liu,S.B.;Deng,F.J.Phys. Chem.B 2008,112,4496.

    (32)Zheng,A.M.;Huang,S.;Chen,W.;Wu,P.;Zhang,H.;Lee,H.; Ménorval,L.;Deng,F.;Liu,S.B.J.Phys.Chem.A 2008,112, 7337.

    (33) Brand,H.V.;Curtiss,L.A.;Iton,L.E.J.Phys.Chem.1993,97, 12773.

    (34) Datka,J.;Boczar,M.;Rymarowicz,P.J.Catal.1988,114,368.

    (35) Dunne,J.A.;Rao,M.;Sircar,S.;Corte,R.J.;Myers,A.L. Langmuir 1996,12,5896.

    (36) Savitz,S.;Siperstein,F.;Rorte,R.J.;Myers,A.L.J.Phys. Chem.B 1998,102,6865.

    (37) Lee,C.;Parrillo,D.J.;Gorte,R.J.;Farneth,W.E.J.Am.Chem. Soc.1996,118,3262.

    (38) Fuchs,A.H.;Adamo,C.J.Phys.Chem.Lett.2010,1,763.

    October 9,2011;Revised:November 17,2011;Published on Web:November 23,2011.

    Adsorption Structure and Energy of Pyridine Confined inside Zeolite Pores

    HAN Bing1,2CHU Yue-Ying1,2ZHENGAn-Min1,*DENG Feng1,*
    (1Wuhan Center for Magnetic Resonance,State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics,Chinese Academy of Sciences,Wuhan 430071,P.R.China;2Graduate University of Chinese Academy of Sciences,Beijing 100049,P.R.China)

    The performance of different exchange-correlation functionals for the description of the interaction of pyridine with different cluster models of ZSM-5 zeolite has been assessed.Theoretical calculations show that upon increasing the cluster model from 8T to 128T,the adsorption energy of pyridine in ZSM-5 zeolite increases gradually and reaches convergence by the 72T cluster model.On the basis of the 72T cluster model,the pyridine adsorption energy calculated with different functionals is further examined.Compared to the conventional functionals(B3LYP and M06-2X),the B97D functional which takes into account the dispersion correction provides calculated results that agree well with experimental data.The present results indicate that the B97D functional is suitable for studying long-range interactions in weakly interacting systems.

    Quantum chemical calculation;Zeolite;Confinement effect;Adsorption energy

    10.3866/PKU.WHXB201111232

    *Corresponding authors.ZHENGAn-Min,Email:zhenganm@wipm.ac.cn;Tel:+86-27-87197359.DENG Feng,Email:dengf@wipm.ac.cn;

    Tel:+86-27-87198820.

    The project was supported by the National Natural Science Foundation of China(21073228,20933009,20921004)and National Key Basic Research Program of China(973)(2009CB918600).

    國(guó)家自然科學(xué)基金(21073228,20933009,20921004)及國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃項(xiàng)目(973)(2009CB918600)資助

    O641

    猜你喜歡
    鍵長(zhǎng)色散孔道
    陰離子調(diào)控錳基鈣鈦礦中Mn─O的鍵長(zhǎng)和磁性
    “光的折射”“光的色散”知識(shí)鞏固
    “光的折射”“光的色散”知識(shí)鞏固
    “光的折射”“光的色散”知識(shí)鞏固
    基于ANSYS的液壓集成塊內(nèi)部孔道受力分析
    接觸壓力非均勻分布下彎曲孔道摩阻損失分析
    『光的折射』『光的色散』隨堂練
    密度泛函理論研究鎘的二鹵化合物分子的結(jié)構(gòu)和振動(dòng)頻率
    淺議鍵能與鍵長(zhǎng)的關(guān)系
    離子對(duì)SBA-15形貌與孔道結(jié)構(gòu)的影響
    少妇人妻久久综合中文| 男女边吃奶边做爰视频| 精品一区二区免费观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲综合色网址| 一级片'在线观看视频| 亚洲欧洲精品一区二区精品久久久 | 欧美xxxx性猛交bbbb| 精品亚洲成国产av| 插逼视频在线观看| 99久久中文字幕三级久久日本| 国精品久久久久久国模美| 韩国av在线不卡| 男女无遮挡免费网站观看| 亚洲人成77777在线视频| 欧美日韩视频精品一区| 最近中文字幕高清免费大全6| 99热全是精品| 日韩在线高清观看一区二区三区| 午夜日本视频在线| 有码 亚洲区| 久久国内精品自在自线图片| 久久久国产欧美日韩av| 波野结衣二区三区在线| 日本vs欧美在线观看视频| 啦啦啦视频在线资源免费观看| 免费看av在线观看网站| 久久久久精品人妻al黑| 成人黄色视频免费在线看| 91久久精品国产一区二区三区| 热99国产精品久久久久久7| 伊人久久国产一区二区| 夫妻性生交免费视频一级片| av免费在线看不卡| 高清欧美精品videossex| 十八禁高潮呻吟视频| videosex国产| 亚洲综合色惰| 9191精品国产免费久久| 午夜福利视频在线观看免费| 插逼视频在线观看| 2018国产大陆天天弄谢| 少妇人妻精品综合一区二区| 国产伦理片在线播放av一区| 日本av免费视频播放| 青青草视频在线视频观看| 欧美成人午夜精品| 久久鲁丝午夜福利片| 日本与韩国留学比较| 丰满迷人的少妇在线观看| 欧美成人午夜精品| 捣出白浆h1v1| 亚洲国产精品专区欧美| 国产亚洲精品久久久com| 国产日韩欧美亚洲二区| 亚洲欧美一区二区三区黑人 | 最后的刺客免费高清国语| 老熟女久久久| 亚洲成人av在线免费| 视频区图区小说| 啦啦啦视频在线资源免费观看| 亚洲精品国产av蜜桃| 免费大片黄手机在线观看| 国产熟女午夜一区二区三区| 精品少妇内射三级| 涩涩av久久男人的天堂| 成人毛片60女人毛片免费| 精品久久蜜臀av无| 国产精品久久久久成人av| 丝袜喷水一区| 91午夜精品亚洲一区二区三区| 777米奇影视久久| 天美传媒精品一区二区| 欧美日韩视频高清一区二区三区二| 看免费成人av毛片| 狠狠婷婷综合久久久久久88av| 欧美精品一区二区大全| 男女边摸边吃奶| 日韩av在线免费看完整版不卡| 哪个播放器可以免费观看大片| 欧美精品av麻豆av| 日韩精品有码人妻一区| 成人18禁高潮啪啪吃奶动态图| 久久韩国三级中文字幕| 国产精品久久久久久精品古装| 国产成人午夜福利电影在线观看| 欧美精品人与动牲交sv欧美| 99热国产这里只有精品6| 国产视频首页在线观看| 搡老乐熟女国产| 夜夜骑夜夜射夜夜干| 麻豆精品久久久久久蜜桃| 女人被躁到高潮嗷嗷叫费观| 午夜免费观看性视频| 国产毛片在线视频| 国产一区二区激情短视频 | 婷婷色综合www| 久久久精品区二区三区| 亚洲国产毛片av蜜桃av| 在线精品无人区一区二区三| 纵有疾风起免费观看全集完整版| 一区在线观看完整版| 亚洲综合色网址| 久久韩国三级中文字幕| 在线观看www视频免费| 国产色婷婷99| 欧美日韩综合久久久久久| 草草在线视频免费看| 久久ye,这里只有精品| 少妇的逼好多水| 丰满饥渴人妻一区二区三| 国产片特级美女逼逼视频| 久久人人97超碰香蕉20202| 免费看不卡的av| 黑丝袜美女国产一区| 91国产中文字幕| 亚洲av福利一区| 大片免费播放器 马上看| 国产乱来视频区| 天天影视国产精品| 成人亚洲精品一区在线观看| 蜜臀久久99精品久久宅男| 成人国产av品久久久| 亚洲色图 男人天堂 中文字幕 | 亚洲欧美成人精品一区二区| 欧美精品一区二区免费开放| 精品国产国语对白av| 国产亚洲av片在线观看秒播厂| 午夜福利在线观看免费完整高清在| 久久人人爽人人爽人人片va| 欧美3d第一页| 国产成人精品在线电影| 高清av免费在线| 黄网站色视频无遮挡免费观看| 亚洲在久久综合| 又大又黄又爽视频免费| 欧美精品亚洲一区二区| 久久鲁丝午夜福利片| 久久久久久伊人网av| 男人爽女人下面视频在线观看| 亚洲欧美成人精品一区二区| 桃花免费在线播放| 久久人人爽人人片av| 免费大片黄手机在线观看| 国产成人免费观看mmmm| 秋霞在线观看毛片| 成人亚洲精品一区在线观看| 一级毛片电影观看| 亚洲美女搞黄在线观看| 久久99一区二区三区| 91久久精品国产一区二区三区| 欧美成人精品欧美一级黄| 精品国产露脸久久av麻豆| 国产国拍精品亚洲av在线观看| 亚洲欧美一区二区三区国产| 波野结衣二区三区在线| 在线观看国产h片| 内地一区二区视频在线| 纵有疾风起免费观看全集完整版| 一二三四中文在线观看免费高清| 夜夜爽夜夜爽视频| 欧美精品一区二区免费开放| 在现免费观看毛片| 深夜精品福利| 啦啦啦在线观看免费高清www| 久久人人爽人人片av| 岛国毛片在线播放| 久久久久视频综合| 欧美丝袜亚洲另类| 国产免费一级a男人的天堂| 国产高清国产精品国产三级| 久久久国产欧美日韩av| 免费高清在线观看视频在线观看| 美国免费a级毛片| 少妇被粗大的猛进出69影院 | 51国产日韩欧美| 99re6热这里在线精品视频| 美女国产视频在线观看| 国产精品一二三区在线看| 免费高清在线观看日韩| 亚洲欧洲精品一区二区精品久久久 | 熟女电影av网| 亚洲精品色激情综合| 免费日韩欧美在线观看| 欧美激情国产日韩精品一区| 午夜激情久久久久久久| 国产深夜福利视频在线观看| 免费人妻精品一区二区三区视频| 高清欧美精品videossex| 18禁国产床啪视频网站| 色婷婷av一区二区三区视频| 亚洲欧美中文字幕日韩二区| 成人亚洲精品一区在线观看| 久久热在线av| 乱人伦中国视频| 亚洲av欧美aⅴ国产| 在线观看美女被高潮喷水网站| 男女下面插进去视频免费观看 | 日韩熟女老妇一区二区性免费视频| 少妇的丰满在线观看| 国产亚洲一区二区精品| 满18在线观看网站| 国产精品久久久久久av不卡| 国产精品不卡视频一区二区| 我的女老师完整版在线观看| 亚洲精品456在线播放app| 女性生殖器流出的白浆| 精品一区二区三区四区五区乱码 | 热99久久久久精品小说推荐| 一区二区三区精品91| 好男人视频免费观看在线| 丝瓜视频免费看黄片| 欧美激情极品国产一区二区三区 | 不卡视频在线观看欧美| 男女高潮啪啪啪动态图| 99久久人妻综合| 精品人妻偷拍中文字幕| 国产1区2区3区精品| 秋霞伦理黄片| 各种免费的搞黄视频| 欧美成人午夜精品| 在线天堂最新版资源| 国产视频首页在线观看| 精品一区二区三卡| av免费观看日本| 美国免费a级毛片| 尾随美女入室| 一边亲一边摸免费视频| 精品亚洲成国产av| 又大又黄又爽视频免费| 久久久久国产精品人妻一区二区| 精品国产露脸久久av麻豆| 两个人看的免费小视频| 成年动漫av网址| 人人妻人人澡人人看| 国产 精品1| 女性生殖器流出的白浆| 色视频在线一区二区三区| 夜夜爽夜夜爽视频| 我的女老师完整版在线观看| 狂野欧美激情性bbbbbb| 国产精品国产av在线观看| 日韩人妻精品一区2区三区| 一级片免费观看大全| 亚洲av男天堂| 婷婷色综合www| 午夜福利视频在线观看免费| 90打野战视频偷拍视频| 亚洲美女视频黄频| av卡一久久| 99热6这里只有精品| 少妇的逼水好多| 日产精品乱码卡一卡2卡三| 青春草亚洲视频在线观看| 久久人人爽人人片av| 亚洲精品一二三| 街头女战士在线观看网站| 国产精品一区二区在线不卡| 麻豆精品久久久久久蜜桃| 久久久久久久久久人人人人人人| 亚洲av电影在线观看一区二区三区| av在线app专区| 免费人妻精品一区二区三区视频| 中文字幕精品免费在线观看视频 | 午夜91福利影院| 另类亚洲欧美激情| 欧美精品国产亚洲| 一个人免费看片子| 亚洲欧美成人综合另类久久久| 99热6这里只有精品| 免费观看av网站的网址| 啦啦啦视频在线资源免费观看| 国产又色又爽无遮挡免| 国产激情久久老熟女| 亚洲图色成人| 69精品国产乱码久久久| 久久人人爽人人片av| 亚洲欧美精品自产自拍| videossex国产| 亚洲久久久国产精品| 2018国产大陆天天弄谢| 日本欧美视频一区| 日日撸夜夜添| 久久精品国产鲁丝片午夜精品| 只有这里有精品99| 秋霞在线观看毛片| 亚洲精品乱码久久久久久按摩| 高清欧美精品videossex| 伦理电影大哥的女人| 精品一区在线观看国产| 久久久久久久亚洲中文字幕| 色吧在线观看| 最近2019中文字幕mv第一页| 欧美人与善性xxx| 伊人亚洲综合成人网| 午夜免费男女啪啪视频观看| 汤姆久久久久久久影院中文字幕| 国产免费一区二区三区四区乱码| 另类亚洲欧美激情| 久久精品久久久久久噜噜老黄| 色婷婷久久久亚洲欧美| 久久久久精品人妻al黑| 国产成人91sexporn| 97人妻天天添夜夜摸| 亚洲成av片中文字幕在线观看 | 国产毛片在线视频| 国产成人精品婷婷| 欧美 日韩 精品 国产| a级毛片在线看网站| 国产成人免费无遮挡视频| 一个人免费看片子| 男的添女的下面高潮视频| av国产久精品久网站免费入址| 精品一品国产午夜福利视频| 最新中文字幕久久久久| 婷婷色麻豆天堂久久| 国产成人a∨麻豆精品| 老女人水多毛片| av播播在线观看一区| 在现免费观看毛片| 久久鲁丝午夜福利片| 色哟哟·www| 久久精品国产亚洲av天美| 欧美精品一区二区大全| 看免费av毛片| 精品一区二区三卡| 只有这里有精品99| 欧美人与性动交α欧美软件 | 国产国拍精品亚洲av在线观看| 久久久久久久大尺度免费视频| 欧美xxⅹ黑人| 久久久久人妻精品一区果冻| av线在线观看网站| 日本黄大片高清| 两个人看的免费小视频| 老女人水多毛片| 中文字幕人妻熟女乱码| 丰满少妇做爰视频| 国产不卡av网站在线观看| 桃花免费在线播放| 超碰97精品在线观看| 国产探花极品一区二区| 久久久久视频综合| 婷婷成人精品国产| 亚洲国产精品999| 亚洲国产精品专区欧美| 51国产日韩欧美| 国产av码专区亚洲av| 亚洲国产精品一区三区| 中文字幕亚洲精品专区| 中文天堂在线官网| 久久 成人 亚洲| 亚洲精华国产精华液的使用体验| 亚洲国产av新网站| 亚洲精品成人av观看孕妇| 90打野战视频偷拍视频| 亚洲精品国产色婷婷电影| 狠狠精品人妻久久久久久综合| 国产成人精品久久久久久| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品国产精品| 汤姆久久久久久久影院中文字幕| 精品国产国语对白av| 另类亚洲欧美激情| 亚洲美女视频黄频| 51国产日韩欧美| 免费av不卡在线播放| 亚洲欧洲精品一区二区精品久久久 | 91成人精品电影| 国产精品国产三级专区第一集| 精品久久久精品久久久| 人人妻人人爽人人添夜夜欢视频| 日韩三级伦理在线观看| 香蕉精品网在线| 欧美精品一区二区大全| 51国产日韩欧美| 亚洲精品乱码久久久久久按摩| 欧美xxxx性猛交bbbb| 1024视频免费在线观看| 久久精品久久久久久久性| 黄片无遮挡物在线观看| 美女视频免费永久观看网站| 99香蕉大伊视频| 久热久热在线精品观看| 免费在线观看黄色视频的| 如何舔出高潮| 菩萨蛮人人尽说江南好唐韦庄| 伊人久久国产一区二区| av又黄又爽大尺度在线免费看| 亚洲欧美色中文字幕在线| 在线观看三级黄色| 赤兔流量卡办理| 青春草亚洲视频在线观看| 高清视频免费观看一区二区| 看免费av毛片| 国产精品一区二区在线观看99| 精品国产国语对白av| 国产av码专区亚洲av| 欧美 亚洲 国产 日韩一| 在线看a的网站| 在线 av 中文字幕| 母亲3免费完整高清在线观看 | 欧美亚洲日本最大视频资源| 中文字幕人妻熟女乱码| 十八禁高潮呻吟视频| 18禁裸乳无遮挡动漫免费视频| 在线观看一区二区三区激情| 日韩av免费高清视频| 午夜视频国产福利| 国产成人av激情在线播放| 久久久久精品人妻al黑| 精品国产国语对白av| 国产精品久久久久久久电影| 99久国产av精品国产电影| 男女免费视频国产| 99热6这里只有精品| 久久99蜜桃精品久久| 蜜桃国产av成人99| 国产一区二区在线观看av| 老女人水多毛片| 久久人人爽人人爽人人片va| 亚洲欧美成人精品一区二区| 中文字幕av电影在线播放| 亚洲人成77777在线视频| 春色校园在线视频观看| 纵有疾风起免费观看全集完整版| 日韩一区二区三区影片| 久久精品久久久久久噜噜老黄| 亚洲伊人色综图| 黑人高潮一二区| 精品一区二区三卡| 中国国产av一级| 国产一区二区三区综合在线观看 | 最近最新中文字幕大全免费视频 | 久久这里只有精品19| 精品人妻一区二区三区麻豆| 大片电影免费在线观看免费| 久久亚洲国产成人精品v| 综合色丁香网| 90打野战视频偷拍视频| 免费日韩欧美在线观看| 蜜臀久久99精品久久宅男| 如何舔出高潮| 亚洲一级一片aⅴ在线观看| 久久国产精品男人的天堂亚洲 | 免费大片黄手机在线观看| 王馨瑶露胸无遮挡在线观看| 国产女主播在线喷水免费视频网站| 亚洲精品第二区| 一区在线观看完整版| 乱人伦中国视频| 亚洲av福利一区| 18+在线观看网站| 最近中文字幕高清免费大全6| 伦理电影大哥的女人| 又黄又粗又硬又大视频| 美女国产高潮福利片在线看| 亚洲av中文av极速乱| 亚洲国产欧美日韩在线播放| 成人国产av品久久久| 啦啦啦啦在线视频资源| 五月伊人婷婷丁香| 国产日韩欧美亚洲二区| 国产精品欧美亚洲77777| 亚洲成av片中文字幕在线观看 | 水蜜桃什么品种好| 国产熟女午夜一区二区三区| 久久久a久久爽久久v久久| 一边亲一边摸免费视频| 在线精品无人区一区二区三| 最后的刺客免费高清国语| 男的添女的下面高潮视频| 在线观看www视频免费| 精品国产一区二区久久| 夜夜爽夜夜爽视频| 黄片播放在线免费| 最近中文字幕2019免费版| 在线 av 中文字幕| 香蕉国产在线看| 十八禁高潮呻吟视频| 精品久久国产蜜桃| 成年人免费黄色播放视频| av电影中文网址| 亚洲精品乱码久久久久久按摩| 欧美 日韩 精品 国产| 成人国语在线视频| 赤兔流量卡办理| 亚洲精品av麻豆狂野| 亚洲第一区二区三区不卡| 狂野欧美激情性bbbbbb| 免费在线观看黄色视频的| 国产成人精品一,二区| 丰满迷人的少妇在线观看| videosex国产| 午夜激情av网站| 美女主播在线视频| 久久久a久久爽久久v久久| 亚洲精品国产色婷婷电影| 国产一区二区激情短视频 | 一级毛片黄色毛片免费观看视频| 美女国产高潮福利片在线看| 亚洲天堂av无毛| 少妇精品久久久久久久| 9色porny在线观看| 一级毛片黄色毛片免费观看视频| 三级国产精品片| 人人妻人人澡人人看| 久久精品熟女亚洲av麻豆精品| 两个人看的免费小视频| 香蕉丝袜av| 青春草国产在线视频| 人人妻人人澡人人看| 精品视频人人做人人爽| 国精品久久久久久国模美| 在线看a的网站| 欧美亚洲日本最大视频资源| 国产一区亚洲一区在线观看| a级片在线免费高清观看视频| 一二三四中文在线观看免费高清| 高清黄色对白视频在线免费看| 在线免费观看不下载黄p国产| 26uuu在线亚洲综合色| 亚洲国产看品久久| 校园人妻丝袜中文字幕| 日本黄大片高清| 亚洲精品中文字幕在线视频| 国产一区亚洲一区在线观看| 十八禁网站网址无遮挡| 老司机影院成人| 日本午夜av视频| 欧美97在线视频| 观看美女的网站| 日本-黄色视频高清免费观看| 国产69精品久久久久777片| www.色视频.com| 赤兔流量卡办理| 国产免费一区二区三区四区乱码| 国产黄色免费在线视频| 亚洲少妇的诱惑av| 亚洲欧洲精品一区二区精品久久久 | 一区二区三区乱码不卡18| 男女边摸边吃奶| 亚洲国产精品国产精品| 久久影院123| 香蕉丝袜av| 婷婷色av中文字幕| 少妇的丰满在线观看| 99re6热这里在线精品视频| 亚洲四区av| 观看av在线不卡| 国产色婷婷99| 国产免费视频播放在线视频| 国产免费福利视频在线观看| 亚洲国产色片| 91午夜精品亚洲一区二区三区| 欧美bdsm另类| 亚洲人与动物交配视频| 美女大奶头黄色视频| 丝袜脚勾引网站| 久久99一区二区三区| 亚洲综合精品二区| 国产熟女欧美一区二区| 一区二区三区乱码不卡18| 夫妻性生交免费视频一级片| 久久青草综合色| 岛国毛片在线播放| 嫩草影院入口| 全区人妻精品视频| a级毛片黄视频| 久久99精品国语久久久| 欧美日韩亚洲高清精品| 亚洲美女视频黄频| tube8黄色片| 黑人猛操日本美女一级片| 人人妻人人爽人人添夜夜欢视频| 最近中文字幕高清免费大全6| 免费高清在线观看日韩| 老司机亚洲免费影院| 亚洲av中文av极速乱| 美女视频免费永久观看网站| 草草在线视频免费看| 午夜影院在线不卡| 一级a做视频免费观看| 午夜福利,免费看| 毛片一级片免费看久久久久| 免费女性裸体啪啪无遮挡网站| av女优亚洲男人天堂| av有码第一页| 9热在线视频观看99| 亚洲国产精品专区欧美| 亚洲,欧美,日韩| 亚洲欧美中文字幕日韩二区| www.色视频.com| 99re6热这里在线精品视频| 人妻少妇偷人精品九色| 国产麻豆69| 人人妻人人爽人人添夜夜欢视频| 曰老女人黄片| 少妇人妻 视频| 亚洲性久久影院| 丝袜脚勾引网站| av免费在线看不卡| 久久久久国产网址| 欧美精品亚洲一区二区| 十分钟在线观看高清视频www| 亚洲国产日韩一区二区| 少妇人妻精品综合一区二区| 日韩三级伦理在线观看| 女性生殖器流出的白浆| 久久99热6这里只有精品| 啦啦啦在线观看免费高清www| 人人妻人人澡人人爽人人夜夜| 国产精品一二三区在线看| 高清av免费在线| 九九在线视频观看精品|