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

    Smith-Purcell radiation improved by multi-grating structure

    2022-04-12 03:46:12JingShu舒靖PingZhang張平ManLiang梁滿ShengPengYang楊生鵬ShaoMengWang王少萌andYuBinGong宮玉彬
    Chinese Physics B 2022年4期
    關鍵詞:張平

    Jing Shu(舒靖), Ping Zhang(張平), Man Liang(梁滿), Sheng-Peng Yang(楊生鵬),Shao-Meng Wang(王少萌), and Yu-Bin Gong(宮玉彬)

    School of Electronic Science and Engineering,University of Electronic Science and Technology of China,Chengdu 610054,China

    Keywords: Smith-Purcell radiation,multi-grating,resonance mode

    1. Introduction

    Modern nanofabrication techniques have promoted the development of photonics greatly and have brought new opportunities for fundamental effects study involving the interactions among free electrons, matter and light.[1,2]The interaction between an electron and photonic crystal structure gives rise to transient induced charges and currents that result in the emission of Smith-Purcell radiation (SPR).[3,4]The SPR is a wide spectrum electromagnetic radiation,[5,6]which has been obtained from the metal chains,[7]sub-wavelength hole array,[8,9]metamaterials[10-12]and photonic crystal,[13-18]etc.The physical mechanisms of superradiate SPR,[19,20]surface plasmon enhanced SPR[8-10]and compression SPR[8]have been found and analyzed by researchers.

    From the view of image charge,the SPR radiation mechanism can be explained as below. When the electron beam is moving above the grating, the image charges are induced.If the grating groove and tooth change periodically, the distance between the electron and its image charge has a periodical variation like a dipole oscillation,then it emits radiation.

    Recently, the SPR in photonic crystal attracted many interests, and an important reason is that there are still some novel physical phenomena and mechanism unclear or undiscovered in photonic crystal for SPR.Ochiaiet al. and Hoanget al.showed that the SPR could be realized by using the quasi-guided modes in the photonic crystal slab[13]and the planar defect mode of the woodpile photonic crystal.[18]The radiation is due to the coupling of a free electron and the electromagnetic modes of the photonic crystal,or the coupling between the electron and defect mode in the band gap. In addition,Zenget al. reported that more intensive scatter fields can be obtained in the photonic time crystals.[21]

    Compared with the dielectric photonic crystal,[22]the metal photonic crystal will support stronger dipole oscillations and have stronger resonances. In this paper,we show that the improved SPR happens in a metal photonic crystal consisted of multi-grating structure,as shown in Fig.1.This structure leads to two advantages:(1)The radiation intensity is improved due to the increase of induced image charges,and(2)the spectrum distribution is improved due to the resonance characteristics.By optimization, the coherent directional THz radiation can be realized.

    Fig.1. The schematic diagram of the multi-grating structure. When an electron bunch passed through this structure,the SPR emitted.The SPR is enhanced and has an obvious frequency peak.

    2. Method

    Precise understanding of the radiation excited by free moving electrons requires, in general, recourse to numerical simulations. The electrons behavior as well as the emitted radiation is well investigated by using the particle-in-cell(PIC)numerical method. Its processes involve interaction between the electrons and electromagnetic fields. Beginning from a specified initial state,the method simulates a physical process as it evolves in time. This approach provides self-consistency and is applicable to broad classes of free electron radiation problems.

    In Fig.2,our metallic structure is set in the center of the simulation box. The radiation is calculated in three dimensions. The boundary of the simulation box is open (the electromagnetic fields are absorbed without reflection).It is equivalent to placing our structure in an infinite space and observing the radiation field in free space. We use a Gaussian shape line electron bunch with beam velocity 0.65c(cis the light velocity in vacuum),and place it at the left side of the box.[23]The length and charge quantity of the electron bunch is 160μm and 1×10-17C,respectively. In order to facilitate the analysis of the radiation field, the electromagnetic fields as functions of space or time variables are recorded.

    Fig.2. The PIC simulation model to calculate the SPR.

    In order to analyze the phenomenon that the increasement of the radiation intensity due to more dipole oscillations occurring at the same time,we carry out researches on the radiations of one-row grating, double-row grating, and four-row grating structure. In addition, when the number of rows of grating increases, a strong resonance effect will happen in the multigrating structure. We study the resonance frequencies and resonance modes by plane wave excitation method,and verify the phenomenon that the resonance effect regulates the radiation frequency peak.

    3. Results

    3.1. Radiation from the one-row grating structure

    Figure 3 is a schematic diagram of a one-row grating structure.It is composed by periodic metal pillars placed along thezdirection. The side lengthsaandbof the metal pillars are both 100μm,and the perioddin thezdirection is 200μm.The free electrons move above the metal grating along thezdirection,and the SPR will be excited.

    The radiation mechanism of the SPR in a one-row grating structure is shown in Fig.3(c). When an electron moves to the top of the metal pillar,an image charge will be induced in the metal pillar. When the electron moves to the gap between two metal cylinders, an image charge will be induced at infinity.With periodic changes of metal pillar and gap, the free electrons and their corresponding image charges oscillate in the form of dipoles,and then the electromagnetic fields generate.

    Fig. 3. (a) The schematic diagram of a one-row grating, d =200 μm.(b) Parameters of the cross section of the metal pillar, a=100 μm,b=100 μm. (c) Radiation mechanism of the SPR in a one-row grating. When an electron is moving above the grating, the image charge is induced. The grating tooth and gap change periodically,the distance between the electron and its image charge has a periodical variation like a dipole oscillation,then it emits radiation.

    The time domain waveform of the radiation field at the position of 2000μm distance away from the grating surface is shown in Fig.4(b). From the time domain of the probe signal,the higher frequency part arrives at the probe first,and then the lower frequency part arrives. The contour map of the electric fieldEzat 20 ps is given in Fig.4(a),where red represents the positive electric field strength,and blue represents the negative electric field strength. For the convenience of comparison,the following simulation results will use the same scale.

    Fig.4. (a)The contour map of the electric field Ez at the time of 20 ps.(b) Time domain waveform of the electric field in a one-row grating structure. The inset is the spectrum of the radiation from one-row grating structure.

    3.2. Radiation from the double-row grating structure

    The radiation intensity of Smith-Purcell radiation decreases as the grating size decreases,so the radiation of a onerow grating structure is weak in THz. To achieve the purpose of enhancing radiation,the double-row grating structure is proposed as shown in Fig. 5. The distance between the gratings inydirection is 100μm. The free electrons move between the two rows of gratings at a constant speed along thezdirection to generate SPR.

    Fig.5. The schematic diagram of a double-row grating structure,where the distance between the gratings is 100μm.

    Figure 6 shows the time domain waveform of the radiation field in the double-row grating structure. Compared with the one-row grating structure, the time domain waveform of the double-row grating structure is steeper,and the wave crest connection is smooth.It has a similar radiation frequency band with the one-row grating structure,but the field intensity is enhanced. The contour map of the electric fieldEzin the doublerow grating structure is given in Fig.6(a). It can be obviously noticed that the field of radiation from the double-row grating is stronger than that from the one-row grating structure.

    Fig.6. (a)The contour map of Ez. (b)Time domain waveform diagram of the electric field in a double-row grating structure. The inset is the spectrum of the radiation from double-row grating structure.

    3.3. Radiation from the multi-grating structure

    Since the double-row grating structure can induce more image charges to achieve the purpose of enhancing radiation,we subsequently suppose that the multi-grating structure could obtain stronger radiation than the double-row grating structure.Figure 1 gives the schematic diagram of the multi-grating structure which consists of four gratings.In this structure,four gratings are periodically arranged along theydirection with a period ofl=200μm. Free electrons move at 0.65calong thezdirection.

    Fig.7. (a)The contour map of Ez in a multi-grating structure. (b)The relation between the radiation frequency and the radiation angle. (c)The time domain waveform of the radiation field in a multi-grating structure. (d)The corresponding spectrum of the radiation field.

    The contour map of the electric field of radiation is given in Fig. 7(a). It can be seen that the radiation of the multigrating structure is the greatest among that of the one-row structure and that of the double-row structure grating. Observing the spatial distribution of the SPR in multi-grating structure from Fig. 7(a), the contrast between red and blue is sharp, indicating the radiation in the multi-grating structure has a strong intensity. There is a radiation direction in the contour map. It indicates that the radiation spectrum may has an obvious peak. In addition,there are two anglesθandαon the contour map. The angleθindicates the radiation angle when the radiation field spot just moves away from the grating surface. The angleαreflects the Doppler effect because of the movement of the electrons along the gratings. Based on the relation between the radiation frequency and radiation angle shown in Fig. 7(b), the radiation angleθ=96°corresponds to the frequency 0.9 THz. Figure 7(c) is the time domain of the radiation field,and its corresponding spectrum is given in Fig. 7(d). In Fig. 7(d), an obvious radiation frequency peak can be found,and it corresponds to the dot shown in Fig.7(b).

    4. Discussion

    The time domain waveforms of radiation field under different grating structures are shown in Fig.8. In order to facilitate the comparison, the observed probes are set at the same position for the three cases. It shows that the maximum values of the radiation fields for the one-row grating structure,double-row grating structure and four-row grating structure excited by the same electron bunch are 0.2 V/m,0.3 V/m and 0.45 V/m,respectively. With the increase of the row number,the radiation field intensity gradually increases.

    Fig.8. The radiation fields for different grating structures.

    Figure 9 are the spectrums of SPR for different grating structures. The spectrum has no obvious radiation peak in the one-row grating structure and the double-row grating structure, but the radiation of the double-row grating structure is stronger than that of the one-row grating structure. Compared with the one-row and double-row grating structures,the fourrow grating structure is prone to have a stronger resonance,its radiation spectrum becomes narrow, and a sharp radiation frequency peak appears.[24]

    Fig. 9. The spectrum of the radiation for different grating structures.The black curve,red curve,and blue curve represent the SPR spectrum curves of a one-row grating structure, a double-row grating structure,and a four-row grating structure excited by the same electron bunch,respectively.

    The increasement of the radiation intensity can be explained from the view of the dipole oscillations. When free electrons pass over the periodic metal structure, the image charges are induced, and the increase in the number of rows causes to more induced image charges. The radiation mechanism of the SPR in a double-row grating structure is shown in Fig. 10. Due to an extra row of grating, the more image charges are induced based on the complex image theory. For simplicity, only the two main mirror charges are considered here. It is the physical mechanism of the SPR enhancement in the double-row grating structure.[25,26]

    Fig.10. Radiation mechanism of SPR in a double-row grating structure.

    For the four-row grating structure, more dipole oscillations occur as shown in Fig.11.The increase of image charges results in increase of the number of dipole oscillations,which makes the intensity of SPR from the four-row grating stronger than that from the double-row grating structure.

    In order to explain the generation of peak in the spectrum for four-row grating structure,we firstly give the Brillouin diagram of the radiation and the position of the peak,as shown in Fig. 12(a). The figure corresponds to the situation where electrons with a velocity of 0.65cpass over a structure with a period of 200μm. The electron beam line is expanded into harmonics in the periodical structure.The red line corresponds to the-1stspace harmonic of the electron. The gray shaded area shows the intersection of the electron line and the fast wave area, indicating the area corresponding to the radiation.The yellow point represents the position where the radiation frequency peak appears in the simulation, and it is located in the radiation region (the gray shaded area). In general, the appearance of this peak is closely related to the resonance in the four-row grating structure. It can be seen from Fig.12(b)that the structure has five resonance peaks,and the 0.907 THz resonance peak is located in the SPR radiation region. The resonance peaks below this frequency are all below the SPR radiation threshold and cannot be radiated. From the inset of Fig.12(b),it can be seen that the 0.907 THz resonance mode is symmetrical. In the electron channel,it has a strong longitudinal field and is easily excited by moving electrons. Comparing the curves shown in Fig.12(b),the resonance point ofS11is in good agreement with the peak of radiation spectrum.Through the change of the grating structure parameters, the resonance characteristics will be changed, and then the control of the grating resonant frequency can be realized.

    Fig.11. Radiation mechanism of SPR in a multi-grating structure.

    Only the period of the structure is adjusted, and the PIC simulations andS-parameter simulations are carried out for the structure periods of 160 μm and 240 μm. TheS-parameters are shown as solid black lines in Fig. 13, and the resonant frequency varies with the period. The shaded areas represent the radiation spectrum of the structure under electron excitation, which also changes with the period. By comparing theS-parameter and the radiation spectrum,it can be seen that the corresponding relationship between the radiation frequency peak and the resonance point of the structure does not change.

    Fig.12. (a)The SPR diagram. The red line is the-1st harmonic of the electron beam with the velocity 0.65c. The gray area is for the SPR region. (b)The S11 curve and spectrum of the SPR.The inset figure is the Ez field distribution for the mode at 0.907 THz in the four-row grating structure.

    Fig.13. The S11 curve and spectrum of the SPR for(a)d=160μm and(b)d=240μm in the four-row grating structure.

    We increase the width of the grating teeth(a)to reduce the gap between the grating teeth, so that the resonance effect of the structure can be enhanced. The radiation spectrums corresponding to differentaare shown in Fig.14. When the parameteraincreases, the intensity of the radiation frequency peak is greatly enhanced, while the frequency of the peak position is almost unchanged. The intensity of the radiation peak has a linear relationship witha,as shown on the left in the figure.Therefore,by enhancing the resonance effect of the structure,it is expected to obtain a stronger radiation frequency peak,and thereby a coherent SPR is obtained.

    Fig.14. The spectrums of the radiation in a multi-grating structure for a=80μm,100μm,120μm and 140μm. The maximum values of the spectrums under different a are shown by the dots on the left face.

    While the structural parameters remain unchanged, it means that the resonant frequency of the resonant mode does not change. We changed the frequency band range of SPR by changing the electron beam energy. When the electron velocity varies in the range of 0.55c-0.85c, although the range of SPR changes greatly, the resonant frequency of 0.907 THz is always located in the SPR region. Therefore,a radiation peak about 0.907 THz can be seen in the radiation spectrums of different electron velocities in Fig.15.

    Fig.15. The spectrum of the SPR for 0.55c,0.65c,0.75c and 0.85c in the four-row grating structure.

    5. Conclusion

    In this paper, the physical phenomenon and mechanism of SPR in multi-grating are studied. The intensity enhancement of the radiation is explained by the occurrence of more dipole oscillations. The number of dipole oscillations of onerow, double-row, and multi-row grating structures gradually increases. And the radiation intensity obtained by the simulation is obviously gradually increased. The frequency tunability is considered to be caused by the resonance characteristics,and the multi-grating structure has a distinct radiation peak at the resonance point in the SPR radiation frequency domain. It will tailor the spectrum of the dipole oscillation energy.Therefore,the multi-grating structure not only offers more energy by increasing the image charges,but also tailors the SPR through the resonance characteristics. This mechanism is proposed to develop a method to improve the SPR,and it can be used for developing THz radiation source combined with micro-nano fabrication technology.

    Acknowledgement

    Project supported by the National Natural Science Foundation of China (Grant Nos. 61921002, 62171098, and 92163204).

    猜你喜歡
    張平
    Topological properties of Sb(111)surface: A first-principles study
    嘰嘰喳喳的小喜鵲
    Magnetic phase diagram of single-layer CrBr3?
    High-pressure elastic anisotropy and superconductivity of hafnium:A first-principles calculation*
    這是你爺倆
    金秋(2020年16期)2020-12-09 01:41:50
    高中生英語學習焦慮現(xiàn)狀研究
    Enhanced spin-dependent thermopower in a double-quantum-dot sandwiched between two-dimensional electron gases?
    千年煙云大西安
    金秋(2019年2期)2019-04-16 09:14:42
    張平書法作品選
    Mechanical Behavior of Bistable Bump Surface for Morphing Inlet
    免费观看人在逋| 久久伊人香网站| 亚洲精品av麻豆狂野| 亚洲五月色婷婷综合| 无遮挡黄片免费观看| 美女高潮到喷水免费观看| 丝袜人妻中文字幕| 黄色丝袜av网址大全| 精品久久久精品久久久| 亚洲精品在线观看二区| 久久久国产精品麻豆| 国产欧美日韩一区二区精品| 久久青草综合色| 日韩中文字幕欧美一区二区| 天堂动漫精品| 男人操女人黄网站| 窝窝影院91人妻| 亚洲精品国产区一区二| 精品一品国产午夜福利视频| 18美女黄网站色大片免费观看| 99久久综合精品五月天人人| 日韩精品中文字幕看吧| 男女做爰动态图高潮gif福利片 | 在线av久久热| 韩国av一区二区三区四区| 国产欧美日韩一区二区三区在线| 亚洲无线在线观看| 女警被强在线播放| 精品久久久久久久人妻蜜臀av | 亚洲精华国产精华精| 色综合欧美亚洲国产小说| 欧美日本视频| 在线视频色国产色| 黑丝袜美女国产一区| 18禁黄网站禁片午夜丰满| 啦啦啦韩国在线观看视频| 久久久久久免费高清国产稀缺| 国产蜜桃级精品一区二区三区| 大型av网站在线播放| 精品欧美一区二区三区在线| 国产精品国产高清国产av| 久久久国产欧美日韩av| 久久精品国产亚洲av高清一级| 成年女人毛片免费观看观看9| 亚洲一区二区三区不卡视频| 99re在线观看精品视频| 国产精品久久久久久人妻精品电影| 成年女人毛片免费观看观看9| 亚洲专区字幕在线| 亚洲片人在线观看| 嫁个100分男人电影在线观看| 亚洲欧洲精品一区二区精品久久久| 一级a爱视频在线免费观看| 亚洲精品久久成人aⅴ小说| 国产成人精品无人区| 成人三级做爰电影| 亚洲欧美激情在线| 久久久水蜜桃国产精品网| 后天国语完整版免费观看| 精品一区二区三区av网在线观看| 天天躁夜夜躁狠狠躁躁| 香蕉丝袜av| 成人国语在线视频| 在线十欧美十亚洲十日本专区| 在线视频色国产色| 视频在线观看一区二区三区| 夜夜躁狠狠躁天天躁| 欧美日韩黄片免| АⅤ资源中文在线天堂| 欧美色欧美亚洲另类二区 | 一级,二级,三级黄色视频| 美女扒开内裤让男人捅视频| 亚洲国产高清在线一区二区三 | 成在线人永久免费视频| 老司机福利观看| 美女国产高潮福利片在线看| 他把我摸到了高潮在线观看| 99在线视频只有这里精品首页| 99精品在免费线老司机午夜| 午夜福利18| 久久中文字幕一级| 国产激情久久老熟女| 国产精品98久久久久久宅男小说| АⅤ资源中文在线天堂| 国产区一区二久久| 亚洲av五月六月丁香网| 欧美激情 高清一区二区三区| 欧美成狂野欧美在线观看| 啦啦啦观看免费观看视频高清 | 久久久久国产一级毛片高清牌| 日韩国内少妇激情av| 在线观看一区二区三区| 午夜精品在线福利| 91九色精品人成在线观看| 亚洲美女黄片视频| 色精品久久人妻99蜜桃| 级片在线观看| 国产99白浆流出| 久久久水蜜桃国产精品网| 在线观看免费午夜福利视频| 中文字幕精品免费在线观看视频| 他把我摸到了高潮在线观看| 人人妻人人澡人人看| 国产熟女xx| 在线视频色国产色| 高清在线国产一区| 国产精品1区2区在线观看.| 日韩欧美一区二区三区在线观看| 99精品久久久久人妻精品| 看片在线看免费视频| 嫩草影视91久久| 黄色视频,在线免费观看| 老司机午夜福利在线观看视频| 色老头精品视频在线观看| 日韩欧美一区视频在线观看| 麻豆国产av国片精品| 久久国产亚洲av麻豆专区| 亚洲av成人不卡在线观看播放网| 天堂√8在线中文| 久久久久九九精品影院| 欧美日韩亚洲国产一区二区在线观看| 美国免费a级毛片| 人人妻,人人澡人人爽秒播| 老汉色av国产亚洲站长工具| 欧美成人一区二区免费高清观看 | 一级a爱视频在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 99国产极品粉嫩在线观看| 在线观看免费视频网站a站| 欧美日本中文国产一区发布| 嫩草影院精品99| 日本免费a在线| 色播亚洲综合网| 人人妻人人澡人人看| 国产精品 欧美亚洲| 欧美性长视频在线观看| 国产精品久久久av美女十八| 大陆偷拍与自拍| 好男人在线观看高清免费视频 | 国产成人精品久久二区二区免费| 天天躁夜夜躁狠狠躁躁| АⅤ资源中文在线天堂| 国产成人影院久久av| 中出人妻视频一区二区| 韩国av一区二区三区四区| 一进一出抽搐gif免费好疼| 国产欧美日韩一区二区三| 搡老妇女老女人老熟妇| 老司机午夜福利在线观看视频| 黄色成人免费大全| 亚洲熟妇熟女久久| 亚洲成a人片在线一区二区| 757午夜福利合集在线观看| 国产精品免费一区二区三区在线| 桃色一区二区三区在线观看| 国产精品99久久99久久久不卡| 岛国视频午夜一区免费看| 国产国语露脸激情在线看| 欧美绝顶高潮抽搐喷水| 91精品三级在线观看| 桃色一区二区三区在线观看| 国产一级毛片七仙女欲春2 | 欧美成人免费av一区二区三区| 国产精品98久久久久久宅男小说| 久久久精品欧美日韩精品| 欧美另类亚洲清纯唯美| 日韩欧美国产一区二区入口| 一区二区三区精品91| 看黄色毛片网站| 美女午夜性视频免费| 狠狠狠狠99中文字幕| 久久精品国产99精品国产亚洲性色 | 黄色丝袜av网址大全| 一级a爱视频在线免费观看| av在线播放免费不卡| 日韩视频一区二区在线观看| 午夜精品在线福利| 欧洲精品卡2卡3卡4卡5卡区| 国产精品亚洲美女久久久| 精品电影一区二区在线| 18禁国产床啪视频网站| 日本a在线网址| 久久这里只有精品19| 亚洲在线自拍视频| 黄色成人免费大全| 在线观看一区二区三区| 久久人妻熟女aⅴ| 精品国产一区二区久久| 九色国产91popny在线| 精品少妇一区二区三区视频日本电影| 国产不卡一卡二| 精品免费久久久久久久清纯| 亚洲成国产人片在线观看| 麻豆一二三区av精品| 亚洲国产精品sss在线观看| 中国美女看黄片| 色综合婷婷激情| 一区二区日韩欧美中文字幕| 欧美色欧美亚洲另类二区 | 中文字幕人妻熟女乱码| 亚洲av五月六月丁香网| 99国产综合亚洲精品| 亚洲精品国产色婷婷电影| 伊人久久大香线蕉亚洲五| 大型av网站在线播放| 精品国内亚洲2022精品成人| 精品人妻在线不人妻| 99精品久久久久人妻精品| 岛国在线观看网站| 美女高潮到喷水免费观看| 亚洲国产中文字幕在线视频| 午夜免费成人在线视频| 如日韩欧美国产精品一区二区三区| 好看av亚洲va欧美ⅴa在| 欧美午夜高清在线| 亚洲欧美精品综合久久99| 亚洲精品中文字幕在线视频| 久久精品91无色码中文字幕| 一本大道久久a久久精品| 岛国视频午夜一区免费看| 欧美黄色淫秽网站| 国产精品综合久久久久久久免费 | 搡老熟女国产l中国老女人| 一进一出好大好爽视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产99白浆流出| 亚洲一码二码三码区别大吗| 咕卡用的链子| 久久久久久久午夜电影| 国产在线精品亚洲第一网站| 免费少妇av软件| 男人操女人黄网站| 一边摸一边做爽爽视频免费| 国产成年人精品一区二区| 可以在线观看的亚洲视频| 黑人操中国人逼视频| 婷婷六月久久综合丁香| 亚洲成国产人片在线观看| 亚洲第一欧美日韩一区二区三区| 嫩草影视91久久| 视频在线观看一区二区三区| 又黄又粗又硬又大视频| 午夜福利免费观看在线| 中文字幕久久专区| 亚洲视频免费观看视频| 人妻久久中文字幕网| 午夜视频精品福利| 给我免费播放毛片高清在线观看| 曰老女人黄片| 欧美激情高清一区二区三区| 男男h啪啪无遮挡| 日韩欧美一区二区三区在线观看| 亚洲人成伊人成综合网2020| 少妇熟女aⅴ在线视频| 亚洲av美国av| www国产在线视频色| 成人永久免费在线观看视频| 色综合站精品国产| 制服人妻中文乱码| 男女下面进入的视频免费午夜 | 日韩欧美在线二视频| 97碰自拍视频| 精品熟女少妇八av免费久了| 精品国产亚洲在线| 1024视频免费在线观看| 色老头精品视频在线观看| 欧美老熟妇乱子伦牲交| 一级作爱视频免费观看| 一级毛片女人18水好多| 夜夜看夜夜爽夜夜摸| 亚洲第一青青草原| 女警被强在线播放| 国产乱人伦免费视频| 国产精品二区激情视频| 国产区一区二久久| 成人特级黄色片久久久久久久| 亚洲无线在线观看| www国产在线视频色| 精品少妇一区二区三区视频日本电影| 国产一区二区激情短视频| 日本vs欧美在线观看视频| 亚洲五月色婷婷综合| 久久久精品欧美日韩精品| 日韩欧美在线二视频| 91字幕亚洲| a级毛片在线看网站| 国产1区2区3区精品| av在线播放免费不卡| 国产高清有码在线观看视频 | av欧美777| 99在线视频只有这里精品首页| 国产99久久九九免费精品| 每晚都被弄得嗷嗷叫到高潮| 叶爱在线成人免费视频播放| 精品国产一区二区久久| 欧美乱色亚洲激情| 啦啦啦 在线观看视频| 国产人伦9x9x在线观看| 免费av毛片视频| 久久久久国产一级毛片高清牌| 乱人伦中国视频| 国产av又大| 亚洲精品一卡2卡三卡4卡5卡| 午夜亚洲福利在线播放| 精品久久久精品久久久| 久久中文字幕人妻熟女| 日韩大尺度精品在线看网址 | 欧美色欧美亚洲另类二区 | 国内毛片毛片毛片毛片毛片| 久久午夜综合久久蜜桃| 日韩欧美一区二区三区在线观看| 国内久久婷婷六月综合欲色啪| 欧美av亚洲av综合av国产av| 可以免费在线观看a视频的电影网站| www.自偷自拍.com| 精品国产一区二区久久| 亚洲av成人不卡在线观看播放网| 久久久精品国产亚洲av高清涩受| 日韩欧美国产一区二区入口| 日韩一卡2卡3卡4卡2021年| 每晚都被弄得嗷嗷叫到高潮| 中文字幕人妻熟女乱码| 日韩免费av在线播放| 久久午夜亚洲精品久久| 国产一区二区三区视频了| 一边摸一边抽搐一进一小说| 少妇 在线观看| 高清黄色对白视频在线免费看| 国产伦一二天堂av在线观看| 久久人妻福利社区极品人妻图片| 中文字幕精品免费在线观看视频| 51午夜福利影视在线观看| 亚洲一区二区三区不卡视频| 免费高清视频大片| 大码成人一级视频| 国产精品一区二区精品视频观看| 精品乱码久久久久久99久播| 亚洲人成伊人成综合网2020| 亚洲专区字幕在线| 亚洲第一青青草原| 99久久精品国产亚洲精品| 在线观看66精品国产| 久久久国产精品麻豆| 欧美成人免费av一区二区三区| 亚洲熟女毛片儿| 久久香蕉国产精品| 欧美成人免费av一区二区三区| 亚洲自拍偷在线| x7x7x7水蜜桃| 最好的美女福利视频网| 日本撒尿小便嘘嘘汇集6| 日韩免费av在线播放| 午夜精品久久久久久毛片777| 亚洲三区欧美一区| 欧美乱妇无乱码| 日本黄色视频三级网站网址| 91国产中文字幕| 中文亚洲av片在线观看爽| 这个男人来自地球电影免费观看| 亚洲国产高清在线一区二区三 | 精品国产一区二区久久| 亚洲精品一区av在线观看| 精品无人区乱码1区二区| 自线自在国产av| 国产亚洲av嫩草精品影院| 亚洲国产精品合色在线| 长腿黑丝高跟| 多毛熟女@视频| 国产精品久久久久久亚洲av鲁大| 亚洲aⅴ乱码一区二区在线播放 | 色精品久久人妻99蜜桃| 啦啦啦观看免费观看视频高清 | av天堂久久9| 国产精品久久久人人做人人爽| 日韩欧美在线二视频| 嫁个100分男人电影在线观看| 可以免费在线观看a视频的电影网站| 国产亚洲精品第一综合不卡| 操美女的视频在线观看| 丁香欧美五月| 精品久久久久久成人av| 免费av毛片视频| 中文字幕最新亚洲高清| 亚洲成av片中文字幕在线观看| 亚洲欧美日韩高清在线视频| 国产黄a三级三级三级人| 99精品在免费线老司机午夜| 久久久久久国产a免费观看| 窝窝影院91人妻| 亚洲av日韩精品久久久久久密| 在线国产一区二区在线| 少妇裸体淫交视频免费看高清 | 成年女人毛片免费观看观看9| 国产精品久久久av美女十八| 午夜成年电影在线免费观看| 美女高潮到喷水免费观看| 日日爽夜夜爽网站| 国产精品精品国产色婷婷| 欧美日本视频| 男女下面插进去视频免费观看| 欧美不卡视频在线免费观看 | 国产成人精品在线电影| 欧美日本视频| 欧美日韩乱码在线| 国产精品一区二区精品视频观看| 国产亚洲av嫩草精品影院| 久久婷婷成人综合色麻豆| 免费在线观看黄色视频的| 欧美+亚洲+日韩+国产| 美女 人体艺术 gogo| 人人妻,人人澡人人爽秒播| 亚洲一区高清亚洲精品| 欧美黑人精品巨大| 两个人看的免费小视频| 日韩三级视频一区二区三区| 国产熟女xx| 欧美成人性av电影在线观看| 国产视频一区二区在线看| 给我免费播放毛片高清在线观看| 成人国产一区最新在线观看| 在线av久久热| 久久婷婷人人爽人人干人人爱 | 91成年电影在线观看| 韩国av一区二区三区四区| 日本a在线网址| av免费在线观看网站| 日韩大尺度精品在线看网址 | 久久精品亚洲熟妇少妇任你| 久久久久久久久久久久大奶| 国产一级毛片七仙女欲春2 | 国产精品亚洲美女久久久| 高清黄色对白视频在线免费看| 一级a爱视频在线免费观看| 禁无遮挡网站| 此物有八面人人有两片| 首页视频小说图片口味搜索| 午夜a级毛片| 午夜久久久久精精品| 91精品三级在线观看| 精品高清国产在线一区| 国产高清有码在线观看视频 | 级片在线观看| 天堂动漫精品| 丁香六月欧美| 嫩草影院精品99| 国产成人精品久久二区二区免费| 亚洲,欧美精品.| 97超级碰碰碰精品色视频在线观看| 国产精品免费视频内射| 久久久久久国产a免费观看| 免费久久久久久久精品成人欧美视频| 久久婷婷成人综合色麻豆| 亚洲狠狠婷婷综合久久图片| 欧美日韩亚洲国产一区二区在线观看| 色哟哟哟哟哟哟| 男女下面插进去视频免费观看| 丝袜美足系列| 在线免费观看的www视频| 性欧美人与动物交配| 麻豆久久精品国产亚洲av| 午夜亚洲福利在线播放| 丁香六月欧美| 在线免费观看的www视频| 首页视频小说图片口味搜索| 999久久久国产精品视频| 又大又爽又粗| 欧美精品啪啪一区二区三区| 国产一区二区三区在线臀色熟女| 久久国产精品人妻蜜桃| 99re在线观看精品视频| 国产麻豆成人av免费视频| 亚洲国产高清在线一区二区三 | 亚洲人成网站在线播放欧美日韩| 亚洲av成人av| 日韩欧美国产在线观看| 国产在线观看jvid| 精品一区二区三区四区五区乱码| 女同久久另类99精品国产91| 久久香蕉国产精品| 操美女的视频在线观看| 色哟哟哟哟哟哟| 母亲3免费完整高清在线观看| 热re99久久国产66热| 亚洲视频免费观看视频| 精品午夜福利视频在线观看一区| 不卡一级毛片| 国产精品 国内视频| 日韩国内少妇激情av| 国产精品爽爽va在线观看网站 | 久久精品国产亚洲av高清一级| 看片在线看免费视频| 亚洲av美国av| 免费久久久久久久精品成人欧美视频| 亚洲自拍偷在线| 成人亚洲精品av一区二区| 91麻豆精品激情在线观看国产| 午夜久久久久精精品| 日本精品一区二区三区蜜桃| 久久精品国产综合久久久| 久久中文字幕一级| 国产精品一区二区精品视频观看| 欧美成狂野欧美在线观看| 99久久久亚洲精品蜜臀av| 亚洲人成伊人成综合网2020| 搡老岳熟女国产| 成人18禁在线播放| 国产精品一区二区三区四区久久 | 国产精品99久久99久久久不卡| 久久久国产成人免费| 国产成人系列免费观看| 人人澡人人妻人| 欧美日韩亚洲综合一区二区三区_| 麻豆成人av在线观看| 国产欧美日韩综合在线一区二区| 国产精品香港三级国产av潘金莲| 脱女人内裤的视频| 日本a在线网址| 久久精品91蜜桃| 国产97色在线日韩免费| 日韩欧美三级三区| 好看av亚洲va欧美ⅴa在| 丝袜人妻中文字幕| 免费看a级黄色片| 国产精品一区二区在线不卡| 国产av精品麻豆| 一级黄色大片毛片| 老司机午夜十八禁免费视频| 国产免费av片在线观看野外av| 中文字幕高清在线视频| 1024视频免费在线观看| av天堂在线播放| 国产片内射在线| 亚洲狠狠婷婷综合久久图片| 国产视频一区二区在线看| 老司机午夜十八禁免费视频| 国产私拍福利视频在线观看| 女人爽到高潮嗷嗷叫在线视频| a在线观看视频网站| 人人妻人人澡欧美一区二区 | 亚洲人成77777在线视频| 亚洲专区国产一区二区| aaaaa片日本免费| 久久亚洲精品不卡| 午夜日韩欧美国产| 亚洲人成电影观看| 波多野结衣av一区二区av| 国产精品一区二区免费欧美| 十八禁网站免费在线| 一区二区三区国产精品乱码| 大陆偷拍与自拍| 757午夜福利合集在线观看| 夜夜爽天天搞| 欧美成人午夜精品| tocl精华| 男人舔女人的私密视频| 桃色一区二区三区在线观看| 亚洲国产毛片av蜜桃av| 免费不卡黄色视频| 国产成人系列免费观看| 国产人伦9x9x在线观看| 精品国产乱码久久久久久男人| 妹子高潮喷水视频| 可以免费在线观看a视频的电影网站| 一二三四社区在线视频社区8| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩精品免费视频一区二区三区| 黄色毛片三级朝国网站| 1024视频免费在线观看| 91成人精品电影| 国产成人啪精品午夜网站| 免费在线观看完整版高清| 窝窝影院91人妻| 大陆偷拍与自拍| 又黄又爽又免费观看的视频| 麻豆av在线久日| 在线十欧美十亚洲十日本专区| 少妇裸体淫交视频免费看高清 | 国产欧美日韩综合在线一区二区| 啦啦啦免费观看视频1| 国产精品野战在线观看| 日本 欧美在线| 最好的美女福利视频网| 怎么达到女性高潮| 日本 欧美在线| 国产三级在线视频| 国产成人av激情在线播放| 日韩av在线大香蕉| 久久人人97超碰香蕉20202| 欧美大码av| 久久精品国产亚洲av高清一级| 日韩免费av在线播放| 级片在线观看| 亚洲精品中文字幕一二三四区| 精品电影一区二区在线| 亚洲欧洲精品一区二区精品久久久| 免费观看人在逋| 欧美一区二区精品小视频在线| 国产av一区在线观看免费| 久久人妻福利社区极品人妻图片| 亚洲国产欧美日韩在线播放| 日本在线视频免费播放| 免费观看人在逋| 日本免费一区二区三区高清不卡 | 久久久水蜜桃国产精品网| 美女国产高潮福利片在线看| aaaaa片日本免费| 长腿黑丝高跟| 一级黄色大片毛片| 中出人妻视频一区二区| 亚洲无线在线观看| 国产午夜福利久久久久久| 国产成人啪精品午夜网站| 免费看美女性在线毛片视频| 久久久久久人人人人人| 久久影院123|