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

    堿性介質(zhì)中CO對(duì)甲醇在PtAu(111)和Pt(111)表面氧化生成甲酸的影響的第一性原理研究

    2016-09-05 13:04:09劉建紅1呂存琴1春1王貴昌31山西大同大學(xué)化學(xué)與環(huán)境工程學(xué)院山西大同0370092先進(jìn)能源材料化學(xué)教育部重點(diǎn)實(shí)驗(yàn)室南開大學(xué)天津3000713南開大學(xué)化學(xué)學(xué)院天津市金屬與分子基材料化學(xué)重點(diǎn)實(shí)驗(yàn)室南開大學(xué)天津300071
    物理化學(xué)學(xué)報(bào) 2016年4期
    關(guān)鍵詞:山西大同南開大學(xué)堿性

    劉建紅1 呂存琴1,2 晉 春1,* 王貴昌3,*(1山西大同大學(xué)化學(xué)與環(huán)境工程學(xué)院,山西大同037009;2先進(jìn)能源材料化學(xué)教育部重點(diǎn)實(shí)驗(yàn)室(南開大學(xué)),天津300071;3南開大學(xué)化學(xué)學(xué)院,天津市金屬與分子基材料化學(xué)重點(diǎn)實(shí)驗(yàn)室(南開大學(xué)),天津300071)

    堿性介質(zhì)中CO對(duì)甲醇在PtAu(111)和Pt(111)表面氧化生成甲酸的影響的第一性原理研究

    劉建紅1呂存琴1,2晉春1,*王貴昌3,*
    (1山西大同大學(xué)化學(xué)與環(huán)境工程學(xué)院,山西大同037009;2先進(jìn)能源材料化學(xué)教育部重點(diǎn)實(shí)驗(yàn)室(南開大學(xué)),天津300071;3南開大學(xué)化學(xué)學(xué)院,天津市金屬與分子基材料化學(xué)重點(diǎn)實(shí)驗(yàn)室(南開大學(xué)),天津300071)

    采用密度泛函理論計(jì)算研究了堿性介質(zhì)中甲醇在清潔的PtAu(111)和Pt(111)表面、及有CO存在的PtAu(111)和Pt(111)表面的氧化。計(jì)算結(jié)果表明,在堿性介質(zhì)中,預(yù)吸附的CO促進(jìn)了甲醇在PtAu(111)和Pt(111)表面氧化的每一步反應(yīng),這與其在Au(111)表面的作用相似。究其原因,是由于CO的吸附增強(qiáng)了OH的穩(wěn)定性和堿性,從而增強(qiáng)了OH奪取氫原子的能力。

    甲醇氧化;M(111)(M=PtAu,Pt);堿性介質(zhì);CO促進(jìn)作用;密度泛函理論計(jì)算

    [Article]

    www.whxb.pku.edu.cn

    In recent years,there are several successful attempts to prove the ability of theoretical modeling based on density functional theory(DFT)calculations in providing detailed mechanism and kinetics for complex surface reactions.However,although some theoretical studies about methanol oxidation have been performed14,15,the physical origin is unclear.The minimal energy pathways for COOH formation from CO and OH co-adsorbed on Au(111)and Au(110)surfaces have been investigated.The calculation results show that the ability of CO to enhance the adsorption of OH is related to the surface structure of gold and changes in the work function as induced by the adsorption of CO14. Wang et al.15investigated all elementary steps of methanol dehydrogenation and oxidation on atomic-oxygen-covered or OH-covered Au(111)surfaces,and they found that the reaction pathway is related to the oxygen coverage.On low oxygen coverage,methanol decomposition starts from α-H elimination and β-H elimination to produce CH2O and CO.On high oxygen coverage,the CH2O formed through the elimination of α-H,and one β-H is suitable for the cooperative interaction of two nearby atomic oxygen.Recently Lv et al.16investigated the kinetics of methanol oxidation into formic acid catalyzed by Au(111)surface with and without CO in alkaline media.The calculation results show that the pre-adsorbed CO promotes almost every step of methanol oxidation both thermochemically and kinetically.Yuan et al.17have studied the methanol oxidation on PtAu(111)surface via CO pathway,and proposed that proper arrangement of Au and Pt sites is suitable for non-CO pathway.Inspired by Yuan et al., methanol oxidation on the PtAu(111)surface along the non-CO pathway has been investigated,and a comparison has been made with that along the CO pathway by Zhong et al.18They found that the non-CO pathway of methanol oxidation on PtAu(111)is energetically more favorable than the CO pathway.

    The catalytic properties of thin metal overlayers deposited on a foreign substrate could be enhanced with respect to those of the parent metals19-23.PtAu bimetallic surfaces show unique catalytic properties for methanol oxidation23,which could be attributed to the synergistic activity of Au and Pt.Analogous to the promoting effect of adsorbed carbon monoxide on the oxidation of methanol on a gold catalyst verified by experiments6and theoretical calculation16,is there the same effect for methanol oxidation to formic acid on platinum and PtAu alloy?To the best of our knowledge, there are no relevant reports about the effect of CO on the electrocatalytic activity of metal platinum and PtAu alloy in alkaline media obtained for methanol oxidation to formic acid.Therefore, we replaced the top slab of Au atoms on Au(111)surface with Pt atoms,referred as PtAu(111)surface hereafter,and investigated the methanol oxidation to formic acid on this surface and Pt(111) surface with and without CO by using density functional theory methods.Detailed comparisons for methanol oxidation to formic acid in alkaline media have been made among PtAu(111),Pt(111), and Au(111)surfaces,and similar promotion effect of CO on methanol oxidation has been found.

    2 Calculation methods and models

    All DFT calculations were implemented by using the Vienna ab initio simulation package(VASP)24,25.The exchange-correlation energy and electron-ion interaction were described at the generalized gradient approximation level(Perdew-Wang 9126)and the projector-augment wave(PAW)scheme27,28.The energy cutoff applied to the plane wave basis has been set as 400 eV.The calculations were done on a(3×3)four-layer slab separated by a vacuum region around 1.5 nm.A 3×3×1 Monkhorst-Pack29kpoint sampling was implemented in all calculations.The geometry optimizations were performed by using the conjugated gradient method,and the top two metal layers as well as the adsorbed molecule were allowed to relax until the maximal forces on each relaxed atom were less than 0.5 eV?nm-1.The molecules in the gas phase have been calculated using a box with the size of 1.7 nm× 1.8 nm×1.9 nm,and only a Г-point was used during the calculation.Spin-polarized calculations were performed when needed. To find the saddle points along the adiabatic minimal energy path (MEP),which connected each initial state(IS)and final state(FS) of a given elementary step,the nudged elastic band(NEB)with quasi-Newton optimizer30implemented in the VASP has been systematically utilized.Eventually,the transition states were identified by the existence of only one normal vibrational mode associated with a pure imaginary frequency.The optimized lattice constants of Au and Pt are 0.418 and 0.399 nm,which agrees well with the experimentally measured values of 0.408 and 0.392 nm31for Au and Pt,respectively.Lattice mismatch resulted in the contraction or expansion of lattice of the bimetal surface compared with that of bulk pure metal32.For the PtAu(111)system,the calculations were performed with the DFT optimized clean Au lattice constant of 0.418 nm.In addition,we have also performed test calculations for the geometry convergence criterion(0.5 eV?nm-1)used in the present work.When the force tolerance is improved to 0.3 eV?nm-1,the calculated adsorption energy of OH is only changed by 0.004 eV,confirming the enough accuracy of the conclusion obtained using the more loose convergence criterion of 0.5 eV?nm-1.

    There are four kinds of adsorption sites on the surface:top, bridge,hcp,and fcc sites.The adsorption energy(Eads)is defined as follows:Eads=Eadsorbate/substrate-Eadsorbate-Esubstrate,where Eadsorbate/substrateis the total energy of the adsorbate-substrate system,Eadsorbateexpresses the energy of the adsorbate in the gas phase,and Esubstrateexpresses the total energy of the substrate.In this work,for a reaction like AB=A+B,the calculated total energy change or reaction enthalpy(?H)is defined by the formula of?H=E(A+ B)/M-E(AB)/M,where E(A+B)/M and E(AB)/M are the total energies for the coadsorption system of the product and of the reactant,respectively.Activation energy,Ea,is calculated based on the equation of Ea=ETS-E(AB)/M,where ETSis the total energy of the transition state(TS).

    Owing to the complexity of real-world direct methanol fuel cells systems,gas phase methanol oxidation investigations are often been used as prototype reactions to probe the mechanism of low-temperature methanol fuel cells18.Therefore,our present calculations were performed with gas phase models to investigate the effect of CO on the methanol oxidation to formic acid in alkaline media.

    3 Results and discussion

    Our previous study show that the CO adsorbed at the top site on Au(111)surface can promote the oxidation of methanol to formic acid16.We consider whether CO adsorbed on the other kinds of metals can have the same or the opposite effect on methanol oxidation to formic acid.For analyzing the effect of CO on methanol oxidation,we calculated the transition states for relevant reactions involved in methanol oxidation to formic acid on PtAu (111)and Pt(111)surfaces with and without CO.In order to compare the calculated results with that onAu(111)surface better, the following possible reaction mechanism were considered in our present work.

    3.1Effect of CO to oxidize the methanol to formic acid on PtAu(111)surface

    To make the comparison more reliable,we give the same initial state on PtAu(111)surface with that on Au(111)surface,i.e.,the molecular CO adsorbs on the top site of Pt for all the initial states, and the optimized configurations are shown in Figs.1 and 2.The corresponding energetic data are reported in Table 1 and Fig.3.

    3.1.1Reaction of CH3OH and OH(CH3OH+OH→CH3O+H2O)

    For the initial oxidative dehydrogenation step of methanol, methanol is adsorbed on the top site via oxygen atom on PtAu(111)and CO/PtAu(111)surfaces.The distance between the O atom of O―H and the H atom of the CH3O―H is 0.157 nm on CO/PtAu(111)surface,which is 0.007 nm smaller than that on PtAu(111)surface.In fact,on the PtAu(111)surface,the reaction of dehydrogenation for methanol is almost a neutral reaction with a ΔH equal to-0.01 eV,while on the CO/PtAu(111)surface,it becomes an exothermic reaction releasing only 0.11 eV in energy. Moreover,the energy barrier only increased by 0.10 eV with the presence of adsorbed CO.The slight difference of reaction energy and activation barrier on PtAu(111)surface with and without CO are the reasons that we donot take this step as the example to talk about the effect of CO.

    3.1.2Reaction of CH3O and OH(CH3O+OH→CH2O+ H2O)

    In relation to the second step,from Table 1,we can find that the barrier of this step on PtAu(111)surface is higher than that of 0.56 eV on Au(111)surface16,and the reason why this step has so high barrier to happen is that Pt has stronger adsorption ability thanAu, which makes CH3O more stable on the PtAu(111)surface and harder to dissociate from the surface to react with OH.In the presence of CO on PtAu(111)surface,the energy barrier is reduced to 0.64 eV from 0.99 eV on PtAu(111)surface.As we cannot find any obvious differences between the initial states on PtAu(111)surface with and without CO,we try to find where the differences come from with the analysis of the reaction steps. Through the transition states analysis,shown in Figs.1 and 2,we can find that this reaction goes along similar paths on the different surfaces.On PtAu(111)and CO/PtAu(111)surfaces,the CH3O first dissociates from the surface and the abstracted H atom which will react with OH tries to approach the OH.Then after the C―H bond cleavage and O―H formation,the C―O bond rotates until the C,O,and H atoms of CH2O on a same plane.This induces the formation of the hydrogen bond between the CH2O and H2O to lower the energy of this system.From Figs.1 and 2,it can be seen that one H bonding of C―H…O is formed between C―H bond of CH3O and O of OH without and with co-adsorbed CO.However,besides the C―H…O bond,it seems another H bonding is formed with the presence of CO,i.e.,the O―H…O bond between H in OH and O in CO,which may lead to the increased stability in the transition state and final state and therefore the change in reaction heat and activation barrier.

    3.1.3Reaction of CH2O and OH(CH2O+OH→CHO+ H2O)

    The third step of methanol oxidation to formic acid has a little difference on the two different surfaces in the enthalpy by 0.05 eV, which can be explained by the data listed in Figs.1 and 2.We can point out that there is no big difference between the final stateson the two different surfaces,while the initial state on the CO/ PtAu(111)surface is less stable than that on the PtAu(111)surface, with a longer distance of 0.013 nm between the H and O which combines to react.However,with the existence of CO,it is obviously advantageous for the initial state to overcome the energy barrier of 0.38 eV compared with that of 0.85 eV on PtAu(111) surface.From the configurations shown in Figs.1 and 2,we can find that there are many differences between the two TSs.First, the C―H bond,which will break in this step,is elongated to 0.130 nm from 0.113 nm on the CO/PtAu(111)surface.In comparison, it is elongated to 0.136 nm from 0.112 nm on the PtAu(111) surface.Therefore,it needs more energy for the surface with CO to catalyze the reaction.Besides,there exists a big difference of almost 0.07 nm for the distance between C and Pt on the surface, which lowers the energy of the transition state on the surface with CO.Through the figures of the transition states,we can figure out that this step also goes in different paths on the two different surfaces.On the surface with CO existence,CH2O comes to the surface firstly and then the C―H bond breaks to react with OH. And finally,the O―H bond forms.However,on the surface without CO,the C―H bond breaks down first and at the same time,CHO goes down to the surface to form a C―Pt bond.Therefore,the energy barrier on the PtAu(111)surface is much higher than that on the surface with the existence of CO by 0.47 eV.

    Fig.1 Optimized adsorption configurations of the initial states,transition states,and final states involved in the methanol oxidation to formic acid on the PtAu(111)surface

    Fig.2 Optimized adsorption configurations of the initial states,transition states,and final states involved in the methanol oxidation to formic acid on the CO/PtAu(111)surface

    3.1.4Reaction of CHO and OH(CHO+OH→HCOOH)

    For the fourth step,from Table 1,we can find that the reaction on the surface with CO existence releases more energy than that on the clean PtAu(111)surface by 0.29 eV,which means that it is more beneficial for this step on the PtAu(111)surface with the existence of CO.From Figs.1 and 2,it can be found that the distance of OH and CHO is 0.011 nm closer with the existence of CO,which may explain the difference of the enthalpy for this reaction.

    From the computed energetic data on PtAu(111)surface with and without CO shown above,we find that the existence of CO is not beneficial for the first step,but it is beneficial for the other three steps.From Figs.1 and 2,we can find that there is a big difference between the adsorption sites on the surface for CO in the four steps.In the first and second steps,CO adsorbs on the fcc site and its C―O bond is about 0.120 nm.While in the third and fourth steps,CO sits on the top site with a bond length of 0.116nm.When CO is adsorbed on the fcc site,it has a larger adsorption energy which makes the PtAu(111)surface poisonous and it is hard to catalyze the reaction.However,when CO adsorbs on the top site,it has the same effect as that on the Au(111)surface,that is,CO is beneficial for OH to oxidize the methanol to HCOOH16. To check out the reason why CO adsorbed on the fcc site will hinder the reaction,we calculated the third step on the PtAu(111) surface with CO on the fcc site.According to the results,we found that the activation barrier is 0.48 eV,which is 0.10 eV higher than that when CO adsorbs on the top site,but it is still lower than that on PtAu(111)surface without CO.This result can prove that CO adsorbed on the fcc site will also promote the reaction,and the promotion effect is smaller than that on the top site.Through the analysis of the adsorption energy of OH,we found that the adsorption energy of OH with the presence of CO at the fcc site is also smaller than that on the clean surface by 0.16 eV,which may be the reason why it is beneficial for the reaction.And,CH2O adsorbs on the top site of one Pt atom,which is not affected by CO.So the reason why CO is not beneficial for the oxidation of methanol to HCOOH on the PtAu(111)surface is that CO adsorbed on the fcc site will prevent the second step,which makes it impossible to react to become CH2O and thus impossible to produce HCOOH.

    Table 1 Energetic data(in eV)for the reactions involved in methanol oxidation to formic acid on OH/PtAu(111)and CO/OH/PtAu(111)surfaces

    Fig.3 Schematic diagram of the production of HCOOH by methanol oxidation on OH/PtAu(111)and CO/OH/PtAu(111)surfaces

    Besides,we calculated the binding energy of CO/OH and coadsorption on PtAu(111)surface(Table 2),and the optimized adsorption configurations of CO and OH co-adsorbed on the PtAu(111)surface with CO at the top site,bridge site and hollow site(fcc site)are shown in Fig.4.Then,it can be found that CO prefers to adsorb on the fcc site,which will obviously reduce the adsorption energy of OH at fcc site.This is disadvantageous for the production of HCOOH by methanol oxidation,especially for the second step.However,when CO adsorbs at the top site,it will not affect the adsorption energy of OH obviously.In fact,it may conclude that the promotion effect may come from the weaken binding energy of OH due to the presence of CO species.This can be explained as following:as the R―H+OH→R+H2O belongs to the association reaction type,and the weaken binding energy of these two species,or at least one of the species,would lead to the high activity for the respective reaction.In our present study,the effect of CO on the binding energy of OH depends on its adsorption sites,either top site,bridge site,or the three-hollow site. From Table 2 we can find that the most significance case occurred when CO adsorbed at the fcc site,followed by bridge site,and the least one is the top site.So structure dependence of CO promotion is available.Moreover,to provide insight into the Pt―CO interaction with OH at the fcc site and CO at the top site,and OH at fcc site and CO at the fcc site,we first plot the projected density of states on d orbitals(d-PDOS)of Pt atoms(before and after CO adsorption on top and fcc site)in PtAu(111)configurations(shown in Fig.5),and then the d band center has been calculated.The dband centers are-1.70,-1.96,and-3.86 eV for the PtAu(111), CO(top)/PtAu(111)andCO(fcc)/PtAu(111)surfaces,respectively.It is consistent well with the order of d band center for the Au(111),CO(top)/Au(111),and CO(fcc)/Au(111)surfaces (shown in Fig.6).In general,the closer the d-band center to the fermi level,the more reactive the metal is.Therefore,the calculated d-band centers suggests that the adsorption energy of OH on PtAu(111)surface is slightly decreased with the presence of CO at the top site,whereas the OH adsorption is relatively large decreased with CO on the fcc site,compared with that on the clean PtAu(111)surface.This is consistent well with the adsorption energy of OH at fcc site(OH moves from initial top site to fcc site) with and without CO at the top and fcc sites on the PtAu(111) surface listed in Table 2.From the Bader charge analysis shown is Table 2,we find there is a repel interaction between CO and OH species owing to the same signs of charges for CO and OH in all investigated co-adsorption cases.

    Table 2 DFT results of binding energy of CO,OH and CO/OH co-adsorption on PtAu(111)surface

    3.2Effect of CO to oxidize the methanol to formic acid on Pt(111)surface

    Fig.4 Optimized adsorption configurations of CO and OH co-adsorbed on the PtAu(111)surface with CO at the top site(a), bridge site(b),and hollow site(fcc site)(c)

    In order to investigate the effect of CO,we also calculated methanol oxidation to formic acid on Pt(111)surface.The calculated binding energies of CO/OH and co-adsorption on PtAu (111)surface were calculated first,and the corresponding energetic data are listed in Table 3.From the calculation results,we found that OH is adsorbed stably at the fcc site with the adsorption energy of-2.46 eV and CO also prefers the fcc site with the adsorption energy being-1.76 eV,which agrees well with that on PtAu(111)surface.To investigate the methanol oxidation with production of formic acid better,the co-adsorption of OH and CO at different adsorption sites were investigated.The calculation results show that the adsorption energies of OH at the fcc site all decreased,with CO coadsorbed at the top,bridge,and fcc sites. Whereas the adsorption energy of OH in the initial top adsorption is not affected with the CO coadsorbed at the fcc site.It is valuable to point out that the OH moves to the fcc site from the initial top site in the co-adsorption case of OH(top)+CO(fcc).To provide insight into the Pt―CO interaction,the d-PDOS of Pt atoms (before and after CO adsorption on fcc and top site)was plotted (shown in Fig.7)and the d band center has been calculated.The d-band center is-2.19 eV for Pt(111),and it is-2.33 and-3.64 eV for Pt(111)surface with CO on the fcc and top sites.The calculated d-band centers indicate that the adsorption energy of OH at the fcc site of Pt(111)surface is slightly decreased with CO at fcc site,whereas it is relatively large decreased with CO on the top site,compared with that on the clean Pt(111)surface.This is agreement well with the adsorption energy of OH at the fcc site with and without CO at the fcc and top site on the Pt(111)surface shown in Table 3.As shown in Table 3,the Bader charge analysis indicates there is a repel interaction between CO and OH species because of the same signs of charges for CO and OH in all investigated co-adsorption cases.

    Fig.5 d-PDOS of the Pt atom in clean PtAu(111),CO(top)/ PtAu(111),and CO(fcc)/PtAu(111)configurations

    After determining the most stable co-adsorption of OH and CO, the relevant reactions involved in methanol oxidation to formic acid on Pt(111)surface were calculated.The optimized configurations are shown in Figs.8 and 9,and the corresponding energetic data are listed in Table 4 and Fig.10.

    Fig.6 d-PDOS of theAu atom in cleanAu(111),CO(top)/Au(111), and CO(fcc)/Au(111)configurations

    Table 3 DFT results of binding energy of CO,OH and CO/OH co-adsorption on Pt(111)surface

    For the first dehydrogenation step(CH3OH+OH→CH3O+H2O),on the surface of CO/Pt(111),it only releases 0.02 eV energy more than that on the Pt(111)surface.And from the activation energies shown in Table 3,we can see that this reaction is nonactivated with and without the presence of adsorbed CO.Therefore,this step is not chosen as an example to talk about the effect of CO.

    Fig.7 d-PDOS of the Pt atom in clean Pt(111),CO(fcc)/Pt(111), and CO(top)/Pt(111)configurations

    For methoxy dehydrogenation into CH2O and H,we find that this reaction is likely to occur on Pt(111),Au(111),and PtAu(111) surfaces with the barriers of 0.53,0.5616,and 0.99 eV,respectively. On Pt(111)surface,the barrier reduced to 0.15 eV from 0.53 eV and the reaction energy is more exothermic by 0.07 eV in the presence of CO,which indicates that the oxidative dehydrogenation of CH3O into CH2O by OH can be promoted by the preadsorbed CO both thermochemically and kinetically.On Pt(111) and CO/Pt(111)surfaces,there is not any obvious difference between the initial states and final states.We try to analyze the transition states and find the differences that bring the different barriers.At the TS,the breaking C―H bond is elongated to 0.129 and 0.127 nm from 0.110 nm on Pt(111)and CO/Pt(111)surfaces, respectively.The distances between the abstracted hydrogen and the O of OH are 0.130 nm on Pt(111)surface and 0.134 nm on CO/ Pt(111)surface.These both lower the barrier of methoxy dehydrogenation on Pt(111)surface in the presence of CO.

    Fig.8 Optimized adsorption configurations of the initial states,transition states,and final states involved in the methanol oxidation on the Pt(111)surface

    After CH2O is formed,it can be further oxidized by OH to CHO.As one part of initial state,CH2O is adsorbed on the surface with the C―O bond tilting at 59.4°and 59.6°from the surface normal on Pt(111)and CO/Pt(111)surfaces.We also computed the geometry reported by Greeley and Mavrikakis33that CH2O is adsorbed on Pt(111)surface in a“top-bridge-top”geometry through oxygen and carbon with the C―O bond almost paralleling to the surface,and found that it is only 0.09 eV more stable than it is adsorbed via oxygen atom at the top site of Pt.Herein we describe the CH2O dehydrogenation reaction with CH2O adsorbed at the top site of Pt through O atom as the initial state.Similar to CH3O dehydrogenation,there is also not any obvious difference between the initial states and final states on Pt(111)and CO/Pt (111)surfaces.The different barriers come from the configurationsof TSs.On Pt(111)surface,CHO adsorbs at the top site by carbon atom with the distance between the C atom and the abstracted H atom being 0.129 nm,which is longer than that of 0.114 nm on CO/Pt(111)surface.This needs more energy to catalyze the dehydrogenation of CH2O.Besides,the distance between the abstracted H and the O of OH is 0.173 nm on the CO/Pt(111)surface,which becomes longer than that of 0.138 nm on the Pt(111) surface.This also decreases the barrier of CH2O dehydrogenation on Pt(111)surface with the presence of CO.From the data listed in Table 3,we can see that the pre-adsorbed CO promotes the CH2O oxidation both thermochemically and kinetically:the barrier for CH2O dehydrogenation to CHO is 0.1 eV on CO/Pt(111) surface and it is 0.41 eV on the Pt(111)surface,and the reaction energy changes from-1.70 eV on Pt(111)to-1.76 eV on CO/ Pt(111)surface.

    For the fourth step,the barriers are both 0.2 eV on Pt(111)and CO/Pt(111)surfaces,which indicates that the CO has little effect on the reaction of CHO and OH.The energy of this reaction is-0.66 eVon Pt(111)surface and-0.70 eVon CO/Pt(111)surface, which means that it is also beneficial for this step with the existence of CO.

    Based on the computed energetic data listed above,we find that the existence of CO has little effect on the first and fourth steps, but it is beneficial for the second and the third steps.In the second and third steps,CO adsorbs on the fcc site,and it has the same effect as that on PtAu(111)and Au(111)surfaces,which is beneficial for OH to oxidize the methanol producing HCOOH. However,when CO adsorbs on the top site,it has not obvious promotion effect on methanol oxidation to formic acid.We calculated the second step on the Pt(111)surface with CO on the top site,and found that the barrier is 0.52 eV,which is 0.37 eV higher than that when CO adsorbs on the fcc site,but it is almost the same as that on the surface without the presence of CO.This result indicates that CO adsorbed on the fcc site will promote the reaction more than it on the top site.On this Pt(111)surface,the ratedetermining step is the second step,which is the same as that on PtAu(111)surface and different from that onAu(111)surface.This indicates that it is hard for the methanol to be oxidized to CH2O, which can deter the reaction to continue to form CO2on Pt(111) surface.

    Fig.9 Optimized adsorption configurations of the initial states,transition states,and final states involved in the methanol oxidation on the CO/Pt(111)surface

    3.3Comparison about effect of CO to oxidize the methanol to formic acid on Au(111)16,PtAu(111), and Pt(111)surfaces

    On Au(111)surface16,the activation energy barriers for the first three steps considered in this work generally decrease in the presence of CO adsorbed on the top site,particularly for the second and third steps.The third step or the conversion of CH2O into CHO species is the rate-controlling step owing to its relatively high energy barriers of 0.72 and 0.61 eV without and with the presence of CO.Therefore,methanol can just be oxidized to CH2O without the product of HCOOH on Au(111)surface at the low potential.On PtAu(111)surface,the energy barriers are decreased for the second and third steps with the CO adsorbed on the top site,which is the same as that on Au(111)surface.On this sur-face,the methoxy dehydrogenation is the rate-determining step, which means that the oxidation of methoxy to CH2O is difficult. On Pt(111)surface,the presence of CO is beneficial for the four steps,especially for the second and third steps,which is similar to that onAu(111)surface and PtAu(111)surface.Whereas the adsorption site of CO is fcc site on Pt(111)surface,which is different from the top adsorption of CO onAu(111)and PtAu(111)surfaces. Moreover,the rate-determining step is the oxidation of methoxy, which is the same as that on PtAu(111)surface and different from that onAu(111)surface.

    Based on the analysis above,we find that the adsorbed CO species can promote methanol oxidation to formic acid in aqueous alkaline media.Moreover,the promotion effect of CO is related to the adsorption site on different surfaces.It is the top adsorption on Au(111)and PtAu(111)surfaces,and the fcc adsorption on Pt(111)surface,which is owing to the higher stability of OH with CO adsorbed on the top sites of Au(111)and PtAu(111) surfaces,and on the fcc site of Pt(111)surface.For the PtAu(111) surface,the adsorption energy of OH at the fcc is not obviously decreased with the CO at the top site compared with other coadsorption cases.Whereas on the Pt(111)surface,OH adsorbed at the top site with the CO co-adsorbed at fcc site moves to the fcc site,and the corresponding adsorption energy is not decreased obviously in all the investigated co-adsorption cases.

    3.4Physical original of CO promotion

    Our present work indicates that CO can promotes the OH oxidative dehydrogenation of methanol on Au(111),PtAu(111),and Pt(111)surfaces,and the reaction mechanism should be determined.The ability of striping hydrogen atoms for OH species is related to its stability and basicity character34.In general,the more negative adsorption energy,the more stable adsorption.For Au (111)and PtAu(111)surfaces,the OH adsorption at the fcc site with the CO co-adsorbed at the top site is most stable,whereas OH adsorbed at the fcc site stably with CO at fcc site is most stable on Pt(111)surface,which promotes the methanol oxidation to formic acid.Moreover,the stronger basicity of OH species,the stronger promotion of X―H(X=CH3O,CH2O,CHO)bond breaking.CO2/ SO2is a typical electrophilic agent,and it is usually used to measure the basicity of surface O/OH species35,36.For methanol oxidation,hydrogen atom is one of the products.So,hydrogen atom is used as the probe atom to measure the OH basicity directly.And the basicity of OH species can be predicted by the reaction energy of X―H+OH→X+H2O,in which a more negative reaction energy indicates stronger basicity of OH species, and a stronger basicity means higher activity for the X―H bond breaking.More negative reaction energy of reaction,the stronger basicity of OH species.The reaction energy of each dehydrogenation step with the CO is generally more exothermic than that without CO,which indicates that the basicity of OH species is stronger with CO than that without CO.

    Table 4 Energetic data(in eV)for the reactions involved in methanol oxidation to formic acid on OH/Pt(111)and CO/OH/Pt(111)surfaces

    Fig.10 Schematic diagram of the production of HCOOH by methanol oxidation on OH/Pt(111)and CO/OH/Pt(111)surfaces

    4 Conclusions

    Our present DFT calculations demonstrate that the promotion effect of CO was acted on the reactions involved in methanol oxidation to formic acid on the PtAu(111)and Pt(111)surfaces, which is similar to that on Au(111)surface.The ability that CO could promote the methanol oxidation to formic acid is owing to the higher stability and stronger basicity of OH species induced by the co-adsorption of CO.The stability of OH and the basicity of OH species can be predicted by the adsorption energy and by the reaction energy of X―H+OH→X+H2O,respectively.More negative adsorption energy of OH,higher stability of OH.More negative reaction energy of reaction,the stronger basicity of OH species,and the higher activity for the X―H bond breaking. Considering the effect of CO on the different adsorption sites, such as top,bridge,and hollow,we found that the effect of CO on the methanol oxidation producing formic acid is different.On PtAu(111)surface,the CO adsorbed on the top site can promote the oxidation of methanol,which is the same as that on Au(111) surface,whereas the promotion effect only occurred with the fcc adsorption for CO on Pt(111)surface.

    References

    (1)Zhang,K.;Yang,W.;Ma,C.;Wang,Y.;Sun,C.W.;Chen,Y. J.;Duchesne,P.;Zhou,J.G.;Wang,J.;Hu,Y.F.;Banis,M.N.; Zhang,P.;Li,F.;Li,J.Q.;Chen,L.Q.NPG Asia Mater.2015, 7,e153.doi:10.1038/am.2014.122

    (2)Sun,S.;Zhang,G.;Gauquelin,N.;Chen,N.;Zhou,J.;Yang, S.;Chen,W.;Meng,X.;Geng,D.;Banis,M.N.;Li,R.;Ye,S.; Knights,S.;Botton,G.A.;Sham,T.;Sun,X.Sci.Rep.2013,3, 1775.doi:10.1038/srep01775

    (3)Ganesan,R.;Lee,J.S.Angew.Chem.Int.Edit.2005,44,6557. doi:10.1002/anie.200501272

    (4)Wasmus,S.;Küver,A.J.Electroanal.Chem.1999,461,14. doi:10.1016/S0022-0728(98)00197-1

    (5)Roberts,J.L.,Jr.;Sawyer,D.T.Electrochim.Acta 1965,10, 989.doi:10.1016/0013-4686(65)80011-1

    (6)Rodriguez,P.;Kwon,Y.;Koper,M.T.M.Nat.Chem.2012,4, 177.doi:10.1038/nchem.1221

    (7)Zope,B.N.;Hibbits,D.D.;Neurock,M.;Davis,R.J.Science 2010,330,74.doi:10.1126/science.1195055

    (8)Ketchie,W.C.;Fang,Y.L.;Wong,M.S.;Murayama,M.; Davis,R.J.J.Catal.2007,250,94.doi:10.1016/j. jcat.2007.06.001

    (9)Burke,L.D.;Nugent,P.F.Gold Bull.1998,31,39.doi: 10.1007/BF03214760

    (10)Kita,H.;Nakajima,H.;Hayashi,K.J.Electroanal.Chem. 1985,190,141.doi:10.1016/0022-0728(85)80083-8

    (11)Zhang,T.F.;Liu,Z.P.;Driver,S.M.;Pratt,S.J.;Jenkins,S. J.;King,D.A.Phys.Rev.Lett.2005,95,266102.doi:10.1103/ PhysRevLett.95.266102

    (12)Gan,L.Y.;Zhao,Y.J.J.Chem.Phys.2010,133,094703.doi: 10.1063/1.3483235

    (13)Liu,Y.H.;Wei,L.;Hu,Y.Z.;Huang,X.Y.;Wang,J.Q.;Li,J. Q.;Hu,X.L.;Zhuang,N.F.J.Alloy.Compd.2016,656,452. doi:10.1016/j.jallcom.2015.10.004

    (14)Shubina,T.E.;Hartnig,C.;Koper,M.T.M.Phys.Chem. Chem.Phys.2004,6,4215.doi:10.1039/B407669A

    (15)Wang,L.;He,C.Z.;Zhang,W.H.;Li,Z.Y.;Yang,J.L. J.Phys.Chem.C 2014,118,17511.doi:10.1021/jp501620h

    (16)Lv,C.Q.;Liu,J.H.;Wang,H.;Wang,G.C.Catal.Commun. 2015,60,60.doi:10.1016/j.catcom.2014.11.013

    (17)Yuan,D.W.;Gong,X.G.;Wu,R.Q.J.Chem.Phys.2008, 128,064706.doi:10.1063/1.2835545

    (18)Zhong,W.H.;Liu,Y.X.;Zhang,D.J.J.Phys.Chem.C 2012, 116,2994.doi:10.1021/jp210304z

    (19)Montero,M.A.;Gennero de Chialvo,M.R.;Chialvo,A.C. Int.J.Hydrog.Energy 2011,36,3811.doi:10.1016/j. ijhydene.2010.12.115

    (20)Ren,H.;Humbert,M.P.;Menning,C.A.;Chen,J.G.;Shu,Y. Y.;Singh,U.G.;Cheng,W.C.Appl.Catal.A 2010,375,303. doi:10.1016/j.apcata.2010.01.018

    (21)Qiu,C.C.;Guo,Y.G.;Zhang,J.T.;Ma,H.Y.;Cai,Y.Q. Mater.Chem.Phys.2011,127,484.doi:10.1016/j. matchemphys.2011.02.041

    (22)Xu,J.B.;Zhao,T.S.;Yang,W.W.;Shen,S.Y.Int.J.Hydrog. Energy 2010,35,8699.doi:10.1016/j.ijhydene.2010.05.008

    (23)Yang,L.;Chen,J.H.;Zhong,X.X.;Cui,K.Z.;Xu,Y.;Kuang, Y.F.Colloids Surf.A 2007,295,21.doi:10.1016/j. colsurfa.2006.08.023

    (24)Kresse,G.;Hafner,J.Phys.Rev.B 1994,49,14251.doi: 10.1103/PhysRevB.49.14251

    (25)Kresse,G.;Furthmüller,J.Comp.Mater.Sci.1996,6,15.doi: 10.1016/0927-0256(96)00008-0

    (26)Perdew,J.P.;Chevary,J.A.;Vosko,S.H.;Jackson,K.A.; Pederson,M.R.;Singh,D.J.;Fiolhais,C.Phys.Rev.B 1992, 46,6671.doi:10.1103/PhysRevB.46.6671

    (27)Kresse,G.;Joubert,D.Phys.Rev.B 1999,59,1758.doi: 10.1103/PhysRevB.59.1758

    (28)Bl?hl,P.E.Phys.Rev.B 1994,50,17953.doi:10.1103/ PhysRevB.50.17953

    (29)Monkhorst,H.J.;Pack,J.D.Phys.Rev.B 1976,13,5188.doi: 10.1103/PhysRevB.13.5188

    (30)Henkelman,G.;Uberuaga,B.P.;Jónsson,H.J.Chem.Phys. 2000,113,9901.doi:10.1063/1.1329672

    (31)Kittel,C.Introduction to Solid State Physics,8th ed.;Wiley: New York,2004

    (32)Hammer,B.Surf.Sci.2000,459,323.doi:10.1016/S0039-6028(00)00467-2

    (33)Greeley,J.;Mavrikakis,M.J.Am.Chem.Soc.2004,126, 3910.doi:10.1021/ja037700z

    (34)Pacchioni,G.;Ricart,J.M.;Illas,F.J.Am.Chem.Soc.1994, 116,10152.doi:10.1021/ja00101a038

    (35)Torres,D.;Lopez,N.;Illas,F.;Lambert,R.M.Angew.Chem. Int.Edit.2007,46,2055.doi:10.1002/ange.200603803

    (36)Zhao,S.;Ma,X.D.;Pang,Q.;Sun,H.W.;Wang,G.C.Phys. Chem.Chem.Phys.2014,16,5553.doi:10.1039/C3CP55048F

    First-Principles Study of Effect of CO to Oxidize Methanol to Formic Acid in Alkaline Media on PtAu(111)and Pt(111)Surfaces

    LIU Jian-Hong1Lü Cun-Qin1,2JIN Chun1,*WANG Gui-Chang3,*
    (1College of Chemistry and Enviromental Engineering,Shanxi Datong University,Datong 037009,Shanxi Province,P.R.China;2Key LaboratoryofAdvancedEnergyMaterialsChemistry(Ministry ofEducation),Nankai University,Tianjin300071,P.R.China;3Department of Chemistry and the Tianjin Key Laboratory of Metal and Molecule-based Material Chemistry, Nankai University,Tianjin 300071,P.R.China)

    Density functional theory calculations have been performed to investigate methanol oxidation to formic acid on PtAu(111)and Pt(111)surfaces with and without CO in alkaline media.The calculated results show that the pre-adsorbed CO species promotes almost every step involved in the oxidation of methanol on PtAu(111)and Pt(111)surfaces,which is similar to that observed on aAu(111)surface.These findings may be attributed to the relatively high stability and strong basicity of the OH species induced by the adsorption of CO, and the enhanced ability to strip the H atoms.

    Methanol oxidation;M(111)(M=PtAu,Pt);Alkaline medium;CO promotion effect; Density functional theory calculation

    1 Introduction

    Given that methanol exhibits their distinctive advantages of high energy density and the ready availability,inexpensiveness, and the safe storage and transportation potential as a liquid fuel, direct methanol fuel cells have attracted much attention for portable power applications1-4.As an intermediate species producedin methanol oxidation,carbon monoxide has a strong tendency to adsorb on the catalyst surface to hinder the sites needed by reactants in electrocatalysis,which makes the adsorbed CO act evidently as a poison or poisoning intermediate in methanol oxidation.Gold is known to be an excellent electrocatalyst for carbon monoxide oxidation5,and alkaline conditions are necessary for efficient aqueous-phase electrocatalysis6-10.Rodriguez et al.6investigated the oxidation of several alcohols,including methanol, on cleanAu(111)and CO-modifiedAu(111)electrodes in aqueous alkaline media,and found that the CO can promote alcohol oxidation in a solution.This can be attributed to the strong bonding of OH to the surface induced by adsorbed CO.Moreover,the adsorbed OH may act as the oxidant for CO oxidation,leading to a significant rate enhancement.Besides the co-adsorption of CO and OH,the effect of other co-adsorption systems on the reaction has been reported.Zhang et al.11observed an interesting stabilizing co-adsorption system of CO and electronegative NO2on Au(111) surface in ultra-high vacuum(UHV);they proposed that such stabilizing system may be responsible for catalytic promotion, particularly for the association reaction in IB metals.Gan and Zhao12investigated the interaction of the co-adsorption of CO and different adatoms(such as Na,S,O,and Cl)on Cu(111),Ag(111), Au(111),and Pd(111)surfaces,and found that the enhancement of CO adsorption by S originates from S-induced positive polarization of Au and Ag surfaces.The recent experimental results show that the CO has an anti-poisoning effect on graphene supported PtPd alloy for the ethanol oxidation reaction13.

    September 28,2015;Revised:January 18,2016;Published on Web:January 19,2016.*Corresponding authors.WANG Gui-Chang,Email:wangguichang@nankai.edu.cn;Tel:+86-22-23503824. JIN Chun,Email:jinchun0828@126.com;Tel:+86-352-7563093. The project was supported by the National Natural Science Foundation of China(21503122,21346002),Natural Science Foundation of Shanxi for Youths,China(2014021016-2),Scientific and Technological Programs in Shanxi Province,China(2015031017),and Foundation of Key Laboratory ofAdvanced Energy Materials Chemistry(Ministry of Education),China.

    O641;O647

    10.3866/PKU.WHXB201601191

    國(guó)家自然科學(xué)基金(21503122,21346002),山西省青年科技研究基金(2014021016-2),山西省科技攻關(guān)項(xiàng)目(2015031017)及先進(jìn)能源材料化學(xué)教育部重點(diǎn)實(shí)驗(yàn)室開放基金資助

    猜你喜歡
    山西大同南開大學(xué)堿性
    山西大同 黃花菜豐收在望
    《山西大同大學(xué)學(xué)報(bào)(自然科學(xué)版)》征稿簡(jiǎn)則
    山西大同大學(xué)采礦研究所簡(jiǎn)介
    山西大同邀客共賞“小黃花大產(chǎn)業(yè)”
    南開大學(xué)制備新型超強(qiáng)韌人造蜘蛛絲
    一道南開大學(xué)自主招生試題的推廣
    堿性磷酸酶鈣-鈷法染色的不同包埋方法比較
    堿性土壤有效磷測(cè)定的影響因素及其控制
    堿性溶液中鉑、鈀和金析氧性能比較
    堿性介質(zhì)中甲醇在PdMo/MWCNT上的電化學(xué)氧化
    欧美激情久久久久久爽电影 | 久久精品91无色码中文字幕| 日本vs欧美在线观看视频| 高清在线国产一区| 久久久久国产精品人妻aⅴ院| 成人av一区二区三区在线看| 午夜福利,免费看| 亚洲七黄色美女视频| aaaaa片日本免费| 日日夜夜操网爽| 久久精品91蜜桃| 国产一区在线观看成人免费| or卡值多少钱| 亚洲人成电影免费在线| 真人做人爱边吃奶动态| ponron亚洲| 久久精品国产清高在天天线| 嫁个100分男人电影在线观看| 免费久久久久久久精品成人欧美视频| 777久久人妻少妇嫩草av网站| 日韩成人在线观看一区二区三区| 国产亚洲精品久久久久5区| 热99re8久久精品国产| 亚洲欧美日韩高清在线视频| 国产精品综合久久久久久久免费 | 成人三级做爰电影| 成人国产综合亚洲| 久久狼人影院| 国产精品久久视频播放| xxx96com| 丝袜在线中文字幕| 亚洲av日韩精品久久久久久密| 国产欧美日韩一区二区精品| 中文字幕高清在线视频| 亚洲少妇的诱惑av| 男女之事视频高清在线观看| 亚洲男人的天堂狠狠| 久久精品亚洲精品国产色婷小说| 国产在线观看jvid| 久久亚洲真实| 久久国产精品人妻蜜桃| 国产成人欧美| 欧美日本中文国产一区发布| 久久久精品国产亚洲av高清涩受| 国产伦一二天堂av在线观看| 精品国产一区二区三区四区第35| 最近最新免费中文字幕在线| 两个人看的免费小视频| 少妇粗大呻吟视频| or卡值多少钱| 色av中文字幕| 手机成人av网站| 午夜福利欧美成人| 丝袜人妻中文字幕| 人妻丰满熟妇av一区二区三区| 国产一级毛片七仙女欲春2 | 国产成人av教育| 淫妇啪啪啪对白视频| 可以在线观看毛片的网站| 亚洲三区欧美一区| 国产精品久久久久久亚洲av鲁大| 久久久久久国产a免费观看| 久9热在线精品视频| 91老司机精品| 乱人伦中国视频| 露出奶头的视频| 美女 人体艺术 gogo| 麻豆一二三区av精品| 大型黄色视频在线免费观看| 极品教师在线免费播放| 亚洲av电影不卡..在线观看| 亚洲人成网站在线播放欧美日韩| 91麻豆精品激情在线观看国产| 老鸭窝网址在线观看| 日本三级黄在线观看| 国产精品1区2区在线观看.| 日韩欧美一区视频在线观看| 国产精华一区二区三区| 国产激情欧美一区二区| 国产视频一区二区在线看| 久久精品91蜜桃| 久久久久久久久久久久大奶| 电影成人av| 99国产精品免费福利视频| 久久中文看片网| 亚洲人成77777在线视频| 久久香蕉国产精品| 国产精品av久久久久免费| 免费观看人在逋| 免费高清在线观看日韩| www.www免费av| 非洲黑人性xxxx精品又粗又长| 国产精品秋霞免费鲁丝片| 可以在线观看毛片的网站| 国产精品乱码一区二三区的特点 | 亚洲三区欧美一区| 国产视频一区二区在线看| 亚洲精品中文字幕一二三四区| 狠狠狠狠99中文字幕| 国产熟女午夜一区二区三区| 欧美成人免费av一区二区三区| 精品国内亚洲2022精品成人| www日本在线高清视频| 亚洲国产欧美日韩在线播放| 精品人妻在线不人妻| 欧美av亚洲av综合av国产av| 国产免费av片在线观看野外av| cao死你这个sao货| 男人操女人黄网站| 国产日韩一区二区三区精品不卡| 老司机靠b影院| 美女高潮到喷水免费观看| 三级毛片av免费| 久久精品国产综合久久久| 亚洲国产精品999在线| 色播亚洲综合网| 丰满人妻熟妇乱又伦精品不卡| 亚洲熟女毛片儿| 一二三四社区在线视频社区8| 可以在线观看的亚洲视频| 亚洲午夜理论影院| 国产xxxxx性猛交| 色精品久久人妻99蜜桃| 亚洲在线自拍视频| 国产精品综合久久久久久久免费 | 午夜免费鲁丝| 亚洲国产精品999在线| 亚洲国产欧美一区二区综合| 成人欧美大片| 免费少妇av软件| 日本 欧美在线| 国产在线精品亚洲第一网站| av天堂在线播放| 狠狠狠狠99中文字幕| 亚洲第一电影网av| 日韩欧美国产一区二区入口| 中亚洲国语对白在线视频| 欧美久久黑人一区二区| 久久国产乱子伦精品免费另类| 色播亚洲综合网| xxx96com| 丁香欧美五月| 真人一进一出gif抽搐免费| 看黄色毛片网站| 美女大奶头视频| 国产精品久久久久久人妻精品电影| 天堂影院成人在线观看| 国产精品影院久久| 嫁个100分男人电影在线观看| 日日摸夜夜添夜夜添小说| 最好的美女福利视频网| 一本久久中文字幕| 国产成+人综合+亚洲专区| 国产在线观看jvid| 中文字幕人妻熟女乱码| 老熟妇乱子伦视频在线观看| 琪琪午夜伦伦电影理论片6080| 欧美中文日本在线观看视频| 丝袜在线中文字幕| 一个人观看的视频www高清免费观看 | 欧美不卡视频在线免费观看 | 久久精品aⅴ一区二区三区四区| 午夜免费激情av| 男男h啪啪无遮挡| 搡老岳熟女国产| 久久精品成人免费网站| 91九色精品人成在线观看| 男女下面进入的视频免费午夜 | x7x7x7水蜜桃| av免费在线观看网站| 亚洲精品美女久久av网站| 亚洲成a人片在线一区二区| 女性被躁到高潮视频| 欧美激情 高清一区二区三区| bbb黄色大片| 欧美性长视频在线观看| 国产午夜精品久久久久久| 女人精品久久久久毛片| 欧美日韩黄片免| 国内毛片毛片毛片毛片毛片| 免费在线观看日本一区| 久热这里只有精品99| 波多野结衣一区麻豆| av电影中文网址| 国产av一区二区精品久久| 午夜精品久久久久久毛片777| 熟妇人妻久久中文字幕3abv| www.999成人在线观看| 亚洲成a人片在线一区二区| 国产精品1区2区在线观看.| 国产精品九九99| 精品欧美一区二区三区在线| 国产99白浆流出| 国内毛片毛片毛片毛片毛片| 成人欧美大片| 宅男免费午夜| 欧美一级毛片孕妇| 亚洲专区中文字幕在线| 99riav亚洲国产免费| 99久久99久久久精品蜜桃| 午夜福利一区二区在线看| 国产欧美日韩综合在线一区二区| 欧美激情久久久久久爽电影 | 欧美 亚洲 国产 日韩一| 精品乱码久久久久久99久播| 色综合欧美亚洲国产小说| 老汉色∧v一级毛片| 免费女性裸体啪啪无遮挡网站| 久久热在线av| 一级a爱片免费观看的视频| 禁无遮挡网站| 又黄又爽又免费观看的视频| 国产精品爽爽va在线观看网站 | 国产1区2区3区精品| 午夜a级毛片| 日本一区二区免费在线视频| 日韩中文字幕欧美一区二区| 欧美激情高清一区二区三区| 看片在线看免费视频| 一区二区日韩欧美中文字幕| 两人在一起打扑克的视频| 女同久久另类99精品国产91| 男女下面插进去视频免费观看| 我的亚洲天堂| 婷婷六月久久综合丁香| 可以免费在线观看a视频的电影网站| 国产99白浆流出| 亚洲精品中文字幕一二三四区| 日韩大尺度精品在线看网址 | 午夜免费鲁丝| 久久天堂一区二区三区四区| 成人国语在线视频| 自线自在国产av| 免费看十八禁软件| 女性被躁到高潮视频| 久久香蕉国产精品| 淫妇啪啪啪对白视频| 久久国产精品男人的天堂亚洲| 69精品国产乱码久久久| 国产免费男女视频| 一二三四社区在线视频社区8| 婷婷丁香在线五月| 国产精品影院久久| 欧美日韩一级在线毛片| 国产精品1区2区在线观看.| 嫩草影院精品99| 制服人妻中文乱码| 男人舔女人下体高潮全视频| 国产单亲对白刺激| 亚洲专区中文字幕在线| 真人一进一出gif抽搐免费| 欧美不卡视频在线免费观看 | 十八禁网站免费在线| 人妻久久中文字幕网| 嫩草影视91久久| 丝袜人妻中文字幕| 女性被躁到高潮视频| 久久久国产精品麻豆| 91大片在线观看| 亚洲成人精品中文字幕电影| 日本一区二区免费在线视频| 悠悠久久av| 电影成人av| 美国免费a级毛片| 亚洲精品美女久久久久99蜜臀| 亚洲国产中文字幕在线视频| 欧美在线黄色| 国产单亲对白刺激| 人人妻人人澡欧美一区二区 | 最好的美女福利视频网| 黄色片一级片一级黄色片| 亚洲五月色婷婷综合| 亚洲九九香蕉| 日日爽夜夜爽网站| 男女床上黄色一级片免费看| 在线播放国产精品三级| 在线观看免费视频网站a站| 非洲黑人性xxxx精品又粗又长| 日本黄色视频三级网站网址| 一级,二级,三级黄色视频| 国产成+人综合+亚洲专区| 国产一区二区在线av高清观看| 国产精品 欧美亚洲| 亚洲欧美日韩无卡精品| 最新在线观看一区二区三区| 18禁观看日本| 后天国语完整版免费观看| 熟妇人妻久久中文字幕3abv| 宅男免费午夜| 一卡2卡三卡四卡精品乱码亚洲| 午夜视频精品福利| 51午夜福利影视在线观看| 久久人妻av系列| 午夜a级毛片| 亚洲一区高清亚洲精品| 国产亚洲精品第一综合不卡| 国产欧美日韩综合在线一区二区| 一区二区三区高清视频在线| 欧美人与性动交α欧美精品济南到| av免费在线观看网站| 51午夜福利影视在线观看| 99久久99久久久精品蜜桃| 久久精品影院6| 制服人妻中文乱码| 18禁国产床啪视频网站| 级片在线观看| 国产不卡一卡二| 两人在一起打扑克的视频| 纯流量卡能插随身wifi吗| 黄片大片在线免费观看| 免费无遮挡裸体视频| 日韩欧美免费精品| 国产亚洲精品一区二区www| 亚洲精品粉嫩美女一区| 日日爽夜夜爽网站| 国产一区二区激情短视频| 久久人人爽av亚洲精品天堂| 黄片小视频在线播放| 日本三级黄在线观看| 久久国产精品男人的天堂亚洲| 99久久精品国产亚洲精品| 悠悠久久av| 给我免费播放毛片高清在线观看| 在线观看www视频免费| 国产精品久久视频播放| 侵犯人妻中文字幕一二三四区| 狠狠狠狠99中文字幕| av电影中文网址| 亚洲av电影在线进入| 看免费av毛片| 欧美日本中文国产一区发布| 亚洲精品av麻豆狂野| 99国产精品免费福利视频| 最近最新中文字幕大全电影3 | 国产亚洲欧美在线一区二区| 中文字幕高清在线视频| or卡值多少钱| 美国免费a级毛片| 久久久精品欧美日韩精品| 日日夜夜操网爽| 国产精品一区二区三区四区久久 | 麻豆一二三区av精品| 国产一区二区三区综合在线观看| 国产激情欧美一区二区| 国产三级黄色录像| 精品高清国产在线一区| 91成年电影在线观看| 日日爽夜夜爽网站| 亚洲无线在线观看| 在线观看舔阴道视频| 色综合站精品国产| 精品国产国语对白av| 90打野战视频偷拍视频| netflix在线观看网站| 久久精品影院6| 精品国产乱子伦一区二区三区| 免费在线观看完整版高清| 久久久久精品国产欧美久久久| 国产主播在线观看一区二区| 国产精品 欧美亚洲| 日本 av在线| 国产精品久久电影中文字幕| 他把我摸到了高潮在线观看| 国产精品一区二区在线不卡| 亚洲成av片中文字幕在线观看| 久久久久久久久免费视频了| 长腿黑丝高跟| 日日摸夜夜添夜夜添小说| 日韩三级视频一区二区三区| 一级毛片高清免费大全| 天天添夜夜摸| 欧美 亚洲 国产 日韩一| 日本免费a在线| 国产亚洲精品综合一区在线观看 | 一二三四社区在线视频社区8| 50天的宝宝边吃奶边哭怎么回事| 久久香蕉国产精品| 欧美日本视频| 亚洲av第一区精品v没综合| 欧美乱色亚洲激情| а√天堂www在线а√下载| 不卡av一区二区三区| 男女之事视频高清在线观看| 老司机在亚洲福利影院| 亚洲人成电影观看| 亚洲精品在线美女| 在线观看免费午夜福利视频| 久久久久亚洲av毛片大全| 97超级碰碰碰精品色视频在线观看| 国产av在哪里看| 免费在线观看黄色视频的| 亚洲 国产 在线| 精品一区二区三区四区五区乱码| 亚洲美女黄片视频| 校园春色视频在线观看| 在线永久观看黄色视频| 国产在线精品亚洲第一网站| 国产三级黄色录像| 国产亚洲精品一区二区www| 国产亚洲精品久久久久5区| 亚洲最大成人中文| 欧美成人一区二区免费高清观看 | 性少妇av在线| av中文乱码字幕在线| 首页视频小说图片口味搜索| 日韩精品免费视频一区二区三区| 正在播放国产对白刺激| 老汉色av国产亚洲站长工具| 色播在线永久视频| 亚洲欧美一区二区三区黑人| 亚洲av五月六月丁香网| 91国产中文字幕| 免费在线观看影片大全网站| x7x7x7水蜜桃| 啦啦啦免费观看视频1| 午夜成年电影在线免费观看| 99热只有精品国产| 麻豆一二三区av精品| 亚洲性夜色夜夜综合| 老司机深夜福利视频在线观看| 亚洲av第一区精品v没综合| 亚洲av五月六月丁香网| 91国产中文字幕| 国产高清视频在线播放一区| 波多野结衣av一区二区av| 精品日产1卡2卡| 久久婷婷人人爽人人干人人爱 | 国产一区二区三区在线臀色熟女| 国产成+人综合+亚洲专区| 大型av网站在线播放| 变态另类成人亚洲欧美熟女 | 国产野战对白在线观看| 日韩高清综合在线| 国产午夜福利久久久久久| 一区在线观看完整版| 精品电影一区二区在线| 国产97色在线日韩免费| 精品国产一区二区三区四区第35| 九色国产91popny在线| 啦啦啦 在线观看视频| 国产成人精品久久二区二区免费| а√天堂www在线а√下载| 国产精品98久久久久久宅男小说| 亚洲成av人片免费观看| 嫩草影视91久久| 中文字幕色久视频| 国产精品久久久人人做人人爽| 我的亚洲天堂| 久久中文字幕人妻熟女| 777久久人妻少妇嫩草av网站| 午夜免费激情av| 国产欧美日韩一区二区三区在线| 少妇的丰满在线观看| 日韩欧美一区视频在线观看| 人妻丰满熟妇av一区二区三区| 国产成人精品在线电影| 岛国在线观看网站| 一区二区三区精品91| 亚洲全国av大片| 最新在线观看一区二区三区| 少妇被粗大的猛进出69影院| 啦啦啦观看免费观看视频高清 | 777久久人妻少妇嫩草av网站| 午夜免费激情av| 18禁观看日本| 人妻久久中文字幕网| 久久中文看片网| 免费一级毛片在线播放高清视频 | 亚洲成人免费电影在线观看| 一进一出抽搐动态| 精品国产一区二区三区四区第35| 两个人看的免费小视频| 免费高清在线观看日韩| 国语自产精品视频在线第100页| 禁无遮挡网站| 最好的美女福利视频网| 黄色丝袜av网址大全| 在线天堂中文资源库| 亚洲专区字幕在线| 欧美成人一区二区免费高清观看 | 大码成人一级视频| 99热只有精品国产| 在线观看舔阴道视频| 久久亚洲真实| ponron亚洲| 久久久久国产一级毛片高清牌| 少妇被粗大的猛进出69影院| 看片在线看免费视频| 亚洲七黄色美女视频| 亚洲国产高清在线一区二区三 | √禁漫天堂资源中文www| 亚洲精品美女久久久久99蜜臀| 女人爽到高潮嗷嗷叫在线视频| 老鸭窝网址在线观看| 黄色 视频免费看| 亚洲成a人片在线一区二区| 精品人妻在线不人妻| 亚洲 欧美一区二区三区| 91av网站免费观看| 日韩中文字幕欧美一区二区| 97超级碰碰碰精品色视频在线观看| 日韩欧美一区视频在线观看| 国产精品野战在线观看| 两个人视频免费观看高清| 亚洲精品在线美女| 久久久国产成人精品二区| 欧美午夜高清在线| 精品国产国语对白av| 色精品久久人妻99蜜桃| 好男人电影高清在线观看| 午夜福利成人在线免费观看| 欧美日本中文国产一区发布| www日本在线高清视频| tocl精华| av天堂在线播放| 在线播放国产精品三级| 露出奶头的视频| 国产精品1区2区在线观看.| 操出白浆在线播放| 日韩精品免费视频一区二区三区| 可以在线观看毛片的网站| 最近最新中文字幕大全免费视频| 国产一区二区三区视频了| 性欧美人与动物交配| 久久中文看片网| 亚洲国产日韩欧美精品在线观看 | 91老司机精品| 性色av乱码一区二区三区2| 亚洲电影在线观看av| 女人精品久久久久毛片| 日韩欧美三级三区| 搞女人的毛片| 亚洲色图 男人天堂 中文字幕| 欧美日本亚洲视频在线播放| xxx96com| 18禁美女被吸乳视频| 999精品在线视频| 中出人妻视频一区二区| 久久精品影院6| 女性被躁到高潮视频| 别揉我奶头~嗯~啊~动态视频| 精品国产一区二区三区四区第35| 日日爽夜夜爽网站| 亚洲精华国产精华精| 日韩欧美国产一区二区入口| 亚洲精品国产色婷婷电影| 国产一区二区三区在线臀色熟女| 精品少妇一区二区三区视频日本电影| 亚洲成a人片在线一区二区| 国产精品免费一区二区三区在线| 国产成+人综合+亚洲专区| 欧美乱妇无乱码| 巨乳人妻的诱惑在线观看| 久久青草综合色| cao死你这个sao货| 伊人久久大香线蕉亚洲五| √禁漫天堂资源中文www| 91成年电影在线观看| 国产亚洲av高清不卡| 97超级碰碰碰精品色视频在线观看| 欧美老熟妇乱子伦牲交| 一级毛片高清免费大全| 一边摸一边抽搐一进一小说| 久久人人精品亚洲av| 亚洲在线自拍视频| 99国产精品一区二区蜜桃av| 一进一出抽搐gif免费好疼| 校园春色视频在线观看| 久久久国产精品麻豆| 一边摸一边做爽爽视频免费| 首页视频小说图片口味搜索| 亚洲av成人av| 亚洲最大成人中文| 女警被强在线播放| 久久精品国产清高在天天线| 熟妇人妻久久中文字幕3abv| 99香蕉大伊视频| 午夜福利一区二区在线看| 欧美乱码精品一区二区三区| 亚洲成av人片免费观看| 在线观看一区二区三区| 丝袜美足系列| 一区二区三区精品91| 国产xxxxx性猛交| 成人精品一区二区免费| 日韩大码丰满熟妇| 搡老岳熟女国产| 精品久久久久久久人妻蜜臀av | 又黄又爽又免费观看的视频| 国产成人av教育| 国产精品免费一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久| 日韩有码中文字幕| 国产欧美日韩一区二区三| 久久精品亚洲精品国产色婷小说| 午夜精品久久久久久毛片777| 久久精品国产亚洲av香蕉五月| 香蕉久久夜色| 国产欧美日韩一区二区精品| 久久婷婷人人爽人人干人人爱 | 午夜免费激情av| 国产xxxxx性猛交| 久久久久久久久中文| 久久久久九九精品影院| 丁香欧美五月| 男男h啪啪无遮挡| 91国产中文字幕| 一区福利在线观看| 国产成人精品久久二区二区免费| 青草久久国产| 国产高清激情床上av| 午夜福利影视在线免费观看| 中文字幕人妻熟女乱码| 国产91精品成人一区二区三区| 女性被躁到高潮视频|