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

    In-Au(111)和Ir-Au(111)合金表面的性質(zhì)及其對(duì)巴豆醛的吸附比較

    2016-12-29 08:20:16蔣軍輝錢夢(mèng)丹薛繼龍夏盛杰倪哲明邵蒙蒙
    物理化學(xué)學(xué)報(bào) 2016年12期
    關(guān)鍵詞:巴豆鍵長(zhǎng)構(gòu)型

    蔣軍輝 錢夢(mèng)丹 薛繼龍 夏盛杰 倪哲明 邵蒙蒙

    (浙江工業(yè)大學(xué)化學(xué)工程學(xué)院,杭州310014)

    In-Au(111)和Ir-Au(111)合金表面的性質(zhì)及其對(duì)巴豆醛的吸附比較

    蔣軍輝 錢夢(mèng)丹 薛繼龍 夏盛杰 倪哲明*邵蒙蒙

    (浙江工業(yè)大學(xué)化學(xué)工程學(xué)院,杭州310014)

    采用密度泛函理論研究了M(M=In,Ir)原子修飾的M-Au(111)合金表面的穩(wěn)定性,并選其最優(yōu)模型探討了合金表面的活性及其對(duì)巴豆醛的吸附。合金的幾何構(gòu)型、形成能和結(jié)合能等性質(zhì)表明,In-Au(111)面的穩(wěn)定性隨In原子的間距增大而提高,Ir-Au(111)面的穩(wěn)定性隨Ir原子的間距增大而降低。對(duì)于巴豆醛在MAu(111)面上的吸附,當(dāng)其通過(guò)C=O吸附于合金表面的TopM位時(shí),吸附能最大,吸附構(gòu)型最穩(wěn)定。從巴豆醛的結(jié)構(gòu)變化、態(tài)密度、差分電荷密度以及Mulliken電荷布居等分析可以看出,穩(wěn)定吸附構(gòu)型的巴豆醛分子形變較大,電荷轉(zhuǎn)移明顯。其中,位于-7.04 eV至費(fèi)米能級(jí)處的p、d軌道雜化,對(duì)體系的吸附具有重要貢獻(xiàn)。分析比較In-Au(111)面與Ir-Au(111)面,發(fā)現(xiàn)后者的配體效應(yīng)更佳,不僅具有更高的穩(wěn)定性和活性,而且對(duì)于巴豆醛具有更強(qiáng)的吸附力。此外,相比于改性前的Au(111)面,M原子的修飾明顯提升了金屬表面的穩(wěn)定性及吸附能力。

    In-Au(111)面;Ir-Au(111)面;修飾;巴豆醛;吸附;密度泛函理論

    1 引言

    1987年,Haruta等1發(fā)現(xiàn)負(fù)載在氧化物上的納米Au對(duì)低溫下的CO氧化反應(yīng)具有很好的催化活性,自此打破了Au催化惰性的觀念,打開了Au催化應(yīng)用的新篇章。Bus等2研究了Au/Al2O3以及Pt/ Al2O3催化劑對(duì)肉桂醛的催化加氫反應(yīng),發(fā)現(xiàn)Au/ Al2O3的催化活性明顯高于Pt/Al2O3。Bailie和Hutching3通過(guò)Au/ZrO2催化劑首次完成了巴豆醛的C=O選擇性氫化。蔣軍輝等4探討了巴豆醛在Au (111)上的選擇性催化加氫機(jī)理,結(jié)果表明Au(111)面對(duì)于巴豆醛的C=O具有較好的選擇性氫化能力。雖然純Au催化劑具有較好的低溫氧化能力以及較高的選擇性,但是純Au催化劑存在活性低且易燒結(jié)等缺陷5,6。因此,對(duì)于Au催化劑的改良,將有助于其未來(lái)的發(fā)展和應(yīng)用范圍的開闊。

    目前,實(shí)驗(yàn)和理論研究上改善單金屬催化劑的方法有載體和前驅(qū)體的替換7、催化劑活性粒子的改變以及第二元素的修飾8,9。其中第二元素的修飾可優(yōu)化催化劑的局域電子效應(yīng)和幾何構(gòu)型,生成雙金屬混合位,從而增強(qiáng)催化劑的活性,進(jìn)而提高其催化性能。張連陽(yáng)等10探討了Au修飾的Au/ Pd(111)面的穩(wěn)定性,同時(shí)研究了噻吩在其表面的吸附和加氫脫硫機(jī)理,發(fā)現(xiàn)改性后的催化劑具有更高的穩(wěn)定性以及催化活性。Cárdenas-Lizana等11研究了Au/Al2O3以及PdAu/Al2O3催化劑在氯代硝基苯加氫中的催化性能,發(fā)現(xiàn)Au/Al2O3催化劑的選擇性較好,但活性較低。Pd修飾后的PdAu/Al2O3催化劑的催化活性明顯提高,且對(duì)氯代硝基苯的選擇性為100%。

    巴豆醛(CAL)是α,β-不飽和醛的典型代表之一,其一系列加氫產(chǎn)物作為重要的化工原料在精細(xì)化學(xué)品的合成當(dāng)中具有廣泛應(yīng)用12-14。例如,C=O選擇性氫化得到的巴豆醇是醫(yī)藥及香精的中間體,C=C選擇性氫化得到的丁醛是樹脂以及殺蟲劑的中間體。近年來(lái),第二元素修飾的Au催化劑被普遍應(yīng)用于α,β-不飽和醛的加氫反應(yīng)中。Yang等15探究了SBA-15負(fù)載的Au和In-Au對(duì)于CAL的催化加氫,結(jié)果顯示In-Au/SBA-15催化劑的C=O選擇性催化能力更強(qiáng),不僅對(duì)于巴豆醇的選擇性和產(chǎn)率分別高達(dá)78%和71%,而且催化劑的循環(huán)使用并不影響其選擇性。He等16計(jì)算表明,In對(duì)于Au(110)面的修飾,顯著增強(qiáng)了其與丙烯醛的相互作用,同時(shí)也提高了對(duì)于丙烯醇的選擇性。Rojas等17分析了Ir-Au/TiO2的性質(zhì)及其對(duì)于檸檬醛的催化加氫,一系列的表征表明Ir-Au/TiO2催化劑具有較好的穩(wěn)定性,并且相對(duì)于Au/TiO2,合金在檸檬醛的選擇性催化加氫中具有更高的催化活性。Zhao等18研究表明,對(duì)于肉桂醛的催化加氫反應(yīng),Ir修飾的Ir-Au催化劑相比較于Au催化劑,既有較高的C=O選擇性,又有較強(qiáng)的H2活化能力。

    綜上,在Au催化劑中摻雜M(M=In,Ir)原子所形成的修飾M-Au合金,具有優(yōu)異的α,β-不飽和醛選擇性催化加氫性能。然而,有關(guān)M-Au合金的性質(zhì)及其對(duì)于CAL吸附的理論研究較少。因此,本研究采用密度泛函理論(DFT)構(gòu)建了M原子修飾的M-Au(111)面,對(duì)可能的合金模型的穩(wěn)定性和活性進(jìn)行了分析,同時(shí)研究了CAL在M-Au(111)面上的吸附行為及其優(yōu)勢(shì)吸附構(gòu)型的電子性質(zhì),以期為M-Au合金催化劑的設(shè)計(jì)及CAL在其表面的選擇性加氫提供理論依據(jù)。

    2 計(jì)算方法和模型

    本研究的DFT計(jì)算利用Materials Studio 5.5軟件包中的DMol3程序包進(jìn)行19,20。為了減小范德華力引起的吸附誤差,本研究運(yùn)用了DFT-D的方法。電子交換相關(guān)泛函基于廣義梯度積分(GGA)的Perdew-Wang-91(PW91)泛函21,價(jià)電子波函數(shù)采用雙數(shù)值極化函數(shù)(DNP)展開,且未限制電子自旋。計(jì)算精度設(shè)為medium,smearing設(shè)為13.10 kJ·mol-1,k點(diǎn)設(shè)為3×3×1。結(jié)構(gòu)優(yōu)化合理性的判據(jù)為:能量差異小于4.80×10-2kJ·mol-1,原子位移小于5.00×10-4nm以及各個(gè)原子上的力小于9.70 kJ·mol-1·nm-1。

    本文首先選用低覆蓋率(1/16)p(4×4)的周期性四層平板模型來(lái)模擬Au(111)面,且將真空層厚度選為1.20 nm,避免了平板間發(fā)生鏡像作用。在此基礎(chǔ)上,通過(guò)2個(gè)修飾原子取代Au(111)表面的2個(gè)Au原子后構(gòu)成M-Au(111)面18,22,并根據(jù)修飾原子取代位置的相鄰距離不同,構(gòu)建了如圖1所示由近及遠(yuǎn)的3種模型分別記為(a)、(b)、(c)。計(jì)算過(guò)程中考慮到表面弛豫的影響,對(duì)底下2層原子進(jìn)行了固定,上2層原子可自由移動(dòng)。圖1同時(shí)列出了CAL的結(jié)構(gòu)。

    3 結(jié)果與討論

    3.1 M-Au(111)面性質(zhì)分析

    圖1 M-Au(111)面的不同模型和CAL的結(jié)構(gòu)Fig.1 Different models of M-Au(111)surfaces and the structure of CAL

    3.1.1 M-Au(111)面穩(wěn)定性分析

    M-Au(111)面是以Au為溶劑,在其晶格點(diǎn)陣中置換溶入溶質(zhì)原子M所形成的置換合金表面,其中原子尺寸因素(ΔR)大小是影響其形成的關(guān)鍵23?,F(xiàn)根據(jù)計(jì)算模型將ΔR定義為:ΔR=|RMRAu|/RAu,式中:RM為溶質(zhì)原子M的原子半徑(RIn= 0.166 nm,RIr=0.136 nm),RAu為溶劑原子Au的原子半徑(RAu=0.146 nm)。當(dāng)ΔR<0.150時(shí),溶質(zhì)M與溶劑Au之間可以形成穩(wěn)定的合金,且ΔR越小,引起的晶格畸變?cè)叫 1?列出了M-Au(111)面的ΔR,由表1可見,對(duì)于In-Au(111)面和Ir-Au(111)面,其ΔR都小于0.150,分別為0.137和0.068。因此,In和Ir原子都能夠較好地與Au面形成合金,并且Ir原子置換所引起的晶格畸變較In原子小。從而可以推測(cè),Ir-Au(111)面較In-Au(111)面的穩(wěn)定性更高。

    表1同時(shí)列出了優(yōu)化后合金的M―Au鍵平均鍵長(zhǎng)d以及置換原子M在M-Au(111)面上的凸起高度h,進(jìn)而從幾何構(gòu)型角度展現(xiàn)了模型的變化情況。在In-Au合金中,模型(a)、(b)、(c)的In―Au鍵平均鍵長(zhǎng)分別為0.300、0.299以及0.297 nm,并且In原子分別凸出表面0.092、0.091以及0.089 nm。這可能是由于RIn大于RAu,以及Au、In原子晶體結(jié)構(gòu)的差異(Au屬于立方晶系,In屬于四方晶系),使得In―Au鍵平均鍵長(zhǎng)大于Au―Au鍵(0.288 nm),從而引起修飾表面凸起。同時(shí)由于模型(c)的幾何形變最小,因此可以推測(cè)In以(c)模型修飾Au(111)面的可能性最大。在Ir-Au合金中,模型(a)、(b)、(c)的Ir―Au鍵平均鍵長(zhǎng)分別為0.283、0.287、0.287 nm,并且Ir原子較表面分別下凹0.036、0.033、0.033 nm。引起表面下凹的可能原因有兩個(gè),其一是RIr小于RAu,其二是Ir和Au之間的電負(fù)性差值大于Au和Au之間的電負(fù)性差值,引起的Ir―Au的作用力大于Au―Au的作用力。對(duì)比分析三種模型可得,在Ir-Au(111)面中模型(a)的穩(wěn)定性最強(qiáng),模型(c)的穩(wěn)定性最弱。綜上,相比于In原子,Ir原子的結(jié)構(gòu)性質(zhì)更加接近Au原子,從而使得Ir-Au(111)面的穩(wěn)定性較In-Au(111)面的穩(wěn)定性更高。

    為了進(jìn)一步研究M-Au合金的穩(wěn)定性,本文還探討了合金的形成能(ΔEf)和結(jié)合能(ΔEc)等熱力學(xué)性質(zhì)。其中ΔEf定義為:ΔEf=EAu(111)+2EM-(EM-Au(111)+2EAu),式中EM-Au(111)和EAu(111)分別為M-Au (111)面和Au(111)面的能量,EAu和EM分別為孤立Au和M原子的能量。ΔEc定義為:ΔEc=(∑Ei-EM-Au(111)/Au(111))/64,式中Ei為孤立原子(Au,In,Ir)的能量,EM-Au(111)/Au(111)為M-Au(111)面或Au(111)面的能量。表1顯示In-Au(111)面不同模型的形成能隨著表面In原子間距的增加近似遞增,而對(duì)于Ir-Au (111)面恰好相反。根據(jù)形成能的定義,形成能越大,體系越穩(wěn)定。In-Au(111)面的模型穩(wěn)定性順序?yàn)?c)>(b)>(a),Ir-Au(111)面的模型穩(wěn)定性順序?yàn)?a)>(b)>(c),且Ir-Au(111)面的穩(wěn)定性高于In-Au (111)面。對(duì)于結(jié)合能,In-Au(111)面和Ir-Au(111)面的模型(a)、(b)、(c)能量分別為 391.15、391.54、391.71以及425.78、425.52、425.44 kJ·mol-1,且都較Au(111)面的能量(265.66 kJ·mol-1)大。根據(jù)結(jié)合能定義,結(jié)合能是由原子自由態(tài)形成化合物所釋放的能量。因此,Ir-Au(111)面的穩(wěn)定性高于In-Au(111)面,同時(shí)都比Au(111)面的穩(wěn)定性要好。In-Au(111)面的構(gòu)型穩(wěn)定性順序?yàn)?c)>(b)>(a),Ir-Au(111)面的構(gòu)型穩(wěn)定性順序?yàn)?a)>(b)>(c),此結(jié)論與形成能以及幾何結(jié)構(gòu)分析結(jié)果一致。

    表1 M-Au(111)面不同模型的平均凸起高度、M―Au鍵平均鍵長(zhǎng)、形成能以及結(jié)合能Table 1 Geometrical parameters of the M-Au(111)surfaces,including the average embossment height of doped atoms, average bond length of M―Au,formation and cohesive energies

    綜上所述,In原子和Ir原子都能夠較好地通過(guò)置換取代的方法與Au(111)面形成相應(yīng)的穩(wěn)定合金,并且Ir-Au(111)面的穩(wěn)定性高于In-Au(111)面。In-Au(111)面和Ir-Au(111)面上的最穩(wěn)定模型分別為(c)和(a),因此這兩類模型存在的可能性最大。下文研究將以In-Au(111)面的模型(c)以及Ir-Au(111)面的模型(a)為基礎(chǔ),對(duì)其進(jìn)行化學(xué)活性以及吸附性能的探討。圖2給出了其優(yōu)化后的M-Au (111)面平板模型。

    3.1.2 M-Au(111)面化學(xué)活性分析

    d帶中心可以較好地反映出金屬表面的化學(xué)活性,且其與費(fèi)米能級(jí)的間距越小,化學(xué)活性越高24。對(duì)于M-Au(111)面,d帶電子對(duì)總態(tài)密度的貢獻(xiàn)最大。第二金屬的修飾使得合金表面的電荷密度發(fā)生變化,從而導(dǎo)致d帶電子的變化,進(jìn)而可用d帶中心的大小來(lái)預(yù)測(cè)和驗(yàn)證不同金屬表面的吸附能力。計(jì)算表明,In-Au(111)面和Ir-Au(111)面的d帶中心分別位于-3.54和-3.48 eV,都比Au(111)面(-3.63 eV)更靠近費(fèi)米能級(jí)。因此In和Ir原子對(duì)Au(111)面的修飾,使得合金的d帶中心向費(fèi)米能級(jí)遷移,并且Ir-Au(111)面的d帶中心向右遷移的距離較In-Au(111)面更遠(yuǎn),這可能是由于Ir和Au原子間的相互作用大于In原子,即Ir原子修飾引起的配體效應(yīng)大于In原子。據(jù)此可以推測(cè)金屬表面吸附能力的大小為Ir-Au(111)面>In-Au(111)面>Au(111)面。另外,Mulliken電荷布居也是一種基于軌道分布定量描述體系電荷分布情況的分析方法25。In和Ir原子加入后,金屬表面上Au原子的電荷由-0.01e分別降低為-0.02e和-0.06e,而In和Ir原子上的電荷分別增加到0.18e和0.47e。所以可以推斷,在M-Au(111)面上,吸附將圍繞M原子展開,且在Ir原子上的吸附力大于In原子,此結(jié)論與d帶中心預(yù)測(cè)結(jié)果一致。

    3.2 M-Au(111)面上CAL的吸附

    3.2.1 吸附構(gòu)型及其吸附能

    CAL催化加氫過(guò)程首先是CAL吸附于催化劑M-Au(111)面,進(jìn)而CAL與解離的H原子發(fā)生選擇性氫化反應(yīng),最終得到的產(chǎn)物脫離M-Au(111)表面,因此研究CAL在合金催化劑表面的吸附是探究催化加氫的關(guān)鍵步驟。根據(jù)本課題組之前的研究以及相關(guān)文獻(xiàn)報(bào)道4,26,CAL存在單點(diǎn)(O,C=O及C=C)和雙點(diǎn)(O協(xié)同C=C或C=O協(xié)同C=C)兩種吸附點(diǎn)位。在In-Au(111)面上,CAL可通過(guò)單點(diǎn)吸附于TopAu、TopIn、BriAu-Au、BriAu-In、Hcp以及Fcc六種吸附位置,共計(jì)18種吸附構(gòu)型,通過(guò)雙點(diǎn)吸附于上述六種吸附位置,共計(jì)36種吸附構(gòu)型,總計(jì)優(yōu)化了54種吸附構(gòu)型。在Ir-Au(111)面上,CAL可通過(guò)單點(diǎn)吸附于 TopAu、TopIr、BriAu-Au、BriAu-Ir、BriIr-Ir、Hcp以及Fcc七種吸附位置,共計(jì)21種吸附構(gòu)型,通過(guò)雙點(diǎn)吸附于上述七種吸附位置,共計(jì)49種吸附構(gòu)型,總計(jì)優(yōu)化了70種吸附構(gòu)型。為了確定其吸附的程度,本研究引入了吸附能(Eads)一值,并定義為:Eads=ECAL+EM-Au(111)-ECAL/M-Au(111),式中ECAL/M-Au(111)表示CAL吸附在M-Au (111)面時(shí)體系的能量。Eads越大表示吸附體系越穩(wěn)定。

    圖2 M-Au(111)面穩(wěn)定平板模型的俯視和側(cè)視圖Fig.2 Top and side views of the stable M-Au(111)plate models

    表2和表3分別列出了優(yōu)化后CAL在In-Au (111)面和Ir-Au(111)面上能夠穩(wěn)定存在的吸附構(gòu)型及其吸附能值和O原子距離表面的高度(dO-layer)。優(yōu)化后CAL在In-Au(111)面和Ir-Au(111)面上分別只存在17種和22種穩(wěn)定吸附構(gòu)型。引起吸附構(gòu)型減少的可能原因是M原子的修飾使得M-Au(111)面對(duì)于CAL具有較大的吸附力,從而引起構(gòu)型的普遍變化。同時(shí),M-Au(111)面上不均勻分布的電荷對(duì)于吸附構(gòu)型的減少也有一定影響。特別是對(duì)于雙點(diǎn)吸附構(gòu)型,分別縮減到9種和12種。這可能是由于甲基引起的空間阻礙效應(yīng)較M-Au(111)面對(duì)于CAL的C=C吸附力更強(qiáng),使得雙點(diǎn)吸附較難存在。對(duì)于CAL/In-Au(111)體系,其吸附能在79.88-145.45 kJ·mol-1之間,O原子距離表面為0.242-0.595 nm之間。當(dāng)CAL以C=O吸附于TopIn位時(shí),即以η2-diσCO(TopIn)構(gòu)型存在時(shí),其吸附能最大,為145.45 kJ·mol-1,O原子距離表面為0.242 nm。且吸附前后,C=O以及C=C的鍵長(zhǎng)分別被拉長(zhǎng)0.003和0.002 nm,連接C=O和C=C的C―C鍵長(zhǎng)縮短0.002 nm,說(shuō)明該種吸附方式導(dǎo)致了分子的鍵長(zhǎng)平均化。對(duì)于CAL/Ir-Au(111)體系,其吸附能在85.18-152.21 kJ·mol-1之間,O原子距離表面為0.173-0.471 nm之間。當(dāng)CAL以C=O吸附于TopIr位時(shí),即以η2-diσCO(TopIr)構(gòu)型存在時(shí),其吸附能最大,為152.21 kJ·mol-1,較Ir(111)面得到提高27,O原子距離表面為0.233 nm。且吸附前后,C=O以及C=C的鍵長(zhǎng)分別被拉長(zhǎng)0.004和0.002 nm,連接C=O和C=C的C―C鍵長(zhǎng)縮短0.003 nm,同樣說(shuō)明該類吸附方式導(dǎo)致了分子的鍵長(zhǎng)平均化。圖3分別為CAL在In-Au(111)面和Ir-Au (111)面上的最穩(wěn)定吸附構(gòu)型。

    表2 CAL在In-Au(111)面上不同吸附位的計(jì)算結(jié)果Table 2 Calculation results of CALadsorption at different sites of In-Au(111)surface

    表3 CAL在Ir-Au(111)面上不同吸附位的計(jì)算結(jié)果Table 3 Calculation results of CALadsorption at different sites of Ir-Au(111)surface

    圖3 CAL在M-Au(111)面上的最穩(wěn)定吸附構(gòu)型Fig.3 The most stable adsorption configurations of CALon M-Au(111)surfaces

    綜上,當(dāng)CAL的C=O以σ鍵吸附于In-Au (111)面和Ir-Au(111)面的TopM位時(shí)最穩(wěn)定,其對(duì)應(yīng)的吸附能最大,表現(xiàn)為化學(xué)吸附。相比較該兩種穩(wěn)定吸附構(gòu)型,對(duì)于吸附能,CAL在Ir-Au(111)面上的吸附能略高于In-Au(111)面,同時(shí)都高于Au (111)面(93.70 kJ·mol-1)。對(duì)于鍵長(zhǎng)變化,CAL在Ir-Au(111)面上的鍵長(zhǎng)變化較In-Au(111)面更加明顯,且C=O的鍵長(zhǎng)變化也相對(duì)較大,從而可以推測(cè)CAL的C=O最容易受到破壞。因此,Ir-Au (111)面對(duì)CAL的吸附作用比In-Au(111)面對(duì)CAL的吸附作用更強(qiáng)。此外,相對(duì)于Pt(111)面26,28、Sn/ Pt(111)面29和Au團(tuán)簇30,M-Au(111)面表現(xiàn)出更強(qiáng)的CAL吸附能力,并且由于其獨(dú)特的C=O吸附構(gòu)型,使之更有利于C=O的活化。

    3.2.2 吸附構(gòu)型的電子性質(zhì)分析

    由上述分析可知,M原子的修飾改善了Au (111)面的吸附性能。為進(jìn)一步了解CAL與M-Au (111)面的電子交互作用,本文對(duì)最穩(wěn)定吸附構(gòu)型的態(tài)密度進(jìn)行了分析,如圖4所示。對(duì)于CAL/MAu(111)體系,在-17.18至-7.04 eV處,CAL與M-Au(111)面的電子雜化作用較弱,但CAL分子內(nèi)的s和p軌道雜化作用明顯,因此可以推測(cè)其主要構(gòu)鍵方式為共價(jià)鍵。在-7.04 eV至費(fèi)米能級(jí)處,主要是M-Au(111)的d軌道與CAL的p軌道發(fā)生雜化作用,這是CAL能夠穩(wěn)定吸附于M-Au(111)面的主要原因,M-Au(111)面的d軌道與CAL的s軌道雜化作用對(duì)此也有一定的貢獻(xiàn)。縱觀全圖,兩吸附體系的態(tài)密度圖總體類似,從而可以推斷其吸附及成鍵方式的一致性,此結(jié)論與吸附構(gòu)型分析的結(jié)果一致。

    圖4 CAL在M-Au(111)面上吸附的分波態(tài)密度(PDOS)和總態(tài)密度(TDOS)Fig.4 Partial densities of states(PDOS)and total densities of states(TDOS)of CALadsorption on M-Au(111)surfaces

    差分電荷密度是一種能夠直觀描述吸附后體系電子重排的分析方法31。為了對(duì)CAL與M-Au (111)面的電子交互作用有更加鮮明的理解,本文計(jì)算了最穩(wěn)定吸附構(gòu)型的差分電荷密度,如圖5所示。差分電荷密度定義為CAL吸附于M-Au(111)面前后電荷密度之差,即:Δρ=ρCAL/M-Au(111)-ρCAL-ρM-Au(111),其中ρCAL/M-Au(111)、ρCAL和ρM-Au(111)分別為體系的總電荷密度、CAL的電荷密度以及M-Au(111)面的電荷密度。圖中藍(lán)色區(qū)域代表失去電子,紅色區(qū)域代表得到電子,且顏色越深、區(qū)域越大代表得失電子數(shù)越多。由圖5可知,CAL在M-Au(111)面上的電子得失區(qū)域較為一致,從而進(jìn)一步驗(yàn)證了CAL在M-Au(111)面上吸附方式的一致性。同時(shí)可以看出,CAL/Ir-Au(111)體系中的顏色區(qū)域及深度都略大于CAL/In-Au(111)體系,因此可知Ir-Au (111)面對(duì)于CAL的吸附力略大于In-Au(111)面,該結(jié)論與吸附能的分析結(jié)果相吻合。

    圖5 CAL在M-Au(111)面上的差分電荷密度Fig.5 Deformation density of CALadsorption on M-Au(111)surfaces

    3.2.3 CAL的電子性質(zhì)分析

    為了能夠更加深入了解CAL吸附前后的電子性質(zhì)以及不同催化劑對(duì)其的影響,本文對(duì)吸附前以及吸附于M-Au(111)面上的CAL的分波態(tài)密度和Mulliken電荷布居進(jìn)行了研究。

    根據(jù)圖4的分析可知,CAL的電子貢獻(xiàn)主要來(lái)自于其p軌道,圖6作出了CAL吸附前以及吸附于M-Au(111)面后的分波態(tài)密度圖。從該圖可知,當(dāng)CAL吸附于M-Au(111)面后,CAL的分波態(tài)密度費(fèi)明顯降低,說(shuō)明吸附后CAL向M-Au(111)面轉(zhuǎn)移了大量電子。費(fèi)米能級(jí)右側(cè)的峰基本消失,說(shuō)明吸附后CAL的反鍵π*軌道受到破壞。費(fèi)米能級(jí)左側(cè)的峰明顯變小且向低能級(jí)移動(dòng),說(shuō)明吸附后CAL的p軌道與金屬表面發(fā)生強(qiáng)相互作用,形成σ鍵。相比較兩類吸附體系,CAL在Ir-Au(111)面上的p軌道偏移和形變程度都略大于In-Au(111)面,進(jìn)一步說(shuō)明了CAL與Ir-Au(111)面的相互作用更強(qiáng)。其中,位于-6.53及-3.39 eV處的高峰主要來(lái)自于C=O的貢獻(xiàn)4。

    圖6 CAL吸附前后p軌道的分波態(tài)密度Fig.6 p oribital partial density of states ofCALbefore and after adsorption

    結(jié)合文獻(xiàn)4報(bào)道以及表4可知,CAL的C=O相比較于C=C具有更強(qiáng)的極性,即C=O更有利于吸附于金屬表面。對(duì)于In-Au(111)面和Ir-Au(111)面上吸附的CAL,其O1、C2、C3、C4以及分子整體所帶電荷量分別為-0.34e、0.29e、-0.08e、0.02e、0.26e和-0.27e、0.29e、-0.08e、0.03e、0.36e。由此可以看出,吸附后CAL向金屬表面轉(zhuǎn)移電子,使得CAL帶有一定量的正電荷,并且Ir-Au(111)面上CAL的轉(zhuǎn)移電荷數(shù)相對(duì)更多,造成了CAL與該面更強(qiáng)的作用力。這可能是由于In和Ir原子都具有較大的有效核電荷,同時(shí)Ir原子的半徑較In原子小,使得In-Au(111)面上的空軌道具有更強(qiáng)的電子對(duì)接受能力。相比較分子鏈上的電荷轉(zhuǎn)移數(shù)可得,C=O的電荷轉(zhuǎn)移數(shù)較C=C多,說(shuō)明CAL的C=O與金屬的作用更強(qiáng)。其中,C2原子上的電荷變化量相對(duì)細(xì)微,這可能是由于金屬反饋部分電子到C=O上,并通過(guò)π鍵傳遞電子至C2原子上。

    綜上所述,在M-Au(111)面上,CAL最有可能通過(guò)C=O垂直吸附于M-Au(111)面的TopM位,且在Ir-Au(111)面上的吸附相對(duì)更穩(wěn)定。同時(shí)對(duì)于穩(wěn)定吸附構(gòu)型的態(tài)密度、差分電荷密度以及Mulliken電荷布居等分析表明:當(dāng)CAL以C=O穩(wěn)定吸附于M-Au(111)面時(shí),C=O能與金屬表面發(fā)生強(qiáng)相互作用,并且CAL/Ir-Au(111)體系的電荷轉(zhuǎn)移量較CAL/In-Au(111)體系更加顯著,進(jìn)一步說(shuō)明了Ir-Au(111)面對(duì)于CAL的吸附作用強(qiáng)于In-Au(111)面。

    表4 CAL在M-Au(111)面上的Mulliken電荷布居數(shù)Table 4 Mulliken charges of CALadsorption on M-Au(111)surfaces

    4 結(jié)論

    本研究基于DFT,首先計(jì)算了M原子修飾Au (111)面形成的3種模型合金的穩(wěn)定性,然后探討了最穩(wěn)定模型的化學(xué)活性以及CAL在其表面吸附的相關(guān)性質(zhì),結(jié)果如下:

    (1)In原子修飾的In-Au(111)面的模型穩(wěn)定性為(c)>(b)>(a)。對(duì)于以間位態(tài)存在的模型(c),其金屬表面幾何形變最小,形成能為8068.22 kJ·mol-1,結(jié)合能為391.71 kJ·mol-1,d帶中心為-3.54 eV。Ir原子修飾的Ir-Au(111)面的模型穩(wěn)定性為(a)>(b)>(c)。對(duì)于以聚集態(tài)存在的模型(a),其金屬表面幾何形變最小,形成能為10252.83 kJ·mol-1,結(jié)合能為425.78 kJ·mol-1,d帶中心為-3.48 eV。

    (2)CAL在In-Au(111)面和Ir-Au(111)面上的吸附,分別具有54種和70種初始吸附構(gòu)型,經(jīng)計(jì)算優(yōu)化后能穩(wěn)定存在的構(gòu)型各自只有17種和22種,并且當(dāng)CAL以C=O吸附于TopM位時(shí),CAL形變最大,吸附構(gòu)型最穩(wěn)定,吸附能分別為145.45和152.21 kJ·mol-1,轉(zhuǎn)移電荷量分別為0.26e和0.36e。其中,位于-7.04 eV至費(fèi)米能級(jí)處的p、d軌道雜化,對(duì)于CAL在M-Au(111)面上的吸附具有重要貢獻(xiàn)。

    (3)分析比較In-Au(111)面和Ir-Au(111)面,發(fā)現(xiàn)后者不僅幾何形變更小、形成能及結(jié)合能更大、構(gòu)型更穩(wěn)定以及d帶中心更低,而且對(duì)于CAL的吸附,表現(xiàn)出更加優(yōu)良的性能,即具有更大的吸附能、更低的態(tài)密度以及更明顯的電子轉(zhuǎn)移。因此,Ir-Au(111)面較In-Au(111)面的配體效應(yīng)更佳。同時(shí),相比較于Au(111)面,M原子的修飾改善了金屬表面的性質(zhì)及其對(duì)于CAL的吸附能力。

    (1) Haruta,M.;Kobayashi,T.;Sano,H.Chem.Lett.1987,16,405. doi:10.1246/cl.1987.405

    (2) Bus,E.;Prins,R.;Van Bokhoven,J.A.Catal.Commun.2007, 8,1397.doi:10.1016/j.catcom.2006.11.040

    (3) Bailie,J.E.;Hutching,G.J.Chem.Commun.1999,21,2151. doi:10.1039/a906538e

    (4) Jiang,J.H.;Xia,S.J.;Ni,Z.M.;Zhang,L.Y.Chem.J.Chin. Univ.2016,37,693.[蔣軍輝,夏盛杰,倪哲明,張連陽(yáng).高等學(xué)?;瘜W(xué)學(xué)報(bào),2016,37,693.]doi:10.7503/cjcu20150844

    (5) Aguirre,A.;Barrios,C.E.;Aguilar-Tapia,A.;Zanella,R.; Baltanás,M.A.;Collins,S.E.Top.Catal.2016,59,347. doi:10.1007/s11244-015-0425-6

    (6) Ngoepe,P.N.M.;Meyer,W.E.;Auret,F.D.;Omotoso,E.; Diale,M.;Swart,H.C.;Duvenhage,M.M.;Coetsee,E.Physica B 2016,480,209.doi:10.1016/j.physb.2015.08.051

    (7) Torres,C.;Campos,C.;Fierro,J.L.G.;Oportus,M.;Reyes,P. Catal.Lett.2013,143,763.doi:10.1007/s10562-013-1034-2

    (8) Lyu,J.H.;Wang,J.G.;Lu,C.S.;Ma,L.;Zhang,Q.F.;He,X. B.;Li,X.N.J.Phys.Chem.C 2014,118,2594.doi:10.1021/ jp411442f

    (9) Han,X.X.;Zhou,R.X.;Lai,G.H.;Zheng,X.M.Catal.Today 2004,93,433.doi:10.1016/j.cattod.2004.06.053

    (10) Zhang,L.Y.;Shi,W.;Xia,S.J.;Ni,Z.M.Acta Phys.-Chim. Sin.2014,30,1847.[張連陽(yáng),施 煒,夏盛杰,倪哲明.物理化學(xué)學(xué)報(bào),2014,30,1847.]doi:10.3866/PKU.WHXB201407141

    (11) Cárdenas-Lizana,F.;Gómez-Quero,S.;Hugon,A.;Delannoy, L.;Louis,C.;Keane,M.A.J.Catal.2009,262,235. doi:10.1016/j.jcat.2008.12.019

    (12) Huang,T.H.;Zhao,Y.J.;Tian,Z.F.;Li,X.L.;Liu,Q.;Zhao, D.Y.Acta Phys.-Chim.Sin.2014,30,2307.[黃天輝,趙玉娟,田兆福,李小蘭,劉 茜,趙東元.物理化學(xué)學(xué)報(bào),2014,30, 2307.]doi:10.3866/PKU.WHXB201410142

    (13) Li,X.H.;Zheng,W.L.;Pan,H.Y.J.Catal.2013,300,9. doi:10.1016/j.jcat.2012.12.007

    (14) Delbecq,F.;Li,Y.;Loffreda,D.J.Catal.2016,334,68. doi:10.1016/j.jcat.2015.10.028

    (15) Yang,Q.Y.;Zhu,Y.;Tian,L.;Xie,S.H.;Pei,Y.;Li,H.;Li,H. X.;Qiao,M.H.;Fan,K.N.Appl.Catal.A-Gen.2009,369,67. doi:10.1016/j.apcata.2009.08.032

    (16) He,X.;Chen,Z.X.;Kang,G.J.J.Phys.Chem.C 2009,113, 12325.doi:10.1021/jp9006729

    (17) Rojas,H.A.;Martínez,J.J.;Díaz,G.;Gómez-Cortés,A.Appl. Catal.A-Gen.2015,503,196.doi:10.1016/j.apcata.2015.07.023 (18) Zhao,J.;Ni,J.;Xu,J.H.;Xu,J.T.;Cen,J.;Li,X.N.Catal. Commun.2014,54,72.doi:10.1016/j.catcom.2014.05.012

    (19) Obot,I.B.;MacDonald,D.D.;Gasem,Z.M.Corrosion Sci. 2015,99,1.doi:10.1016/j.corsci.2015.01.037

    (20) Yu,Y.X.ACS Appl.Mater.Interfaces 2014,6,16267. doi:10.1021/am504452a

    (21) Manyer,H.G.;Yang,B.;Daly,H.;Moor,H.;McMonagle,S.; Tao,Y.;Yadav,G.D.;Goguet,A.;Hu,P.;Hardacre,C. ChemCatChem 2013,5,506.doi:10.1002/cctc.201200447

    (22) Mohr,C.;Hofmeister,H.;Radnik,J.;Claus,P.J.Am.Chem. Soc.2003,125,1905.doi:10.1021/ja027321q

    (23) Tao,J.;Yao,Z.J.;Xue,F.Fundamentals of Material Science;Chemical Industry Press:Beijing,2006;pp 50-51.[陶 杰,姚正軍,薛 烽.材料科學(xué)基礎(chǔ).北京:化學(xué)工業(yè)出版社,2006: 50-51.]

    (24) Chang,C.R.;Long,B.;Yang,X.F.;Li,J.J.Phys.Chem.C 2015,119,16072.doi:10.1021/acs.jpcc.5b03965

    (25) Amaya-Roncancio,S.;Linares,D.H.;Sapag,K.;Rojas,M.I. Appl.Surf.Sci.2015,346,438.doi:10.1016/j. apsusc.2015.03.215

    (26) Pirillo,S.;López-Corral,I.;Germán,E.;Juan,A.Vacuum 2014, 99,259.doi:10.1016/j.vacuum.2013.06.013

    (27) Yang,Y.J.;Teng,B.T.;Liu,Y.;Wen,X.D.Appl.Surf.Sci. 2015,357,369.doi:10.1016/j.apsusc.2015.09.039

    (28) Haubrich,J.;Loffreda,D.;Delbecq,F.;Jugnet,Y.;Sautet,P.; Krupski,A.;Becker,C.;Wandelt,K.Chem.Phys.Lett.2006, 433,188.doi:10.1016/j.cplett.2006.10.123

    (29) Delbecq,F.;Sautet,P.J.Catal.2003,220,115.doi:10.1016/ S0021-9517(03)00249-5

    (30) Zeinalipour-Yazdi,C.D.;Willock,D.J.;Machado,A.;Wilson, K.;Lee,A.F.Phys.Chem.Chem.Phys.2014,16,11236. doi:10.1039/c3cp53691b

    (31) Shi,W.;Zhang,L.Y.;Xia,S.J.;Ni,Z.M.Acta Phys.-Chim. Sin.2014,30,2249.[施 煒,張連陽(yáng),夏盛杰,倪哲明.物理化學(xué)學(xué)報(bào),2014,30,2249.]doi:10.3866/PKU.WHXB201408283

    Comparison of Properties of In-Au(111)and Ir-Au(111)Alloy Surfaces, and Their Adsorption to Crotonaldehyde

    JIANG Jun-Hui QIAN Meng-Dan XUE Ji-Long XIASheng-Jie NI Zhe-Ming*SHAO Meng-Meng
    (College of Chemical Engineering,Zhejiang University of Technology,Hangzhou 310014,P.R.China)

    The stability of M-Au(111)surfaces modified by M(M=In,Ir)was investigated using density functional theory.The most favorable model was selected to explore the chemical reactivity and adsorption of crotonaldehyde.The stability of the M-Au(111)surfaces was calculated using geometric configuration,and formation and cohesive energies.The calculations showed that the stability of the In-Au(111)surface increased as the atomic spacing of In was increased.Conversely,the Ir-Au(111)showed the opposite trend.The adsorption at the TopMsite was most stable when the crotonaldehyde on the M-Au(111)surfaces interacting via the C=O. Additionally,the adsorption energies were at their maximum.A combination of structural changes,density of state,deformation density and Mulliken charge analysis showed that the deformation of crotonaldehyde was larger than other adsorption modes,with an obvious charge transfer.Additionally,p and d orbital hybridization between-7.04 eV to Fermi level was found to have an important contribution to the adsorption process. Compared with the Au(111)surface,the stability and adsorption capacity of the M-Au(111)surfaces were significantly improved following modification by M atoms.Importantly,the Ir-Au(111)surface was found to have higher stability and activity,and stronger adsorption of crotonaldehyde than the In-Au(111)surface.

    In-Au(111)surface;Ir-Au(111)surface;Modification;Crotonaldehyde;Adsorption; Density functional theory

    O641

    10.3866/PKU.WHXB201609302

    Received:July 13,2016;Revised:September 30,2016;Published online:September 30,2016.

    *Corresponding author.Email:jchx@zjut.edu.cn;Tel:+86-13858123256.

    The project was supported by the National Natural Science Foundation of China(21503188).

    國(guó)家自然科學(xué)基金(21503188)資助項(xiàng)目

    猜你喜歡
    巴豆鍵長(zhǎng)構(gòu)型
    陰離子調(diào)控錳基鈣鈦礦中Mn─O的鍵長(zhǎng)和磁性
    分子和離子立體構(gòu)型的判定
    航天器受迫繞飛構(gòu)型設(shè)計(jì)與控制
    密度泛函理論研究鎘的二鹵化合物分子的結(jié)構(gòu)和振動(dòng)頻率
    淺議鍵能與鍵長(zhǎng)的關(guān)系
    遙感衛(wèi)星平臺(tái)與載荷一體化構(gòu)型
    巴豆養(yǎng)生漫畫
    健康必讀(2014年2期)2014-06-23 09:12:26
    兩個(gè)具stp三維拓?fù)錁?gòu)型的稀土配位聚合物{[Ln2(pda)3(H2O)2]·2H2O}n(Ln=Nd,La)
    風(fēng)吹來(lái)一封信
    意林注音版(2013年9期)2013-04-29 00:44:03
    巴豆“抓”小偷
    国产亚洲欧美精品永久| 天堂√8在线中文| 日本 av在线| 久热这里只有精品99| 国产欧美日韩精品亚洲av| 久久中文字幕人妻熟女| 一边摸一边抽搐一进一小说| 好男人在线观看高清免费视频 | or卡值多少钱| 香蕉丝袜av| 男女下面插进去视频免费观看| 大型av网站在线播放| 午夜精品久久久久久毛片777| 后天国语完整版免费观看| 69精品国产乱码久久久| 国内精品久久久久精免费| 精品少妇一区二区三区视频日本电影| 国产一卡二卡三卡精品| 日韩高清综合在线| 国产片内射在线| 天天一区二区日本电影三级 | 午夜久久久在线观看| 亚洲专区字幕在线| 人人妻,人人澡人人爽秒播| 黄色视频不卡| 亚洲人成伊人成综合网2020| 香蕉久久夜色| 日本欧美视频一区| 国产精品久久久久久人妻精品电影| 国产单亲对白刺激| 老熟妇乱子伦视频在线观看| 久久狼人影院| 18美女黄网站色大片免费观看| 校园春色视频在线观看| 久久国产精品影院| 亚洲成av片中文字幕在线观看| 日韩高清综合在线| 久久青草综合色| 久久久久久免费高清国产稀缺| 久久中文看片网| 一边摸一边抽搐一进一小说| 久久国产精品人妻蜜桃| 成人av一区二区三区在线看| 啦啦啦 在线观看视频| 亚洲色图 男人天堂 中文字幕| 国产日韩一区二区三区精品不卡| 久久婷婷成人综合色麻豆| 精品乱码久久久久久99久播| 一区在线观看完整版| 国产在线精品亚洲第一网站| 纯流量卡能插随身wifi吗| 天堂√8在线中文| 波多野结衣一区麻豆| 欧美一级a爱片免费观看看 | 妹子高潮喷水视频| 国产精品日韩av在线免费观看 | 99re在线观看精品视频| av网站免费在线观看视频| 欧美日韩亚洲国产一区二区在线观看| 国产精品98久久久久久宅男小说| 免费av毛片视频| 免费看十八禁软件| 亚洲欧美日韩高清在线视频| 女人被躁到高潮嗷嗷叫费观| 亚洲性夜色夜夜综合| 日韩一卡2卡3卡4卡2021年| 亚洲国产欧美网| 成人国产一区最新在线观看| 国产三级黄色录像| 无人区码免费观看不卡| 中文字幕精品免费在线观看视频| avwww免费| 久久久久久久精品吃奶| 久久影院123| 一区二区三区激情视频| 国产成年人精品一区二区| 给我免费播放毛片高清在线观看| 国产99久久九九免费精品| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久久毛片微露脸| 极品人妻少妇av视频| 亚洲国产日韩欧美精品在线观看 | 女性被躁到高潮视频| 91av网站免费观看| 久久久久久久久久久久大奶| 天天添夜夜摸| 久久久久久久久久久久大奶| 国产色视频综合| 99香蕉大伊视频| 中文字幕av电影在线播放| 日韩精品青青久久久久久| 久久久国产精品麻豆| 国产精品秋霞免费鲁丝片| 看黄色毛片网站| 亚洲精品一区av在线观看| 一个人免费在线观看的高清视频| 亚洲av电影不卡..在线观看| 成人国产综合亚洲| 中文字幕人妻丝袜一区二区| 无人区码免费观看不卡| 欧美日韩精品网址| 亚洲激情在线av| 桃红色精品国产亚洲av| 国产99久久九九免费精品| 97碰自拍视频| av网站免费在线观看视频| 级片在线观看| 波多野结衣高清无吗| av在线天堂中文字幕| 久久久国产成人精品二区| 久久精品亚洲熟妇少妇任你| 国产一卡二卡三卡精品| x7x7x7水蜜桃| 老司机在亚洲福利影院| 国产私拍福利视频在线观看| 无人区码免费观看不卡| 欧美绝顶高潮抽搐喷水| 国产主播在线观看一区二区| 久9热在线精品视频| 看片在线看免费视频| 午夜免费成人在线视频| 国产精品久久视频播放| 中文字幕av电影在线播放| 免费看a级黄色片| 免费看美女性在线毛片视频| 欧美国产日韩亚洲一区| 在线观看www视频免费| 在线天堂中文资源库| 非洲黑人性xxxx精品又粗又长| 又黄又粗又硬又大视频| 精品国产美女av久久久久小说| 成年人黄色毛片网站| 国产1区2区3区精品| 最好的美女福利视频网| 又大又爽又粗| 91在线观看av| 精品国产美女av久久久久小说| 如日韩欧美国产精品一区二区三区| 波多野结衣巨乳人妻| 波多野结衣av一区二区av| 免费搜索国产男女视频| 久久 成人 亚洲| 午夜福利18| 18美女黄网站色大片免费观看| 免费在线观看日本一区| 国产成人系列免费观看| 一边摸一边抽搐一进一小说| 91九色精品人成在线观看| 午夜福利欧美成人| 免费在线观看黄色视频的| 搡老熟女国产l中国老女人| 国产又爽黄色视频| 精品熟女少妇八av免费久了| 成年版毛片免费区| 琪琪午夜伦伦电影理论片6080| 亚洲成a人片在线一区二区| 色婷婷久久久亚洲欧美| www.999成人在线观看| 亚洲人成伊人成综合网2020| 一边摸一边做爽爽视频免费| 久久人妻福利社区极品人妻图片| 黑人欧美特级aaaaaa片| 亚洲国产精品sss在线观看| 18美女黄网站色大片免费观看| 免费观看精品视频网站| 成人三级黄色视频| 久9热在线精品视频| 国产精品亚洲美女久久久| 国产高清激情床上av| 午夜精品在线福利| 熟妇人妻久久中文字幕3abv| 岛国视频午夜一区免费看| 久久这里只有精品19| 欧美日韩瑟瑟在线播放| 亚洲av日韩精品久久久久久密| 老司机福利观看| 久久精品91蜜桃| 久久久国产精品麻豆| 欧美另类亚洲清纯唯美| 12—13女人毛片做爰片一| 国内毛片毛片毛片毛片毛片| 午夜免费鲁丝| 免费少妇av软件| 午夜免费激情av| 丁香六月欧美| 黄色 视频免费看| 国产精品久久电影中文字幕| 久久久久久亚洲精品国产蜜桃av| 欧美另类亚洲清纯唯美| 精品少妇一区二区三区视频日本电影| 成人三级黄色视频| 国产高清激情床上av| av视频免费观看在线观看| 亚洲无线在线观看| 叶爱在线成人免费视频播放| 精品免费久久久久久久清纯| 日本 欧美在线| 欧美成狂野欧美在线观看| 很黄的视频免费| 一本综合久久免费| 国产三级在线视频| or卡值多少钱| 日韩一卡2卡3卡4卡2021年| 我的亚洲天堂| 悠悠久久av| 久久精品人人爽人人爽视色| 亚洲人成伊人成综合网2020| 亚洲在线自拍视频| 久久婷婷成人综合色麻豆| 久久精品影院6| 久99久视频精品免费| 久久 成人 亚洲| 国产精品影院久久| 国产午夜精品久久久久久| 一进一出抽搐动态| 波多野结衣巨乳人妻| 人人妻人人澡欧美一区二区 | 午夜福利成人在线免费观看| 欧美日韩精品网址| 精品一品国产午夜福利视频| 午夜免费成人在线视频| 九色亚洲精品在线播放| 国产精品久久久av美女十八| 国产成人精品久久二区二区91| 午夜成年电影在线免费观看| 亚洲五月色婷婷综合| 国产精品九九99| av片东京热男人的天堂| 国产精品国产高清国产av| 亚洲av美国av| x7x7x7水蜜桃| 又黄又粗又硬又大视频| 欧美黑人精品巨大| 给我免费播放毛片高清在线观看| 99国产精品99久久久久| 村上凉子中文字幕在线| 丝袜美腿诱惑在线| 日本vs欧美在线观看视频| 国产亚洲精品一区二区www| 男人舔女人的私密视频| 国产一区二区三区综合在线观看| 久久伊人香网站| 人人妻,人人澡人人爽秒播| a在线观看视频网站| 亚洲国产欧美网| 无遮挡黄片免费观看| 美女大奶头视频| 又紧又爽又黄一区二区| 久久久国产欧美日韩av| 男女下面进入的视频免费午夜 | 亚洲伊人色综图| 国产精品九九99| 久久久国产成人免费| 亚洲av成人一区二区三| 国产精品影院久久| 黄色视频不卡| 真人一进一出gif抽搐免费| 黄色丝袜av网址大全| 亚洲精品在线美女| 亚洲熟妇中文字幕五十中出| 一进一出抽搐动态| 中文字幕精品免费在线观看视频| 国产精品香港三级国产av潘金莲| 亚洲国产欧美一区二区综合| 亚洲欧洲精品一区二区精品久久久| 亚洲成人免费电影在线观看| 亚洲熟妇熟女久久| 国产麻豆69| 两性午夜刺激爽爽歪歪视频在线观看 | 91麻豆精品激情在线观看国产| 国产av又大| 精品高清国产在线一区| 国产欧美日韩一区二区三| 热99re8久久精品国产| 精品久久久精品久久久| 亚洲国产毛片av蜜桃av| 成人免费观看视频高清| x7x7x7水蜜桃| netflix在线观看网站| 老汉色∧v一级毛片| 91精品三级在线观看| 一区在线观看完整版| 国产精品野战在线观看| 亚洲熟妇中文字幕五十中出| 9色porny在线观看| 亚洲激情在线av| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av在线| 国产欧美日韩一区二区精品| 国产亚洲av高清不卡| 久久精品国产亚洲av高清一级| 精品卡一卡二卡四卡免费| 亚洲精品美女久久av网站| 亚洲精品国产色婷婷电影| 一级a爱片免费观看的视频| 免费在线观看黄色视频的| 久久精品91蜜桃| 免费人成视频x8x8入口观看| 757午夜福利合集在线观看| 天天添夜夜摸| 中文字幕av电影在线播放| 99热只有精品国产| 黑人欧美特级aaaaaa片| 岛国视频午夜一区免费看| 日本三级黄在线观看| 无遮挡黄片免费观看| 国产成人av教育| 丝袜在线中文字幕| tocl精华| 国产一级毛片七仙女欲春2 | 久久午夜综合久久蜜桃| 一区福利在线观看| 熟女少妇亚洲综合色aaa.| 日韩成人在线观看一区二区三区| 国产区一区二久久| 日本五十路高清| 别揉我奶头~嗯~啊~动态视频| 黑人巨大精品欧美一区二区mp4| cao死你这个sao货| 日韩精品中文字幕看吧| 狠狠狠狠99中文字幕| 在线观看舔阴道视频| 国产黄a三级三级三级人| 亚洲情色 制服丝袜| 久久婷婷人人爽人人干人人爱 | 亚洲av美国av| 91老司机精品| 人人妻人人爽人人添夜夜欢视频| 9色porny在线观看| 国产私拍福利视频在线观看| 成人18禁高潮啪啪吃奶动态图| 日韩欧美一区二区三区在线观看| 在线播放国产精品三级| 国产一卡二卡三卡精品| 精品第一国产精品| 国产成年人精品一区二区| av天堂久久9| 免费av毛片视频| 日韩欧美在线二视频| 巨乳人妻的诱惑在线观看| 桃色一区二区三区在线观看| av中文乱码字幕在线| 国产国语露脸激情在线看| 悠悠久久av| 欧美色视频一区免费| 夜夜看夜夜爽夜夜摸| 久久精品91无色码中文字幕| 一级a爱视频在线免费观看| 久久久久久久精品吃奶| 老司机在亚洲福利影院| 一边摸一边抽搐一进一小说| 亚洲成人免费电影在线观看| 久久午夜亚洲精品久久| 女人被躁到高潮嗷嗷叫费观| 操出白浆在线播放| 久久久久亚洲av毛片大全| 国产精品美女特级片免费视频播放器 | 高清黄色对白视频在线免费看| 国产伦一二天堂av在线观看| 亚洲欧美一区二区三区黑人| 熟妇人妻久久中文字幕3abv| 99国产精品一区二区蜜桃av| 亚洲精品一卡2卡三卡4卡5卡| 欧美在线一区亚洲| 亚洲专区国产一区二区| 午夜免费激情av| 一区二区三区激情视频| 欧美乱妇无乱码| 亚洲精品国产区一区二| 久99久视频精品免费| 日本免费一区二区三区高清不卡 | 免费少妇av软件| 欧美一级毛片孕妇| 亚洲国产看品久久| 日韩欧美在线二视频| 日韩欧美一区视频在线观看| 成人国语在线视频| 欧美在线黄色| 性欧美人与动物交配| 涩涩av久久男人的天堂| www国产在线视频色| 日本五十路高清| 亚洲成a人片在线一区二区| 人人妻人人澡欧美一区二区 | 午夜免费观看网址| 高潮久久久久久久久久久不卡| 午夜福利,免费看| 午夜免费鲁丝| 美女高潮到喷水免费观看| 亚洲精品av麻豆狂野| 亚洲av日韩精品久久久久久密| 免费一级毛片在线播放高清视频 | 亚洲av五月六月丁香网| 亚洲精品久久国产高清桃花| 国产一级毛片七仙女欲春2 | 99久久99久久久精品蜜桃| 成在线人永久免费视频| 色婷婷久久久亚洲欧美| 大型av网站在线播放| 老司机在亚洲福利影院| 亚洲五月婷婷丁香| 久久久久久人人人人人| 国产亚洲欧美98| 亚洲美女黄片视频| 成人三级做爰电影| 午夜免费观看网址| 亚洲色图综合在线观看| 天天添夜夜摸| 久久久久国产精品人妻aⅴ院| 久久久久久久久久久久大奶| 激情视频va一区二区三区| 给我免费播放毛片高清在线观看| 999久久久国产精品视频| 精品无人区乱码1区二区| 91成人精品电影| 欧美一级毛片孕妇| 午夜免费观看网址| 啦啦啦韩国在线观看视频| 天堂动漫精品| 黄色丝袜av网址大全| 悠悠久久av| 国产高清视频在线播放一区| 免费女性裸体啪啪无遮挡网站| 淫妇啪啪啪对白视频| 成人三级黄色视频| 三级毛片av免费| 欧美 亚洲 国产 日韩一| 欧美黑人欧美精品刺激| 三级毛片av免费| 国产午夜福利久久久久久| 99在线视频只有这里精品首页| 久久人人97超碰香蕉20202| 国产成人av教育| 波多野结衣av一区二区av| 精品久久蜜臀av无| 1024香蕉在线观看| 精品久久久久久久人妻蜜臀av | xxx96com| 国产黄a三级三级三级人| 精品卡一卡二卡四卡免费| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲欧美日韩无卡精品| 亚洲精品中文字幕一二三四区| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美日韩无卡精品| 99精品久久久久人妻精品| 女性被躁到高潮视频| 一个人免费在线观看的高清视频| 黄片播放在线免费| 美女免费视频网站| 国产精品美女特级片免费视频播放器 | 可以在线观看毛片的网站| 久久精品国产清高在天天线| 在线av久久热| 亚洲 国产 在线| 又大又爽又粗| 女性被躁到高潮视频| 极品教师在线免费播放| 欧美最黄视频在线播放免费| 欧美激情 高清一区二区三区| 亚洲 欧美一区二区三区| 国产精品av久久久久免费| 亚洲午夜精品一区,二区,三区| 亚洲人成电影观看| 嫩草影视91久久| 首页视频小说图片口味搜索| 在线永久观看黄色视频| 亚洲精品av麻豆狂野| 91av网站免费观看| 亚洲人成网站在线播放欧美日韩| 一进一出抽搐动态| 在线观看一区二区三区| 亚洲欧美日韩高清在线视频| 露出奶头的视频| 欧美日韩瑟瑟在线播放| 午夜精品在线福利| 午夜福利18| 在线播放国产精品三级| 成人免费观看视频高清| 曰老女人黄片| 国产av又大| 亚洲美女黄片视频| 久久影院123| 久久久久久大精品| 欧美精品亚洲一区二区| 99香蕉大伊视频| 久久人人精品亚洲av| 国内久久婷婷六月综合欲色啪| 91精品三级在线观看| 亚洲一区二区三区不卡视频| 99riav亚洲国产免费| 91字幕亚洲| e午夜精品久久久久久久| 少妇被粗大的猛进出69影院| 免费一级毛片在线播放高清视频 | 精品一区二区三区视频在线观看免费| 免费在线观看亚洲国产| 婷婷六月久久综合丁香| 丁香六月欧美| 亚洲国产精品成人综合色| 欧美日韩黄片免| 欧美黄色淫秽网站| 免费在线观看完整版高清| 男男h啪啪无遮挡| 国产精品久久久久久人妻精品电影| 久久婷婷成人综合色麻豆| 国产97色在线日韩免费| 国产欧美日韩一区二区三| 精品高清国产在线一区| av超薄肉色丝袜交足视频| 午夜福利高清视频| 欧美中文综合在线视频| 久久久久久国产a免费观看| 一边摸一边抽搐一进一出视频| 真人一进一出gif抽搐免费| 亚洲欧美日韩无卡精品| 精品一区二区三区av网在线观看| 99国产综合亚洲精品| 国产亚洲av高清不卡| 很黄的视频免费| 男女做爰动态图高潮gif福利片 | 欧美黄色片欧美黄色片| 一本综合久久免费| 国产色视频综合| 波多野结衣av一区二区av| 精品第一国产精品| 国产伦人伦偷精品视频| 久久婷婷成人综合色麻豆| 久久草成人影院| 国产成年人精品一区二区| 成人永久免费在线观看视频| 成人手机av| 日韩成人在线观看一区二区三区| 亚洲在线自拍视频| 日本a在线网址| 国产99白浆流出| 老司机午夜福利在线观看视频| 人人妻人人爽人人添夜夜欢视频| 免费在线观看视频国产中文字幕亚洲| 国产av又大| 国产真人三级小视频在线观看| 久久青草综合色| 男女下面插进去视频免费观看| 一进一出抽搐gif免费好疼| 亚洲avbb在线观看| 视频区欧美日本亚洲| 亚洲色图av天堂| 可以在线观看的亚洲视频| 精品一品国产午夜福利视频| 一区二区三区国产精品乱码| 久久精品91无色码中文字幕| av中文乱码字幕在线| 久久人妻av系列| 国产伦人伦偷精品视频| 亚洲国产欧美日韩在线播放| 成人亚洲精品av一区二区| 国产三级黄色录像| 久久国产精品影院| 久久久久久亚洲精品国产蜜桃av| 免费在线观看完整版高清| 亚洲自偷自拍图片 自拍| 国产熟女午夜一区二区三区| 美女午夜性视频免费| 可以在线观看毛片的网站| 老鸭窝网址在线观看| 日韩精品青青久久久久久| 国产精品永久免费网站| 国产精品av久久久久免费| 国产99白浆流出| 一区二区三区高清视频在线| 日韩欧美在线二视频| 亚洲精品国产精品久久久不卡| 亚洲av五月六月丁香网| 国内精品久久久久久久电影| 亚洲,欧美精品.| 9191精品国产免费久久| 午夜福利在线观看吧| 精品人妻1区二区| 91老司机精品| cao死你这个sao货| 人人妻,人人澡人人爽秒播| 亚洲情色 制服丝袜| 满18在线观看网站| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧洲精品一区二区精品久久久| 日韩欧美免费精品| 老司机靠b影院| 国产成人啪精品午夜网站| 黑丝袜美女国产一区| 国产一级毛片七仙女欲春2 | 黑人操中国人逼视频| 一本久久中文字幕| 99re在线观看精品视频| 色综合婷婷激情| 欧美乱妇无乱码| 久久精品aⅴ一区二区三区四区| 9191精品国产免费久久| 国产av又大| 精品人妻在线不人妻| 国产伦一二天堂av在线观看| 99香蕉大伊视频| 久久精品国产清高在天天线| 欧美国产精品va在线观看不卡| 国产成人av教育| 国产成年人精品一区二区| 欧美日韩精品网址| 精品无人区乱码1区二区| 人人妻人人爽人人添夜夜欢视频| 激情在线观看视频在线高清| 免费不卡黄色视频| 免费高清视频大片| 麻豆成人av在线观看| 欧美一级毛片孕妇| 一个人免费在线观看的高清视频| 欧美日韩黄片免|