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

    Topological resonators based on hexagonal-star valley photonic crystals

    2023-12-02 09:22:36XinWan萬(wàn)鑫ChenyangPeng彭晨陽(yáng)GangLi李港JunhaoYang楊俊豪andXinyuanQi齊新元
    Chinese Physics B 2023年11期
    關(guān)鍵詞:楊俊新元

    Xin Wan(萬(wàn)鑫), Chenyang Peng(彭晨陽(yáng)), Gang Li(李港), Junhao Yang(楊俊豪), and Xinyuan Qi(齊新元)

    School of Physis,Northwest University,Xi’an 710127,China

    Keywords: valley photonic crystals,triangular resonant cavity,topological edge state,unidirectional transmission

    1.Introduction

    The photonic quantum valley Hall effect, an analogy to the quantum valley Hall effect in condensed physics, has attracted intensive attention due to its unique properties in light manipulations.[1]Realizing the photonic quantum valley Hall effect requires the valley degrees of freedom.Generally,electromagnetic duality in valley photonic crystals or all-dielectric valley photonic crystals is used to introduce valley degrees of freedom in order to break the band degeneracy near the valleys (KandK′) by breaking spatial inversion symmetry.[2]These enable the realization of optical pseudospin-dependent,path-entangled transport analogous to the electronic valley Hall effect.Nowadays, valley photonic crystals with different structures based on the all-dielectric materials, e.g.,sublattice-broken honeycomb lattices,[3–5]triangular lattices with different-shaped scatterers,[6–8]stretched or compressed kagome lattices,[9,10]and waveguide arrays with alternating refractive indices,[11]have been proposed to breakC6symmetry while preservingC3symmetry.

    The materials for realizing topological boundary states in optical systems are also quite diverse, ranging from anisotropic magneto-electric medium,[12,13]to dualanisotropic materials,[14,15]hyperbolic metamaterials,[16]and equivalent dual-anisotropic materials constructed with metal–air combinations.[17]Furthermore,isotropic all-dielectric photonic crystals have become widely used in recent years due to their simple fabrication and straightforward design of topological states based on the all-dielectric materials.Besides,such valley photonic crystals do not require external magnetic fields,or specific magneto-electric or dual-anisotropic materials,making them highly versatile for various applications.

    The physical nature of valley topological edge states is mainly related to the following effects, such as the quantum Hall effect,[18,19]Floquet topological insulators,[20,21]quantum spin Hall effect,[18,21,22]or quantum valley Hall effect.[22,23]An external magnetic field was initially employed to break the time-reversal symmetry,[24,25]and then spatial inversion symmetry breaking was proposed to obtain the effects above.[26]The emergence of topological photonics provides possibilities for better control and manipulation of light.The topological edge states generated in such systems exhibit excellent properties,[27–30]including defect immunity,robust transmission around obstacles, and scattering suppression.Therefore, they have been extensively studied and hold great potential for applications in photonics integration[31,32]and optical communication fields.[33,34]

    A microcavity[35]refers to a microresonator with a size in the order of wavelength.The study of microcavities has made it possible to integrate devices such as optical filters,[36]optical modulators/demodulators,and optical logic gates,enabling a wide range of applications in fundamental physics research and device applications.[37,38]However, micro-resonant cavities can suffer from increased losses due to fabrication defects and disorder.The combination of micro-resonant cavities with topological edge states in photonic crystals has been explored to suppress the backscattering, e.g., microcavities based on heterostructure,[39]two-dimensional photonic crystal line defect microcavities,[40]photonic crystal nanobeam microcavities,[41]and photonic crystal micro-ring refractive index sensors.[42]Besides, micro-resonators using traditional photonic crystals with defects[43,44]were also proposed.Nevertheless, all-dielectric topological photonic resonator with high-quality factor is still rarely reported.

    In this study, we design aC3symmetric valley photonic crystal based on hexagrams and hexagons.By reducing the lattice symmetry and breaking the spatial inversion symmetry of the system,the degenerate Dirac points are suppressed,and the band gaps appear.Supercell structures with straight,Z-shaped,and ?-type interface waveguides consisting of two valley photonic crystal units,which exhibit opposite but equal magnitudes of Berry curvature distributions, are constructed.When light propagates in the waveguide,the valley edge states emerge at the topological photonic waveguides, exhibiting unidirectional propagation and the suppression of backscattering.Furthermore, valley photonic crystal resonators are designed, and the quality factor is explored.The influence of resonator lengthLon the quality factorQis also investigated.Our findings of photonic crystal resonators may have potential applications in designing photonic crystal lasers, wavelength division multiplexers,and other photonic devices.

    2.Valley photonic crystals

    In this section,a novel all-dielectric two-dimensional valley photonic crystal is constructed.The photonic crystal is constructed by alternating between hexagrams and hexagons.The schematic diagram of a 2D all-dielectric valley photonic crystal is shown in Fig.1(a).The lattice constant is denoted asa=400 nm, and the dielectric constant of the dielectric rods isεd=14.The dielectric constant of the surrounding environment of the dielectric rods isεa=1.The composition of an individual dielectric rod is depicted in Fig.1(b),where the hexagon has a side length ofa1=0.225a, and the difference between the stretched hexagon and the hexagon is given by.The initial primitive cell without modification is depicted in Fig.1(c).The symmetry at the high-symmetry pointsK(K′) in the reduced Brillouin zone of this structure exhibitsC6symmetry.In previous studies,there have various methods to break theC6symmetry of the photonic lattices.This paper mainly discusses breaking spatial inversion symmetry by adjusting the distancedbetween the six-pointed star and the hexagon to achieve a structurally simple topological valley photonic crystal.

    Due to the photonic lattice being composed of two sublattices, it is necessary to define two different types of dielectric rods in order to break the spatial inversion symmetry.The sizes of the modified dielectric rods are as follows:three hexagonal structure dielectric rods and three six-pointed star structure dielectric rods.After reducing the lattice symmetry,the structure of the primitive cell is shown in Fig.1(d)and labeled as valley photonic crystal 1 (VPhC1).If we exchange the positions of the hexagon and the six-pointed star,we can construct valley photonic crystal 2 (VPhC2), where the primitive cell structure is represented by the red structure in Fig.1(e).By keeping the geometric and material parameters of VPhC1 and VPhC2 unchanged, rotating the center of the VPhC1 primitive cell by 180?yields VPhC2.Therefore,compared to VPhC1,it can be considered that VPhC2 simply exchanges the positions of the hexagon and the six-pointed star.Clearly,VPhC1 and VPhC2 have identical photonic band structures.

    Fig.1.Structure diagram of two-dimensional all-dielectric valley photonic crystal.(a) Schematic diagram of the undistorted valley photonic crystal.(b) Schematic diagram of the dielectric rod structure.(c) Unit cell of the undistorted valley photonic crystal.(d)The green structure is a schematic diagram of VPhC1.(e)The orange structure is a schematic diagram of VPhC2.

    In this study, we first consider the TM polarized mode.Scanning along the wave vector pathΓ–M–K(K′)–Γin the reduced Brillouin zone of the honeycomb lattice,the photonic band structure of the valley photonic crystal can be computed using the COMSOL Multiphysics.The computed results are shown in Fig.2.

    From Fig.2(a), it can be seen that the first band (blue curve) and the second band (green curve) exhibit double degeneracy at the valley pointsKandK′,resulting in the closure of the band gap.This double degeneracy atKandK′is an inevitable consequence of the honeycomb lattice exhibitingC6symmetry when the two dielectric rods in the primitive cell are identical.This degeneracy is also known as the Dirac degeneracy.The linear dispersion relation formed around the Dirac points is referred to as the Dirac cone.The normalized frequency of the Dirac points is 0.278c/a,wherecis the speed of light in vacuum.The Dirac cone is protected by both time-reversal symmetry and spatial inversion symmetry.From Fig.2(b), when we consider VPhC1(VPhC2)as the scanned object,the first and second bands are gapped out,ranging from 0.249c/ato 0.283c/a.The Dirac cone is no longer protected by spatial inversion symmetry,leading to the disappearance of the double degeneracy.Besides,we also simulate with TE polarized mode and find that the photonic band gap could not be opened for the first and second bands,which is different from the work in Ref.[45].Therefore,only the TM polarized mode is considered in the following content.

    Figure 2(c) shows the phase distribution of the first and second bands at theKpoint.The first band exhibits a counterclockwise phase distribution(indicated by the orange arrows)atK, with phase singularities primarily located at the center of the six-pointed star dielectric rods.On the other hand, the second band exhibits a clockwise phase distribution(indicated by the blue arrows) atK, with phase singularities located at the center of the hexagon dielectric rods.The phase variations are represented by left-handed circular polarization(LCP)and right-handed circular polarization (RCP).Similarly, Fig.2(d)analyzing the phase distribution at theK′point reveals that the phase distribution atK′is opposite to that atK.At theK′point,the phase of the first band exhibits a clockwise distribution,with phase singularities located at the center of the hexagon dielectric rods, while the phase of the second band exhibits a counterclockwise distribution,with phase singularities located at the center of the six-pointed star dielectric rods.

    The two valley photonic crystals maintain their timereversal symmetry,and the total Chern number of a Brillouin zone is always zero.However,the Berry curvatures are locally nonzero at each valley and are opposite atKandK′valleys due to the quantum valley Hall effect.[46,47]To characterize the topological property of our system, we introduce valley Chern numberCV,[47]

    which is calculated within half of the Brillouin zone around the valley.

    The Chern number for each energy band can be expressed as

    where the Berry curvatureFn(k)=?×〈un(k)|i?k|un(k)〉,andun(k) represents the eigenstate of then-th energy band.By modifying the integration region in Eq.(2)to half of the Brillouin zone for valleysKandK′, respectively, we can obtain the Chern numbersCKandCK′.

    Fig.2.Band structure of the valley photonic crystals.(a) Band structure of the valley photonic crystal with C6 symmetry, illustrated with the first Brillouin zone diagram.(b) Band structure corresponding to the VPhC1 (VPhC2).(c) Phase distribution map at the K point.(d)Phase distribution map at the K′ point.

    Fig.3.Topological edge states of the valley photonic crystals.(a)Schematic diagram of the structure of a valley photonic crystal.(b)Projected band diagram of the dispersion of edge states in the valley photonic crystal.The blue region represents the bulk bands of the valley photonic crystal,while the purple curve represents the dispersion of the edge states.

    The results of Eq.(2)areCK,K′=±1/2 for VPhC1, andCK,K′=?1/2 for VPhC2.Clearly, the valley Chern numberCVare +1 and-1, respectively, for VPhC1 and VPhC2.Therefore, as shown in Fig.3(a), the VPhC1 and VPhC2 are combined to form two different structures at the boundary, constituting a two-dimensional valley photonic crystal.By extracting a long strip-shaped supercell from the twodimensional valley photonic crystal,covering the width of the primitive unit cell of the photonic quasicrystal, we can calculate its projected band structure along thekxdirection, as shown in Fig.3(b).It displays the projected band structure of the valley photonic crystal, where the blue regions represent the bulk bands of the photonic crystal,and the band gaps exist between the bulk bands.A purple-colored dispersion curve of a boundary state appears within the band gap.The band group velocity at theKvalley is positive,indicating the presence of a forward-propagating topological edge state.On the other hand,the band group velocity at theK′valley is negative,indicating the presence of a backward-propagating topological edge state.

    3.Transport properties of topological edge states in valley photonic crystals

    Due to breaking the spatial inversion symmetry,combining two different valley photonic crystals can generate topologically protected boundary states, leading to various interesting physical phenomena.In this section,several topological photonic waveguides are designed based on valley photonic crystals that can exhibit topological valley boundary states.These waveguides are utilized to investigate their characteristics of one-way transmission and strong robustness.

    3.1.Unidirection transmission of TM polarized mode

    One of the significant advantages of the topological edge states in valley photonic crystals is their characteristic of unidirectional transmission.Valley coupling opens up a complete bandgap,where photons within a specific frequency range are prohibited from entering the photonic crystal.However,light waves of TM polarized mode can propagate along the boundaries of the valley photonic crystals, formed by the combination of VPhC1 and VPhC2, in different directions.This phenomenon gives rise to optical topological properties.

    The valley photonic crystal is formed by combining VPhC1 and VPhC2 (see Fig.4(a)).Firstly, a TM polarized mode is incident from the left edge along the interface(zaxis)to obtain the transmittance.The frequency range of the light source is calculated to be within 175 THz and 215 THz based on the dispersion curve shown in Fig.3(a).The transmittance spectrum curve at the output port (right handside port)is shown in Fig.4(b).It can be seen from the diagram that the valley photonic crystal waveguide exhibits high transmittance across the entire wide spectral range of 183 THz to 208 THz.Secondly, we simulate the transmission of light beams when the light sources are incident at the center of the interface(the yellow pentagram) along thezaxis to excite the edge states(see Figs.4(c)and 4(d)).When a right-handed circularly polarized light(RCP)is used as the light source, the TM polarized mode component is excited and propagates along the positive direction ofzaxis; when a left-handed circularly polarized light(RLCP)is used as the light source,the TM polarized mode component is excited and propagates along the negative direction ofzaxis.For both two cases, the energy is primarily concentrated at the interface between VPhC and VPhC2 and the frequencies of light waves are 198.48 THz(0.257c/a).Circularly polarized states of different handedness will propagate in different directions within the combined topological insulator structure.In contrast, as shown in Figs.4(e) and 4(f),such a polarization-selective excitation phenomenon disappears at a higher or lower frequency, e.g., 188.17 THz or 200.35 THz.Therefore,this phenomenon is quite sensitive to light frequency,and only a single-frequency excitation is supported.We attribute such a polarization-selectively unidirectional transmission to the topological valley boundary states.In other words,the designed valley photonic crystal structure achieves unidirectional light transmission for different polarization states.

    Fig.4.The unidirectional transmission and field distribution characteristics of the valley edge states in straight waveguides.(a) Schematic diagram of the straight waveguide in the valley photonic crystal, composed of VPhC1 (green region) and VPhC2 (orange region), with the position of the light source indicated by a yellow pentagram.(b) Transmittance curve.(c) Electric field intensity distribution diagram of the beam excited by right-handed circularly polarized light in the straight waveguide.(d)Electric field intensity distribution diagram of a beam excited by lefthanded circularly polarized light in a straight waveguide.(e)and(f)The distribution diagram of the electric field intensity in the straight waveguide of the beam excited by right-handed circularly polarized light at different frequencies.

    3.2.Robustness

    Another significant advantage of the valley photonic crystal edge states is their robustness against scattering.This section verifies the topological properties of the valley edge states at corners.Straight waveguides,Z-shaped waveguides,and ?-shaped waveguides are designed based on the valley photonic crystals,as shown in Fig.5(a).A light source is placed at the left entrance port of the waveguide,and the frequency range of the light source remains between 175 THz and 215 THz.The Z-shaped waveguide contains two sharp corners,while the ?-shaped waveguide contains four sharp corners, each with an internal angle of 60?.

    Fig.5.Robust transport and field distribution characteristics of energy valley boundary states in different waveguides.Panels(a)–(c)are straight waveguide, Z-shaped and ?-shaped waveguide energy valley photonic crystal waveguide structure diagram and electric field intensity distribution.(d)Transmittance curves of different waveguides.

    The transmission behavior of the boundary states in Zshaped and ?-shaped waveguides is computed using COMSOL Multiphysics.Figures 5(a)–5(c) show the field distribution of the boundary states along the straight waveguide,Z-shaped bent waveguide, and ?-shaped bent waveguide, respectively, at a light wave frequency of 192.36 THz.The following observations can be made from the electric field plots.(i) There is no scattering of the light wave towards the inner regions at each corner.(ii) The light wave’s energy remains relatively unchanged before and after each corner.This analysis indicates that the boundary states are topologically protected and can be transmitted almost losslessly in straight waveguides.In order to better demonstrate the transmission behavior of light waves in the Z-shaped and ?-shaped waveguides,Fig.5(d)presents the transmission spectra of the straight waveguide,Z-shaped waveguide,and ?-shaped waveguide at the output ports.By comparing the transmission spectra of the Z-shaped and ?-shaped waveguides with those of the straight waveguide, it is observed that the transmission spectra in the Z-shaped and ?-shaped waveguides almost overlap entirely with the transmission spectra in the straight waveguide within the bandgap range, and the transmission rates are also high.The topological boundary states transmitted at the sharp corners are also protected by topology,enabling efficient propagation through the bending corners of the waveguide along the predetermined path, demonstrating strong robustness.

    4.Triangular ring resonator based on valley photonic crystal

    The pursuit of optical devices is towards high integration and low energy loss.However,traditional optical devices suffer from losses caused by backscattering,such as waveguides and resonators.The topological boundary states of valley photonic crystals possess immune transmission properties against defects and impurities.Therefore,combining them in the design of optical devices holds excellent research value.

    Photon crystal resonators are a crucial component in photonics integration, and their performance evaluation relies on two essential parameters: resonator volume and quality factor.The resonator volume reflects the degree of miniaturization of the photonic crystal device, with smaller volumes indicating higher integration levels.The quality factor characterizes the resonator’s ability to confine photons,with higher values indicating stronger photon confinement within the resonator.

    In this section, a topological triangular cavity is constructed using the edge states of the valley topological crystal.The achieved topological cavity differs significantly from traditional photonic crystal nanocavities.The valley boundary states are topologically protected and insensitive to the geometric structure,which provides more possibilities for designing and fabricating resonators on the photonic crystal platform.As an essential branch of topological photonics,the valley photonic crystal exhibits stable transmission characteristics.Its relatively simple structure makes it suitable for designing compact on-chip optical devices.

    The resonator is depicted in Fig.6(a),surrounded by the boundaries of the VPhC1 structure outside the cavity and the VPhC2 structure inside the cavity.The topologically protected valley photonic crystal triangular resonator has three corners.To investigate the resonator’s performance, random point sources were placed near the labeled positions on the boundary of the resonator to excite all resonant modes.The frequency range of excitation was from 175 THz to 215 THz.Five models with different resonator lengths of 15a,18a,23a,26a, and 29awere constructed to investigate the influence of resonator length on the quality factor.As shown in Fig.6(b),it was discovered that light can be confined to the topological edge states of the photonic crystal with energy valleys.As depicted in Fig.6(c), at a frequency of 199.21 THz, light can propagate steadily along the boundary state with minimal losses.However, as shown in Fig.6(d), at a frequency of 207.31 THz,the electric field distribution of light becomes non-uniform,resulting in increased losses.

    Fig.6.Topological valley photonic crystal triangular ring resonator.(a) Schematic diagram of a valley photonic crystal triangular resonator.(b) Electric field distribution of a valley photonic crystal resonator.(c) Light field distribution diagram with a frequency of 199.21 THz.(d)Light field distribution diagram with a frequency of 191.67 THz.(e)Q factor of the topological triangular ring resonator.(f)Highest Q factor among different cavity lengths.

    As shown in Fig.6(e), a resonant cavity with a specific cavity length (side length 23a) is studied, and theQfactor exhibits a trend of initially increasing and then decreasing with the increase of the resonant frequency.At a specific frequency,the resonant cavity reaches its highestQfactor,which falls within the frequency range of the valley photonic crystal boundary state.Within this frequency range,topological edge states in the valley photonic crystal provide topological protection, resulting in reduced energy loss in the resonant cavity.As shown in Fig.6(f),comparing the maximumQfactors of the resonant cavities with five different cavity lengths, the highest value is 12531,which occurs at a cavity length of 23acorresponding to a frequency of 203.61 THz.The frequency from pointAto pointEalso gradually increases.The reasons for the decrease in the quality factor of the resonant cavity with increasing cavity length differ.When the side length of the resonant cavity is too short,the light waves at the cavity boundary will couple with each other when the light waves propagate along the resonant cavity,increasing loss.When the side length of the cavity increases,that is,when the size of the cavity increases, theQfactor of the cavity continues to increase before 24a.Continuing to increase the side length will lead to a decrease in theQfactor.The reason for this phenomenon is mainly due to the fact that when light is transmitted in the cavity,the loss it experiences is related to the propagation length it experiences.When the cavity length becomes longer,the loss will lead to a decrease in the energy stored in the cavity.According to the definition,theQfactor equals the energy stored in the cavity divided by the energy lost per unit of time.As theQfactor decreases,the propagation loss of the boundary mode will increase.

    5.Conclusion

    We have proposed a method to achieve a hexagonal-star valley photonic crystal.By alternating the sizes of hexagram scatters,the degeneracy between the high-symmetry pointsKandK′is broken.A valley photonic crystal is realized by combining two sets ofC3symmetric hexagonal star photonic lattices.Supercell structures with straight,Z-shaped,and ?-type interface waveguides consisting of two valley photonic crystal units, which exhibit opposite but equal magnitudes of Berry curvature distributions,are designed.Simulation results show that the topological edge states can be achieved at the interface of the supercell structure, and the robust and stable transmission, particularly at corner bends, further confirms the highly unidirectional light transmission.Based on these findings, a topological triangular resonant cavity is simulated using supercell structures.The maximumQfactor of the resonant cavity reaches 1.25×104for the cavity length of 23aand frequency of 203.61 THz.Our work riches the content of valley topological photonics and may have potential applications in topological lasers,photonic integrated circuits,and optical modulation devices.

    Acknowledgment

    Project supported by the National Natural Science Foundation of China(Grant No.1217040857).

    猜你喜歡
    楊俊新元
    Topological states switching and group velocity control in two-dimensional non-reciprocal Hermitian photonic lattice
    舊歲千重錦,新元百尺竿
    楊俊德:農(nóng)業(yè)豐收的“守護(hù)神”
    DYNAMICS ANALYSIS OF A DELAYED HIV INFECTION MODEL WITH CTL IMMUNE RESPONSE AND ANTIBODY IMMUNE RESPONSE?
    李泊城 隋邦平 楊俊顯 王綠竹 作品
    大眾文藝(2019年3期)2019-01-24 13:39:44
    《紅蜻蜓》教案
    《新元史·高防傳》勘補(bǔ)
    Eあects of Correlation between Network Structure and Dynamics of Oscillators on Synchronization Transition in a Kuramoto Model on Scale-Free Networks?
    王新元與“萬(wàn)家美”聯(lián)手推出“面子”針織時(shí)尚秀
    流行色(2005年5期)2005-04-29 18:26:58
    丝袜人妻中文字幕| 俺也久久电影网| 欧美日韩亚洲综合一区二区三区_| 中文字幕人成人乱码亚洲影| 亚洲真实伦在线观看| 久久精品人妻少妇| 午夜两性在线视频| 亚洲专区国产一区二区| 国产乱人伦免费视频| 成人永久免费在线观看视频| 看片在线看免费视频| 国产aⅴ精品一区二区三区波| 99热这里只有是精品50| 日韩大尺度精品在线看网址| 久久午夜综合久久蜜桃| 国产成人av教育| 国产精华一区二区三区| 91成年电影在线观看| 欧美一区二区国产精品久久精品 | 女同久久另类99精品国产91| 国产av在哪里看| 黄色女人牲交| 69av精品久久久久久| x7x7x7水蜜桃| 亚洲avbb在线观看| 法律面前人人平等表现在哪些方面| 亚洲av成人av| 欧美日韩中文字幕国产精品一区二区三区| 黄色视频不卡| 久久婷婷人人爽人人干人人爱| 日韩成人在线观看一区二区三区| 欧美三级亚洲精品| 日本黄色视频三级网站网址| 国产三级在线视频| 美女高潮喷水抽搐中文字幕| 91国产中文字幕| 亚洲国产精品久久男人天堂| 国产精品1区2区在线观看.| 亚洲av成人不卡在线观看播放网| 国产高清激情床上av| 午夜免费观看网址| 亚洲国产高清在线一区二区三| 99在线视频只有这里精品首页| 国产野战对白在线观看| 精品电影一区二区在线| 日韩免费av在线播放| 真人做人爱边吃奶动态| 可以在线观看毛片的网站| xxx96com| 国产1区2区3区精品| 亚洲 欧美一区二区三区| 丁香六月欧美| 日韩欧美在线乱码| 日本熟妇午夜| 婷婷亚洲欧美| 好男人电影高清在线观看| 久久久久久九九精品二区国产 | 欧美日韩一级在线毛片| 91大片在线观看| 欧美久久黑人一区二区| 日韩三级视频一区二区三区| 丰满人妻一区二区三区视频av | 亚洲 欧美一区二区三区| 国产成人欧美在线观看| 脱女人内裤的视频| 麻豆成人av在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 三级毛片av免费| 少妇的丰满在线观看| 淫妇啪啪啪对白视频| 好看av亚洲va欧美ⅴa在| 黄色成人免费大全| 我的老师免费观看完整版| 亚洲av片天天在线观看| 人妻丰满熟妇av一区二区三区| 国产激情偷乱视频一区二区| 中文字幕av在线有码专区| 欧美日韩国产亚洲二区| 三级男女做爰猛烈吃奶摸视频| 亚洲专区字幕在线| 悠悠久久av| 国产成人系列免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 一级毛片高清免费大全| 两个人的视频大全免费| 少妇被粗大的猛进出69影院| 极品教师在线免费播放| 精品熟女少妇八av免费久了| 一本精品99久久精品77| 亚洲专区字幕在线| 欧美日韩瑟瑟在线播放| 日本黄大片高清| 啦啦啦免费观看视频1| 精品国产美女av久久久久小说| 91成年电影在线观看| 精品免费久久久久久久清纯| 中国美女看黄片| 日日摸夜夜添夜夜添小说| 精品日产1卡2卡| 别揉我奶头~嗯~啊~动态视频| 啪啪无遮挡十八禁网站| 男人的好看免费观看在线视频 | xxxwww97欧美| 久久午夜综合久久蜜桃| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲专区中文字幕在线| 久久久国产欧美日韩av| 黄色女人牲交| 国产av一区二区精品久久| 在线永久观看黄色视频| 人人妻人人澡欧美一区二区| 日本黄色视频三级网站网址| 亚洲中文字幕一区二区三区有码在线看 | 国产精品,欧美在线| 国产高清视频在线播放一区| 亚洲国产精品合色在线| 亚洲精品一卡2卡三卡4卡5卡| 精品不卡国产一区二区三区| 91字幕亚洲| 国产高清视频在线播放一区| 婷婷丁香在线五月| 成人亚洲精品av一区二区| 精品国内亚洲2022精品成人| 99国产精品一区二区三区| 国产三级在线视频| 99国产精品一区二区蜜桃av| 级片在线观看| 国产精品亚洲一级av第二区| 两个人免费观看高清视频| 亚洲黑人精品在线| 亚洲精品av麻豆狂野| 国产一区二区三区视频了| 国产激情久久老熟女| 亚洲在线自拍视频| 久久精品亚洲精品国产色婷小说| 亚洲国产精品合色在线| 天天一区二区日本电影三级| 久久久久久免费高清国产稀缺| 啦啦啦免费观看视频1| 天天添夜夜摸| 日本一二三区视频观看| 真人做人爱边吃奶动态| 亚洲专区中文字幕在线| 天天添夜夜摸| 两人在一起打扑克的视频| netflix在线观看网站| 一级a爱片免费观看的视频| 一个人免费在线观看电影 | 国产亚洲av嫩草精品影院| 夜夜夜夜夜久久久久| 中文资源天堂在线| 午夜免费成人在线视频| 久久天堂一区二区三区四区| 亚洲,欧美精品.| 亚洲狠狠婷婷综合久久图片| 精品乱码久久久久久99久播| 国产一区二区激情短视频| 久久久国产精品麻豆| 岛国在线免费视频观看| 久久天躁狠狠躁夜夜2o2o| 成人18禁在线播放| 久久香蕉精品热| 他把我摸到了高潮在线观看| 18禁裸乳无遮挡免费网站照片| 国产成人av激情在线播放| 国产又黄又爽又无遮挡在线| 9191精品国产免费久久| 国产精品 国内视频| 俺也久久电影网| 99久久精品国产亚洲精品| 久久国产精品人妻蜜桃| 国产午夜精品论理片| 日本免费一区二区三区高清不卡| 国产精品香港三级国产av潘金莲| 女同久久另类99精品国产91| 天天躁夜夜躁狠狠躁躁| 亚洲中文av在线| 99在线视频只有这里精品首页| 精品国产乱子伦一区二区三区| 国产一级毛片七仙女欲春2| 亚洲乱码一区二区免费版| 亚洲精品在线观看二区| 黄色女人牲交| 嫁个100分男人电影在线观看| 夜夜夜夜夜久久久久| 免费在线观看视频国产中文字幕亚洲| 色综合亚洲欧美另类图片| 这个男人来自地球电影免费观看| 日日爽夜夜爽网站| 男女午夜视频在线观看| 日本a在线网址| 成年女人毛片免费观看观看9| 黄色成人免费大全| 哪里可以看免费的av片| 国产成人欧美在线观看| 麻豆一二三区av精品| av天堂在线播放| 欧美性猛交╳xxx乱大交人| 国产区一区二久久| 国产视频一区二区在线看| 免费在线观看完整版高清| 国产成人一区二区三区免费视频网站| 欧美不卡视频在线免费观看 | 亚洲一区高清亚洲精品| 久久香蕉精品热| 天天添夜夜摸| 天堂影院成人在线观看| 最好的美女福利视频网| 午夜福利视频1000在线观看| 国产伦在线观看视频一区| 国产av一区二区精品久久| 18禁黄网站禁片午夜丰满| 成年版毛片免费区| 国产成年人精品一区二区| 国产精华一区二区三区| 欧美一区二区国产精品久久精品 | 国内毛片毛片毛片毛片毛片| 亚洲av美国av| 国产精品 国内视频| 淫妇啪啪啪对白视频| 在线观看66精品国产| av欧美777| 欧美日韩乱码在线| 女人被狂操c到高潮| 亚洲最大成人中文| 久久国产乱子伦精品免费另类| 女生性感内裤真人,穿戴方法视频| 麻豆成人午夜福利视频| 精品国内亚洲2022精品成人| 999精品在线视频| 久久中文字幕一级| 看片在线看免费视频| 一区福利在线观看| 亚洲精品久久国产高清桃花| 伦理电影免费视频| 中文亚洲av片在线观看爽| 丰满的人妻完整版| 精品久久久久久久末码| 老司机午夜福利在线观看视频| 1024视频免费在线观看| 嫁个100分男人电影在线观看| 欧美大码av| 妹子高潮喷水视频| 欧美日韩中文字幕国产精品一区二区三区| 中国美女看黄片| 亚洲欧美精品综合久久99| 成人永久免费在线观看视频| 亚洲国产精品sss在线观看| av在线天堂中文字幕| 非洲黑人性xxxx精品又粗又长| 色老头精品视频在线观看| 亚洲激情在线av| 最新在线观看一区二区三区| 好男人电影高清在线观看| 丝袜人妻中文字幕| 18禁黄网站禁片免费观看直播| 最近最新中文字幕大全免费视频| 看片在线看免费视频| 国产私拍福利视频在线观看| 少妇的丰满在线观看| 每晚都被弄得嗷嗷叫到高潮| 日韩高清综合在线| 免费在线观看视频国产中文字幕亚洲| 国产精品爽爽va在线观看网站| 制服人妻中文乱码| 日本 av在线| 亚洲欧美日韩高清专用| 国产精品一区二区三区四区免费观看 | 18美女黄网站色大片免费观看| 黄色毛片三级朝国网站| 亚洲国产欧洲综合997久久,| 精品久久蜜臀av无| 亚洲一区中文字幕在线| 淫秽高清视频在线观看| 蜜桃久久精品国产亚洲av| 日韩精品中文字幕看吧| 香蕉丝袜av| 99精品在免费线老司机午夜| 久久久久久久久久黄片| 国产高清videossex| 天天躁狠狠躁夜夜躁狠狠躁| 在线看三级毛片| 色哟哟哟哟哟哟| 日日干狠狠操夜夜爽| 中文字幕高清在线视频| 99在线人妻在线中文字幕| svipshipincom国产片| 2021天堂中文幕一二区在线观| 国产欧美日韩一区二区精品| 国产亚洲精品久久久久久毛片| 麻豆国产97在线/欧美 | 老司机靠b影院| 日日干狠狠操夜夜爽| 久久久久久久久免费视频了| 亚洲国产精品sss在线观看| 久久午夜亚洲精品久久| 欧美 亚洲 国产 日韩一| 久久精品影院6| 国产激情久久老熟女| 91字幕亚洲| 一边摸一边抽搐一进一小说| 在线观看舔阴道视频| 亚洲男人天堂网一区| 亚洲成a人片在线一区二区| 亚洲欧美精品综合一区二区三区| 我要搜黄色片| 国产激情欧美一区二区| 黄片小视频在线播放| 又爽又黄无遮挡网站| 国产爱豆传媒在线观看 | 欧美乱色亚洲激情| 丁香欧美五月| 一进一出抽搐动态| 国产精品一区二区精品视频观看| 国产伦人伦偷精品视频| 男男h啪啪无遮挡| 亚洲av美国av| 国产激情久久老熟女| 国产三级中文精品| 青草久久国产| 精品少妇一区二区三区视频日本电影| 精品国产乱码久久久久久男人| 高潮久久久久久久久久久不卡| www日本黄色视频网| 一本一本综合久久| 国产精品久久久久久久电影 | 亚洲欧美日韩东京热| 变态另类丝袜制服| 国产片内射在线| 国产三级在线视频| 在线观看午夜福利视频| 亚洲色图av天堂| 久久中文字幕人妻熟女| 一区二区三区国产精品乱码| 99在线人妻在线中文字幕| 久久草成人影院| 国产黄片美女视频| 亚洲国产欧洲综合997久久,| 在线看三级毛片| 亚洲午夜精品一区,二区,三区| 宅男免费午夜| 午夜精品在线福利| 成人欧美大片| 大型黄色视频在线免费观看| 久久久久久久午夜电影| 亚洲av电影不卡..在线观看| 久久国产乱子伦精品免费另类| 老司机午夜十八禁免费视频| 国产三级中文精品| 日本撒尿小便嘘嘘汇集6| 999久久久国产精品视频| 麻豆国产av国片精品| 一边摸一边做爽爽视频免费| 男人的好看免费观看在线视频 | 成人av在线播放网站| 国产午夜福利久久久久久| 在线观看午夜福利视频| 日韩欧美在线乱码| 宅男免费午夜| 国产亚洲欧美98| 欧美不卡视频在线免费观看 | 久久草成人影院| 欧美精品亚洲一区二区| 岛国视频午夜一区免费看| 国产成人av教育| 国产精品国产高清国产av| 久久精品影院6| 欧美色视频一区免费| 成人永久免费在线观看视频| 色综合亚洲欧美另类图片| 欧美成狂野欧美在线观看| 黄色a级毛片大全视频| 精品不卡国产一区二区三区| 国产亚洲精品久久久久5区| 高潮久久久久久久久久久不卡| 黄片小视频在线播放| 最近最新中文字幕大全电影3| 1024香蕉在线观看| 黑人欧美特级aaaaaa片| 长腿黑丝高跟| 成人三级做爰电影| 国产爱豆传媒在线观看 | www日本在线高清视频| 一级a爱片免费观看的视频| 亚洲中文字幕一区二区三区有码在线看 | 精品免费久久久久久久清纯| 18美女黄网站色大片免费观看| 久久久精品大字幕| 又粗又爽又猛毛片免费看| 国产亚洲精品久久久久5区| 国产精品 欧美亚洲| 国产视频一区二区在线看| a级毛片a级免费在线| 亚洲精品国产一区二区精华液| 成人国语在线视频| 亚洲av成人不卡在线观看播放网| 精品无人区乱码1区二区| 97碰自拍视频| 亚洲成a人片在线一区二区| 村上凉子中文字幕在线| 叶爱在线成人免费视频播放| 99国产精品一区二区三区| 最近最新中文字幕大全免费视频| av免费在线观看网站| 久久精品成人免费网站| 亚洲中文字幕一区二区三区有码在线看 | 变态另类丝袜制服| 国产精品电影一区二区三区| 99在线人妻在线中文字幕| 精品国内亚洲2022精品成人| 亚洲国产欧美人成| 波多野结衣巨乳人妻| 中文字幕熟女人妻在线| 亚洲一区二区三区不卡视频| 日本a在线网址| 欧美不卡视频在线免费观看 | 日本黄大片高清| 1024手机看黄色片| 男人舔奶头视频| 欧美3d第一页| 欧美日韩瑟瑟在线播放| 男男h啪啪无遮挡| 久久香蕉精品热| 亚洲欧美一区二区三区黑人| 国产在线精品亚洲第一网站| 国产三级中文精品| 色精品久久人妻99蜜桃| 精品高清国产在线一区| 两个人免费观看高清视频| 国产69精品久久久久777片 | 亚洲精品粉嫩美女一区| 观看免费一级毛片| 男女做爰动态图高潮gif福利片| 性色av乱码一区二区三区2| 国产亚洲av嫩草精品影院| 久久国产精品影院| 变态另类丝袜制服| 91大片在线观看| 成人三级黄色视频| 精品久久久久久久久久免费视频| 精品一区二区三区视频在线观看免费| 一级片免费观看大全| 又黄又爽又免费观看的视频| 听说在线观看完整版免费高清| 亚洲 欧美一区二区三区| 精品国产美女av久久久久小说| e午夜精品久久久久久久| 少妇粗大呻吟视频| 午夜免费激情av| 国产精品乱码一区二三区的特点| 国产亚洲欧美在线一区二区| 午夜久久久久精精品| 在线观看一区二区三区| 亚洲成a人片在线一区二区| 久久久久久九九精品二区国产 | 看黄色毛片网站| 女生性感内裤真人,穿戴方法视频| 男女那种视频在线观看| 国产免费av片在线观看野外av| 国产午夜精品久久久久久| 免费在线观看亚洲国产| 精品一区二区三区视频在线观看免费| 免费看十八禁软件| 午夜激情av网站| 亚洲七黄色美女视频| 中文资源天堂在线| 国产一区二区三区在线臀色熟女| 亚洲 欧美一区二区三区| 十八禁网站免费在线| 久久国产精品人妻蜜桃| 久久天堂一区二区三区四区| 蜜桃久久精品国产亚洲av| av在线播放免费不卡| 人妻久久中文字幕网| 无遮挡黄片免费观看| 精品一区二区三区四区五区乱码| 久久久久久九九精品二区国产 | 久久久久国产精品人妻aⅴ院| 制服丝袜大香蕉在线| 国产熟女xx| av在线天堂中文字幕| 午夜日韩欧美国产| 中文字幕熟女人妻在线| 国产精品野战在线观看| 五月玫瑰六月丁香| 制服人妻中文乱码| 国产91精品成人一区二区三区| 久久伊人香网站| 两个人免费观看高清视频| 午夜成年电影在线免费观看| 五月伊人婷婷丁香| 99久久综合精品五月天人人| 午夜影院日韩av| 又爽又黄无遮挡网站| 亚洲成人精品中文字幕电影| 久久久久久亚洲精品国产蜜桃av| 欧美精品亚洲一区二区| 少妇粗大呻吟视频| 亚洲激情在线av| 久久精品国产99精品国产亚洲性色| 国产精品影院久久| 精品无人区乱码1区二区| 又粗又爽又猛毛片免费看| 狂野欧美白嫩少妇大欣赏| 91字幕亚洲| 美女 人体艺术 gogo| 波多野结衣高清无吗| x7x7x7水蜜桃| 国产av在哪里看| 亚洲一区高清亚洲精品| 久久人人精品亚洲av| 久久久久久亚洲精品国产蜜桃av| 国产三级在线视频| 成人国语在线视频| 美女免费视频网站| 亚洲精品国产一区二区精华液| av在线播放免费不卡| 法律面前人人平等表现在哪些方面| 国产主播在线观看一区二区| 在线播放国产精品三级| 黑人操中国人逼视频| 欧美黄色片欧美黄色片| 怎么达到女性高潮| 91成年电影在线观看| 国产伦一二天堂av在线观看| 日本免费a在线| 两性午夜刺激爽爽歪歪视频在线观看 | 天堂影院成人在线观看| av中文乱码字幕在线| 最新在线观看一区二区三区| 欧美日本亚洲视频在线播放| 夜夜夜夜夜久久久久| 亚洲自拍偷在线| 欧美成人免费av一区二区三区| 午夜福利欧美成人| 国产精品 欧美亚洲| 午夜免费观看网址| 国产人伦9x9x在线观看| 日本精品一区二区三区蜜桃| 日韩大尺度精品在线看网址| 国产精品日韩av在线免费观看| 男男h啪啪无遮挡| 99国产综合亚洲精品| 老司机深夜福利视频在线观看| 最近最新中文字幕大全免费视频| 免费无遮挡裸体视频| 国模一区二区三区四区视频 | 成人国产一区最新在线观看| 99热6这里只有精品| xxxwww97欧美| 最近最新中文字幕大全电影3| 美女黄网站色视频| 一级毛片女人18水好多| 亚洲无线在线观看| 免费观看人在逋| 国产精品野战在线观看| 亚洲精品国产精品久久久不卡| 国产人伦9x9x在线观看| 亚洲国产精品999在线| 97超级碰碰碰精品色视频在线观看| 国产在线精品亚洲第一网站| 亚洲精品一区av在线观看| 欧美丝袜亚洲另类 | 美女扒开内裤让男人捅视频| 日韩欧美一区二区三区在线观看| 国产又黄又爽又无遮挡在线| 亚洲男人天堂网一区| 丰满人妻熟妇乱又伦精品不卡| 99久久精品热视频| 成人一区二区视频在线观看| 天天添夜夜摸| 日本熟妇午夜| 亚洲国产精品成人综合色| 亚洲专区国产一区二区| 一夜夜www| 黑人操中国人逼视频| 欧美日韩乱码在线| 手机成人av网站| 69av精品久久久久久| 欧美精品亚洲一区二区| 在线观看免费日韩欧美大片| 日韩精品青青久久久久久| 老司机午夜福利在线观看视频| 一进一出好大好爽视频| 一级毛片高清免费大全| 国产成人精品无人区| 国产一级毛片七仙女欲春2| av福利片在线| 无遮挡黄片免费观看| 精品久久久久久久久久久久久| 女生性感内裤真人,穿戴方法视频| 免费观看精品视频网站| 在线观看美女被高潮喷水网站 | 一本综合久久免费| 亚洲精华国产精华精| 一区二区三区高清视频在线| 精品久久蜜臀av无| 亚洲男人天堂网一区| 中文字幕av在线有码专区| 舔av片在线| 欧美 亚洲 国产 日韩一| 久久久久九九精品影院| www.自偷自拍.com| 亚洲avbb在线观看| 国产av不卡久久| 两性午夜刺激爽爽歪歪视频在线观看 | 日本三级黄在线观看| 美女黄网站色视频| 最近视频中文字幕2019在线8| 日日摸夜夜添夜夜添小说| 国产亚洲精品久久久久5区| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久人妻精品电影| 久久婷婷成人综合色麻豆| 老司机午夜十八禁免费视频|