• <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)的影響
    国产高清不卡午夜福利| 亚洲成人久久性| av专区在线播放| 日本三级黄在线观看| 国产中年淑女户外野战色| 亚洲va在线va天堂va国产| 一级二级三级毛片免费看| 国产成年人精品一区二区| 久久中文看片网| 国产精品,欧美在线| 丰满人妻一区二区三区视频av| 国产大屁股一区二区在线视频| 久久久久久大精品| 国产乱人视频| 人妻制服诱惑在线中文字幕| 免费看日本二区| 亚洲自拍偷在线| 国产成人a区在线观看| 久久九九热精品免费| 欧美激情久久久久久爽电影| 国产探花在线观看一区二区| 久久综合国产亚洲精品| 亚洲欧美日韩无卡精品| 成年女人看的毛片在线观看| 黄色日韩在线| 免费一级毛片在线播放高清视频| 欧美成人精品欧美一级黄| 哪里可以看免费的av片| 久久鲁丝午夜福利片| 又黄又爽又刺激的免费视频.| 精品久久久久久久末码| 国产精品一区二区三区四区久久| 麻豆久久精品国产亚洲av| 久久精品国产99精品国产亚洲性色| 嫩草影院新地址| a级毛片a级免费在线| 亚洲精品色激情综合| 成人永久免费在线观看视频| 亚洲精品456在线播放app| 在线观看美女被高潮喷水网站| 亚洲激情五月婷婷啪啪| 免费看光身美女| av视频在线观看入口| 中文字幕av成人在线电影| 晚上一个人看的免费电影| 国内精品宾馆在线| 国产一区二区在线av高清观看| 自拍偷自拍亚洲精品老妇| 熟女电影av网| 久久亚洲国产成人精品v| 九九爱精品视频在线观看| 精品一区二区三区视频在线| 成人午夜高清在线视频| 成人三级黄色视频| 99视频精品全部免费 在线| 国产 一区 欧美 日韩| 亚洲av免费高清在线观看| 日韩成人av中文字幕在线观看| 少妇被粗大猛烈的视频| 亚洲经典国产精华液单| 国产精品一二三区在线看| 男人的好看免费观看在线视频| 长腿黑丝高跟| 国产精品无大码| www.av在线官网国产| 国产日韩欧美在线精品| 99久国产av精品国产电影| 成熟少妇高潮喷水视频| 在线观看免费视频日本深夜| 久久精品国产亚洲av天美| 男女视频在线观看网站免费| 日本-黄色视频高清免费观看| 级片在线观看| 一区二区三区四区激情视频 | 免费电影在线观看免费观看| 亚洲成av人片在线播放无| 91精品国产九色| 精品久久久久久久末码| 亚洲aⅴ乱码一区二区在线播放| 狂野欧美激情性xxxx在线观看| 久久婷婷人人爽人人干人人爱| 日韩三级伦理在线观看| 你懂的网址亚洲精品在线观看 | 国产 一区 欧美 日韩| 99热全是精品| 最近视频中文字幕2019在线8| 色哟哟·www| 国产精品福利在线免费观看| 欧美区成人在线视频| av专区在线播放| 色视频www国产| 成人av在线播放网站| 日韩强制内射视频| 日韩三级伦理在线观看| 国语自产精品视频在线第100页| 能在线免费观看的黄片| 午夜a级毛片| 一本精品99久久精品77| 亚洲七黄色美女视频| 午夜福利在线在线| 欧美bdsm另类| 少妇人妻一区二区三区视频| 男人舔女人下体高潮全视频| av天堂在线播放| 男的添女的下面高潮视频| 亚洲精品亚洲一区二区| 亚洲欧美精品专区久久| 国产av在哪里看| 亚洲av中文av极速乱| 久久国产乱子免费精品| 欧美日韩在线观看h| 看十八女毛片水多多多| 亚洲在久久综合| 99热这里只有精品一区| 桃色一区二区三区在线观看| 一级毛片我不卡| 欧美成人a在线观看| 99九九线精品视频在线观看视频| 淫秽高清视频在线观看| 中国国产av一级| 久久这里只有精品中国| 亚洲欧洲国产日韩| 26uuu在线亚洲综合色| 3wmmmm亚洲av在线观看| 黄片wwwwww| 久久久国产成人精品二区| 99riav亚洲国产免费| 99久国产av精品国产电影| 亚洲七黄色美女视频| 一夜夜www| 99在线人妻在线中文字幕| 99九九线精品视频在线观看视频| 欧美zozozo另类| 一级黄片播放器| 在线免费十八禁| 男人的好看免费观看在线视频| 最后的刺客免费高清国语| 欧美性感艳星| 午夜福利在线观看吧| 久久久久久久午夜电影| 一夜夜www| 日本爱情动作片www.在线观看| 一本久久精品| 18禁在线无遮挡免费观看视频| 一级毛片电影观看 | 国产亚洲精品久久久com| 国产成人影院久久av| 成人无遮挡网站| 国产精品久久视频播放| 午夜亚洲福利在线播放| 欧美zozozo另类| 岛国在线免费视频观看| 亚洲精品久久国产高清桃花| 丝袜美腿在线中文| 超碰av人人做人人爽久久| 日本一本二区三区精品| 中文字幕制服av| 国产精品一二三区在线看| 特大巨黑吊av在线直播| 99久国产av精品| 精品久久久久久成人av| 蜜桃久久精品国产亚洲av| 青春草国产在线视频 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av中文字字幕乱码综合| 国产综合懂色| 啦啦啦韩国在线观看视频| 久久精品国产鲁丝片午夜精品| 亚洲av成人精品一区久久| 男人狂女人下面高潮的视频| 欧美3d第一页| 晚上一个人看的免费电影| 99热全是精品| 国产单亲对白刺激| 欧美性感艳星| 欧美bdsm另类| 久久精品影院6| 最近中文字幕高清免费大全6| 亚洲综合色惰| 国产老妇女一区| 波野结衣二区三区在线| 日韩一区二区视频免费看| 成人午夜精彩视频在线观看| 国产成人午夜福利电影在线观看| 老熟妇乱子伦视频在线观看| 国产大屁股一区二区在线视频| 国产精品久久电影中文字幕| 亚洲国产精品成人久久小说 | 激情 狠狠 欧美| 午夜福利在线在线| 久久精品国产亚洲av香蕉五月| 久久精品国产亚洲av香蕉五月| 99久久中文字幕三级久久日本| 黑人高潮一二区| 久久久久久伊人网av| 色噜噜av男人的天堂激情| 亚洲成人精品中文字幕电影| avwww免费| 亚洲精品久久久久久婷婷小说 | 麻豆成人午夜福利视频| 国产亚洲欧美98| 别揉我奶头 嗯啊视频| 可以在线观看毛片的网站| 国产成人一区二区在线| 国产欧美日韩精品一区二区| 亚洲,欧美,日韩| 人人妻人人澡人人爽人人夜夜 | 国产在线精品亚洲第一网站| av专区在线播放| 熟女电影av网| 国内精品美女久久久久久| 国产av在哪里看| 国产成人精品婷婷| 亚洲av男天堂| 亚洲内射少妇av| 少妇的逼好多水| 国产黄色小视频在线观看| 一级毛片我不卡| 久久午夜福利片| videossex国产| 青春草国产在线视频 | 亚洲天堂国产精品一区在线| 欧美最黄视频在线播放免费| 国产一级毛片七仙女欲春2| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产一区二区在线av高清观看| 黄色视频,在线免费观看| 日韩成人av中文字幕在线观看| 六月丁香七月| 一边摸一边抽搐一进一小说| 18禁裸乳无遮挡免费网站照片| 亚洲国产精品成人久久小说 | 欧美变态另类bdsm刘玥| av专区在线播放| 女的被弄到高潮叫床怎么办| 色视频www国产| 午夜福利成人在线免费观看| 亚洲不卡免费看| 真实男女啪啪啪动态图| 99久久精品一区二区三区| 日韩成人伦理影院| 国产成人午夜福利电影在线观看| 在线观看av片永久免费下载| 国产av一区在线观看免费| 免费一级毛片在线播放高清视频| 不卡一级毛片| 哪个播放器可以免费观看大片| 国产精品精品国产色婷婷| 国产精华一区二区三区| 男人狂女人下面高潮的视频| 能在线免费观看的黄片| 99久久无色码亚洲精品果冻| 97热精品久久久久久| 色播亚洲综合网| 亚洲精品久久国产高清桃花| 在线免费观看不下载黄p国产| 亚洲欧美中文字幕日韩二区| 欧美3d第一页| 国产精品久久电影中文字幕| 亚洲一区高清亚洲精品| 亚洲成人久久爱视频| 99热全是精品| 国产一级毛片在线| 成人二区视频| 国产单亲对白刺激| 少妇熟女欧美另类| 人妻少妇偷人精品九色| 国产又黄又爽又无遮挡在线| 中出人妻视频一区二区| 级片在线观看| 看免费成人av毛片| 欧美日韩一区二区视频在线观看视频在线 | 乱码一卡2卡4卡精品| 久久久久久久久久久丰满| 国产亚洲av片在线观看秒播厂 | 一级毛片电影观看 | 国产淫片久久久久久久久| 欧美成人a在线观看| 亚洲国产精品久久男人天堂| 国产黄片美女视频| 久久中文看片网| 精品久久国产蜜桃| 如何舔出高潮| 欧美三级亚洲精品| 国产成人精品婷婷| 麻豆精品久久久久久蜜桃| 黄色欧美视频在线观看| 亚洲激情五月婷婷啪啪| av视频在线观看入口| 中文字幕人妻熟人妻熟丝袜美| 国产成人精品久久久久久| a级毛片免费高清观看在线播放| 非洲黑人性xxxx精品又粗又长| av在线老鸭窝| 亚洲激情五月婷婷啪啪| 国产精品爽爽va在线观看网站| 最近的中文字幕免费完整| 午夜精品国产一区二区电影 | 午夜福利视频1000在线观看| 欧美潮喷喷水| 九九在线视频观看精品| 最好的美女福利视频网| 男女视频在线观看网站免费| av黄色大香蕉| 看十八女毛片水多多多| 亚洲熟妇中文字幕五十中出| 好男人在线观看高清免费视频| 在线播放无遮挡| 久久久午夜欧美精品| 国产精品精品国产色婷婷| 久久久国产成人免费| 国产午夜精品久久久久久一区二区三区| 伊人久久精品亚洲午夜| 亚洲第一区二区三区不卡| 国产精品无大码| 亚洲国产欧洲综合997久久,| 能在线免费看毛片的网站| 国产视频内射| 大又大粗又爽又黄少妇毛片口| 在线观看66精品国产| 少妇人妻一区二区三区视频| 乱人视频在线观看| 国产老妇伦熟女老妇高清| 国产午夜精品论理片| 亚洲欧洲日产国产| 亚洲国产精品久久男人天堂| 又粗又爽又猛毛片免费看| 两个人视频免费观看高清| 国产乱人偷精品视频| 久久精品国产自在天天线| 成人av在线播放网站| 看免费成人av毛片| 国产精品爽爽va在线观看网站| 别揉我奶头 嗯啊视频| 亚洲乱码一区二区免费版| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久久成人| 哪里可以看免费的av片| 美女cb高潮喷水在线观看| 亚洲在久久综合| 爱豆传媒免费全集在线观看| 日日摸夜夜添夜夜爱| 国产在线男女| 91久久精品国产一区二区三区| 又爽又黄a免费视频| 日日啪夜夜撸| 一区二区三区免费毛片| 中文字幕av在线有码专区| 免费不卡的大黄色大毛片视频在线观看 | 看黄色毛片网站| 亚洲av男天堂| 欧美丝袜亚洲另类| 亚洲在久久综合| 99热只有精品国产| 高清毛片免费看| 成人毛片a级毛片在线播放| 日本爱情动作片www.在线观看| 黄色欧美视频在线观看| 成人国产麻豆网| 禁无遮挡网站| 不卡视频在线观看欧美| 91久久精品国产一区二区三区| 午夜久久久久精精品| 2022亚洲国产成人精品| 久久人人爽人人片av| 国产精品一区二区性色av| 国产高清有码在线观看视频| 国内精品久久久久精免费| 国产成年人精品一区二区| 欧美另类亚洲清纯唯美| 国产亚洲5aaaaa淫片| 一级二级三级毛片免费看| 又爽又黄a免费视频| 天天躁日日操中文字幕| a级毛片a级免费在线| 亚洲欧洲国产日韩| 成人三级黄色视频| 六月丁香七月| 乱码一卡2卡4卡精品| 国产爱豆传媒在线观看| 两个人视频免费观看高清| 国内久久婷婷六月综合欲色啪| 97人妻精品一区二区三区麻豆| a级毛片免费高清观看在线播放| 99热这里只有精品一区| 一本一本综合久久| 精品久久久噜噜| 国产一区二区三区av在线 | 91午夜精品亚洲一区二区三区| 亚洲综合色惰| 三级毛片av免费| 国产精品不卡视频一区二区| 日韩精品有码人妻一区| 性色avwww在线观看| 欧美潮喷喷水| 精品国产三级普通话版| 久99久视频精品免费| 国产探花极品一区二区| 亚州av有码| 男人狂女人下面高潮的视频| 一区二区三区高清视频在线| 熟女人妻精品中文字幕| 国产三级中文精品| 麻豆乱淫一区二区| 日韩欧美精品v在线| 久久草成人影院| 午夜a级毛片| 秋霞在线观看毛片| 久久精品91蜜桃| 免费不卡的大黄色大毛片视频在线观看 | 国产成年人精品一区二区| 欧美色视频一区免费| 色5月婷婷丁香| 久久这里只有精品中国| 毛片女人毛片| 日韩国内少妇激情av| 最新中文字幕久久久久| 国产亚洲精品久久久com| 18+在线观看网站| 欧美日韩在线观看h| 给我免费播放毛片高清在线观看| 国产伦理片在线播放av一区 | 亚洲人成网站在线播| 日韩中字成人| av黄色大香蕉| 久久久久久伊人网av| 亚洲人成网站高清观看| 国产一区二区激情短视频| 久久久欧美国产精品| 在线观看美女被高潮喷水网站| 两性午夜刺激爽爽歪歪视频在线观看| 色播亚洲综合网| 插逼视频在线观看| 午夜精品在线福利| 99国产极品粉嫩在线观看| 欧美xxxx性猛交bbbb| 此物有八面人人有两片| 成人鲁丝片一二三区免费| 精华霜和精华液先用哪个| 99久久人妻综合| 国产精品一及| a级一级毛片免费在线观看| 一级毛片电影观看 | 欧美日韩综合久久久久久| 国产精华一区二区三区| 卡戴珊不雅视频在线播放| 欧美激情国产日韩精品一区| 69av精品久久久久久| .国产精品久久| 给我免费播放毛片高清在线观看| 夜夜看夜夜爽夜夜摸| 麻豆国产av国片精品| 99国产极品粉嫩在线观看| 国产精品伦人一区二区| 我的女老师完整版在线观看| 精品久久久噜噜| 丰满乱子伦码专区| 不卡视频在线观看欧美| 成人无遮挡网站| 麻豆久久精品国产亚洲av| 国产成人aa在线观看| 亚洲真实伦在线观看| 亚洲人成网站在线播| 99久久中文字幕三级久久日本| 成人欧美大片| 一级毛片久久久久久久久女| 亚洲丝袜综合中文字幕| 国产精品一及| 在线免费十八禁| 非洲黑人性xxxx精品又粗又长| 波野结衣二区三区在线| 在线观看美女被高潮喷水网站| 国产高清三级在线| 国产 一区精品| 亚洲欧美中文字幕日韩二区| 日本成人三级电影网站| 天美传媒精品一区二区| 国产精品野战在线观看| 在线观看一区二区三区| 免费看日本二区| 国产精品av视频在线免费观看| av黄色大香蕉| 久久精品综合一区二区三区| 99久久九九国产精品国产免费| 亚洲精品国产av成人精品| 国产精品综合久久久久久久免费| 免费看日本二区| 美女国产视频在线观看| 午夜免费男女啪啪视频观看| 亚洲av第一区精品v没综合| 变态另类成人亚洲欧美熟女| 最近的中文字幕免费完整| 91av网一区二区| 亚洲天堂国产精品一区在线| 噜噜噜噜噜久久久久久91| 成人高潮视频无遮挡免费网站| 99riav亚洲国产免费| 天天躁夜夜躁狠狠久久av| 干丝袜人妻中文字幕| 亚洲av男天堂| 久久久久久久久中文| 日本av手机在线免费观看| 在线天堂最新版资源| 欧美xxxx性猛交bbbb| 精华霜和精华液先用哪个| 又黄又爽又刺激的免费视频.| 欧美性猛交黑人性爽| 欧美不卡视频在线免费观看| 亚洲最大成人手机在线| 亚洲乱码一区二区免费版| 午夜久久久久精精品| 久久精品综合一区二区三区| 国产白丝娇喘喷水9色精品| 男人的好看免费观看在线视频| 欧美+亚洲+日韩+国产| 中文字幕免费在线视频6| 国产高清视频在线观看网站| 伊人久久精品亚洲午夜| 久久国内精品自在自线图片| 人人妻人人澡欧美一区二区| 亚洲图色成人| 久久精品综合一区二区三区| 国产精品1区2区在线观看.| 亚洲国产色片| 在线观看一区二区三区| 99热精品在线国产| 欧美在线一区亚洲| 亚洲成人久久性| 中出人妻视频一区二区| 国产av不卡久久| 欧美+日韩+精品| 久久国内精品自在自线图片| 久久人人精品亚洲av| 最近最新中文字幕大全电影3| 噜噜噜噜噜久久久久久91| 婷婷精品国产亚洲av| 日韩 亚洲 欧美在线| 男人舔女人下体高潮全视频| 天堂av国产一区二区熟女人妻| 国产黄色小视频在线观看| 国产一级毛片七仙女欲春2| 人妻系列 视频| 国产成人91sexporn| 在线国产一区二区在线| videossex国产| av在线观看视频网站免费| 国产在视频线在精品| 中文字幕熟女人妻在线| 精品不卡国产一区二区三区| 韩国av在线不卡| 国产私拍福利视频在线观看| 中文资源天堂在线| 国产亚洲精品久久久com| 中文字幕av成人在线电影| 亚洲人成网站在线播| 国产国拍精品亚洲av在线观看| 熟妇人妻久久中文字幕3abv| 老司机影院成人| 婷婷色综合大香蕉| 看片在线看免费视频| 欧美成人a在线观看| 99国产极品粉嫩在线观看| 天堂影院成人在线观看| av在线天堂中文字幕| 永久网站在线| 蜜桃亚洲精品一区二区三区| 床上黄色一级片| 看免费成人av毛片| 亚洲人与动物交配视频| 91麻豆精品激情在线观看国产| 国产午夜精品论理片| 老司机影院成人| 国产不卡一卡二| 亚洲在线自拍视频| 亚洲七黄色美女视频| 久久热精品热| 欧洲精品卡2卡3卡4卡5卡区| 狂野欧美激情性xxxx在线观看| 夜夜看夜夜爽夜夜摸| 国产精品蜜桃在线观看 | 国产中年淑女户外野战色| 男人舔女人下体高潮全视频| 国产午夜精品一二区理论片| 午夜福利在线观看吧| 村上凉子中文字幕在线| 99热精品在线国产| 午夜福利在线观看吧| 国产一区二区三区av在线 | 狠狠狠狠99中文字幕| 天天躁夜夜躁狠狠久久av| 精品人妻一区二区三区麻豆| 国产黄a三级三级三级人| 天堂av国产一区二区熟女人妻| 校园人妻丝袜中文字幕| 久久热精品热| 老女人水多毛片| 淫秽高清视频在线观看| 婷婷亚洲欧美| 国产老妇伦熟女老妇高清| 成人国产麻豆网| 91狼人影院| 99在线视频只有这里精品首页| 国产精品av视频在线免费观看| 亚洲成av人片在线播放无| 日韩亚洲欧美综合| 国内揄拍国产精品人妻在线| av在线亚洲专区| 好男人在线观看高清免费视频| 99久久精品热视频| 91狼人影院| 乱系列少妇在线播放| 久久精品国产亚洲av香蕉五月| 久久久久久久久久成人| 国产精品久久久久久精品电影小说 | 亚洲国产欧美在线一区|