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

    Non-Hermitian quasicrystal in dimerized lattices?

    2021-10-28 07:16:10LongwenZhou周龍文andWenqianHan韓雯岍
    Chinese Physics B 2021年10期
    關(guān)鍵詞:龍文

    Longwen Zhou(周龍文) and Wenqian Han(韓雯岍)

    College of Physics and Optoelectronic Engineering,Ocean University of China,Qingdao 266100,China

    Keywords: non-Hermitian system,quasicrystal,localization,topological phase

    1. Introduction

    Quasicrystals have long-range order without spatial periodicity. They form a class of system in between crystals and fully disordered lattices.[1–6]Experimentally, quasicrystals have been realized in a variety of solid state materials and quantum simulators.[7–13]Rich phenomena induced by quasiperiodicity have been revealed, such as topological phases,[14–21]quantized adiabatic pumping,[15,20,21]anomalous transport and localization transitions,[22–51]attracting attention over a broad range of research fields.

    Recently, the study of quasicrystals has been extended to non-Hermitian systems, where the interplay between quasiperiodicity and gain/loss or nonreciprocity could induce exotic dynamical, localization and topological phenomena.[52–69]In particular, in non-Hermitian variants of the Aubry–Andr′e–Harper (AAH) model, complex onsite potential or nonreciprocal hopping could induce aPT-breaking transition and a metal–insulator transition, which can be further characterized by a spectral topological winding number.[52,56]In related studies, non-Hermiticity induced mobility edges in generalized and superconducting AAH models have also been found and described by topological invariants.[54,57–59,62,67]Besides,the investigation of wavepacket spreading in non-Hermitian quasicrystals (NHQCs) has revealed their anomalous dynamical features,[55,60]such as the disordered-enhanced transport.[55]

    Meanwhile, the Su–Schrieffer–Heeger (SSH) model[70]provides another paradigm for the study of localization,[71]topological[72]and non-Hermitian[73–75]physics. Specially,the hopping dimerization allows the SSH model to possess a topological phase characterized by an integer winding number and degenerate edge modes.[72]However, in the context of non-Hermitian quasicrystals, phases and phenomena that could arise due to the interplay between spatial aperiodicity and hopping dimerization have not been revealed. In this paper, we address this issue by introducing a dimerized quasiperiodic lattice in Section 2, which forms a hybridization of the SSH and non-Hermitian AAH models. The system is found to possess rich patterns ofPT-breaking and localization transitions,together with three phases of distinct transport nature. In Section 3,we perform detailed analyses of the spectrum and localization nature of these phases, and construct a pair of topological winding numbers to characterize the transitions between them. Despite an extended and a localized phase, we also find a mobility edge phase, which is absent without the hopping dimerization. These discoveries thus uncover the richness of topological and transport phenomena in dimerized NHQCs. In Section 4, we summarize our results and discuss potential future directions. Further details about the spectrum,Lyapunov exponents and wavepacket dynamics are provided in the Appendix.

    2. Model

    In this section, we introduce the dimerized NHQC that will be investigated in this work. Our model can be viewed as an extension of the AAH model,[5,6]which is prototypical in the study of localization transitions in one-dimensional (1D) quasicrystals. In position representation, the Hamiltonian of the AAH model is ?HAAH= ∑n(J|n〉〈n+1|+H.c.+Vcos(2παn+λ)|n〉〈n|),where{|n〉}represents the eigenbasis of the lattice,Jis the hopping amplitude,Vcontrols the strength of the onsite potential, andλis a phase shift. Whenαis irrational, the potentialVn=Vcos(2παn+λ)is quasiperiodic inn,and ?HAAHdescribes a 1D quasicrystal. WhenV<2J, the spectrum of ?HAAHis continuous and all its eigenstates are extended under the periodic boundary condition(PBC).Comparatively, ?HAAHpossesses a point spectrum with localized eigenstates at all energies whenV>2J. WhenV=2J, the spectrum of ?HAAHis purely singular continuous with critical wavefunctions,and the system undergoes a localization transition whenVchanges fromV<2JtoV>2J.[4]

    Recently, a non-Hermitian variant of ?HAAHwas introduced by settingλ=β+iγ, with (β,γ)∈R.[52]Such a non-Hermitian AAH model possessesPT-breaking and localization transitions atγc=ln(2J/V), which are accompanied by the quantized jump of a spectral winding numberw.Forγ<γc(γ>γc),the spectrum is real(complex)with winding numberw=0(w=?1),and each eigenstate is extended(localized).[52]Since all states subject to the same localization transition atγ=γc, no mobility edges are found. Similar results were reported in a nonreciprocal AAH model,[56]which might be related to the model in Ref.[52]by Fourier transformations.

    Fig.1. Schematic illustration of the dimerized NHQC. Solid balls denote lattice sites with index n ∈Z. Hopping amplitudes are alternated between J ?Δ (solid lines)and J+Δ (dashed lines). Vn represents the strength of the onsite potential,with irrational modulation parameter α and phase shift iγ.

    In this work, we consider another extension of the AAH quasicrystal by introducing hopping dimerizations. We setJ →Jn=J+(?1)nΔas the hopping amplitude between then’s and the(n+1)’s lattice sites,and letλ=iγbe the imaginary phase shift in the onsite potentialVn. Our model Hamiltonian thus takes the form

    A schematic illustration of the model is shown in Fig.1.WhenV=0,Eq.(1)reduces to the SSH model with hopping dimerizationΔ. WithV,γ/= 0, ?Hdescribes an SSH chain with quasiperiodically correlated non-Hermitian disorder. Therefore, our system realizes a hybridization between the AAH and SSH models. Since the realizations of SSH and non-Hermitian AAH models have both been discussed,[52,76]we expect our model to be within reach under current or near term experimental conditions. Furthermore, our construction allows us to explore the interplay between hopping dimerizations and quasiperiodic non-Hermitian potential,which could lead to rich phase diagrams and transport phenomena,as will be shown in the following section.

    3. Results

    3.1. Spectrum and level statistics

    The spectrum of the dimerized NHQC is obtained by solving the eigenvalue equation ?H|ψ〉=E|ψ〉. Projecting the equation to the lattice representation,we find

    Fig.2. (a)Density of states with complex eigenvalues and(b)averaged AGRs of ?H versus the dimerization amplitude Δ and the imaginary part of phase shift γ. System parameters are set as J =1, V =0.5, and. PBC is taken in the calculation and the length of the lattice is L=2584.

    To have a comprehensive view of the spectrum and its connection with localization properties of states, we present the density of states with nonzero imaginary parts of energies and the mean of AGRs versusγandΔin Figs.2(a)and 2(b),respectively. The density of states with complex energies is defined as

    3.2. Inverse participation ratio

    When max(IPR)diverges from zero but min(IPR)remains at zero, localized states start to appear and the system switches from the extended to the mobility edge phase.When min(IPR)also deviates from zero, the last extended state vanishes and the system enters a phase in which all states are localized. The maximal and minimal values of IPR can thus be used to distinguish phases with different localization nature and locate the corresponding phase boundaries.

    Fig.3. (a)The maximum and(b)minimum of IPRs versus the imaginary part of phase shift γ and hopping dimerization Δ. System parameters are set as J=1,V =0.5,and α =. The length of the lattice is L=610. PBC is taken in the diagonalization ?H.

    To summarize, we find that the IPR can be viewed as an“order parameter”to describe the localization transitions in the dimerized NHQC. The behaviors of IPRs in distinct parameter domains could be employed to characterize the three different phases in the system. For completeness, we study the Lyapunov exponents of the system in Appendix B, and find consistent results as predicted by IPRs. We have also checked the spectrum and IPRs of the system under the open boundary condition (OBC), and find consistent results with those obtained under the PBC, excluding possible impact of non-Hermitian skin effects. In the following,we demonstrate that the localization transitions in our system are of topological origin,and can be depicted by a pair of spectral winding numbers.

    3.3. Topological winding number

    In previous studies,spectral winding numbers have been employed to characterize non-Hermitian topological matter in 1D clean and disordered systems.[52,73]For the dimerized NHQC, we can introduce a pair of winding numbers to describe the topological nature of its localization transitions.These numbers can be defined as

    HeteNis the number of dimerized cells of the lattice and?=1,2. The phase shiftβ/Nis introduced into the Hamiltonian ?Hvia settingVn →Vcos(2παn+β/N+iγ). (?1,?2)are two real-valued base energies,and(w1,w2)count the number of times that the spectrum ofH(β/N)winds around these energies whenβsweeps over a cycle from zero to 2π. It is clear thatw1(w2) can be nonzero only if the spectrum of ?Haround?1(?2)takes complex values. These winding numbers are thus closely related to the complex spectrum structure of the system. When the spectrum does not possess a mobility edge,there is only a single base energy that can in principle be chosen arbitrarily,[52]and we would always havew1=w2. In our model this is the case whenΔ=0(uniform hopping). If mobility edges exist in the spectrum,the choice of base energies(?1,?2)should be related to its boundaries.[54]More precisely,in a given range of imaginary phase shift iγor hopping dimerizationΔ,we choose?1to be the real part of energy of the first eigenstate of ?Hwhose IPR starts to deviate from zero,i.e., the first eigenstate that becomes localized. Similarly, we set?2as the real part of energy of the last eigenstate whose profile changes from extended to localized.?1and?2thus decide the lower and upper bounds of the mobility edge on the ReE–γor ReE–Δplane. The winding numbersw1andw2defined with respect to?1and?2are expected to have quantized jumps when the mobility edge appears and vanishes in the spectrum.[54]

    Fig. 4. The winding numbers (w1,w2) versus the imaginary part of phase shift γ. System parameters are set as J =1,V =0.5, Δ =0.4,α =. The size of the lattice is L=610 with PBC. The vertical dotted lines highlight the critical values γc1 and γc2, where transitions between different phases happen.

    In Fig. 5, we present the topological phase diagrams of the dimerized NHQC versus the imaginary part of phase shift and hopping dimerization for two typical sets of system parameters. In the region ofρ ?0 [real spectrum, see also Fig.2(a)]and max(IPR)?0[all bulk states are extended,see also Fig. 3(a)], we findw1=w2=0, implying that thePTinvariant extended phase is topologically trivial according to the spectral winding numbers. Whenρa(bǔ)nd max(IPR) start to deviate from zero, the system undergoes aPT-breaking transition and enters a mobility edge phase. In the meantime,w1takes a quantized jump from 0 to?1, whereasw2remains at zero in Fig. 5(a). The transition from extended to mobility edge phases is thus topological and captured by the quantized change of winding numberw1. The mobility edge phase can also be viewed as a topological phase characterized by (w1,w2)=(?1,0). Whenρ ?1 and min(IPR) starts to deviate from zero,the system enters a localized phase and all bulk states take complex eigenvalues.Meanwhile,the winding numberw2jumps from 0 to?1 whereasw1remains at?1,as shown in Fig.5(a). Therefore,the transition from the mobility edge to localized phases is also topological and accompanied by the quantized change of winding numberw2by?1.The localized phase can then be viewed as a topological phase with(w1,w2)=(?1,?1). These observations are all demonstrated in the phase diagram Fig.5(a). We also considered other possible amplitudes of onsite potentialV[with one example given in Fig.5(b)],and obtain similar kinds of diagrams,which verifies the generality of our approach to the characterization of localization transitions and topological phases in dimerized NHQCs. We summarize the key results of this section in Table 1.

    Fig. 5. Topological phase diagrams of the dimerized NHQC. System parameters are set as J=1,α =,and V =0.5(V =1)for panel(a)[panel(b)]. The length of the lattice is L=610,containing N=305 dimerized cells. In panels(a)and(b),each region with a uniform color corresponds to a topological phase, whose winding numbers are denoted therein. The extended, mobility edge and localized phases have(w1,w2)=(0,0),(?1,0)and(?1,?1). w1 or w2 jumps at a boundary between different phases. PBC is taken in the calculation.

    Table 1. Summary of the results for the dimerized NHQC.

    4. Discussion

    In experiments, the dimerized NHQC might be engineered in photonic systems. The uniform part of hopping amplitude and non-Hermitian quasiperiodic potential could be realized by a frequency-modulated mode-locked laser with gain medium, phase modulator and low-finesse intracavity etalon,as proposed in Ref. [52]. The dimerized hopping amplitude could be realized by engineering the profile of refractive index in the model-locked laser setup.[76]Therefore,our model should be within reach in current or near-term experimental situations. To promote the detection and characterization of different phases in the dimerized NHQC, we also investigate its wavepacket dynamics in Appendix C,and find connections between the dynamical signatures and localization properties of the system.

    In conclusion, we find localization and topological transitions in a dimerized NHQC, which are originated from the cooperation between hopping dimerizations and complex onsite quasiperiodic potential.In the region of weak dimerization and non-Hermiticity, the system is in an extended phase with real spectrum and delocalized eigenstates. With the increase of hopping dimerization and complex potential, the system transforms into a mobility edge phase. When the strength of hopping dimerization and non-Hermitian modulation become stronger,the system enters a third phase in which the spectrum is complex and all eigenstates are localized. Moreover, the transitions between the extended,mobility edge and localized phases are of topological nature. They can be characterized by the quantized jumps of two spectral winding numbers. Our results thus uncover the unique spectrum,topological and transport features of quasicrystals due to the interplay between hopping dimerizations and non-Hermitian onsite potential. The different phases and transitions found in our system further reveal the richness of localization and topological phenomena in non-Hermitian quasicrystals. In future work, it would be interesting to consider the impact of hopping dimerizations in other types of nonreciprocal and non-Hermitian quasicrystals,and explore the effect of nonlinearity,many-body interactions and skin effects on the localization and topological physics in dimerized non-Hermitian systems.

    Appendix A:Details of the spectrum

    In this appendix,we provide more details about the spectrum of the dimerized NHQC.In Fig.A1,we present the spectrum for a typical set of system parameters. When the imaginary part of phase shiftγis small, we observe real spectrum in Figs. A1(a) and A1(d), implying that the system is in thePT-invariant region.

    Fig. A1. The spectrum of the dimerized NHQC model under PBC on the complex plane. System parameters are set as J =1, V =0.5,α =,and Δ =0.2(Δ =0.4)for panels(a)–(c)[(d)–(f)]. The values of the imaginary part of phase shift are shown in the panel captions.The length of the lattice is L=2584.

    With the increase ofγ, the spectrum starts to become complex and develop loops on the complexEplane, which means that the system has undergone aPT-transition and roamed into aPT-broken region. However,parts of the spectrum are still pinned to the real axis in Figs.A1(b)and A1(e),indicating that the eigenstates with real and complex eigenvalues coexist in these cases. Whenγfurther increases, the number of eigenstates with real eigenvalues tends to decrease,and finally almost all states have complex energies,as shown in Figs.A1(c)and A1(f). Besides,the range and shape of the spectrum also depend on the dimerization strengthΔ, which reveals that both the hopping dimerization and complex onsite potential could affect the properties of the spectrum.

    Appendix B:Lypunov exponent

    Fig.B1. IPRs and LEs versus the imaginary part of phase shift γ. System parameters are J =1,V =0.5, and α = The length of the lattice is L=610, and PBC (OBC) is taken in the calculation of IPRs(LEs). In panels (a)–(d), the solid (dashed) lines represent the maximum (minimum) of IPRs and the dotted (dash-dotted) lines show the derivative of the averaged(minimum of)LE.

    Fig. B2. The minimum and derivative of LEs, as defined in Eqs.(B2)and(B3),versus the imaginary part of phase shift γ and hopping dimerization Δ. System parameters are set as J=1,V =0.5,and α=The length of the lattice is L=610 and the OBC is taken in the diagonalization ?H.

    Appendix C:Wavepacket dynamics

    In this appendix, we investigate the dynamical properties of the dimerized NHQC, and suggest to characterize the different phases by the spreading velocity and return probability of wavepackets. The evolution of a state in our sys-

    In Fig.C1, we show the profiles of a wavepacket during the evolution and the averaged spreading velocity for a typical set of system parameters. In Fig. C1(a), the system is set in the extended phase (see also Fig. 4). We observe that the wavepacket indeed performs a ballistic spreading and presents a light cone pattern. In Fig. C1(b), the system is set in the mobility edge phase(see also Fig.4). The wavepacket is initially found to undergo a slow spreading over a finite range of sites. At a later stage,the spreading tends to terminate and the wavepacket retains a finite width around its original site. The initial spreading process may be assisted by the remaining extended states of the mobility edge phase.But during the evolution,the wavepacket develops more overlap with the localized states, which finally help to shut off its transport in the lattice. In Fig.C1(c),the parameters of ?Hare set in the localized phase(see also Fig.4). As expected,the wavepacket could not spread and remains exponentially localized around its original location. The distinctive signatures of wavepacket spreading in these three phases could then provide us with a dynamical way of distinguishing them. In Fig. C1(d), the spreading velocity of a wavepacket is found to decrease with the increase ofγwhen the system is set in the extended phase. At largeγ, the system enters the localized phase and the velocity approaches zero. Whenγgoes across the boundary between the extended and mobility edge phases, the velocity shows a quick drop in a small range ofγ, implying the appearance of a dynamical transition between these two phases. Notably,the critical pointγc1of the transition is different for different hopping dimerizations, and the velocity could also show an anomalous growth with the increase ofγin the mobility edge phase [e.g., see the solid line in Fig. C1(d)]. These observations suggest that the interplay between hopping dimerizations and complex onsite potential could not only create new types of dynamical phases in non-Hermitian quasicrystals, but also cause non-Hermiticity enhanced transport observed previously in systems with nonreciprocal hopping.[55]

    Fig.C1. Profiles [in panels (a)–(c)]and spreading velocities [in panel(d)] of wavepackets initially localized at the central site of the lattice.The imaginary phase shift is set as(a)γ=0.4,(b)γ=1 and(c)γ=1.6.System parameters are J=1,V =0.5, Δ =0.4 in panels(a)–(c). The length of the lattice is L=1000 with PBC.The spreading velocity v is obtained in panel(d)after averaging v(t)over a time span of t=250.

    Fig.C2. Spreading velocities of an initially localized wavepacket versus the imaginary part of phase shift γ and hopping dimerization Δ.v(t)is averaged over a time duration t=250. The strength of the onsite potential is V =0.5 in panel(a)and V =1 in panel(b).System parameters are set as J =1 and α =. The length of the lattice is L=1000 with PBC,and the initial state satisfies〈n|ψ(0)〉=δn0.

    In Fig. C3, we show the return probabilities of the initial states prepared in and quenched to different phases of the dimerized NHQC. In Fig. C3(a), the system is initialized in the extended phase with state profile|ψi|. When the postquench system is in the same phase,P(t) undergoes oscillations without a global decay, as shown by the blue solid line in Fig. C3(d). If the postquench Hamiltonian is in the mobility edge phase,P(t) first subjects to oscillations with a global decay profile,and approaches zero in the long-time domain, as demonstrated by the red dashed line in Fig. C3(d).The initial oscillations ofP(t)can be traced back to the presence of extended states in the mobility edge phase. With the progress of time, the state is trapped by localized states in the mobility edge phase, leading to the decay ofP(t)at larget. If ?Hfis in the localized phase,P(t) decays very fast and quickly approaches zero,as depicted by the yellow dotted line in Fig.C3(d). Since all states are localized in the postquench system, the evolving state can never go back. The distinctive features ofP(t)in three phases can thus help us to discriminate them following the postquench dynamics if the initial state is prepared in the extended phase.

    In Figs. C3(b) and C3(c), we prepare the system into a localized state in the mobility edge and insulator phases. The results in Figs. C3(e) and C3(f) suggest that the return probabilities also behave distinctly when ?Hfis set in phases with different transport nature. Putting together, we conclude that wherever the initial state is prepared,its return probability following quenches to different phases could help us to distinguish them dynamically. In experiments,the averaged spreading velocity and return probability of the wavepackets can thus be employed to detect phases with different localization nature in the dimerized NHQC.

    Fig. C3. Profiles of the initial states [in panels (a)–(c)] and their return probabilities following quenches [in panels (d)–(f)] in the dimerized NHQC.The states in panels(a)–(c)are chosen to be the ones with the largest imaginary parts of energy at the corresponding phase shift γi.System parameters for panels(a)–(c)are set as J=1,V =0.5,Δ =0.4,α =,and the length of the lattice is L=610 with PBC.In panels(d)–(f),the curves correspond to return probabilities of the initial states in panels(a)–(c)following quenches of γ from the corresponding γi to different values of γf applied at t=0.

    猜你喜歡
    龍文
    Relativistic Hartree–Fock model and its recent progress on the description of nuclear structure*
    Free-boundary plasma equilibria with toroidal plasma flows
    典故逸事龍文鞭影
    Experimental study of sheath potential coefficient in the J-TEXT tokamak
    Effect of edge turbulent transport on scrapeoff layer width on HL-2A tokamak
    紅軍第一位飛行員龍文光
    勤上光電收購(gòu)標(biāo)的經(jīng)營(yíng)亂象
    向北 向北 再向北
    明成祖朱棣:成就大業(yè)不忘恩師
    龍文未駕 鞭影猶存
    人人妻人人爽人人添夜夜欢视频| 国产精品国产高清国产av| 精品久久久久久电影网| 久久精品91无色码中文字幕| 午夜福利在线免费观看网站| 精品久久久久久,| 757午夜福利合集在线观看| 老司机午夜福利在线观看视频| 香蕉国产在线看| 亚洲片人在线观看| 久久久久国产精品人妻aⅴ院| 中文字幕最新亚洲高清| 麻豆国产av国片精品| 国产亚洲精品久久久久5区| 免费av毛片视频| 久久国产精品影院| 色播在线永久视频| 一级片'在线观看视频| 在线观看日韩欧美| 在线看a的网站| 成年人免费黄色播放视频| 日本五十路高清| 老熟妇乱子伦视频在线观看| 国产激情欧美一区二区| 黑人猛操日本美女一级片| 美女 人体艺术 gogo| 夜夜爽天天搞| 成人免费观看视频高清| 另类亚洲欧美激情| 国产精品免费一区二区三区在线| 国产xxxxx性猛交| 夜夜看夜夜爽夜夜摸 | 国产精品日韩av在线免费观看 | 欧美日韩亚洲高清精品| 黄色毛片三级朝国网站| 桃色一区二区三区在线观看| 大型av网站在线播放| 久久精品aⅴ一区二区三区四区| 亚洲精品一卡2卡三卡4卡5卡| 久久久久精品国产欧美久久久| 女警被强在线播放| av网站免费在线观看视频| 18禁国产床啪视频网站| 午夜精品国产一区二区电影| 午夜福利影视在线免费观看| 最近最新中文字幕大全电影3 | 国产国语露脸激情在线看| 女生性感内裤真人,穿戴方法视频| 精品一区二区三卡| 女同久久另类99精品国产91| 亚洲精品av麻豆狂野| 久久精品国产99精品国产亚洲性色 | 国产一区二区三区视频了| 国产精品久久电影中文字幕| 国产精品免费一区二区三区在线| 19禁男女啪啪无遮挡网站| 午夜福利欧美成人| 亚洲片人在线观看| 午夜91福利影院| 欧美中文综合在线视频| bbb黄色大片| 中文字幕人妻熟女乱码| 国产精品久久视频播放| 欧美日韩乱码在线| 视频区欧美日本亚洲| 美女高潮喷水抽搐中文字幕| 黄色女人牲交| 亚洲国产毛片av蜜桃av| 黑人巨大精品欧美一区二区蜜桃| 不卡一级毛片| 欧美日韩黄片免| 麻豆成人av在线观看| 女人爽到高潮嗷嗷叫在线视频| 欧美最黄视频在线播放免费 | 岛国在线观看网站| 欧美 亚洲 国产 日韩一| 国产99久久九九免费精品| 久久伊人香网站| 成人亚洲精品av一区二区 | 久久婷婷成人综合色麻豆| 亚洲欧美激情在线| 老熟妇乱子伦视频在线观看| 欧美日韩瑟瑟在线播放| 欧美 亚洲 国产 日韩一| 美国免费a级毛片| 美国免费a级毛片| 搡老熟女国产l中国老女人| 精品高清国产在线一区| 啦啦啦在线免费观看视频4| 黄片播放在线免费| 亚洲欧洲精品一区二区精品久久久| 中出人妻视频一区二区| 国产亚洲欧美98| 桃色一区二区三区在线观看| 99久久精品国产亚洲精品| 99热国产这里只有精品6| 免费搜索国产男女视频| 99国产精品一区二区蜜桃av| 在线看a的网站| 日韩人妻精品一区2区三区| 久久亚洲真实| 满18在线观看网站| 大陆偷拍与自拍| 久久久国产精品麻豆| 男女床上黄色一级片免费看| 两个人免费观看高清视频| 精品久久久久久,| 最好的美女福利视频网| 好男人电影高清在线观看| 麻豆久久精品国产亚洲av | 国产色视频综合| 亚洲色图 男人天堂 中文字幕| 久久精品国产亚洲av高清一级| 亚洲av美国av| 精品乱码久久久久久99久播| 久久天堂一区二区三区四区| 亚洲精品久久成人aⅴ小说| 精品久久久久久,| 巨乳人妻的诱惑在线观看| 免费看a级黄色片| 亚洲伊人色综图| 国产三级在线视频| 亚洲精品一区av在线观看| 神马国产精品三级电影在线观看 | 99国产精品99久久久久| 午夜免费鲁丝| 国产成人精品久久二区二区免费| 一级毛片高清免费大全| 亚洲中文日韩欧美视频| 91在线观看av| 中文字幕高清在线视频| 国产成人精品久久二区二区免费| 国产熟女xx| 免费日韩欧美在线观看| 亚洲欧美精品综合一区二区三区| 亚洲三区欧美一区| 成人国产一区最新在线观看| 亚洲国产毛片av蜜桃av| 男女做爰动态图高潮gif福利片 | 搡老熟女国产l中国老女人| 亚洲欧美一区二区三区久久| 一级a爱视频在线免费观看| 亚洲欧美日韩另类电影网站| 久久久久久免费高清国产稀缺| 少妇粗大呻吟视频| 国产高清视频在线播放一区| 97碰自拍视频| 久久人人97超碰香蕉20202| 亚洲全国av大片| 亚洲成人精品中文字幕电影 | 久久久久久大精品| 黄频高清免费视频| 国产成人精品无人区| 日韩人妻精品一区2区三区| 这个男人来自地球电影免费观看| 亚洲七黄色美女视频| 国产成人精品久久二区二区91| 在线看a的网站| 久久天躁狠狠躁夜夜2o2o| 18禁裸乳无遮挡免费网站照片 | 性少妇av在线| 亚洲一区高清亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 成人18禁在线播放| 黑人欧美特级aaaaaa片| 精品久久久久久久久久免费视频 | 中文亚洲av片在线观看爽| 91字幕亚洲| 精品国产一区二区三区四区第35| av视频免费观看在线观看| 1024香蕉在线观看| 成年人免费黄色播放视频| 久久中文字幕一级| 日韩av在线大香蕉| 久久精品亚洲av国产电影网| 99久久国产精品久久久| 久久精品国产99精品国产亚洲性色 | ponron亚洲| 亚洲伊人色综图| av欧美777| 国产欧美日韩精品亚洲av| 久久中文字幕人妻熟女| 嫁个100分男人电影在线观看| 真人一进一出gif抽搐免费| 亚洲免费av在线视频| 色尼玛亚洲综合影院| 亚洲欧美日韩高清在线视频| 日日摸夜夜添夜夜添小说| 成人永久免费在线观看视频| 中文亚洲av片在线观看爽| 纯流量卡能插随身wifi吗| 神马国产精品三级电影在线观看 | 一区二区三区激情视频| 欧美一区二区精品小视频在线| 亚洲自拍偷在线| 一区二区日韩欧美中文字幕| 国产精品二区激情视频| 国产精品九九99| 性色av乱码一区二区三区2| 日韩人妻精品一区2区三区| 99香蕉大伊视频| 两人在一起打扑克的视频| 高清在线国产一区| 欧美人与性动交α欧美软件| 九色亚洲精品在线播放| 国产一区二区三区视频了| 男女午夜视频在线观看| 久久亚洲精品不卡| 国产午夜精品久久久久久| 精品国产乱码久久久久久男人| 亚洲五月天丁香| 一进一出抽搐动态| 欧美丝袜亚洲另类 | 欧美一级毛片孕妇| 成人影院久久| 午夜91福利影院| 亚洲欧美日韩另类电影网站| 亚洲精品在线观看二区| 亚洲视频免费观看视频| 中文字幕人妻丝袜一区二区| 最近最新中文字幕大全免费视频| 麻豆国产av国片精品| e午夜精品久久久久久久| 久久人妻av系列| 又黄又粗又硬又大视频| 国产视频一区二区在线看| 波多野结衣高清无吗| 在线观看免费视频网站a站| 精品久久久久久久毛片微露脸| 在线看a的网站| 在线观看免费午夜福利视频| 成熟少妇高潮喷水视频| 真人做人爱边吃奶动态| 老司机在亚洲福利影院| 人人妻,人人澡人人爽秒播| 国产精品影院久久| 国产真人三级小视频在线观看| 99国产精品免费福利视频| 久久人妻熟女aⅴ| 免费在线观看黄色视频的| 高清av免费在线| 脱女人内裤的视频| 777久久人妻少妇嫩草av网站| 妹子高潮喷水视频| 亚洲精品一区av在线观看| 国产91精品成人一区二区三区| 国产伦人伦偷精品视频| 国产欧美日韩综合在线一区二区| www.熟女人妻精品国产| 香蕉久久夜色| 色哟哟哟哟哟哟| 精品一品国产午夜福利视频| 国产精品爽爽va在线观看网站 | 性欧美人与动物交配| 夜夜躁狠狠躁天天躁| 亚洲一区中文字幕在线| 一进一出抽搐gif免费好疼 | 可以免费在线观看a视频的电影网站| 一级毛片女人18水好多| 日韩精品中文字幕看吧| 99国产极品粉嫩在线观看| 中文字幕人妻丝袜一区二区| 在线看a的网站| 精品一区二区三区四区五区乱码| 久久热在线av| 在线观看66精品国产| 免费一级毛片在线播放高清视频 | 老熟妇仑乱视频hdxx| 日韩欧美一区二区三区在线观看| 99在线人妻在线中文字幕| 久久亚洲真实| 成人国语在线视频| 法律面前人人平等表现在哪些方面| 女人精品久久久久毛片| 国产亚洲精品久久久久5区| 中文亚洲av片在线观看爽| 亚洲成a人片在线一区二区| av天堂在线播放| 亚洲三区欧美一区| 婷婷精品国产亚洲av在线| 99热只有精品国产| 午夜福利在线观看吧| 色综合欧美亚洲国产小说| 国产成人精品久久二区二区免费| 久久精品91蜜桃| 亚洲国产看品久久| 香蕉丝袜av| 亚洲欧美激情综合另类| 中文字幕av电影在线播放| 亚洲免费av在线视频| 亚洲成国产人片在线观看| 美女福利国产在线| 久久婷婷成人综合色麻豆| a级片在线免费高清观看视频| 母亲3免费完整高清在线观看| 精品日产1卡2卡| 亚洲一区二区三区不卡视频| 亚洲 欧美 日韩 在线 免费| 国产高清激情床上av| 亚洲国产精品一区二区三区在线| 久久精品国产99精品国产亚洲性色 | 国产不卡一卡二| 国产有黄有色有爽视频| 国产区一区二久久| 午夜91福利影院| 乱人伦中国视频| 久久香蕉国产精品| 在线观看日韩欧美| 精品久久久久久久久久免费视频 | 国产野战对白在线观看| 国产欧美日韩精品亚洲av| 这个男人来自地球电影免费观看| 欧美另类亚洲清纯唯美| 久久久久久久久中文| 日本免费一区二区三区高清不卡 | 久久久久亚洲av毛片大全| 男女床上黄色一级片免费看| 精品免费久久久久久久清纯| e午夜精品久久久久久久| 国产激情欧美一区二区| 日韩中文字幕欧美一区二区| 男人舔女人下体高潮全视频| 高清黄色对白视频在线免费看| 国产成人一区二区三区免费视频网站| 怎么达到女性高潮| 精品国产亚洲在线| 岛国视频午夜一区免费看| 日韩人妻精品一区2区三区| 欧美一级毛片孕妇| 亚洲国产精品合色在线| 欧美av亚洲av综合av国产av| 18禁美女被吸乳视频| 女人精品久久久久毛片| 91九色精品人成在线观看| 老熟妇仑乱视频hdxx| 99精品欧美一区二区三区四区| 日韩高清综合在线| av天堂久久9| 97人妻天天添夜夜摸| 久久国产精品男人的天堂亚洲| 中文字幕另类日韩欧美亚洲嫩草| 精品一区二区三区四区五区乱码| 丝袜美足系列| 婷婷精品国产亚洲av在线| 少妇裸体淫交视频免费看高清 | 琪琪午夜伦伦电影理论片6080| 中文欧美无线码| 成人18禁高潮啪啪吃奶动态图| 一级黄色大片毛片| 他把我摸到了高潮在线观看| 18禁裸乳无遮挡免费网站照片 | 中文字幕精品免费在线观看视频| 午夜视频精品福利| 日本精品一区二区三区蜜桃| 一区二区三区精品91| av天堂久久9| 一级黄色大片毛片| 国产91精品成人一区二区三区| 91国产中文字幕| 99精品在免费线老司机午夜| 亚洲国产欧美日韩在线播放| 久久久水蜜桃国产精品网| 亚洲人成电影免费在线| 热re99久久国产66热| 午夜免费成人在线视频| 亚洲av第一区精品v没综合| 黄片大片在线免费观看| 制服人妻中文乱码| 亚洲人成网站在线播放欧美日韩| 久久天堂一区二区三区四区| 桃红色精品国产亚洲av| 最近最新中文字幕大全电影3 | 欧美一区二区精品小视频在线| 99精品欧美一区二区三区四区| 成人手机av| 久久国产精品人妻蜜桃| 动漫黄色视频在线观看| 亚洲熟妇中文字幕五十中出 | 精品少妇一区二区三区视频日本电影| 日韩欧美一区视频在线观看| 一边摸一边抽搐一进一小说| 黄色丝袜av网址大全| 热re99久久精品国产66热6| 亚洲精品美女久久久久99蜜臀| 黄色片一级片一级黄色片| 色综合婷婷激情| 久久国产乱子伦精品免费另类| 中文字幕人妻丝袜制服| 久久久久国产一级毛片高清牌| 黑人操中国人逼视频| 91成人精品电影| 母亲3免费完整高清在线观看| 国产高清国产精品国产三级| 精品国产亚洲在线| videosex国产| 亚洲第一欧美日韩一区二区三区| 手机成人av网站| 1024香蕉在线观看| 国产国语露脸激情在线看| 两性午夜刺激爽爽歪歪视频在线观看 | 最近最新中文字幕大全电影3 | 久久精品国产清高在天天线| 久久香蕉激情| 国产成+人综合+亚洲专区| 国产蜜桃级精品一区二区三区| 国产av又大| 国产亚洲精品综合一区在线观看 | 一区二区日韩欧美中文字幕| a级毛片在线看网站| 亚洲欧美日韩高清在线视频| 国产精品影院久久| cao死你这个sao货| 精品久久久久久久久久免费视频 | 亚洲精品粉嫩美女一区| aaaaa片日本免费| 九色亚洲精品在线播放| 国产一区二区三区综合在线观看| 热99re8久久精品国产| 亚洲一区二区三区不卡视频| 久久精品国产亚洲av香蕉五月| 香蕉国产在线看| 国产熟女xx| 91麻豆精品激情在线观看国产 | 激情视频va一区二区三区| 女生性感内裤真人,穿戴方法视频| 久久久国产成人精品二区 | 亚洲一码二码三码区别大吗| 国产片内射在线| 国产深夜福利视频在线观看| 69av精品久久久久久| 国产日韩一区二区三区精品不卡| 日韩欧美国产一区二区入口| 久久人妻福利社区极品人妻图片| 精品久久久久久成人av| 久久精品91蜜桃| 丝袜美腿诱惑在线| 最新美女视频免费是黄的| 成人av一区二区三区在线看| av国产精品久久久久影院| a级毛片黄视频| 88av欧美| 国产区一区二久久| 日本黄色视频三级网站网址| 亚洲狠狠婷婷综合久久图片| 成人亚洲精品av一区二区 | а√天堂www在线а√下载| 制服人妻中文乱码| 91字幕亚洲| 欧美激情高清一区二区三区| 亚洲精品国产色婷婷电影| 日韩免费av在线播放| 老鸭窝网址在线观看| 国产99久久九九免费精品| 日本五十路高清| 丝袜人妻中文字幕| 亚洲精品av麻豆狂野| 久久人人精品亚洲av| 视频在线观看一区二区三区| 久久 成人 亚洲| 国产真人三级小视频在线观看| 黄片大片在线免费观看| 精品国内亚洲2022精品成人| 天天躁狠狠躁夜夜躁狠狠躁| tocl精华| 亚洲午夜精品一区,二区,三区| 久久久水蜜桃国产精品网| 欧美中文综合在线视频| 一进一出抽搐gif免费好疼 | 丝袜美足系列| 另类亚洲欧美激情| 亚洲一区高清亚洲精品| 亚洲精品国产区一区二| 性色av乱码一区二区三区2| 看免费av毛片| 国产精品香港三级国产av潘金莲| 操美女的视频在线观看| 国产精品1区2区在线观看.| 亚洲精品美女久久av网站| 校园春色视频在线观看| 久久久久久人人人人人| 人成视频在线观看免费观看| 亚洲午夜理论影院| 99国产精品免费福利视频| 久久久国产一区二区| 成人国产一区最新在线观看| 最近最新免费中文字幕在线| 真人做人爱边吃奶动态| 亚洲视频免费观看视频| 亚洲一区高清亚洲精品| 身体一侧抽搐| 中亚洲国语对白在线视频| 国产免费男女视频| 老汉色av国产亚洲站长工具| 国产激情久久老熟女| 日韩成人在线观看一区二区三区| 麻豆国产av国片精品| 欧美日韩av久久| 美女高潮到喷水免费观看| www.熟女人妻精品国产| 午夜免费成人在线视频| 99热国产这里只有精品6| 黑人巨大精品欧美一区二区mp4| 18禁美女被吸乳视频| 9191精品国产免费久久| 免费在线观看完整版高清| 高清av免费在线| 久久欧美精品欧美久久欧美| 欧美日韩视频精品一区| 精品福利观看| 国内毛片毛片毛片毛片毛片| 啦啦啦在线免费观看视频4| 精品乱码久久久久久99久播| 国产激情久久老熟女| 日本五十路高清| 亚洲色图 男人天堂 中文字幕| 亚洲精品成人av观看孕妇| 国产成人精品久久二区二区91| 91在线观看av| 五月开心婷婷网| 亚洲在线自拍视频| 两个人免费观看高清视频| 操美女的视频在线观看| 亚洲精品中文字幕在线视频| 夜夜爽天天搞| 欧美日韩瑟瑟在线播放| 淫秽高清视频在线观看| 国产精品一区二区免费欧美| 天天影视国产精品| 夜夜夜夜夜久久久久| 纯流量卡能插随身wifi吗| 一边摸一边抽搐一进一小说| 一本大道久久a久久精品| 99精国产麻豆久久婷婷| 欧美大码av| 婷婷精品国产亚洲av在线| 精品熟女少妇八av免费久了| 超碰97精品在线观看| 亚洲五月婷婷丁香| 欧美精品亚洲一区二区| 午夜亚洲福利在线播放| 91老司机精品| 亚洲国产看品久久| 久久香蕉激情| 老熟妇乱子伦视频在线观看| www.999成人在线观看| 老司机午夜福利在线观看视频| 亚洲欧美精品综合一区二区三区| 成人亚洲精品一区在线观看| 黑人猛操日本美女一级片| 午夜影院日韩av| 无人区码免费观看不卡| 精品国产超薄肉色丝袜足j| 日日干狠狠操夜夜爽| 人人澡人人妻人| 99久久人妻综合| 久久精品亚洲av国产电影网| 欧美精品啪啪一区二区三区| 久久精品人人爽人人爽视色| 看免费av毛片| 一区二区日韩欧美中文字幕| 狂野欧美激情性xxxx| 97碰自拍视频| 午夜免费激情av| 99国产精品一区二区三区| 久久精品人人爽人人爽视色| 人人妻人人添人人爽欧美一区卜| 男人舔女人下体高潮全视频| 视频区图区小说| 夜夜看夜夜爽夜夜摸 | 精品福利永久在线观看| 咕卡用的链子| 90打野战视频偷拍视频| 亚洲精品av麻豆狂野| 亚洲成av片中文字幕在线观看| 午夜成年电影在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 99久久精品国产亚洲精品| 久久精品91蜜桃| 狠狠狠狠99中文字幕| 日日爽夜夜爽网站| 免费在线观看完整版高清| 国产激情久久老熟女| 日本欧美视频一区| 亚洲免费av在线视频| 亚洲专区中文字幕在线| 久久精品成人免费网站| 9热在线视频观看99| 久久久久久大精品| 黄色a级毛片大全视频| 国产精品免费一区二区三区在线| www.精华液| 在线观看一区二区三区激情| 操出白浆在线播放| 色哟哟哟哟哟哟| 欧美精品一区二区免费开放| 一级a爱视频在线免费观看| 国产成人精品久久二区二区91| 高清欧美精品videossex| 天天添夜夜摸| 日日摸夜夜添夜夜添小说| 国产91精品成人一区二区三区| 日本三级黄在线观看| 一a级毛片在线观看| 一级片'在线观看视频| 国产一区在线观看成人免费| 亚洲欧美激情综合另类| 亚洲专区中文字幕在线| 真人一进一出gif抽搐免费| 亚洲中文字幕日韩| 一个人观看的视频www高清免费观看 | 韩国av一区二区三区四区| 天堂动漫精品| 在线播放国产精品三级| 亚洲美女黄片视频| 中文字幕色久视频|