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

    Quantum plasmons in the hybrid nanostructures of double vacancy defected graphene and metallic nanoarrays*

    2021-01-21 02:08:34RuiTang唐睿YangXu徐陽HongZhang張紅andXinLuCheng程新路
    Chinese Physics B 2021年1期
    關(guān)鍵詞:張紅新路

    Rui Tang(唐睿), Yang Xu(徐陽), Hong Zhang(張紅),3,?, and Xin-Lu Cheng(程新路),3

    1College of Physics,Sichuan University,Chengdu 610065,China

    2Institute of Atomic and Molecular Physics,Sichuan University,Chengdu 610065,China

    3Key Laboratory of High Energy Density Physics and Technology of Ministry of Education,Sichuan University,Chengdu 610065,China

    Keywords: plasmon, double vacancy defected graphene, mixed metallic nanoarrays, multi-beam laser, harmonic spectrum

    1. Introduction

    Surface plasmons (SPs), excitations in metal nanostructure, can break the diffraction limit of nanometer level.[1,2]This excellent property can be used in many aspects such as optical nano-antenna,[3]catalytic reactions,[4]optical waveguides,[5]molecular electronics,[6]and optical biosensors.[7]However, the SPs in metal is liable to decay during propagation,and the SPs in metal devices is relatively simple and can not be effectively modulated. Graphene is one-atom-thick two-dimensional layers of sp2-bonded carbon,which has many advantages, e.g., exceptionally high crystal and electronic quality,[8]ultrahigh carrier mobility,[9]and wide range of light absorption,[10]and can be controlled by the grating.[11,12]In particular, compared with metal nanostructure, the plasmon on the surface of the graphene has some more excellent properties. The most one is a relatively long propagation distance with a light decay.[13,14]Therefore,graphene can make up the drawbacks of noble metals in SPs,for example, it can replace noble metals for optical signal transmission.[15,16]Recently,the research result of Abajo et al.indicated that when the graphene scale is down to the molecular size,the layout of the resonance frequency of plasmon can be extended to visible light and ultraviolet range.[17]The photoelectric devices made of graphene have been experimentally obtained, such as photoelectric detectors, ultrafast laser, and optical modulator.[18–20]

    At present, graphene has been experimentally prepared by a number of methods.[21–23]Due to the limitation of manufacturing process, all sorts of defects inevitably appear in graphene,including Stone–Wales defects,vacancy defects,absorption defects, etc.[24–26]One of the most typical and common defects is double vacancy defects (DVDs) (as shown in Fig. 1).[27,28]The defects would change graphene’s properties just like the enhancement of the activity of graphene substrate,[29]prone to the stable adsorption of different molecules or atoms on its surface,[30]breaking the symmetry of graphene,[31]inducing magnetism,[32]and also destroying the metal properties and enhancing the semiconductor properties. Numerous methods have been used for observing the grain boundaries and defects in graphene sheets, for example,transmission electron microscopy,[33,34]atomic force microscopy, scanning tunneling spectroscopy[35]and scanning tunneling microscopy.[36]Although there are several methods to detect the defects,they have some shortcomings such as the need for complex sample preparation, taking up more time,and a limited observation range. Furthermore, much less attention is paid to DVDGr than defect-free graphene,and there is no research on plasmons of DVDGr neither theoretical nor experimental, while the DVDGr is so common in applications. Therefore, theoretical research for DVDGr plasmons has important guiding significance for practical applications of graphene.

    Theoretical calculations have predicted the presence of collective plasmon modes in a few atom chains of several metallic elements: Na,[37]Ag,[38]and Au;[39]see, for example, Ref.[40]for a review. In a previous study,[41]plasmon–exciton interactions in a monolayer graphene–Ag nanoparticle hybrid system were systemically investigated. The experimental results demonstrated that the surface catalytic reactions can not only be controlled by plasmon–exciton coupling, but also be affected by the gate voltages and electric currents (or bias voltages). With a fabricated electro-optical device based on a hybrid of graphene–Ag nanoparticles, Lin et al. studied plasmon–exciton codriven oxidation reactions. It is found that the most contribution to the oxidation reaction is from hot holes on graphene,controlled by the gate voltage and the bias voltage.[42]Nayyar et al.examined the role of transition-metal(TM) doping in the generation of plasmon modes in a single Au chain. It was found that doping leads to several additional excitations in the absorption spectrum reflecting changes in the potential around the TM atom and collective effects of the“l(fā)ocalized”(TM d-)and delocalized(s)electrons.[43]

    The linear and nonlinear optical absorption properties of TwBLG have been carefully studied. This generalized analysis method of electric field and magnetic field separation shows that Moire pattern can precisely regulate the distribution of TMDM to control the chirality of TwBLG.[44]However, the linear and nonlinear optical absorption properties of the double vacancy defected graphene and metallic nanoarrays have not been fully studied theoretically. In this paper,we present first-principles calculations of quantum plasmons in the DVDGr–gold hybrid nanostructures. When we put the DVDGr flake and gold wire together,the interaction between the light and graphene will become stronger. This can produce stronger plasmons in the visible ranges.[45–47]We address some of the fundamental questions about the nature of collective excitations in nanostructure plasmons by considering an arrangement of mixed metallic nanoarrays. The optical properties of both pure and “mixed” arrays of free-standing homonuclear chains are calculated using the time-dependent density functional theory (TDDFT). We also study the harmonics from the hybrid nanostructures exposed to the combination of lasers ranged from ultraviolet to infrared and a controlling pulse. Our simulation results may provide beneficial information for plasmons in a practical DVDGr-gold nano-device and the detection of defects on graphene.

    2. Computational details

    All calculations were performed using the real-space and real-time TDDFT code OCTOPUS,[45]and the carbon,hydrogen and gold atoms were described by the Troullier–Martins pseudopotentials.[48]Local-density approximation (LDA)[49]for the exchange-correlation was used in both the ground state and excited state calculations as employed by Marinopoulos et al.[50]To obtain the linear optical absorption spectrum,the system was excited by a δ-function impulse in the initial state.Then the Kohn–Sham wave functions were propagated for some time. The simulation zone was defined by assigning a sphere around each atom with a radius of 6 ?A and a uniform mesh grid of 0.3 ?A. In the real-time propagation, the excitation spectrums were extracted by Fourier transform of the dipole strength,and the electronic wave packets were evolved for typically 6000 steps with a time step of Δt=0.003/eV.

    Different types and shapes of time-dependent external perturbations are applied to the system,in the form of

    3. Results and discussion

    3.1. Plasmons in the planar hybrid structures

    Comparing the plasmon in defect-free graphene flake and DVDGr flake, we demonstrate that DVDs have an effect on the light absorption and change the charge distribution in graphene surface. Figures 1(a) and 1(d) indicate that the absorption peak at the high-energy modes shows a large blueshift and reduces markedly. For considering the optical application of DVDs,we just focus on the low-energy modes(less than 5 eV). As shown in Figs. 1(b) and 1(c), charge density distribution of graphene and DVDGr in the low-energy modes performs differently. The tunnel effect changes due to the surface charge distribution changes in the DVDGr for the lack of a few carbon atoms compared with the perfect graphene. When the impulse excitation is polarized along the zigzag edge, the charge is separated from the hole (Figs. 1(e) and 1(f)). The direction of the incident light exerts a crucial influence on the graphene plasmon due to the long-range charge transfer excitation of SPs in low energy region. Along the direction of excitation,electrons can vibrate back and forth over a long distance. Thus,we can detect whether or not there are defects in graphene by comparing the spectrum and the induced charge density with defect-free graphene.

    Fig.1. (a)Optical absorption of DVDGr flake and defect-free graphene flake. The impulse excitation is polarized along the armchair edge(parallel to the x-axis);(b)and(c)Fourier transform of the induced charge density for defect-free graphene and DVDGr, at the energy resonance points 4.19 and 4.02 eV,respectively. (d)Optical absorption of DVDGr flake and defect-free graphene flake. The impulse excitation is polarized along the zigzag edge (parallel to the y-axis). (e) and (f) Fourier transform of the induced charge density for defect-free graphene and DVDGr,at the energy resonance points 4.08 and 4.83 eV,respectively.

    To acquire a more particular knowledge of the electron behavior, we study the plasmon in a series of hybrid systems by placing a DVDGr sheet and the Au wires together named as DVDGr-AuX, where X (X = 2, 3, 4, and 5) represents the number of Au atoms on each side of DVDGr. As shown in Fig. 2, the edges of DVDGr are passivated with hydrogen atoms,and the distance between the nearest gold atom and the hydrogen atom was set to 2.89 ?A according to the experimental results of Nilius et al.[51]

    Figures 2(b)–2(d) show the optical absorptions of four DVDGr-AuX (X = 2, 3, 4, and 5) with resonance absorptions at 3.51,2.95,2.33,and 2.06 eV,respectively. The close absorption frequencies of isolated DVDGr (black curve) and isolated Au nanoarrays(red curve)lead to enhancement light absorption of DVDGr-AuX (blue curve) by resonance. The results illustrate vividly that the collective nature of the excitations is ascertained by the strong absorption peak. The free electrons can be transferred between the rectangular DVDGr flake and the Au cluster, and the local electromagnetic field will also be enhanced. Moreover, the gold antennas convert the incident light into plasmon accompanied by the enhancement of the electromagnetic field, and this strong nearfield can directly excite the hole-electron pair in the defective graphene. As more Au atoms are added, the position of resonance absorption peak is red-shifted about 1.5 eV and the value is raising about 17 eV-1in total. The reason for selected enhancement is the coherent superposition of DVDGr-AuX′SPs and the increase of carrier concentration.These processes lead to an increase in carrier concentration,and further affect the optical response. More importantly, the enhancement modes can be chosen by the resonant frequencies of gold wires. The coherent superposition effect will enhance the absorption peak obviously. In order to confirm the above viewpoint,Figs.2(e)–2(h)display the Fourier transform of the induced charge density for the four hybrid structures at the energy resonance points.In a previous work,Niu et al.suggested that an antiparallel dipole will be formed in the graphene because of the electromagnetic field coupling between graphene and Au nanoparticles.[52]We apply this explanation to our hybrid structures of DVDGr and gold wires.When the resonance frequencies are very close to each other, this electromagnetic coupling will evolve into plasmon coherence resonance.In addition,it can be seen that the induced charge in the DVDGr is mainly localized at the boundaries,and different modes correspond to similar spatial charge distributions. The edge state in graphene layer plays a vital role in optical response and plasmon modes.[53]

    Fig.2. Optical absorption of DVDGr-AuX consisting of the same DVDGr flake and different Au atoms: (a)DVDGr-Au2, (b)DVDGr-Au3,(c) DVDGr-Au4, and (d) DVDGr-Au5. The impulse excitation is polarized along the armchair edge. For comparison, the absorptions of the isolated DVDGr flake and the isolated gold wires are represented by a black curve and a red curve, respectively. Fourier transform of the induced charge density for the DVDGr-AuX: (e) DVDGr-Au2, (f) DVDGr-Au3, (g) DVDGr-Au4, and (h) DVDGr-Au5, at the energy resonance points 3.51,2.95,2.33,and 2.06 eV,respectively.

    3.2. Plasmons of the mixed noble metal-TM arrays with the double vacancy defected graphene in planar hybrid structures

    The above-mentioned studies addressed that increasing the number of Au atoms further intensifies the optical absorption and causes a slowly saturating redshift. In contrast, considering an arrangement of mixed metallic nanoarrays,we put the DVDGr flake and Au cluster doped with TM atom Pd together,named as DVDGr-AuXPdY,where X (X =6,4,4,and 2) and Y (Y =4, 4, 2, and 2) represent the numbers of Au and Pd. The distance of 2.884 ?A is employed between Pd and Au. The assignments of configuration and optical absorption of planar hybrid structures with impulse excitation polarizing along the armchair edge are shown in Fig. 3. Resonance absorption of four DVDGr-AuXPdY (X+Y =4,6,8 and 10)is 4.40,3.31,2.64 and 2.37 eV,respectively.It is well understood that an interplay between spill out and interband polarization control the actual sign of the plasmon shift. Presumably, this redshift is related to electron spill out. The plasmon energy increases with the nanowire amount. However,in the situation of the planar arrays of 4 Au and 4 Pd clusters, the plasmon energy does not increase but moderately decreases. The nonmonotonic behavior for the configuration, which differs from the pure Au cluster,is owing to the quantum mechanical confinement of the average electron density between the mixed noble metal cluster and transition metal cluster impacted by the introduction of Pd. The photoabsorption spectrum of the DVDGr-AuXPdY is significantly suppressed and fragmented compared to the DVDGr-AuX. Despite collective effects enhance, this weak damping slows the increase of the plasmon intensity, indicating that Au cluster doped with TM atom Pd can reduce the plasmon energy on account of the more effective screening by the TM Pd.

    Fig. 3. Optical absorption of DVDGr-AuXPdY consisting of the same DVDGr flake and different Au, Pd clusters. The impulse excitation is polarized along the armchair edge. From bottom to top on the left: DVDGr-Au2Pd2, DVDGr-Au4Pd2, DVDGr-Au4Pd4, and DVDGr-Au6Pd10,respectively.

    Leaving the impulse excitation polarized along zigzag edge,the plasmon blueshifts slightly with the increased number of clusters, as seen in Fig. 4. Resonance absorptions of four DVDGr-AuXPdY (X+Y =4, 6, 8 and 10) at lower energy mode are located at 4.37, 4.34, 4.39 and 4.4 eV, respectively,and at higher energy mode they are 5.85,6.36,6.2 and 6.44 eV. Obviously, the doping of transition-metal generates additional excitations in the absorption spectrum, leading to the occurrence of double strong absorption spectrum. The optical absorption is broadened overwhelmingly that enables graphene-made quantum devices to have a wider range of light absorption spectra.

    Fig.4. Optical absorption of DVDGr-AuXPdY consisting of the same DVDGr flake and different Au, Pd clusters. The impulse excitation is polarized along zigzag edge. From bottom to top on the left: DVDGr-Au2Pd2, DVDGr-Au4Pd2, DVDGr-Au4Pd4, and DVDGr-Au6Pd10,respectively.

    Compared with the mixed Au/Pd nanoarrays, pure Au nanoarrays have stronger plasmonic enhancement. Thus, the gold nanoarrays are expected to take as optical antenna, it presents a strong ability to capture visible lights.

    3.3. Plasmons in non-planar hybrid structures

    We analyze the evolution process of plasmon in the nonplanar asymmetrical DVDGr-AuX which consists of a higher symmetry DVDGr(approximate hexagonal)and the different number of Au atoms–X (X = 2, 3, 4, and 5). The highersymmetric DVDGr has a single shape boundary.[35]The gold wire is perpendicular to the DVDGr plane(yz plane)and it is mapped at the center of the octagon.The polarization direction of the incident light is maintained along the x-axis. Figure 5 displays the the DVDGr-Au5 hybrid structure and absorption spectrum of DVDGr-AuX (X =2,3,4,and 5). The coupling excitation of photon and phonon (graphene plasmon) can result in photon-electron phonon interactions,which can significantly enhance resonance Raman scattering of mid-IR region,which is the so-called graphene plasmon-enhanced resonance Raman scattering of mid-IR region.[54]Focus on the optical absorption in mid-IR (below 1 eV) non-planar asymmetrical DVDGr-AuX structures, see Fig.5. It is found that there is a new absorption peak in the infrared frequency ranges, which not only means that this non-planar hybrid structure has a very wide optical absorption but also makes clear that a new plasmon mode has appeared.

    Fig. 5. Optical absorption of DVDGr-AuX consisting of the same DVDGr and the different numbers of gold atoms X (X =2, 3, 4, and 5). The isolated DVDGr is displayed in the inset graphs.

    We also analyze the plasmon of symmetrical DVDGr-AuX (X =2, 3, 4, and 5). The polarization direction of incident light is along the x-axis. Figures 6(a)–6(d) display the optical absorption of the four symmetrical hybrid structures. Compared to the asymmetric structure, symmetrical hybrid structures present a pronounced peak in infrared frequency ranges, the peak value is rising a little. We can also find that the optical absorption spectra are wide in regard of visible,infrared or terahertz regions from Fig.6(a)–6(d). Such wide ranges bring the non-planar symmetrical DVDGr-AuX a possibility to be used as an euryphotic photosensitive detector.Furthermore,increasing the number of Au atoms further intensifies the optical absorption and causes a new absorption peak in mid-IR(below 1 eV).The photoinduced charge transfer of DVDGr and Au can occur in electronic state transition of mid-IR region,and then the transferred electrons can interact with the phonon of DVDGr at mid-IR.

    Fig.6. Optical absorptions for non-planar symmetrical DVDGr-AuX (X =2,3,4,and 5)shown in(a),(b),(c),(d). Red represents the higher symmetry DVDGr with different Au wires (X =2, 3, 4, and 5), and black represents the asymmetrical DVDGr-AuX (X =2, 3, 4, and 5).Fourier transform of the induced charge density for the DVDGr-Au5,at the energy resonance points 2.0 eV:(e)the xy plane,(f)the yz plane.

    Figures 6(e)–6(f) show the Fourier transform of the induced charge density distribution of the symmetric DVDGr-Au5 in xy plane and yz plane. The polarization direction of incident light is perpendicular to DVDGr. The electrons are apparently clustered on the surface of graphene. Due to the shielding effect of interface charge,the antiparallel dipole appears on both sides of graphene, the electric field becomes weaker,and the peak value decreases.

    Previous studies indicated that graphene boundary has a large effect on the plasmon mode and the light response and acts as a charge transfer channel. However, it is difficult to achieve such a function in the non-planar hybrid structure because it is very far away from the gold wire. Finally, the change of antenna size causes the variety of eigen frequency,which further results in a change of the optical response of the hybrid structure. Thus we can change the gold antenna to precisely regulate the range and intensity of the plasmon.

    3.4. Plasmon-induced nonlinear responses

    The couple’s excitation between photon and phonon can result in the graphene plasmon, and the photon-electronphonon interaction can result in the graphene plasmonenhanced resonance Raman scattering,which can be enhanced up to 107 times.[54]Thus, we study the DVDGr plasmoninduced nonlinear responses to realize the extension and enhancement of the absorption. In this section, we adopt three kinds of incident laser beams to study the interaction of ultrafast lasers with DVDGr-Au5. The first and second for a Gaussian-like pulse and the third case for a totally diverse laser,as shown in Figs.7(a),7(c),and 7(e). Similar harmonic spectra are seen in Figs. 7(b) and 7(d) for the same choice of parameters but different amplitudes for the intensity. It is evident that the strength of the laser pulse constitutes important issues for the generation of harmonic. A dominant feature of the two spectra is the strong peak emitted at nearly the frequency that corresponds precisely to the plasmon frequency of DVDGr-Au5. For the case of the last one,the wavelength and intensity of the laser are chosen as λ =248 nm and I =5.3×1014W/cm2, respectively. An interesting feature of the spectra in Fig.7(f)is that the harmonic gets excited at 15 eV due to nonlinear coupling with external perturbations.

    Fig.7. Laser oscillograph and the corresponding harmonic spectra. Laser beams shown in(a),(c),(e)are for various wavelengths of 466,466 and 248(nm)and various intensities of 1.327×1014,1.327×1015,5.30×1014 W/cm2.

    To gain a deeper insight into DVDGr-Au5 plasmoninduced nonlinear responses, a set of simulations were performed with the purpose of exploring the effects of pulse shaping. The lasers ranging from ultraviolet to infrared with the same intensity corresponding to the center wavelengths of 138,414 and 1379 nm respectively were chosen. Considering the momentously influence of the field intensity, four combined laser fields with peak intensity of 2×1013, 6×1013, 1×1014,1.45×1014W/cm2were selected. Figure 8 shows the laser oscillograph with the intensity of 1.23 eV/?A.The actual effect of ultrafast laser is 6 fs.

    Fig. 8. Lasers with a controlling pulse oscillograph. Lasers with various wavelengths of 138,414 and 1379(nm)in the intensity of 1.23 eV/?A.Controlling pulse with the intensity of 2 eV/?A and wavelength of 620 nm.

    The acceleration of the dipole moment is governed by the strength and the shape of the laser, hence the harmonic will be impacted dramatically by introducing a diverse type and shape laser as a modulation factor in a mixed field. As discussed above, Gaussian-like pulse can effectively boost the intensity of the harmonic spectrum. Thus, Gaussianlike pulses that performed as regulatory fields with the center wavelengths of 620 nm and peak intensity of 5.31×1013,8.97×1013, 1.2×1014, 1.63×1014W/cm2are applied to the system. (Fig.9 pink line).

    Fig.9. Harmonic spectrum of DVDGr flake and DVDGr-Au5 obtained from the pulses in Figs. 8. (a) The DVDGr flake without Gaussianlike pulse. (b) The DVDGr flake with Gaussian-like pulse. (c) The DVDGr-Au5 flake without Gaussian-like pulse. (d) The DVDGr flake with Gaussian-like pulse.

    In Fig. 9 we compare the harmonic spectrum of the DVDGr flake and DVDGr-Au5 obtained from the pulses in Figs.8.The harmonic emission spectra show a double-plateau structure,the width increases mildly and the peak value boosts monotonously and non-negligibly with the augment of field intensity. We observe that,for the case in which the laser pulse is shaped, the harmonic at high energy mode can be excited,the electron energy varies according to the modulation of the laser electric field. Comparing the plasmon-induced nonlinear responses in DVDGr flake and DVDGr-Au, we note that the Au has the effect of enhancing harmonics in low-energy mode while suppressing emission in high-energy mode. A straightforward interpretation of these simulations is that the emission can be tuned by the resonant frequency of gold nanoarrays.We also find that the emission is instructive in comparison with the harmonic generated from lasers with and without the laser electric field shaped by a Gaussian envelope. As shown in Fig. 9(d), the double-plateau of the harmonic spectrum is extraordinarily increased. In particular,the peak value in low energy mode even exceeds that in high energy mode. The latter introduces sharp spectral features characterized by the controlling laser, implying that changing the wavelength of the controlling laser can touch its quantum path, boost the specific harmonic peak. Furthermore,the modulation of the field brings about the spectrum becoming much smoother for the harmonics after the second plateau. This complexity provides an opportunity to tune nonlinear optical properties by combining disparate type and shape of external perturbations.

    4. Conclusions

    This work provides the first detailed theoretical analysis of DVD-Gr with multiple pure and mixed noble/TM nanoarrays using TDDFT. Effects of different structure on absorption spectrum and effects of various electronic field on harmonic spectrum have been investigated. It is noted that by adding metallic nanoarrays on double vacancy defected graphene plane, the plasmons produce coherent superposition, which enhances the absorption. This makes it possible to tune the plasmons in the hybrid nanostructures. The gold nanoarrays are expected to as optical antenna, and it presents a stronger ability to capture visible lights than mixed Au/Pd nanoarrays. Furthermore, DVDGr-Au plasmon-induced nonlinear responses suggest that for moderately intense femtosecond lasers, harmonic emission corresponds approximately to the plasmons frequency. We provide an effective method to realize the extension and enhancement of specific peak for the harmonic. This study is helpful to build a theoretical basis for understanding the microscopic mechanism of plasmoninduced nonlinear response and achieving ultimate control at the nanoscale.

    Acknowledgment

    We thank ShiJia Yuan and Tong Liu for fruitful discussions.

    猜你喜歡
    張紅新路
    Core level excitation spectra of La and Mn ions in LaMnO3
    Polyhedral silver clusters as single molecule ammonia sensor based on charge transfer-induced plasmon enhancement
    Theoretical study of M6X2 and M6XX'structure(M =Au,Ag;X,X'=S,Se): Electronic and optical properties,ability of photocatalytic water splitting,and tunable properties under biaxial strain
    李學奇《新路》
    水土保持探新路 三十九年寫春秋
    Stability,electronic structure,and optical properties of lead-free perovskite monolayer Cs3B2X9(B=Sb,Bi;X =Cl,Br,I)and bilayer vertical heterostructure Cs3B2X9/Cs3X9(B,B′=Sb,Bi;X =Cl,Br,I)
    Linear and nonlinear optical response of g-C3N4-based quantum dots*
    我愛你,中國
    蔬果種植走新路
    城鄉(xiāng)一體化走出的新路
    春色校园在线视频观看| 亚洲av一区综合| 亚洲综合精品二区| 在线亚洲精品国产二区图片欧美 | 亚洲一区二区三区欧美精品 | 又黄又爽又刺激的免费视频.| a级一级毛片免费在线观看| 免费看光身美女| 亚洲精品视频女| 国产欧美另类精品又又久久亚洲欧美| 亚洲综合精品二区| 国产 一区 欧美 日韩| 伊人久久精品亚洲午夜| 国产午夜精品久久久久久一区二区三区| 五月玫瑰六月丁香| 天天躁日日操中文字幕| 91精品国产九色| 精品一区二区三卡| 老司机影院毛片| 中文欧美无线码| 中文精品一卡2卡3卡4更新| 少妇猛男粗大的猛烈进出视频 | 亚洲成人中文字幕在线播放| 在线观看三级黄色| 色婷婷久久久亚洲欧美| 男女那种视频在线观看| 午夜福利在线观看免费完整高清在| 国产一区二区在线观看日韩| 91在线精品国自产拍蜜月| 久久久久精品性色| 99久久精品一区二区三区| 免费观看的影片在线观看| 久久久精品欧美日韩精品| 日韩大片免费观看网站| 乱码一卡2卡4卡精品| 我要看日韩黄色一级片| 久久精品久久久久久噜噜老黄| 日韩制服骚丝袜av| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产a三级三级三级| 美女视频免费永久观看网站| 国产 精品1| 狂野欧美激情性xxxx在线观看| av福利片在线观看| 色网站视频免费| 国产精品一区二区三区四区免费观看| 亚洲三级黄色毛片| 别揉我奶头 嗯啊视频| 建设人人有责人人尽责人人享有的 | 三级国产精品欧美在线观看| 精品少妇黑人巨大在线播放| 日本三级黄在线观看| 免费看av在线观看网站| 一个人观看的视频www高清免费观看| 国产av国产精品国产| 一级毛片电影观看| 国产精品国产三级国产av玫瑰| 国产一区亚洲一区在线观看| 国产一区二区在线观看日韩| 日韩av在线免费看完整版不卡| 国产精品99久久久久久久久| 极品教师在线视频| 特大巨黑吊av在线直播| 最近的中文字幕免费完整| 亚洲国产精品专区欧美| 久热这里只有精品99| av在线播放精品| 亚洲av欧美aⅴ国产| 男人狂女人下面高潮的视频| 日产精品乱码卡一卡2卡三| 免费人成在线观看视频色| 国产免费视频播放在线视频| 搞女人的毛片| 蜜桃亚洲精品一区二区三区| 欧美zozozo另类| 欧美日韩视频精品一区| 国产色婷婷99| 色视频www国产| 美女视频免费永久观看网站| 丰满少妇做爰视频| 久久精品久久精品一区二区三区| 亚洲欧洲国产日韩| 国产一区二区三区av在线| 五月开心婷婷网| 精品久久久久久久人妻蜜臀av| 久久久久久久久久成人| 亚洲国产精品999| 内射极品少妇av片p| 成人欧美大片| 成人一区二区视频在线观看| 亚洲av在线观看美女高潮| 国产91av在线免费观看| 成人鲁丝片一二三区免费| 在线免费十八禁| 国模一区二区三区四区视频| 亚洲国产精品成人久久小说| 久久精品久久久久久久性| 欧美 日韩 精品 国产| 久久久久性生活片| 亚洲av不卡在线观看| 成人亚洲精品av一区二区| 亚洲成人一二三区av| 亚洲欧美一区二区三区国产| 小蜜桃在线观看免费完整版高清| av在线app专区| 美女高潮的动态| 在线观看一区二区三区| 国产乱人视频| 亚洲精品视频女| 亚洲欧美日韩无卡精品| 亚洲精品aⅴ在线观看| 久久久久久久久久人人人人人人| 小蜜桃在线观看免费完整版高清| 日韩,欧美,国产一区二区三区| 伊人久久国产一区二区| 青青草视频在线视频观看| eeuss影院久久| 亚洲欧美日韩东京热| 中文字幕制服av| av国产久精品久网站免费入址| www.色视频.com| 国产白丝娇喘喷水9色精品| 国产乱人偷精品视频| 性色avwww在线观看| 看非洲黑人一级黄片| 国产成人freesex在线| 99热这里只有精品一区| 六月丁香七月| 五月伊人婷婷丁香| 白带黄色成豆腐渣| 国产精品久久久久久久久免| 亚洲国产精品成人综合色| 少妇人妻久久综合中文| 亚洲精品色激情综合| 听说在线观看完整版免费高清| 我的女老师完整版在线观看| 青青草视频在线视频观看| 在线精品无人区一区二区三 | 国产一区有黄有色的免费视频| 欧美潮喷喷水| 亚洲经典国产精华液单| 国产有黄有色有爽视频| 成人毛片a级毛片在线播放| 亚洲色图av天堂| 边亲边吃奶的免费视频| 亚洲精品久久久久久婷婷小说| 亚洲欧洲国产日韩| 日日撸夜夜添| 国产国拍精品亚洲av在线观看| 国国产精品蜜臀av免费| 美女xxoo啪啪120秒动态图| 日日撸夜夜添| av又黄又爽大尺度在线免费看| 女人十人毛片免费观看3o分钟| 国产成人a∨麻豆精品| 国产片特级美女逼逼视频| 熟女人妻精品中文字幕| 国产欧美日韩一区二区三区在线 | 亚洲第一区二区三区不卡| 伦理电影大哥的女人| 国产白丝娇喘喷水9色精品| 免费av观看视频| 国产大屁股一区二区在线视频| 久久精品久久久久久久性| 热99国产精品久久久久久7| 能在线免费看毛片的网站| 欧美高清性xxxxhd video| 免费大片黄手机在线观看| 波野结衣二区三区在线| 亚洲av成人精品一二三区| 久久97久久精品| 三级国产精品片| 国产爽快片一区二区三区| 秋霞伦理黄片| 欧美精品一区二区大全| 久久99热这里只有精品18| 国产高清国产精品国产三级 | 久久久精品欧美日韩精品| 精品久久久精品久久久| 日韩 亚洲 欧美在线| 青春草国产在线视频| 免费黄频网站在线观看国产| 国产精品一区二区在线观看99| av在线蜜桃| 日韩中字成人| 久久久国产一区二区| 久久97久久精品| 亚洲国产色片| 亚洲精品乱码久久久v下载方式| 国产精品麻豆人妻色哟哟久久| 丝瓜视频免费看黄片| 极品教师在线视频| 欧美3d第一页| 女人被狂操c到高潮| 男女那种视频在线观看| 国产熟女欧美一区二区| 国产精品久久久久久精品电影| 久久综合国产亚洲精品| 国产69精品久久久久777片| 欧美97在线视频| 国产成人精品福利久久| 日韩伦理黄色片| 国产一区二区亚洲精品在线观看| 亚洲国产精品成人久久小说| 欧美潮喷喷水| 午夜福利高清视频| 欧美另类一区| 最近2019中文字幕mv第一页| 久久这里有精品视频免费| 99视频精品全部免费 在线| 少妇丰满av| 精品一区二区三区视频在线| 成人亚洲精品一区在线观看 | 亚洲av电影在线观看一区二区三区 | 成年女人看的毛片在线观看| 免费观看无遮挡的男女| 国产精品人妻久久久久久| 国产伦精品一区二区三区视频9| 成年人午夜在线观看视频| 十八禁网站网址无遮挡 | 欧美日韩在线观看h| 最后的刺客免费高清国语| 欧美变态另类bdsm刘玥| 最近最新中文字幕免费大全7| 少妇人妻精品综合一区二区| 你懂的网址亚洲精品在线观看| 亚洲最大成人av| 久久久久久久久久久丰满| 天天躁夜夜躁狠狠久久av| 免费观看无遮挡的男女| 亚洲三级黄色毛片| 日韩中字成人| 国产精品一二三区在线看| 在线天堂最新版资源| 深爱激情五月婷婷| 亚洲丝袜综合中文字幕| 国产精品一及| 久久久久久久大尺度免费视频| 人妻夜夜爽99麻豆av| 视频中文字幕在线观看| 成年女人在线观看亚洲视频 | 亚洲在线观看片| 国产精品无大码| 2021少妇久久久久久久久久久| 特级一级黄色大片| 91久久精品国产一区二区三区| 国产久久久一区二区三区| 中文字幕av成人在线电影| 男男h啪啪无遮挡| 国产成人午夜福利电影在线观看| 国产精品人妻久久久影院| 蜜桃亚洲精品一区二区三区| freevideosex欧美| 老司机影院成人| www.av在线官网国产| 国产乱人视频| 男插女下体视频免费在线播放| 国产精品99久久久久久久久| av一本久久久久| av在线app专区| 91aial.com中文字幕在线观看| 欧美成人一区二区免费高清观看| 一个人看视频在线观看www免费| 18禁裸乳无遮挡动漫免费视频 | 亚洲国产精品专区欧美| 成人二区视频| 亚洲精品日韩在线中文字幕| 在现免费观看毛片| 在线天堂最新版资源| 视频中文字幕在线观看| 亚洲精华国产精华液的使用体验| 欧美精品国产亚洲| 尤物成人国产欧美一区二区三区| 3wmmmm亚洲av在线观看| 午夜爱爱视频在线播放| 寂寞人妻少妇视频99o| 啦啦啦中文免费视频观看日本| 国产免费又黄又爽又色| 久久久久久伊人网av| 少妇熟女欧美另类| 日韩人妻高清精品专区| 一级毛片我不卡| 尾随美女入室| 欧美少妇被猛烈插入视频| 中文在线观看免费www的网站| 成人黄色视频免费在线看| 国产精品一区www在线观看| 一个人看的www免费观看视频| 我要看日韩黄色一级片| 午夜精品一区二区三区免费看| 少妇被粗大猛烈的视频| 久热久热在线精品观看| 亚洲自拍偷在线| 亚洲国产精品成人久久小说| 亚洲欧美一区二区三区黑人 | 91久久精品国产一区二区三区| av在线老鸭窝| 久久久久久久久大av| 免费看av在线观看网站| 自拍偷自拍亚洲精品老妇| 日韩免费高清中文字幕av| 久久鲁丝午夜福利片| 综合色av麻豆| 久热这里只有精品99| 国产精品国产三级国产av玫瑰| 18禁动态无遮挡网站| 日韩一本色道免费dvd| 男插女下体视频免费在线播放| 大香蕉久久网| 午夜亚洲福利在线播放| 国产成人一区二区在线| 噜噜噜噜噜久久久久久91| 亚洲三级黄色毛片| 国产亚洲91精品色在线| 国产日韩欧美亚洲二区| 菩萨蛮人人尽说江南好唐韦庄| 超碰97精品在线观看| 亚洲精品乱码久久久久久按摩| 亚洲欧美清纯卡通| 人体艺术视频欧美日本| 少妇猛男粗大的猛烈进出视频 | 插阴视频在线观看视频| 国产高清不卡午夜福利| 六月丁香七月| 欧美最新免费一区二区三区| 国产成人a∨麻豆精品| 久久精品国产亚洲av涩爱| 久久久久久久午夜电影| 亚洲丝袜综合中文字幕| 一级a做视频免费观看| av在线老鸭窝| 免费不卡的大黄色大毛片视频在线观看| 国产爽快片一区二区三区| 亚洲欧美日韩无卡精品| 亚洲精品中文字幕在线视频 | 另类亚洲欧美激情| 精品一区在线观看国产| 国产精品国产av在线观看| 中文精品一卡2卡3卡4更新| 乱码一卡2卡4卡精品| 一区二区av电影网| 亚洲人成网站在线播| 国产一区二区三区av在线| 秋霞在线观看毛片| 国产极品天堂在线| 99热全是精品| 国产精品福利在线免费观看| 热re99久久精品国产66热6| 久久精品人妻少妇| 最近2019中文字幕mv第一页| 日韩av不卡免费在线播放| 好男人在线观看高清免费视频| 插逼视频在线观看| 亚洲精品国产av成人精品| 久久女婷五月综合色啪小说 | kizo精华| 深爱激情五月婷婷| 99re6热这里在线精品视频| 狂野欧美激情性xxxx在线观看| 久热这里只有精品99| 九色成人免费人妻av| 亚洲在线观看片| 免费不卡的大黄色大毛片视频在线观看| 久久久久久久久久久丰满| 亚洲欧洲国产日韩| 高清在线视频一区二区三区| 色吧在线观看| 国产成人a区在线观看| 国产成人91sexporn| 亚洲精品乱码久久久v下载方式| 久久精品久久久久久噜噜老黄| 欧美成人精品欧美一级黄| 亚洲成人精品中文字幕电影| 免费看日本二区| 国产精品无大码| 干丝袜人妻中文字幕| 中文欧美无线码| 最近2019中文字幕mv第一页| 午夜福利在线在线| 女人被狂操c到高潮| 日韩欧美精品v在线| 国产精品久久久久久精品古装| 国产精品偷伦视频观看了| 国产视频内射| 国产真实伦视频高清在线观看| 国产亚洲一区二区精品| 欧美日韩在线观看h| 国产69精品久久久久777片| 国产精品成人在线| 一级av片app| 亚洲熟女精品中文字幕| 亚洲国产成人一精品久久久| 国产精品av视频在线免费观看| 亚洲国产成人一精品久久久| 国产精品av视频在线免费观看| av线在线观看网站| 亚洲av男天堂| 国产色婷婷99| 91久久精品国产一区二区三区| 久久久欧美国产精品| 国产视频首页在线观看| 国产毛片在线视频| 亚洲,一卡二卡三卡| 国产毛片a区久久久久| 免费观看无遮挡的男女| 最近中文字幕2019免费版| 最后的刺客免费高清国语| 高清午夜精品一区二区三区| 国产片特级美女逼逼视频| 欧美精品一区二区大全| 久久久久性生活片| 日本午夜av视频| 在线亚洲精品国产二区图片欧美 | 免费少妇av软件| 免费电影在线观看免费观看| 如何舔出高潮| 欧美最新免费一区二区三区| 在线观看一区二区三区激情| 99精国产麻豆久久婷婷| 97超视频在线观看视频| 女的被弄到高潮叫床怎么办| 蜜臀久久99精品久久宅男| 老司机影院成人| 久久久久久久久久久免费av| 亚洲最大成人av| 国产69精品久久久久777片| 久久鲁丝午夜福利片| 又粗又硬又长又爽又黄的视频| 国产人妻一区二区三区在| av国产久精品久网站免费入址| 亚洲丝袜综合中文字幕| 国内揄拍国产精品人妻在线| 欧美xxxx性猛交bbbb| 久久久久久久久久人人人人人人| 国产精品精品国产色婷婷| 亚洲天堂国产精品一区在线| 国产黄色免费在线视频| 亚洲最大成人手机在线| 久久久久久久久久成人| 日韩一区二区视频免费看| 日韩av免费高清视频| 80岁老熟妇乱子伦牲交| 国产精品爽爽va在线观看网站| 91狼人影院| 综合色丁香网| av天堂中文字幕网| 国产老妇伦熟女老妇高清| 国产爱豆传媒在线观看| 免费看av在线观看网站| 99久久精品热视频| 国产探花在线观看一区二区| 久久影院123| 精品少妇黑人巨大在线播放| 久久韩国三级中文字幕| 肉色欧美久久久久久久蜜桃| 一级a爱视频在线免费观看| 久久性视频一级片| 熟妇人妻不卡中文字幕| av又黄又爽大尺度在线免费看| 久热这里只有精品99| 久久综合国产亚洲精品| 欧美最新免费一区二区三区| 男人爽女人下面视频在线观看| 婷婷色麻豆天堂久久| 妹子高潮喷水视频| 国产成人91sexporn| 美女午夜性视频免费| 欧美国产精品一级二级三级| 男女床上黄色一级片免费看| 亚洲婷婷狠狠爱综合网| 欧美精品一区二区大全| 国产精品久久久久久精品古装| 亚洲美女视频黄频| 极品少妇高潮喷水抽搐| 日本猛色少妇xxxxx猛交久久| 欧美黄色片欧美黄色片| 午夜日韩欧美国产| 国产成人精品久久二区二区91 | 另类亚洲欧美激情| 日日爽夜夜爽网站| 亚洲欧美一区二区三区黑人| 麻豆精品久久久久久蜜桃| 日韩伦理黄色片| 亚洲国产毛片av蜜桃av| 亚洲欧洲精品一区二区精品久久久 | 国产精品久久久久久人妻精品电影 | 99精国产麻豆久久婷婷| 国产精品99久久99久久久不卡 | 国产乱人偷精品视频| 成人三级做爰电影| 欧美人与善性xxx| 王馨瑶露胸无遮挡在线观看| 午夜福利网站1000一区二区三区| 亚洲av综合色区一区| 精品国产一区二区久久| 国产淫语在线视频| 午夜免费鲁丝| 你懂的网址亚洲精品在线观看| 一级a爱视频在线免费观看| 国产熟女午夜一区二区三区| 欧美av亚洲av综合av国产av | 亚洲中文av在线| 99久国产av精品国产电影| 又大又爽又粗| 久久久久久人妻| 午夜激情久久久久久久| 亚洲久久久国产精品| 女人爽到高潮嗷嗷叫在线视频| 国产精品 国内视频| 免费看av在线观看网站| 少妇人妻久久综合中文| 赤兔流量卡办理| 久久 成人 亚洲| 亚洲七黄色美女视频| 丰满饥渴人妻一区二区三| 精品国产一区二区久久| 久久久精品免费免费高清| 人人妻人人添人人爽欧美一区卜| 你懂的网址亚洲精品在线观看| 精品少妇一区二区三区视频日本电影 | 久久热在线av| 精品国产一区二区三区久久久樱花| 欧美精品亚洲一区二区| 日本一区二区免费在线视频| 美女中出高潮动态图| 天堂8中文在线网| 欧美日韩亚洲综合一区二区三区_| 久久久久国产精品人妻一区二区| 人人妻人人澡人人看| 日韩 欧美 亚洲 中文字幕| 日韩成人av中文字幕在线观看| 一级毛片黄色毛片免费观看视频| 黄网站色视频无遮挡免费观看| 亚洲伊人久久精品综合| 国产精品久久久人人做人人爽| 丁香六月欧美| 国产精品久久久久久人妻精品电影 | 免费看不卡的av| 精品第一国产精品| 精品国产露脸久久av麻豆| 国产乱来视频区| 国产成人免费观看mmmm| 欧美黄色片欧美黄色片| 欧美人与性动交α欧美软件| 国产精品一国产av| 欧美在线黄色| 中文字幕av电影在线播放| 午夜福利网站1000一区二区三区| 中文字幕亚洲精品专区| 免费久久久久久久精品成人欧美视频| 大香蕉久久网| 国产精品免费视频内射| 精品久久久久久电影网| 宅男免费午夜| 最近手机中文字幕大全| 日本av免费视频播放| 国产日韩欧美在线精品| 精品第一国产精品| 久久人妻熟女aⅴ| 黄色 视频免费看| 大片电影免费在线观看免费| 精品酒店卫生间| 国产 精品1| 国产麻豆69| a级片在线免费高清观看视频| 亚洲国产最新在线播放| 汤姆久久久久久久影院中文字幕| 91成人精品电影| 欧美日韩视频高清一区二区三区二| 成年美女黄网站色视频大全免费| 亚洲国产欧美日韩在线播放| 成年女人毛片免费观看观看9 | 国产av精品麻豆| 麻豆精品久久久久久蜜桃| 18禁裸乳无遮挡动漫免费视频| 99re6热这里在线精品视频| 五月天丁香电影| 在线观看免费午夜福利视频| 亚洲精品久久成人aⅴ小说| 精品少妇内射三级| 欧美亚洲 丝袜 人妻 在线| 日本91视频免费播放| 亚洲情色 制服丝袜| 亚洲七黄色美女视频| 午夜精品国产一区二区电影| 久久久精品94久久精品| 波野结衣二区三区在线| 国产亚洲一区二区精品| 少妇的丰满在线观看| 捣出白浆h1v1| 国产黄色视频一区二区在线观看| 99久久精品国产亚洲精品| 深夜精品福利| 9热在线视频观看99| 极品人妻少妇av视频| 深夜精品福利| 下体分泌物呈黄色| 少妇人妻 视频| 黄片小视频在线播放| 又黄又粗又硬又大视频| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看免费视频网站a站| av国产久精品久网站免费入址| 黄片无遮挡物在线观看| 日韩精品免费视频一区二区三区| 精品卡一卡二卡四卡免费| 只有这里有精品99| 青春草视频在线免费观看| 午夜日韩欧美国产| 两性夫妻黄色片| 秋霞伦理黄片| 亚洲综合精品二区| 99香蕉大伊视频| 久久久欧美国产精品| 国产精品熟女久久久久浪|