• <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免费观看日本| 国产亚洲最大av| 婷婷色综合www| 午夜激情福利司机影院| 啦啦啦啦在线视频资源| 国产伦理片在线播放av一区| 夫妻性生交免费视频一级片| 午夜视频国产福利| 91精品国产九色| 亚洲av欧美aⅴ国产| 欧美xxⅹ黑人| 国产午夜福利久久久久久| 成年免费大片在线观看| 国精品久久久久久国模美| 男女国产视频网站| 国内少妇人妻偷人精品xxx网站| 日韩一本色道免费dvd| 欧美日韩视频精品一区| 国产精品一及| 国产精品女同一区二区软件| 免费黄色在线免费观看| av黄色大香蕉| 亚洲成人av在线免费| 亚洲一区二区三区欧美精品 | 免费在线观看成人毛片| 最近手机中文字幕大全| 国产av不卡久久| 国产精品一及| 国产国拍精品亚洲av在线观看| 爱豆传媒免费全集在线观看| 午夜福利在线在线| 亚洲欧美精品自产自拍| 国产黄色免费在线视频| 欧美日韩视频精品一区| 久久久久久伊人网av| 久久精品久久久久久噜噜老黄| 亚洲av福利一区| 亚洲欧洲国产日韩| 天堂中文最新版在线下载 | 亚洲精品乱码久久久v下载方式| 97在线人人人人妻| 国产在线男女| 精品人妻一区二区三区麻豆| 在线看a的网站| 看非洲黑人一级黄片| 亚洲一级一片aⅴ在线观看| 亚洲欧洲国产日韩| 看免费成人av毛片| 99久久人妻综合| 中文字幕制服av| 老师上课跳d突然被开到最大视频| 国产亚洲av片在线观看秒播厂| 97超碰精品成人国产| av网站免费在线观看视频| 一本色道久久久久久精品综合| 国产视频首页在线观看| 99久久精品国产国产毛片| 久久99热这里只频精品6学生| 又粗又硬又长又爽又黄的视频| 亚洲自偷自拍三级| 在线观看一区二区三区| www.色视频.com| 欧美日韩视频高清一区二区三区二| 好男人视频免费观看在线| 久久精品国产a三级三级三级| 国产在线男女| 黄色配什么色好看| 亚洲精品久久久久久婷婷小说| 久久久成人免费电影| 国产成人免费观看mmmm| 成人特级av手机在线观看| 国产人妻一区二区三区在| 亚洲丝袜综合中文字幕| 国产黄频视频在线观看| 日本爱情动作片www.在线观看| 日本黄色片子视频| 伦精品一区二区三区| 亚洲精品日韩av片在线观看| 一级毛片我不卡| 新久久久久国产一级毛片| 国产老妇女一区| 爱豆传媒免费全集在线观看| 人妻少妇偷人精品九色| 99re6热这里在线精品视频| 日本午夜av视频| av国产精品久久久久影院| 晚上一个人看的免费电影| 日韩,欧美,国产一区二区三区| 国产极品天堂在线| 一二三四中文在线观看免费高清| 日日啪夜夜爽| 午夜福利网站1000一区二区三区| 日韩强制内射视频| 成年免费大片在线观看| 亚洲成人av在线免费| 免费不卡的大黄色大毛片视频在线观看| 男插女下体视频免费在线播放| 又黄又爽又刺激的免费视频.| 亚洲欧美精品自产自拍| 18禁在线播放成人免费| 成人二区视频| 夜夜看夜夜爽夜夜摸| 欧美精品国产亚洲| 男人爽女人下面视频在线观看| 联通29元200g的流量卡| 亚洲精品,欧美精品| 97在线视频观看| 天天躁日日操中文字幕| a级一级毛片免费在线观看| 免费高清在线观看视频在线观看| 一二三四中文在线观看免费高清| tube8黄色片| 18+在线观看网站| 亚洲丝袜综合中文字幕| 亚洲aⅴ乱码一区二区在线播放| 在线观看美女被高潮喷水网站| 一本久久精品| 国产免费一区二区三区四区乱码| 99久国产av精品国产电影| 成年女人看的毛片在线观看| 成人免费观看视频高清| 日韩在线高清观看一区二区三区| 亚洲欧美成人精品一区二区| 中文资源天堂在线| 男女边吃奶边做爰视频| 久久久国产一区二区| 亚洲国产精品成人综合色| 午夜福利视频1000在线观看| 久久久久久久国产电影| 国产成人午夜福利电影在线观看| 久久久久久久久久久免费av| 成人亚洲精品一区在线观看 | 免费大片18禁| 性色avwww在线观看| 尤物成人国产欧美一区二区三区| 五月开心婷婷网| 精品国产乱码久久久久久小说| 国产精品久久久久久av不卡| 久久精品国产亚洲av天美| 内射极品少妇av片p| 蜜臀久久99精品久久宅男| 久久久国产一区二区| 男插女下体视频免费在线播放| 午夜精品国产一区二区电影 | 亚洲欧美精品自产自拍| 成人黄色视频免费在线看| 97超视频在线观看视频| 视频中文字幕在线观看| 欧美成人a在线观看| 最后的刺客免费高清国语| 亚洲婷婷狠狠爱综合网| 国产片特级美女逼逼视频| 麻豆成人午夜福利视频| 国产精品伦人一区二区| 乱码一卡2卡4卡精品| 18+在线观看网站| 亚洲av二区三区四区| 亚洲av日韩在线播放| 国产极品天堂在线| 91aial.com中文字幕在线观看| 一边亲一边摸免费视频| 久久久久九九精品影院| 欧美成人精品欧美一级黄| 国产69精品久久久久777片| 搞女人的毛片| 男人舔奶头视频| 狂野欧美激情性xxxx在线观看| 最近2019中文字幕mv第一页| 一本一本综合久久| 国产精品蜜桃在线观看| 波多野结衣巨乳人妻| 成年女人看的毛片在线观看| 日韩制服骚丝袜av| 亚洲欧美日韩卡通动漫| 精品人妻熟女av久视频| 成人美女网站在线观看视频| av福利片在线观看| 国产一区二区在线观看日韩| 久久人人爽av亚洲精品天堂 | 国产精品一区www在线观看| 十八禁网站网址无遮挡 | 日韩av不卡免费在线播放| 最近中文字幕高清免费大全6| 国产亚洲av片在线观看秒播厂| 亚洲av日韩在线播放| 视频中文字幕在线观看| 一级毛片aaaaaa免费看小| 久久国产乱子免费精品| av卡一久久| 国产精品久久久久久精品古装| 人人妻人人看人人澡| 久久久午夜欧美精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 狂野欧美白嫩少妇大欣赏| 欧美+日韩+精品| 97超视频在线观看视频| 精品99又大又爽又粗少妇毛片| 两个人的视频大全免费| 国产精品久久久久久av不卡| 国产毛片在线视频| 亚洲av男天堂| 啦啦啦啦在线视频资源| 日韩一本色道免费dvd| 国产成人a∨麻豆精品| 成人欧美大片| 国产伦精品一区二区三区视频9| 亚洲精品日本国产第一区| 日本熟妇午夜| 国产精品爽爽va在线观看网站| 国产爽快片一区二区三区| 国产高潮美女av| 精品一区在线观看国产| av黄色大香蕉| 男女下面进入的视频免费午夜| 午夜爱爱视频在线播放| 国产黄片美女视频| 噜噜噜噜噜久久久久久91| 亚洲精品乱码久久久v下载方式| 男女啪啪激烈高潮av片| 日本欧美国产在线视频| 久久女婷五月综合色啪小说 | 黄色欧美视频在线观看| 99久久精品热视频| 丝袜脚勾引网站| 99久久九九国产精品国产免费| 麻豆国产97在线/欧美| 内地一区二区视频在线| 精品人妻视频免费看| 99久国产av精品国产电影| 亚洲最大成人中文| 免费av毛片视频| 日本一二三区视频观看| 秋霞伦理黄片| 一个人观看的视频www高清免费观看| 久久久色成人| 99久久精品一区二区三区| 在线观看av片永久免费下载| 啦啦啦在线观看免费高清www| 一级二级三级毛片免费看| 久久精品久久久久久久性| 精品少妇黑人巨大在线播放| 性色avwww在线观看| 2021少妇久久久久久久久久久| 在线免费观看不下载黄p国产| a级一级毛片免费在线观看| 汤姆久久久久久久影院中文字幕| 成年人午夜在线观看视频| 国产视频内射| a级一级毛片免费在线观看| 国产伦精品一区二区三区四那| 国产av码专区亚洲av| 欧美另类一区| 免费看av在线观看网站| 久久女婷五月综合色啪小说 | 亚洲av福利一区| 99久久精品国产国产毛片| 国产精品麻豆人妻色哟哟久久| 午夜福利高清视频| 久久久精品94久久精品| 一级av片app| 欧美日韩综合久久久久久| 最新中文字幕久久久久| 亚洲成人一二三区av| 中文天堂在线官网| 日韩强制内射视频| 欧美丝袜亚洲另类| 日韩欧美一区视频在线观看 | 亚洲国产精品999| 国产一区亚洲一区在线观看| av线在线观看网站| 看非洲黑人一级黄片| 麻豆乱淫一区二区| 1000部很黄的大片| 99久久精品国产国产毛片| 女人久久www免费人成看片| 亚洲最大成人av| 亚洲三级黄色毛片| 黑人高潮一二区| 欧美高清性xxxxhd video| 免费高清在线观看视频在线观看| 免费黄色在线免费观看| 永久网站在线| 听说在线观看完整版免费高清| 国产精品久久久久久久久免| 国产高清不卡午夜福利| 欧美精品国产亚洲| 超碰97精品在线观看| av网站免费在线观看视频| 亚洲精品国产av成人精品| 人妻制服诱惑在线中文字幕| 日韩成人av中文字幕在线观看| 人妻 亚洲 视频| 亚洲美女视频黄频| 极品教师在线视频| 亚洲成人精品中文字幕电影| 色视频www国产| 22中文网久久字幕| 久久国内精品自在自线图片| 国产黄片美女视频| 91精品伊人久久大香线蕉| 伦精品一区二区三区| 97在线视频观看| 午夜精品一区二区三区免费看| 视频区图区小说| 色视频www国产| 男女边摸边吃奶| 亚洲av免费高清在线观看| 特级一级黄色大片| 亚洲一级一片aⅴ在线观看| 日韩大片免费观看网站| 午夜日本视频在线| 亚洲av电影在线观看一区二区三区 | 免费看a级黄色片| 婷婷色综合www| .国产精品久久| 日日摸夜夜添夜夜爱| 亚洲丝袜综合中文字幕| 一级a做视频免费观看| 一级爰片在线观看| 日韩制服骚丝袜av| 又爽又黄a免费视频| 国产一区二区三区av在线| 午夜激情福利司机影院| 夜夜看夜夜爽夜夜摸| 丰满少妇做爰视频| 99视频精品全部免费 在线| 一二三四中文在线观看免费高清| 一个人看视频在线观看www免费| 欧美性猛交╳xxx乱大交人| 99久久精品一区二区三区| 欧美精品国产亚洲| 日韩一区二区视频免费看| 国产成年人精品一区二区| 国产久久久一区二区三区| 精品国产乱码久久久久久小说| 亚洲精品456在线播放app| 蜜臀久久99精品久久宅男| 国产中年淑女户外野战色| 精品国产露脸久久av麻豆| 免费不卡的大黄色大毛片视频在线观看| 免费看a级黄色片| 嫩草影院新地址| 日本一本二区三区精品| 97人妻精品一区二区三区麻豆| 干丝袜人妻中文字幕| 国产精品偷伦视频观看了| 人妻系列 视频| 麻豆久久精品国产亚洲av| 久久久色成人| 久久热精品热| 日韩一区二区三区影片| 日韩av不卡免费在线播放| 国产精品99久久99久久久不卡 | 国产亚洲精品久久久com| 日韩,欧美,国产一区二区三区| 黄色怎么调成土黄色| 欧美+日韩+精品| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 在线观看美女被高潮喷水网站| 国产毛片在线视频| 啦啦啦中文免费视频观看日本| 99re6热这里在线精品视频| 少妇裸体淫交视频免费看高清| av在线亚洲专区| 免费大片18禁| 少妇高潮的动态图| 王馨瑶露胸无遮挡在线观看| 高清午夜精品一区二区三区| 亚洲精品日本国产第一区| 国产免费一级a男人的天堂| 在现免费观看毛片| 日韩人妻高清精品专区| 国产精品国产三级国产专区5o| av在线观看视频网站免费| 久久精品久久久久久久性| 国产午夜精品久久久久久一区二区三区| 免费黄色在线免费观看| 91aial.com中文字幕在线观看| 又爽又黄无遮挡网站| 日韩欧美精品v在线| 少妇人妻精品综合一区二区| 国产午夜福利久久久久久| 欧美zozozo另类| kizo精华| 最近最新中文字幕免费大全7| 成人特级av手机在线观看| 久久影院123| 久热久热在线精品观看| 最近中文字幕2019免费版| 天美传媒精品一区二区| 午夜福利在线在线| 日韩人妻高清精品专区| 久久韩国三级中文字幕| 黄色配什么色好看| 亚洲精品乱久久久久久| 新久久久久国产一级毛片| 九九久久精品国产亚洲av麻豆| 国产69精品久久久久777片| 寂寞人妻少妇视频99o| 99精国产麻豆久久婷婷| 国产成年人精品一区二区| 亚洲精品国产成人久久av| 大香蕉97超碰在线| 人体艺术视频欧美日本| av国产精品久久久久影院| 在线播放无遮挡| 精品久久国产蜜桃| 极品少妇高潮喷水抽搐| 国产真实伦视频高清在线观看| 亚洲最大成人手机在线| 成人午夜精彩视频在线观看| 国产成人精品婷婷| 97在线人人人人妻| 亚洲成人精品中文字幕电影| 国产精品99久久99久久久不卡 | 网址你懂的国产日韩在线| 亚洲av电影在线观看一区二区三区 | 欧美三级亚洲精品| 男男h啪啪无遮挡| 欧美+日韩+精品| 狂野欧美激情性xxxx在线观看| 国产淫语在线视频| 美女国产视频在线观看| 一级二级三级毛片免费看| 91久久精品国产一区二区成人| 久久这里有精品视频免费| 在线天堂最新版资源| 亚洲精品久久午夜乱码| 青春草视频在线免费观看| 日日撸夜夜添| 女的被弄到高潮叫床怎么办| 寂寞人妻少妇视频99o| 国产成人免费无遮挡视频| 七月丁香在线播放| 中文字幕免费在线视频6| 国产高清不卡午夜福利| 精品少妇久久久久久888优播| 日韩视频在线欧美| kizo精华| 美女被艹到高潮喷水动态| 欧美老熟妇乱子伦牲交| 少妇高潮的动态图| 中文欧美无线码| videossex国产| 秋霞伦理黄片| 日韩欧美精品v在线| 美女xxoo啪啪120秒动态图| 国产一级毛片在线| 欧美日韩一区二区视频在线观看视频在线 | 久久99热6这里只有精品| 久久精品国产亚洲av天美| 天堂中文最新版在线下载 | 99久久精品国产国产毛片| 久久精品熟女亚洲av麻豆精品| 一区二区三区精品91| 性插视频无遮挡在线免费观看| 亚洲在线观看片| 国产黄频视频在线观看| 美女内射精品一级片tv| 国产探花在线观看一区二区| 神马国产精品三级电影在线观看| 一二三四中文在线观看免费高清| 国产日韩欧美在线精品| 国产精品偷伦视频观看了| 国产亚洲av嫩草精品影院| 80岁老熟妇乱子伦牲交| 国产精品秋霞免费鲁丝片| 51国产日韩欧美| 亚洲色图av天堂| 亚洲精品国产成人久久av| 老女人水多毛片| 少妇人妻一区二区三区视频| 成人无遮挡网站| 日韩成人伦理影院| 少妇人妻久久综合中文| 在线免费观看不下载黄p国产| 亚洲av二区三区四区| 一级毛片电影观看| 免费黄频网站在线观看国产| 国产一区二区三区av在线| 老司机影院成人| 岛国毛片在线播放| 午夜激情久久久久久久| 免费黄网站久久成人精品| 亚洲激情五月婷婷啪啪| 97热精品久久久久久| 亚洲欧美日韩东京热| 男女边吃奶边做爰视频| 精品99又大又爽又粗少妇毛片| 久久热精品热| 国产在视频线精品| 欧美+日韩+精品| 丰满人妻一区二区三区视频av| 欧美丝袜亚洲另类| 欧美三级亚洲精品| 国产高清有码在线观看视频| 日韩欧美一区视频在线观看 | av线在线观看网站| 汤姆久久久久久久影院中文字幕| 婷婷色综合www| 黄色一级大片看看| 国产亚洲午夜精品一区二区久久 | 久久精品久久久久久噜噜老黄| 亚洲精品成人av观看孕妇| 美女脱内裤让男人舔精品视频| 18禁在线播放成人免费| 国产日韩欧美在线精品| 日韩伦理黄色片| 最近2019中文字幕mv第一页| 熟女电影av网| 亚洲欧美日韩卡通动漫| av在线亚洲专区| 久久久国产一区二区| 尤物成人国产欧美一区二区三区| 91aial.com中文字幕在线观看| 少妇 在线观看| 国产有黄有色有爽视频| 日韩成人av中文字幕在线观看| 看黄色毛片网站| 久久人人爽人人爽人人片va| 亚洲天堂国产精品一区在线| 婷婷色综合大香蕉| 少妇熟女欧美另类| 永久免费av网站大全| 欧美日韩国产mv在线观看视频 | 国产伦在线观看视频一区| 国产av不卡久久| 老师上课跳d突然被开到最大视频| 看非洲黑人一级黄片| 建设人人有责人人尽责人人享有的 | 91精品伊人久久大香线蕉| 亚洲国产精品999| 九色成人免费人妻av| 国产日韩欧美在线精品| 亚洲av.av天堂| 亚洲图色成人| 一区二区av电影网| 久久久久久久国产电影| 黄片wwwwww| 成人亚洲精品av一区二区| 亚洲天堂av无毛| 久久精品久久久久久久性| 久久久午夜欧美精品| 国产精品一区www在线观看| 国产av码专区亚洲av| 亚洲国产最新在线播放| 岛国毛片在线播放| 亚洲精品久久午夜乱码| 97超碰精品成人国产| 国产一区二区亚洲精品在线观看| 国产成人a区在线观看| 久久久成人免费电影| 国产在线一区二区三区精| 搡女人真爽免费视频火全软件| 女的被弄到高潮叫床怎么办| 亚洲精品,欧美精品| 美女高潮的动态| 久久精品国产自在天天线| 亚洲人成网站在线播| av免费在线看不卡| 国产色婷婷99| 国国产精品蜜臀av免费| 日韩成人伦理影院| 国产高清有码在线观看视频| 国产白丝娇喘喷水9色精品| 亚洲国产最新在线播放| 国产免费福利视频在线观看| 国产精品99久久99久久久不卡 | 两个人的视频大全免费| 下体分泌物呈黄色| av在线观看视频网站免费| 嘟嘟电影网在线观看| 欧美bdsm另类| 国产 一区 欧美 日韩| 亚洲国产精品专区欧美| 亚洲av成人精品一区久久| 国产成人福利小说| 小蜜桃在线观看免费完整版高清| 99九九线精品视频在线观看视频| 成人高潮视频无遮挡免费网站| 18禁在线无遮挡免费观看视频| 亚洲欧美一区二区三区黑人 | 日本黄大片高清| 亚洲国产色片| 国产黄频视频在线观看| 欧美 日韩 精品 国产| av国产免费在线观看| 国产久久久一区二区三区| 成人国产av品久久久| 亚洲真实伦在线观看| 国产精品蜜桃在线观看| 精华霜和精华液先用哪个| 晚上一个人看的免费电影| 亚洲国产精品专区欧美| 久久久久久久亚洲中文字幕| 麻豆国产97在线/欧美| 中国美白少妇内射xxxbb| 国产av国产精品国产| 大陆偷拍与自拍| 丝袜脚勾引网站| 国产精品福利在线免费观看| 国产爽快片一区二区三区| 爱豆传媒免费全集在线观看| 精品久久久久久久人妻蜜臀av| 综合色丁香网| 亚洲成人av在线免费| 国产精品久久久久久精品电影小说 | 午夜视频国产福利| 日本av手机在线免费观看| 亚洲av二区三区四区| 五月天丁香电影| 国产亚洲5aaaaa淫片|