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

    Identify fine microstructure of multifarious iron oxides via O K-edge EELS spectra

    2022-09-15 03:11:20JunnanChenYujieQiMingLuYimingNiuBingsenZhang
    Chinese Chemical Letters 2022年9期

    Junnan Chen, Yujie Qi, Ming Lu, Yiming Niu, Bingsen Zhang,?

    a Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016, China

    bSchool of Materials Science and Engineering, University of Science and Technology of China, Shenyang 110016, China

    c Key Laboratory of Functional Materials Physics and Chemistry of the Ministry of Education, Jilin Normal University, Changchun 130103, China

    ABSTRACT Relying on the electron energy loss spectrum (EELS) of metallic elements to obtain microstructure analysis is an investigation method of the reaction mechanisms of transition metal oxides (TMOs) in catalysis,energy storage and conversion.However, the low signal from K shell owing to insufficient electron beam energy, and the complicated electronic structure in L shell of the metal element restrict the analysis of the coordination environment of the TMOs.Herein, density functional theory (DFT) calculation, Fourier transform (FT) and wavelet transform (WT) were employed to probe the relationship between the four individual peaks in O K-edge spectra of iron oxides and the microstructure information (chemical bonds and atomic coordination).The findings show that the peak amplitude ration is in a linear correlation with the valence state of Fe element, and that the coordination number obtained by radial distribution function (RDF) is favorably linearly correlative with that from the standard coordination structure model.As a result, the quantitative analysis on the change of valence state and atomic coordination in microstructure can be realized by EELS O K-edge spectra.This study establishes EELS O K-edge spectrum as a promising pathway to quantitatively analyze the valence state and atomic coordination information of TMOs,and offers an effective method to conduct microstructure analysis via the EELS spectra of the non-metal element.

    Keywords:EELS ELNES EXELFS Transition metal oxides DFT Wavelet analyze

    Transition metal oxides (TMOs) have been received widespread attention as catalysts or electrode materials with the merits of various valance states and multiphase transition [1–4].To explore the reaction kinetics of a catalytic and electrode process, the analysis on microstructural transformation of TMOs is absolutely essential [5,6].In general, XRD andin-situXRD are utilized to obtain the phase transformation information.X-ray photoelectron spectroscopy (XPS) and X-ray absorption fine structure (XAFS) are employed to investigate the evolution of valence state and coordination environment for characteristic elements of TMOs.But these approaches cannot satisfy the collection of the local structure signal of TMOs [7,8].Electron energy-loss spectroscopy (EELS) based on transmission electron microscope (TEM) platform could realize the local acquisition of the chemical information on distribution,valence state and coordination environment of characteristic elements of TMOs with a high spatial resolution [9].For instance,the phase transformation behavior of V2O5underin-situvacuum heating and the strong metal-support interaction in Au/TiO2catalysts have been studied by the EELS spectrum to reveal the process mechanisms [10–12].However, the combination of low signal from K shell and complicated electron structure in L shell of the metal element of TMOs restrict the analysis of coordination environmentviacollecting the EELS spectra of metal element [13].Further, a limitation appears on the exploration of the microstructural change of TMOs.Due to the strong bonding between the transition metal atom and the O atom in TMOs, the O EELS spectrum can be collected to realize a comprehensive investigation of the local microstructure of TMOs.

    Herein, based on O K near-edge structure, the O K-edge spectra of iron oxides (FeO, Fe3O4,γ-Fe2O3,α-Fe2O3) have been collected as examples to model the reaction process at different stages.The Fourier transform (FT) and wavelet transform (WT) methods have also been used to analyze O K-edge extended energy-loss fine structures (EXELFS), to construct the relationship between characteristic peaks and microstructure, which provides an effective route to reveal the fine microstructure of TMOsviathe EELS spectra of the non-metal element.

    The O K-edge and Fe L-edge spectra of iron oxides (FeO, Fe3O4,γ-Fe2O3,α-Fe2O3) are shown in Fig.1.There are two peaks in L-edge spectra of iron oxides in Fig.1b, which are located in the ranges of 710?720 eV and 725?735 eV, respectively.This phenomenon originates from Fe electrons transition from 2p3/2and 2p1/2orbitals to unoccupied 3d states.With the increase of valence state for Fe in iron oxides, the L2and L3peaks move toward higher energy and the white-line ratio of Fe L2to Fe L3also increase.In the Fig.1a, there is no distinct difference between the O K-edge spectra of iron oxides with different valence states.There are four individual peaks (a–d) in O K-edge spectra, the relative intensity of peak a (prepeak), located in the range below 540 eV, increases significantly when iron oxide is changing from FeO toα-Fe2O3.It corresponds to the transition from the O 1s core state to the unoccupied state of O 2p hybridized with Fe 3d states.As reported in the literature [14,15], the intensity of peak a is closely related to the occupancy number of the Fe 3d orbital, the coordination geometry around the O atom, and the bond length between O atom and nearest neighbor Fe atom.The Peak b is emerged at around 543 eV,generally speaking, resulted from the O 2p unoccupied state hybridized with Fe 4s and 4p states.The relative weak peak c located at around 550 eV and the broad peak d centered at around 565 eV can also be observed in Fig.1a.However, there exists a remarkable controversy on whether peak c and peak d belong to EXELFS or energy-loss near-edge structure (ELNES) [15,16].There is one point of view, peak d stems from the backscattering in the first shell of Fe atoms and peak c is attributed to the backscattering in the second shell of O atoms [15].

    Fig.1.The O K-edge spectra (a) and Fe L-edge spectra (b) for iron oxides (FeO,Fe3O4, γ-Fe2O3, α-Fe2O3).

    Fig.2.The crystal structures and corresponding simulations of O K-edge ELNES spectra for FeO (a, d), Fe3O4 (b, e) and α-Fe2O3 (c, f).

    Due to random arranged oxygen vacancies, the accurate crystal model ofγ-Fe2O3is hard to be set up.In consideration of similar crystal structure and electronic structure betweenγ-Fe2O3and Fe3O4, the O K-edge ELNES simulation of Fe3O4can be substituted for that ofγ-Fe2O3.In order to explore the origin of these characteristic peaks, the crystal models of FeO, Fe3O4andα-Fe2O3have been established, as shown in Figs.2a–c.Based on these structures of common iron oxides, the O K-edge ELNES spectra for these three iron oxides have been simulated through WIEN2k, as shown in Figs.2d–f.The electron configurations of the associated oxygen polyhedral ligands are shown in Fig.S1 (Supporting information) [17].The FeO is an antiferromagnetic material with NaCl type structure (Fig.2a).There are three peaks centered at 532, 536 and 539 eV in a simulation of FeO ELNES spectra (Fig.2d), and the peaks located at 532 and 539 eV correspond to peak a and peak b in O K-edge spectra (Fig.1a).However, there is no peak in Fig.1a corresponding to the peak at 536 eV observed in the O K-edge spectra in Fig.2d.Moreover, the position of peak b from the simulation is lower than that from the experimental spectra.This is due to the presence of cationic defect sites within FeO as a nonstoichiometric compound, which leads to a slightly higher average valence of Fe than 2+.The Fe3O4is a typical ferrimagnetic material with an inverse spinel structure, in which the oxygen ions form a face-centered cubic lattice and the Fe atoms occupy tetrahedral and octahedral sites.Two different sublattices are formed by these Fe atoms in different oxygen polyhedral center sites, which are coupled in antiferromagnetic forms (Fig.2b).In the simulation of Fe3O4, there are two peaks centered at 530 eV and 540 eV (Fig.2e), corresponding to peak a and peak b in experimental spectra of Fe3O4, respectively.The spin polarization effect on the spectral structure is also shown in Fig.2e.Compared to the spin-up spectra, the peak a of spin-down spectra appears splitting, which is similar to the feature ofα-Fe2O3simulated spectra.Meanwhile,theα-Fe2O3has a hexagonal structure with oxygen close stacking,and exhibits antiferromagnetism and relatively high thermal stability (Fig.2c).Similar to the case of Fe3O4, the simulated spectrum ofα-Fe2O3also has two peaks centered at 530 and 540 eV, corresponding to peak a and peak b.The difference is that the simulated peak a ofα-Fe2O3appears splitting (Fig.2f).Although this phenomenon could not be observed in the experimental spectra due to low energy resolution, it is feasible to observe split peak by continuously improving the resolution.

    In the Fig.S1, the diagrams of molecular orbital for octahedral ligands of (FeO6)10?, (FeO4)5?and (FeO6)9?are shown.Combined with this figure and density of states calculations, the influence of electronic structure on the O K-edge spectra can be further understood.The spin-up states of spin-up polarized Fe atoms and spindown states of spin-down polarized Fe atoms are both under Fermi energy, it is in accord with the character of high-spin electronic configurations.In other words, the suborbitals of d electron are occupied by at least one electron in iron oxides.The O atoms around iron atom not only split d electron orbital into egand t2gstates,but also hybridize with partial energy-level (3d, 4s and 4p states)of Fe atom.The energy states of lone-pair electrons in O atoms are denoted asσ.In the oxygen octahedral ligands ((FeO6)10?and(FeO6)9?), six undifferentiatedσstates are hybridized with the eg,4s and 4p states of Fe atom, to form bonding orbital a1g, t1u, eg,and antibonding orbital a1g?, t1u?, eg?.Without hybridization, Fe t2gstate forms a non-bonding orbital.As for oxygen tetrahedral ligand ((FeO4)5?), four undifferentiatedσstates are hybridized with the 4s and 4p states of Fe atom, to form bonding orbital a1g, t1u,and antibonding orbital a1g?, t1u?.Meanwhile, Fe t2gand egstates form two non-bonding orbital without hybridization.From these diagrams, the different ligand structures of iron oxides can result in different unoccupied energy state distribution.As shown in Fig.S2 (Supporting information), the O 2p orbitals and Fe 3d orbitals of FeO, Fe3O4andα-Fe2O3with different spin directions are clearly presented.By analyzing the electronic structures of these three iron oxides, it can be found that unoccupied states of O 2p electron, which hybridized with Fe 3d electron, are mainly distributed above the Fermi energy 0—5 eV.This indicates that the K-shell ionization edge from threshold energy 0—5 eV (peak a structure) could be caused by electron transferring from 1s orbital to the antibonding orbital formed by Fe and O hybridization.It is confirmed that peak a structure in oxygen ELNES is related to the Fe-O bond.

    In FeO, Fe atoms have four electrons in the t2gorbital, which means that electron numbers of dxydyzdxzorbital are not equal,causing the t2gorbital to be split (known as the Jahn-Teller effect).This can be seen in Fig.S2a, the spin-up Fe atoms (red line) have two relatively sharp peaks in spin-down states, which appear at occupied states (near Fermi energy) and unoccupied states (about 1.2 eV), respectively.These two peaks both belong to t2gorbital.Some suborbitals, which could be fully occupied, are pulled down in energy until they are near the Fermi energy.Others increase in energy and lead to a decrease in the splitting energy between eg?and t2gorbitals.With only three electrons in t2gorbit, there is no Jahn-Teller effect inα-Fe2O3.Therefore, the phenomenon that two relatively sharp peaks appear above Fermi energy can be observed in Fig.S2c, corresponding to eg?and t2gorbits, respectively.The structure of Fe3O4is relatively complex, including not only two kinds of oxygen octahedral ligand structures (FeO type andα-Fe2O3type), but also oxygen tetrahedral ligand structure.In Fig.S2b, we can observe that the unoccupied energy state of spin-up Fe atom (located at the center of octahedral ligand) has both the features ofα-Fe2O3and FeO.The partial t2gstates reduce energy and form occupied states (FeO feature), t2gand egstates are discrete and easy to identify (α-Fe2O3feature).In the spin-down Fe atom (located at the center of tetrahedral ligand), the t2gorbital energy is higher than egorbital energy, however, the splitting energy is smaller than that in octahedral ligand.Therefore, the splitting degree of peak a for Fe3O4is between the two cases of FeO andα-Fe2O3due to the overlapping electron structures from these different ligands.

    The near edge fine structure simulation results of iron oxides do not reveal the origin of peaks b, c, and d, suggesting that these peaks may belong to extended edge fine structure.This part of the spectral structures originates from interference between the ejected electron of the excited atom and the backscattered electron from the neighbor atom.It is beneficial to analyze the extended edge fine structure for revealing the coordination of neighboring atoms around oxygen atoms.In order to show the coordination of oxygen atom, the radial distribution function (RDF) obtained by Fourier transform (FT) is used to describe the distribution of neighbor atoms.Firstly, the oscillating part of EXELFSχ(k) is obtained by the piecewise spline function fitting method with 5thorder [13,18].Then, before Fourier transform,χ(k) is reduced noise,and the de-noising effect is shown in Fig.S3 (Supporting information).As the backscattering amplitude factor presents a strong dependence on the atomic number Z, different atomic backscatters from the heavy atom and light atom are localized in the different k-space.Therefore, the type of atoms can be distinguished byR~ktwo-dimensional relationship diagram of kn-weighted EXELFS [19–21].To achieve this purpose, wavelet transform (WT) is applied for obtaining ofR~ktwo-dimensional relationship diagram.

    Fig.3.The radial distribution function (a), wavelet analysis contour (b) and oscillating part of EXELFS (c) of O atoms in α-Fe2O3.

    The Fourier transform (FT) and WT analysis results ofα-Fe2O3are displayed in Fig.3, and corresponding results for FeO, Fe3O4andγ-Fe2O3are presented in Figs.S4-S6 (Supporting information),respectively.Takingα-Fe2O3as an example, the RDF (Fig.3a) could be obtained by the oscillation signal which isχ(k) spectrum multiplied by window function W(k) and amplitude modulation factork(Fig.3c).Although the distance of neighbor atoms can be determined by the peaks in RDF, some fake peaks derived from the data truncation, would appear in the curve.Therefore, the wavelet analysis diagram that representing the relation between interatomic distanceRand wave vectork, would be helpful in identifying the actual neighboring atoms.In the wavelet analysis results (Fig.3b),a lobe is localized in theΔk~(0–6.5) ?A?1andR~(1.5–3.5) ?A range.By changing window function W(k)range, wavelet analysis results of two regions (with or without b, c, and d structure) were obtained.In the region including b, c, and d structure (Fig.S7b in Supporting information), a lobe centered onR=2.5 ?A could be found in the wavelet analysis contour.The result does not suggest that there are some neighbor atoms at 2.5 ?A from the central O atom.In fact, this situation is due to the common effect of neighbor atoms backscattering at 2 ?A and 3 ?A.However, the broadening effect due to the extreme small range, would cause an overlap of two lobes overlapping.In the region without b, c, and d structure (Fig.S7c in Supporting information), only a lobe centered onR=2 ?A could be found in the wavelet analysis contour.This indicates that lobe in Fig.3b is made of the overlapping of two small lobes.By comparing with the center of two lobes, peaks located at 2 ?A and 3 ?A can be determined that they are corresponding to the nearest neighbor Fe atom and O atom, respectively.The similar results could also be obtained in the cases of FeO, Fe3O4andγ-Fe2O3, as shown in Figs.S4-S6.Overall, the O K edge extended fine structure is mainly caused by the backscattering of the nearest Fe atom.In addition, the b, c and d peaks are affected by the backscattering of the second neighbor O atoms.

    Fig.4.The fitting diagram of O K-edge structures (a) and radial distribution function (b) in α-Fe2O3 nanorods.The coordination structures of oxygen atom in FeO(c), Fe3O4 (d) and α-Fe2O3 (e).

    Based on the analysis above, an effective method could be proposed in order to analyze oxides microstructure quantitatively.It is easy to obtain the valence state from near edge structure, which reflects the electronic structure of the material.The basic edge shapes of the O K spectrum are saw tooth, and with additional fine structures.Since the peak a reflects the chemical bonding between oxygen and transition metal, if more d electrons of transition metals are involved in bonding, more antibonding orbitals will be produced accordingly.In short, the intensity of peak a relative to the ionization edge main structure is positively correlated with the valence state of transition metal.For convenience of use, instead of using the hydrogenic model or Hartree-Slater method to calculate the basic shape of O K-edge, it can be used to estimate valence state by simply measuring the peak amplitude ratio of peak a to peak b.Although, the peak b is affected by the electronic structures of oxide, the body is mainly influenced by the inelastic scattering of incident electrons.The relative amplitude between peak a and peak b has been measured by Gaussian fitting method, shown in Table S1 (Supporting information).In Fig.4a, the fitting diagram ofα-Fe2O3is displayed.The analytic result indicates that there is a linear correlation between the relative amplitude and the valence state of Fe, and that the relative amplitude can be considered as the calibration for a valence state.The extended edge structure reflects the coordination information of atoms around oxygen atom,consequently, the ratio of the nearest neighbor atoms to the second nearest neighbor atoms around the oxygen atom can be obtainedviaGaussian fitting (Fig.4b), the analysis result is listed in Table S2 (Supporting information).Fig.4b shows the fitting diagram of the radial distribution function inα-Fe2O3.The coordination structure models of oxygen in different iron oxides have been established, as shown in Figs.4c–e, to investigate the relationship between the result and oxygen coordination.Therefore, the atoms near oxygen atoms from these iron oxides are focused on during the research process.The coordination structure of FeO with 6 nearest neighbor iron atoms and 12 second nearest neighbor oxygen atoms distributed around the central oxygen atom is presented in Fig.4c, and the valence state of all iron atoms is +2.Fig.4d shows the coordination structure of Fe3O4, in which the central oxygen atom is surrounded with 4 nearest neighbor iron atoms and 12 second nearest neighbor oxygen atoms.The distances from four iron atoms to the central oxygen atom are different, one iron atom is farther away from the central oxygen atom than the other three iron atoms.This iron atom is in the center of the tetrahedral ligand and the other three iron atoms are in the center of the octahedral ligand.As for the coordination structure ofα-Fe2O3, which has 4 nearest neighbor iron atoms and 12 second nearest neighbor oxygen atoms surrounded the central oxygen atom, the distances between four iron atoms to the central oxygen atoms are identical.Considering that FeO as a nonstoichiometric compound, its actual structure can be expressed as Fe0.83O~Fe0.95O due to a large number of cation defects.So this means that the ratio of coordination number needs to be corrected to 1.349~1.544 for FeO model.Although the ratio of coordination number obtained by RDF function is not equal to that from the coordination structure model, there exists a linear correlation between them.When these experimental values multiplied by about 1.176, it is that these experimental values are in good agreement with the values in the model.Although the nature of the coefficient 1.176 is not yet clear, however, the value of the coordination number can be accurately determined by using standard samples in the process of practical application,thus avoiding the determination of this coefficient.Therefore, the quantitative analysis on the change of atomic coordination in microstructure can be realized by extended edge fine structure.In order to reveal the coordination structure of oxygen atoms more accurately, the future research and investigation would be required.

    In summary, an approach of resolving EELS spectra to analyze fine microstructures of TMOs by the O-K than Fe-L was proposed.Relied on the corporation of DFT calculation, Fourier transform (FT)and wavelet transform (WT), the relationship between electron energy loss O K spectra and microstructures for iron oxides was investigated in detail.The origins of four individual peaks in O Kedge spectra have been interpreted by the comparison between experimental spectra and simulated spectra, and the contrasts among wavelet analysis contours of O atoms with different wave vector ranges.Consequently, the quantitative analysis on the change of valence state and atomic coordination in microstructure can be realized by the near edge structure and extended edge fine structure,respectively.This analytical method based on EELS O K-edge spectra could obtain valuable structural information for TMOs, which encourages the quantitative structure-activity relationship building.Furthermore, this methodviaEELS spectra of non-metal element can expand the research on the microstructures of TMOs.

    Declaration of competing interest

    The authors declare no conflict of interest.

    Acknowledgments

    The authors gratefully acknowledge the financial support provided by the National Natural Science Foundation of China (Nos.22072164, 51932005, 21773269, 52161145403), Liao Ning Revitalization Talents Program (No.XLYC1807175), and the Research Fund of SYNL.

    Supplementary materials

    Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2021.12.027.

    又爽又黄a免费视频| av福利片在线| 亚洲天堂av无毛| 亚洲久久久国产精品| 国产白丝娇喘喷水9色精品| 国产高清不卡午夜福利| 精品国产一区二区久久| av卡一久久| 欧美丝袜亚洲另类| 久久精品国产亚洲av涩爱| 男女边摸边吃奶| 99久久精品热视频| 色哟哟·www| 校园人妻丝袜中文字幕| 国产亚洲精品久久久com| 成人18禁高潮啪啪吃奶动态图 | 欧美激情极品国产一区二区三区 | 久久久国产一区二区| 一级a做视频免费观看| 国产高清国产精品国产三级| 一区二区三区免费毛片| 在线观看www视频免费| 女的被弄到高潮叫床怎么办| 男女无遮挡免费网站观看| 久久毛片免费看一区二区三区| 亚洲三级黄色毛片| 日韩精品免费视频一区二区三区 | 亚洲精华国产精华液的使用体验| 欧美最新免费一区二区三区| 日韩人妻高清精品专区| 精品久久久久久电影网| 久久久精品94久久精品| 亚洲成人手机| 免费黄色在线免费观看| 亚洲精品国产av蜜桃| 国产一区有黄有色的免费视频| 蜜桃在线观看..| 欧美日本中文国产一区发布| 男女边摸边吃奶| 国产 一区精品| 视频区图区小说| 自拍欧美九色日韩亚洲蝌蚪91 | 国产在线一区二区三区精| 国产亚洲av片在线观看秒播厂| 国产av一区二区精品久久| 高清视频免费观看一区二区| 日韩,欧美,国产一区二区三区| 九九在线视频观看精品| av福利片在线| 男女边吃奶边做爰视频| 丰满少妇做爰视频| 啦啦啦啦在线视频资源| 人妻 亚洲 视频| 国产精品99久久久久久久久| 天天操日日干夜夜撸| 国产 一区精品| 国产极品天堂在线| 女人精品久久久久毛片| 久久青草综合色| 免费在线观看视频国产中文字幕亚洲 | 制服诱惑二区| 精品第一国产精品| 精品久久久精品久久久| 国产一级毛片在线| 搡老岳熟女国产| 国产伦理片在线播放av一区| 亚洲黑人精品在线| 国产一区二区三区av在线| 一区二区av电影网| 黑人欧美特级aaaaaa片| 精品久久久精品久久久| 久久久久久久国产电影| 亚洲第一青青草原| 亚洲性夜色夜夜综合| 1024视频免费在线观看| 国产97色在线日韩免费| 国精品久久久久久国模美| 一边摸一边抽搐一进一出视频| 黄色 视频免费看| 精品亚洲成a人片在线观看| 一区二区日韩欧美中文字幕| 国产亚洲精品久久久久5区| 欧美日韩成人在线一区二区| 少妇的丰满在线观看| 9色porny在线观看| 天天躁日日躁夜夜躁夜夜| av免费在线观看网站| 久久人人97超碰香蕉20202| 精品人妻在线不人妻| 亚洲精品粉嫩美女一区| 在线观看免费日韩欧美大片| 一级片'在线观看视频| 国精品久久久久久国模美| 在线 av 中文字幕| 亚洲伊人色综图| 欧美黑人精品巨大| 欧美日韩亚洲综合一区二区三区_| 黄片大片在线免费观看| 爱豆传媒免费全集在线观看| 成人影院久久| 国产精品久久久久久精品电影小说| 久久久久久人人人人人| 欧美黑人欧美精品刺激| 交换朋友夫妻互换小说| 中文字幕色久视频| 老司机午夜福利在线观看视频 | 中文字幕人妻丝袜一区二区| 18禁黄网站禁片午夜丰满| 日韩制服骚丝袜av| 日韩中文字幕视频在线看片| 日日摸夜夜添夜夜添小说| 老熟妇乱子伦视频在线观看 | 久久久国产成人免费| 久久 成人 亚洲| 伊人久久大香线蕉亚洲五| 日韩视频一区二区在线观看| 精品一区在线观看国产| 一区二区三区激情视频| 啦啦啦 在线观看视频| 99热国产这里只有精品6| 国产成人精品久久二区二区免费| 美女午夜性视频免费| 母亲3免费完整高清在线观看| 国产人伦9x9x在线观看| 亚洲精品中文字幕在线视频| 天堂中文最新版在线下载| 欧美国产精品一级二级三级| 纯流量卡能插随身wifi吗| 午夜福利免费观看在线| 亚洲国产中文字幕在线视频| 999久久久精品免费观看国产| 亚洲精品美女久久av网站| 亚洲国产av新网站| 久热爱精品视频在线9| 人妻一区二区av| 在线观看免费午夜福利视频| 久久久久网色| 欧美老熟妇乱子伦牲交| 精品一区二区三区四区五区乱码| 久久国产精品影院| 国产成人精品在线电影| 2018国产大陆天天弄谢| 18禁观看日本| 电影成人av| 国产成人av激情在线播放| 亚洲精品国产av蜜桃| 肉色欧美久久久久久久蜜桃| 欧美国产精品一级二级三级| 久久精品国产a三级三级三级| 中文字幕人妻丝袜制服| 久久久久精品国产欧美久久久 | av超薄肉色丝袜交足视频| 久久热在线av| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久久久久婷婷小说| 色视频在线一区二区三区| 韩国精品一区二区三区| 一级,二级,三级黄色视频| 久久99热这里只频精品6学生| 亚洲av成人不卡在线观看播放网 | 丰满饥渴人妻一区二区三| 国产三级黄色录像| 亚洲精品国产色婷婷电影| 国产精品1区2区在线观看. | 成人国产一区最新在线观看| 伊人亚洲综合成人网| 国产精品一区二区免费欧美 | 午夜成年电影在线免费观看| 一边摸一边做爽爽视频免费| 91av网站免费观看| 国产成人av激情在线播放| 精品人妻1区二区| 大型av网站在线播放| 亚洲精品国产色婷婷电影| 999久久久国产精品视频| 黄频高清免费视频| 色视频在线一区二区三区| 亚洲伊人色综图| 欧美精品人与动牲交sv欧美| 久久久久久久久免费视频了| 在线观看人妻少妇| 波多野结衣一区麻豆| 国产在视频线精品| 亚洲久久久国产精品| 国产成人免费观看mmmm| 老汉色∧v一级毛片| 国产精品av久久久久免费| 欧美少妇被猛烈插入视频| 亚洲成人免费电影在线观看| 久久中文字幕一级| 50天的宝宝边吃奶边哭怎么回事| 久久狼人影院| 超色免费av| 亚洲中文日韩欧美视频| 国产伦理片在线播放av一区| av网站在线播放免费| 久久久久久免费高清国产稀缺| 老司机影院毛片| 国产激情久久老熟女| 欧美亚洲日本最大视频资源| 欧美在线一区亚洲| 国产亚洲精品一区二区www | 国产又色又爽无遮挡免| 久久久久久亚洲精品国产蜜桃av| 超碰97精品在线观看| 亚洲性夜色夜夜综合| 高清在线国产一区| 免费黄频网站在线观看国产| 如日韩欧美国产精品一区二区三区| 精品少妇内射三级| 人人妻人人爽人人添夜夜欢视频| 日本猛色少妇xxxxx猛交久久| 午夜精品国产一区二区电影| 亚洲五月婷婷丁香| 女性被躁到高潮视频| 男女床上黄色一级片免费看| 亚洲欧洲日产国产| 亚洲成人免费av在线播放| 桃红色精品国产亚洲av| 欧美精品啪啪一区二区三区 | 成人三级做爰电影| av线在线观看网站| 国产精品成人在线| 男人爽女人下面视频在线观看| 欧美黑人欧美精品刺激| 丰满人妻熟妇乱又伦精品不卡| 每晚都被弄得嗷嗷叫到高潮| 成年美女黄网站色视频大全免费| 久久青草综合色| 女人高潮潮喷娇喘18禁视频| 两性夫妻黄色片| 女人爽到高潮嗷嗷叫在线视频| 欧美午夜高清在线| 777米奇影视久久| 蜜桃国产av成人99| 亚洲一区二区三区欧美精品| 国产精品久久久久成人av| 99re6热这里在线精品视频| 一本色道久久久久久精品综合| 丰满饥渴人妻一区二区三| 黄频高清免费视频| netflix在线观看网站| 国产不卡av网站在线观看| 亚洲av男天堂| 国产欧美日韩一区二区三区在线| 中亚洲国语对白在线视频| 美女福利国产在线| 亚洲成av片中文字幕在线观看| av网站在线播放免费| 亚洲 国产 在线| 最新的欧美精品一区二区| 午夜福利影视在线免费观看| 午夜福利免费观看在线| 国产av国产精品国产| 一区福利在线观看| 成年女人毛片免费观看观看9 | 日日摸夜夜添夜夜添小说| 国产精品久久久人人做人人爽| 日本撒尿小便嘘嘘汇集6| 日韩一区二区三区影片| 欧美精品一区二区大全| 欧美 亚洲 国产 日韩一| 亚洲国产日韩一区二区| 成人18禁高潮啪啪吃奶动态图| 少妇精品久久久久久久| 两人在一起打扑克的视频| 999久久久精品免费观看国产| 69av精品久久久久久 | 国产深夜福利视频在线观看| 欧美日韩视频精品一区| 国产精品国产三级国产专区5o| 欧美激情高清一区二区三区| 少妇猛男粗大的猛烈进出视频| 日韩欧美一区二区三区在线观看 | 久久人妻福利社区极品人妻图片| 亚洲国产欧美网| 亚洲男人天堂网一区| 中文精品一卡2卡3卡4更新| 老熟妇仑乱视频hdxx| 亚洲激情五月婷婷啪啪| 老司机在亚洲福利影院| 成年美女黄网站色视频大全免费| 中文字幕高清在线视频| 亚洲视频免费观看视频| 人妻久久中文字幕网| 自线自在国产av| 国产一区二区三区av在线| 美女国产高潮福利片在线看| 天天添夜夜摸| 91精品伊人久久大香线蕉| 人人妻人人爽人人添夜夜欢视频| 国产精品一区二区在线不卡| 一区二区日韩欧美中文字幕| 亚洲国产中文字幕在线视频| 精品视频人人做人人爽| 久久久久国产一级毛片高清牌| 午夜精品国产一区二区电影| 99香蕉大伊视频| 亚洲精品第二区| 男女免费视频国产| 高清欧美精品videossex| 大香蕉久久网| 日日夜夜操网爽| 国产不卡av网站在线观看| 人妻久久中文字幕网| 国产日韩欧美亚洲二区| av片东京热男人的天堂| 日韩欧美国产一区二区入口| 99热国产这里只有精品6| 另类精品久久| 国产精品偷伦视频观看了| 国产精品99久久99久久久不卡| 久久人妻福利社区极品人妻图片| 精品久久久精品久久久| 美女脱内裤让男人舔精品视频| 看免费av毛片| 中文字幕人妻丝袜一区二区| 我的亚洲天堂| 十八禁网站网址无遮挡| 午夜免费鲁丝| 男女床上黄色一级片免费看| 黄色a级毛片大全视频| 国产一卡二卡三卡精品| 国产日韩一区二区三区精品不卡| 他把我摸到了高潮在线观看 | 成人国产一区最新在线观看| 国产成人av教育| 婷婷成人精品国产| 香蕉国产在线看| 久久久久国产精品人妻一区二区| 中国美女看黄片| 老熟妇乱子伦视频在线观看 | 脱女人内裤的视频| 亚洲国产日韩一区二区| 亚洲精品av麻豆狂野| 久久av网站| svipshipincom国产片| 午夜影院在线不卡| 狂野欧美激情性xxxx| 乱人伦中国视频| 亚洲人成电影免费在线| 久久人人97超碰香蕉20202| 99香蕉大伊视频| 婷婷色av中文字幕| tocl精华| 久久久精品区二区三区| 在线观看免费日韩欧美大片| 久久久久久免费高清国产稀缺| 嫁个100分男人电影在线观看| 久久国产精品影院| 黄片大片在线免费观看| 亚洲精品国产区一区二| 男人操女人黄网站| 狠狠精品人妻久久久久久综合| 欧美日韩视频精品一区| 纯流量卡能插随身wifi吗| 欧美人与性动交α欧美精品济南到| 九色亚洲精品在线播放| 久久天躁狠狠躁夜夜2o2o| 老司机在亚洲福利影院| 亚洲一区二区三区欧美精品| 丰满饥渴人妻一区二区三| 亚洲av国产av综合av卡| 午夜两性在线视频| 欧美精品一区二区免费开放| 欧美日韩亚洲高清精品| 一进一出抽搐动态| 亚洲av成人不卡在线观看播放网 | 巨乳人妻的诱惑在线观看| 最新在线观看一区二区三区| 亚洲第一av免费看| www.精华液| 久热这里只有精品99| √禁漫天堂资源中文www| 精品视频人人做人人爽| 男女国产视频网站| 亚洲成人国产一区在线观看| 亚洲精品久久午夜乱码| 精品人妻熟女毛片av久久网站| 亚洲精品久久午夜乱码| 女人久久www免费人成看片| 在线观看免费高清a一片| 大香蕉久久网| 大片免费播放器 马上看| 日韩视频一区二区在线观看| 不卡av一区二区三区| 国产伦人伦偷精品视频| 岛国在线观看网站| 国产黄色免费在线视频| 在线观看www视频免费| 国产av又大| 丝袜在线中文字幕| 亚洲人成电影免费在线| 极品少妇高潮喷水抽搐| 国产男人的电影天堂91| 欧美在线黄色| 2018国产大陆天天弄谢| 亚洲精品一二三| 亚洲情色 制服丝袜| 国产精品久久久人人做人人爽| 欧美精品av麻豆av| 日本精品一区二区三区蜜桃| 成人18禁高潮啪啪吃奶动态图| 人妻一区二区av| 一级毛片女人18水好多| 久久久久精品人妻al黑| 国产免费一区二区三区四区乱码| 美女脱内裤让男人舔精品视频| 一本大道久久a久久精品| 久久精品国产亚洲av高清一级| 一边摸一边抽搐一进一出视频| 亚洲国产欧美一区二区综合| 激情视频va一区二区三区| 国产精品成人在线| 好男人电影高清在线观看| 日本一区二区免费在线视频| 日本91视频免费播放| 色综合欧美亚洲国产小说| av欧美777| 久久ye,这里只有精品| 亚洲,欧美精品.| 天堂8中文在线网| 男女床上黄色一级片免费看| 精品少妇内射三级| 亚洲性夜色夜夜综合| 欧美性长视频在线观看| 国产精品 国内视频| 老司机在亚洲福利影院| 十八禁人妻一区二区| 久久久精品区二区三区| 午夜老司机福利片| 免费黄频网站在线观看国产| 午夜成年电影在线免费观看| 久久亚洲国产成人精品v| 淫妇啪啪啪对白视频 | 丰满迷人的少妇在线观看| 欧美黑人欧美精品刺激| 色播在线永久视频| 亚洲欧美一区二区三区黑人| av免费在线观看网站| 国产精品成人在线| 乱人伦中国视频| 大码成人一级视频| 国产av一区二区精品久久| 在线看a的网站| 成年人午夜在线观看视频| 纵有疾风起免费观看全集完整版| 国产精品99久久99久久久不卡| 十八禁人妻一区二区| 午夜视频精品福利| 女人精品久久久久毛片| h视频一区二区三区| 欧美久久黑人一区二区| 午夜福利在线免费观看网站| 成年人午夜在线观看视频| 日韩大码丰满熟妇| 狠狠精品人妻久久久久久综合| 最新在线观看一区二区三区| 少妇被粗大的猛进出69影院| 看免费av毛片| 中文字幕人妻丝袜制服| 精品少妇一区二区三区视频日本电影| 一区在线观看完整版| av天堂在线播放| 大码成人一级视频| 欧美黄色片欧美黄色片| 真人做人爱边吃奶动态| 亚洲精品中文字幕一二三四区 | 91麻豆精品激情在线观看国产 | 久久精品国产亚洲av高清一级| 一二三四在线观看免费中文在| 亚洲色图 男人天堂 中文字幕| 老熟女久久久| 亚洲国产中文字幕在线视频| 久久久久精品人妻al黑| 无限看片的www在线观看| 男女边摸边吃奶| 飞空精品影院首页| 可以免费在线观看a视频的电影网站| 国产精品久久久人人做人人爽| 看免费av毛片| 男女之事视频高清在线观看| 韩国高清视频一区二区三区| 色婷婷av一区二区三区视频| 欧美变态另类bdsm刘玥| 真人做人爱边吃奶动态| 人成视频在线观看免费观看| 老汉色∧v一级毛片| 十分钟在线观看高清视频www| 国产av一区二区精品久久| 亚洲 欧美一区二区三区| 成在线人永久免费视频| 狂野欧美激情性bbbbbb| 18禁黄网站禁片午夜丰满| 国产xxxxx性猛交| 搡老乐熟女国产| 国产免费现黄频在线看| 欧美激情久久久久久爽电影 | 亚洲精品久久午夜乱码| 狠狠婷婷综合久久久久久88av| 秋霞在线观看毛片| 极品人妻少妇av视频| 午夜福利乱码中文字幕| 亚洲精品久久午夜乱码| 亚洲欧美一区二区三区黑人| 国产无遮挡羞羞视频在线观看| 日韩制服骚丝袜av| 国产精品久久久人人做人人爽| 久久 成人 亚洲| 永久免费av网站大全| 美女视频免费永久观看网站| 亚洲国产av新网站| 91成人精品电影| 国产精品亚洲av一区麻豆| 国产又色又爽无遮挡免| 精品少妇内射三级| 俄罗斯特黄特色一大片| 欧美性长视频在线观看| 日韩大码丰满熟妇| 国产日韩欧美亚洲二区| 久久久久精品国产欧美久久久 | 欧美国产精品一级二级三级| 女性生殖器流出的白浆| 黄色视频在线播放观看不卡| 大型av网站在线播放| 黄色 视频免费看| 亚洲第一欧美日韩一区二区三区 | 天天添夜夜摸| 一边摸一边做爽爽视频免费| 日本撒尿小便嘘嘘汇集6| 国产又色又爽无遮挡免| 久久精品国产综合久久久| av网站免费在线观看视频| 在线观看www视频免费| 成人国语在线视频| 宅男免费午夜| 大片免费播放器 马上看| 一级,二级,三级黄色视频| 韩国精品一区二区三区| 国产片内射在线| 老司机影院毛片| 亚洲专区国产一区二区| 极品人妻少妇av视频| 青春草视频在线免费观看| 成人18禁高潮啪啪吃奶动态图| 俄罗斯特黄特色一大片| 动漫黄色视频在线观看| 老司机午夜福利在线观看视频 | 大陆偷拍与自拍| 秋霞在线观看毛片| 欧美黄色淫秽网站| 午夜福利一区二区在线看| 国产1区2区3区精品| 久久免费观看电影| av天堂在线播放| 久久影院123| 少妇裸体淫交视频免费看高清 | 国产精品一区二区在线观看99| 久久久久国产精品人妻一区二区| 日韩一卡2卡3卡4卡2021年| 两个人免费观看高清视频| 国产免费视频播放在线视频| 久久免费观看电影| 97在线人人人人妻| 91大片在线观看| 狠狠婷婷综合久久久久久88av| 蜜桃国产av成人99| 亚洲avbb在线观看| 国产精品九九99| 午夜激情av网站| 人人妻人人澡人人看| 日韩欧美一区视频在线观看| 99国产精品99久久久久| 18在线观看网站| 国产精品影院久久| 男女国产视频网站| 91av网站免费观看| 蜜桃国产av成人99| 亚洲国产中文字幕在线视频| 黄色视频不卡| 97人妻天天添夜夜摸| 亚洲人成电影免费在线| 另类亚洲欧美激情| 无限看片的www在线观看| 69精品国产乱码久久久| 久久精品国产亚洲av香蕉五月 | 久久热在线av| 在线看a的网站| 自线自在国产av| 国产日韩一区二区三区精品不卡| 91av网站免费观看| 国产成人欧美| 亚洲精品国产一区二区精华液| 黄色视频不卡| 亚洲精品成人av观看孕妇| 各种免费的搞黄视频| 久久久久精品人妻al黑| 亚洲精品一卡2卡三卡4卡5卡 | 欧美日韩黄片免| www.精华液| 美女脱内裤让男人舔精品视频| 国产一级毛片在线| 国产精品 欧美亚洲| 免费高清在线观看日韩| 建设人人有责人人尽责人人享有的| 日韩大片免费观看网站| 岛国毛片在线播放| 日韩一区二区三区影片| 色综合欧美亚洲国产小说| 亚洲av电影在线进入| 不卡一级毛片| 国产激情久久老熟女| 亚洲九九香蕉| 纯流量卡能插随身wifi吗| 制服人妻中文乱码|