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

    Nonlinear dynamical stability of gap solitons in Bose–Einstein condensate loaded in a deformed honeycomb optical lattice?

    2021-12-22 06:50:06HongjuanMeng蒙紅娟YushanZhou周玉珊XuepingRen任雪平XiaohuanWan萬曉歡JuanZhang張娟JingWang王靜XiaobeiFan樊小貝WenyuanWang王文元andYurenShi石玉仁
    Chinese Physics B 2021年12期
    關(guān)鍵詞:張娟王靜小貝

    Hongjuan Meng(蒙紅娟) Yushan Zhou(周玉珊) Xueping Ren(任雪平)Xiaohuan Wan(萬曉歡) Juan Zhang(張娟) Jing Wang(王靜)Xiaobei Fan(樊小貝) Wenyuan Wang(王文元) and Yuren Shi(石玉仁)

    1College of Physics and Electronic Engineering,Northwest Normal University,Lanzhou 730070,China

    2Key Laboratory of Atomic and Molecular Physics&Functional Materials of Gansu Province,Northwest Normal University,Lanzhou 730070,China

    Keywords: gap soliton,Bose–Einstein condensate,deformed honeycomb optical lattice

    1. Introduction

    Bose–Einstein condensate (BEC) in optical lattices has generated considerable interesting research for studying many quantum phenomena in periodic potentials.[1–4]It arises in many applications ranging such as quantum simulators[5–7]and quantum computation.[8]The cubic lattice is a simple representative example of an optical lattice. Experiments with ultracold atoms have utilized exotic lattices such as triangular,[9]honeycomb,[10]and Kagome lattices.[11]The honeycomb optical lattices can lead to significantly different results due to the underlying symmetries. A special feature of a perfect honeycomb lattice is the two-lowest-band-formed Dirac points.These unique characteristics could provide new opportunities to investigate many interesting physical phenomena in many systems such as monolayer graphene,[12,13]electron gases in semiconductors,[14]coupled photonic lattices,[15,16]and photonic waveguides.[17]

    On the other hand,the interplay between lattice periodicity and atomic interaction brings a large number of interesting nonlinear effects.[18–20]In many cases, periodicity relates to the band structure of the dispersion relation and Bloch waves.Meanwhile,the nonlinearity can induce self-phase modulation and can also couple waves between different modes or Bloch bands. The existence of the nonlinearly localized structures that reside in the linear energy gaps has attracted considerable interest in recent years.[21–23]They are sometimes referred to as gap solitons,which are often considered theoretically as bifurcations from the Bloch-band edges into the band gaps.[24,25]The general properties of gap solitons have been intensively studied in many different nonlinear periodic systems, examples include but are not limited to optical waveguide arrays and photonic lattices,[26,27]semiconductor microcavities,[28]and BEC in optical lattices.[1,23,29,30]

    Gap solitons in honeycomb optical lattice have been studied both theoretically and experimentally.[17]The authors of Ref. [31] demonstrated the gap solitons residing in the gap between the second and the third bands in honeycomb photonic lattices. In Ref. [32], the authors showed solitons and necklaces in the first band gap of a two-dimensional optically induced honeycomb defocusing photonic lattice. In Ref.[33],the authors studied gap solitons in honeycomb dynamical lattices with the cubic nonlinearity. However,there is still great significance to study gap solitons of BEC in honeycomb optical lattice,especially for the interplay between the honeycomb optical lattice and the nonlinearity affects on the existence and dynamical stability of multipole gap solitons.

    Here, we investigate gap solitons of BEC in deformed honeycomb optical lattices with atomic interaction and study their nonlinear dynamical stability. Deformed honeycomb lattices are promising platforms. With the deformation of the honeycomb lattice,the symmetries can be broken,then a gap opens between the first two bands,and the lattice can support a gap soliton.[34,35]The existence of the gap soliton provides a way to detect the symmetry breaking. The deformed honeycomb lattices have been studied extensively,such as nonlinear wave packets,[36]Klein tunneling,[37]fractional topological quantum states,[38]nondegenerate chiral phonons,[39]and topological valley transport.[40]

    In this paper, we show the linear band-gap structure of the deformed honeycomb optical lattices and the existence of the gap solitons. We find that the system admits dipole gap solitons both in the first gap and in the semi-infinite gap, depending on the repulsive or attractive properties of atomic interactions. The dipole gap soliton can have their bright solitary structures be in-phase or out-of-phase. Then,we analysis the nonlinear dynamical stability of these dipole gap solitons by using direct simulations of the Gross–Pitaevskii(GP)equation.Finally,we derive the envelope equations of Bloch waves and analyze these gap solitons.

    The paper is organized as follows: In Section 2,we introduce the model of BEC in deformed honeycomb optical lattices. In Section 3,we show the linear band-gap structure and the existence of the gap solitons. Following this,the analysis of nonlinear dynamical stability of the gap solitons is given in Section 4.Finally,we give summary and conclusion in Section 5.

    2. Model of BEC in a honeycomb optical lattice

    We consider that BEC is loaded in a honeycomb optical lattice. The honeycomb optical lattice can be experimentally generated by superposing three coplanar traveling laser beams.[41,42]The three coplanar traveling laser beams have the same angular frequencyωL=ck0withk0being the wave vector of the radiation. Then, the lattice potential is denoted byV(r)with[35]

    Fig. 1. Contour plot of the potential energy for (a) η =1 (the optical lattice has a perfect standard graphene-like honeycomb),and for(b)η=0.6(the optical lattice becomes a deformed honeycomb).The white dotted line in the figure is a schematic diagram of the corresponding lattice sites.

    3. Gap solitons and theoretical analysis

    3.1. Band-gap structure

    In this subsection, we show the band-gap structure and the existence of the gap solitons. We first investigate the linear band-gap structure of the deformed honeycomb optical lattices. Let us consider quantum mechanical of BEC in such periodic potentials. Gap solitons are stationary solutions to Eq. (2), we can apply Bloch’s theorem, which suggests that the Bloch waves are the nonlinear eigenstates of the following formΨ(r,t)=φk(r)exp(i(k·r ?μt)).[44]Hereμis a real chemical potential, andφkhas the same period as the potential.Then we can obtain the following equation for each Bloch wave stateφk:The set of eigenvaluesμ(k)then forms Bloch bands.

    Since gap solitons only exist in the linear energy gaps between the Bloch bands,it is important to first identify the positions of these gaps. To do this,we can neglect the nonlinear terms, i.e.,g=0, and calculate the linear Bloch spectrum by numerical diagonalization of Eq.(3). The band gap spectrum generated by the solution of the linearized version of the equations is shown in Fig. 2. Due to the expected periodicity of the solution to Eq. (3), there are several numerical methods to find the Bloch waveφk.[45–47]We use Fourier collocation method[48]for numerically computing these Bloch bands.

    In Fig.2,we only show the lowest five Bloch bands. One can see that a band gap[49]is separated by the two Bloch bands.The lowest gap below the lowest Bloch band is called the semiinfinite gap and the higher gaps are called the first band gap,the second band gap,and so on.

    Fig.2. (a)The linear Bloch band gap structure of BEC in a deformed honeycomb optical lattice as a function of k(left). (b)Axial side view of the linear Bloch band gap structure as a function of kx,shaded areas correspond to the linear bands.

    The band gap structure of BEC in honeycomb optical lattices depends not only on the depth of the optical lattices but also on the degree of deformed honeycomb. Figure 2 shows how the band gap structure of BEC in honeycomb optical lattices is affected by lattice depthV0and deformation parameterη.

    Figure 3(a) shows the energies of different Bloch bands overlaping for a suffciiently small lattice depthV0,in this case no other band gap appears except the semi-infniite gap. With increasing lattice depthV0, the window of band gap will become more and more obvious.Figure 3(b)shows that the band gap structure depends on the deformation of the honeycomb.Forηchanging in the range<η ≤1,the smaller theη,the more deformed the honeycomb lattice. With increasingη,the frist gap becomes narrower and narrower and eventually disappears. Figure 3(b) shows an interesting phenomenon, that is,the more deformed the honeycomb lattice,the more prominent the frist gap. For a standard honeycomb lattice, the frist gap does not exist.Therefore,in order to obtain gap solitons in the frist gap,one must choose a deformed non-standard honeycomb lattice.In this article,we choose a deformed honeycomb lattice forη=0.6 as an example.

    Fig. 3. The band gap structure of BEC in honeycomb optical lattices with increasing lattice depth(a)and deformation parameter η (b).

    3.2. The gap solitons

    The interplay between spatial periodicity and nonlinear dynamical evolution harbors a large number of interesting physical many-body effects. Gap solitons is the existence of unique soliton solutions that reside in the linear energy gaps between the Bloch bands. Gap solitons spatially localized in gap modes may appear as a result of the balance between the dispersion and the nonlinearity. For BECs in an optical lattice,the hoppings between different lattice sites play the role of dispersion,and the atom-atom interaction brings the nonlinearity.

    In what follows, we consider the nonlinearity brought from the atom-atom interaction and stationary solutions in the linear Bloch bands gaps. We look for gap solitons by numerical methods such as the Newton conjugate-gradient (NCG)method,[50]which is based on iterating certain time evolution equations associated with the linear stability eigenvalue problem and is one of the powerful numerical techniques to find the solitary wave solutions of a nonlinear evolution equation.The methods are applied to compute both the ground states and excited states in a large number of physical systems with and without periodic potentials. We take the Gaussian profile as the initial condition for iterating. One can adjust the amplitude,width and center position of the initial Gaussian profiles to get the wanted gap solitons. For example, in order to obtain unipole gap solitons, the initial Gaussian wave packet is mainly distributed in a certain potential well.In order to obtain the dipole gap solitons,the initial two Gaussian wave packets are mainly distributed in two adjacent potential wells.

    After numerical techniques mentioned above,we find the unipole gap solitons in the semi-infinite gap. Figure 4(a)shows gap solitons in the semi-infinite gap withμ=4.8 andg=?1 as an example. Figure 4(b)shows gap solitons in the semi-infinite gap withμ=5.01 andg=?1 as an example.The gap solitons in the first gap can be obtained withg>0.

    Fig.4.Gap solitons in the semi-infinite gap with attractive atomic interaction g=?1.0:(a)forμ=4.8,(b)forμ=5.01.The white dotted line in the figure is a schematic diagram of the corresponding lattice sites.

    In order to study the properties of gap solitons more clearly,we pay attention to effect of nonlinear eigenenergyμon the amplitude of the gap solitons and the total power of gap solitons. In order to see this, we define the amplitudeAand the total powerPof the gap solitons given as follows:

    In Fig. 5, we show amplitudeAand powerPof the gap solitons for nonlinear eigenenergyμwith respect to on-site interactiongin the range of the band gap. Here, we modulate the nonlinear eigenenergyμfrom semi-infinite band gap to first band gap. Here,we takeg=?1 andg=1 as examples for gap solitons in semi-infinite band gap and first band gap,respectively. When the nonlinear eigenenergyμmoves toward the band,numerical results indicate that the amplitude of gap solitons both in semi-infinite gap and in first gap approach to infinitesimal. However, when the nonlinear eigenenergyμmoves toward the band,the power of gap solitons both in semiinfinite gap and in first gap approach to infinity.

    There is a very interesting phenomenon shown in Fig. 5 for the powerPof the gap solitons versus eigenenergyμin the range of the band gap. When the nonlinear eigenenergyμmoves toward the band, the powerPof gap solitons both in semi-infinite gap and in first gap decreases slowly first, then increases gradually. This means that the relationship between the amplitudeAand powerPof the gap solitons is inconsistent with the eigenenergyμ,they move in step whenμfarther away from the band,and they move out of step whenμcloser to the band. Thus, they are bifurcated from the Bloch band edges.

    Fig.5.The amplitude A and power P of the gap solitons versus eigenenergy μ in the range of the band gap. Shaded areas correspond to the linear bands.

    It is worth emphasizing that the appearance of gap solitons both in the first gap and in the semi-infinite gap depends on the repulsive or attractive properties of atomic interactions.The first gap solitons appear only in repulsive atomic interactions regime (g> 0). On the other hand, semi-infinite gap solitons appear only in attractive atomic interaction regime(g<0). The gap solitons are intrinsically related to the extended Bloch waves,[51–54]which can be understood from this perspective of the composition relation between Bloch waves and gap solitons.[52,55]With this composition relation, many conclusions can be drawn for gap solitons from Bloch waves in nonlinear periodic systems. The reason why the appearance of gap solitons in different gaps depends on the properties of atomic interactions is that the atomic interactions of different properties have different effects on the Bloch waves.[52,55]

    Next,we show dipole gap solitons in the semi-infinite gap withμ= 4.8 andg=?1 in Fig. 6 as an example, (a) for the out-of-phase gap solitons, (b) for the in-phase gap solitons. The dipole gap solitons in the first gap which are similar to solitons of the semi-infinite gap can also be obtained withg>0. In order to study the properties of dipole gap solitons more clearly,we pay attention to effect of nonlinear eigenenergyμon the total power of dipole gap solitons in Fig.7. The amplitudes of the dipole gap solitons in the first gap and in the semi-infinite gap are similar to unipole gap solitons. In Fig.7(a),when the nonlinear eigenenergyμmoves toward the band, the powerPof gap solitons of in-phase in the semiinfinite gap and out-of-phase in the first gap decreases slowly first, then increases gradually. We call gap solitons such that they are bifurcated from the Bloch band edges. In Fig. 7(b),when the nonlinear eigenenergyμmoves toward the band,the powerPof gap solitons of out-of-phase in the semi-infinite gap and in-phase in the first gap decreases slowly. We call gap solitons such that they are not bifurcated from the Bloch band edges.

    Fig. 6. Dipole gap solitons in the semi-infinite gap with attractive atomic interaction g=?1.0 andμ =4.8(a)for out-of-phase gap solitons and(b)for in-phase gap solitons.

    Fig.7. The power P of the dipole gap solitons versus eigenenergyμ in the range of the band gap. Shaded areas correspond to the linear bands.

    3.3. Theoretical analysis of gap solitons

    In this subsection, we derive the envelope equations of small-amplitude Bloch-wave packets from edges of Bloch bands. Considering the stationary solutionsψ(r) of Eq. (2)with the formΨ(r,t)=ψ(r)e?iμt. When the solutionψ(r)is infinitesimal,this solution is a linear superposition of these two Bloch modes. Whenψ(r) is small but not infinitesimal,this solution then becomes a combination of these two Blochwave packets and can be expanded into a multiscale perturbation series:

    Atε3,the equation forψ2is

    The solution to Eq.(7)has the following form:

    Equation (10) is the leading-order two-Bloch-wave packets,whereX=εxandY=εyare slow spatial variables of envelope functionsA1andA2.

    The solution to Eq.(8)is

    Substituting theψ0,ψ1into Eq. (9) and using the Fredholm condition (which requires being orthogonal to the homogeneous solutionsp1(x)p2(y)andp1(y)p2(x)),we can obtain

    Here,the superscript?represents complex conjugation,Assume that the modulation wave is a Gaussian profile one,then it can be expressed as

    whereAandW(W>0) are the amplitude and width of the solitons, respectively. Note that Eq. (16) is not an exact solution to Eq. (15). Substituting Eq. (16) into Eq. (15) yields the residual error. We minimize the magnitude of the residual error and then have the equation. We solve the equation in a semi-numerical manner,which can be performed easily in Mathematica software by the NSolve command. The results read

    It can be seen from Eq.(17)that the amplitudeAof the soliton depends on the parametersμ1,g,andα,whereas the widthWof the soliton depends on the parametersD1andμ1. Equation(17)also gives the condition for the existence of gap solitons,that is,μ1andgmust have the same sign,butμ1andD1have the opposite signs.In the semi-infinite gap,ifD1>0,one will haveμ1<0 andg<0,with the increase of|μ1|,the value ofμdecreases, and the value ofμis farther away from the band. The amplitude of the soliton increases,which is consistent with the previous numerical results shown in Fig.5. For the first gap,ifD1<0,one will haveμ1>0 andg>0. With the increase of|μ1|,the value ofμincreases,and the value ofμis also farther away from the first band. The amplitude of the soliton increases, which is also consistent with the previous numerical results shown in Fig.5.

    4. Linear stabilities and nonlinear dynamical stabilities of gap solitons

    Until now, we have shown the existence of unipole and dipole gap solitons of BEC in a deformed honeycomb optical lattice. One unique feature of BEC in optical lattices is the occurrence of instability. To do this,one can analyze the nonlinear and the linear stability of these gap solitons using direct simulations of the GP Eq.(2)and its linearized equation.

    4.1. Linear stabilities of gap solitons

    The linear stability (so-called Bogolyubov-de Gennes or BdG) analysis is employed in order to consider the fate of small amplitude perturbations and the potential robustness of the solutions. The linear stability of these gap solitons can be modified state and inserted into the time-dependent GP equation. By keeping only the linear term in the perturbation, one can describe the time evolution of the small perturbation,[46,56–58]which yields the linear eigenvalue problem.If the corresponding eigenvalues are real,the perturbation will grow exponentially,which leads to dynamic instability.

    Now we study the linear stability of these gap solitons by numerical methods. We perturb these solitons as

    The above linear eigenvalue problem (19) can be solved numerically by the finite difference method or the Fourier collocation method.[48]To illustrate,we consider unipole solitons in the semi-infinite gap and in the first gap,as shown in Fig.8.Also,we consider dipole gap solitons in the semi-infinite gap as an example,as shown in Fig.9.Stability spectra for unipole solitons are shown in Fig. 8. It is seen that unipole gap solitons far away from the band have stable eigenvalues and thus are linearly stable in Figs. 8(a) and 8(d). However, unipole gap solitons near the band have unstable eigenvalue and thus are linearly unstable in Figs. 8(b) and 8(c). Figure 9 shows dipole gap solitons of stability spectra in the semi-infinite gap.It is seen that in-phase dipole gap solitons are all linearly unstable in Figs.9(a)and 9(b). Figures 9(c)and 9(d)show that out-of-phase dipole gap solitons far away from the band are linearly stable,and near the band are linearly unstable.

    Fig.8. The linear stability analysis of unipole solitons: (a)and(b)for unipole solitons in the semi-infinite gap;(c)and(d)for unipole solitons in the first gap.

    Fig. 9. The linear stability analysis of the dipole gap solitons in the semi-infinite gap:(a)and(b)for in-phase solitons bifurcated from band edges; (c) and (d) for out-of-phase solitons not bifurcated from band edges.

    4.2. Nonlinear dynamical stabilities of unipole gap solitons

    In the above section, we have analyzed the linear stability problem for unipole and dipole gap solitons. However,the important issue is the nonlinear dynamical stability of these unipole gap solitons and dipole gap solitons,which can be induced by the presence of nonlinear interactions and periodic potential. Dynamical instability implies that small deviations from the stationary solution grow exponentially in time. In order to examine the nonlinear dynamical stability of these unipole gap solitons and dipole gap solitons found above,we directly compute Eq. (2) over a long time. The initial conditions are taken to be the gap soliton with 1%perturbation in the amplitude and phase. We consider a stable gap soliton satisfying the two conditions as follows: it should preserve its peak amplitude,and it should preserve its position on the lattices.

    Fig.10. Nonlinear evolution of unipole gap solitons in the semi-infinite gap with μ =4.8 (corresponding to the point a in Fig. 5). The initial conditions were taken with 1%perturbation in the amplitude and phase to trigger potential instabilities: (a)t=0,(b)t=100.

    We first consider the nonlinear dynamical stabilities of the unipole gap solitons, and takeμ=4.8 in semi-finite gap as an example (corresponding to the pointain Fig. 5). In Fig.10,we show time evolution of unipole gap solitons in the semi-infinite gap fort=0 and 100. It can be seen from this figure that the amplitude and position on the lattices of the unipole gap solitons are preserved. The gap soliton is nonlinear dynamical stable. This is confirmed by the linear stability analysis in Fig.8(a).

    Then, the nonlinear dynamical stability properties of the unipole gap solitons withμ=5.01 is examined, as shown in Figs.11(a)and 11(b)fort=0 and 100,respectively. We cam see the maximum amplitude and the location of the center of mass versus the evolution timet. During the nonlinear dynamical evolution, the peak amplitude undergoes finite-time collapse, its position on the lattice drift-unstable solitons is typically characterized by“humps”that drift from lattice maxima toward nearby minima. Obviously, the unipole gap soliton shows nonlinearly dynamical unstable,which is also confirmed by the linear stability analysis in Fig.8(b).

    Fig.11. Nonlinear evolution of unipole gap solitons in the semi-infinite gap with μ =5.01(corresponding to the point b in Fig.5). The initial conditions are taken with 1% perturbation in the amplitude and phase to trigger potential instabilities: (a)t=0 and(b)t=100.

    In order to examine the nonlinear dynamical stability of unipole gap solitons in the first gap,we also directly compute Eq. (2) over a long time. The initial conditions were taken to be a unipole gap soliton with 1%perturbation in amplitude and phase in Fig.12. It can be seen that the unipole gap soliton shows nonlinearly dynamical unstable and (b) nonlinear dynamical stable. These are also confirmed by the linear stability analysis in Figs.8(c)and 8(d),respectively.

    Fig. 12. Nonlinear evolution of the unipole gap solitons in the first gap with repulsive atomic interactions g=1. The initial conditions are taken with 1%perturbation in the amplitude and phase to trigger potential instabilities:(a)μ=5.36 and(b)μ=5.56.The unipole gap soliton corresponds to the points c and d in Fig.5,respectively.

    It is well known that the power versusμ, i.e., dP/dμ,plays an important role in determining the stability properties of the solitons. This connection between the slope of the power curve and the linear stability is the so-called Vakhitov–Kolokolov stability criterion.[59]This criterion states that with the continuous change ofμ, a pair of real eigenvalues move towards the origin, collide there,and then bifurcate along the imaginary axis. The power minimum point turns out to be precisely the bifurcation point,which indeed suggests that linear stability of these solitary waves is intimately related to the sign of the slope of the power curve.[48]In a large number of simulations,we have ever plotted the total power versusμfor unipole gap solitons corresponding to the semi-infinite gap and the first gap,and have observed two recurring features:

    (a)In the semi-infinite gap,the unipole gap soliton is nonlinear dynamical stable only if its power decreases with increasingμ, i.e., dP/dμ<0(as the nonlinear dynamical evolution of the pointain Fig.5).On the contrary,the unipole gap soliton will be nonlinear dynamical unstable only if its power increases with increasingμ,i.e.,dP/dμ>0(as the nonlinear dynamical evolution of the pointbin Fig.5).

    (b) In the first gap, the unipole gap soliton will be nonlinearly dynamically stable only if its power increases with increasingμ,i.e.,dP/dμ>0(as the nonlinear dynamical evolution of the pointdin Fig.5). Meanwhile,the unipole gap soliton will be unstable only if its power decreases with increasingμ,i.e.,dP/dμ>0(as the nonlinear dynamical evolution of the pointcin Fig.5).

    The nonlinear dynamical stabilities of unipole gap solitons depend on dP/dμ,which are opposite in semi-infinite gap and first gap. The reason for this phenomenon is that there are different requirements for the properties of atomic interaction when solitons occur in two gaps.

    4.3. Nonlinear dynamical stabilities of dipole gap solitons

    In the semi-infinite gap, the out-of-phase and in-phase dipole gap solitons exist (shown in Fig. 6) for the attractive atomic interactions regime. Now,we investigate the nonlinear stability properties of the lattice soliton. The in-phase dipole gap soliton corresponds to Fig.6(b). The dynamics of nonlinear evolution of the in-phase dipole gap soliton withμ=4.8 is shown in Fig. 13. Figures 13(a) and 13(b) depict the density profiles of the in-phase dipole gap soliton fort=0 and 100, respectively. In fact, in the semi-infinite gap, for all inphase dipole gap solitons their nonlinear dynamic evolution will show a similar phenomenon nonlinearly dynamical unstable,as seen in Fig.7(a).In the first gap,the in-phase dipole gap solitons satisfy the condition dP/dμ>0. Solutions near the band have unstable eigenvalues and thus are linearly unstable.However, solitons far away from the band have no unstable eigenvalues and thus are linearly stable.

    Fig. 13. Nonlinear evolution of the in-phase dipole gap soliton in the semi-infinite gap with repulsive atomic interactions g=?1. The inphase dipole gap soliton corresponds to Fig.6(b). The initial conditions are taken with 1% perturbation in the amplitude and phase to trigger potential instabilities.

    Fig. 14. Nonlinear evolution of the out-of-phase dipole gap soliton in the semi-infinite gap with repulsive atomic interactions g=?1 andμ =4.8. The out-of-phase dipole gap soliton corresponds to Fig.6(a).The initial conditions are taken with 1% perturbation in the amplitude and phase to trigger potential instabilities: (a)t=0 and(b)t=100.

    Figure 14 shows nonlinear evolution of the out-of-phase dipole gap soliton in the semi-infinite gap atμ=4.8, corresponds to the soliton in Fig.6(a). Figures 14(a)and 14(b)fort=0 and 100 respectively. It can be seen from this figure that the amplitude and position on the lattices of the out-of-phase dipole gap solitons are preserved.The out-of-phase dipole gap soliton is nonlinear dynamical stable. Figure 15 shows nonlinear evolution of the out-of-phase dipole gap soliton corresponding to the semi-infinite gap atμ=5.0.The dipole gap soliton also corresponds to the Fig. 6(a). Figures 15(a) and 15(b)fort=0,100 respectively. It can be seen from this figure that the maximum amplitude and the location of the center of mass versus the evolution timet. During the nonlinear dynamical evolution, the peak amplitude undergoing finite-time collapse, its position on the lattice drift-unstable solitons are typically characterized by“humps”that drift from lattice maxima toward nearby minima. Obviously the dipole gap soliton shows nonlinearly dynamical unstable. That is, in the semiinfinite gap,when solutions near the band have unstable eigenvalue and thus are linearly unstable.Solitons far away from the band have no unstable eigenvalues and thus are linearly stable.However,for all out-of-phase dipole gap solitons show nonlinearly dynamical unstable in the first gap. It can be seen from Fig.7.

    Fig. 15. Time evolution of the out-of-phase dipole gap solitons in the semi-infinite gap withμ=5.0,g=?1. The initial conditions are taken with 1% perturbation in the amplitude and phase to trigger potential instabilities: (a)t=0 and(b)t=100.

    5. Summary

    We have investigated the existence and dynamical stability of multipole gap solitons in BEC loaded in a deformed honeycomb optical lattice. Honeycomb lattices possess a unique band structure, the first and second bands intersect at a set of so-called Dirac points. We find that,when the lattice is gradually deformed, at some strong enough deformation, a gap forms between the first and second bands. We show that the gap in such a deformed lattice can form gap soliton. Both the unipole gap soliton and dipole gap soliton are found. The dipoles can have their bright solitary structures being in-phase or out-of-phase. The linear stabilities and nonlinear stabilities of these gap solitons are investigated.

    The stabilities of unipole gap solitons depend on dP/dμ.In the semi-infinite gap,the unipole gap soliton will be stable if its power and eigenenergy satisfy dP/dμ<0, whereas the unipole gap soliton will be unstable if its power and eigenenergy satisfy dP/dμ>0.In the first gap,the unipole gap soliton will be stable if its power and eigenenergy satisfy dP/dμ>0,

    whereas the unipole gap soliton will be unstable if its power and eigenenergy satisfy dP/dμ>0. The stabilities of dipole gap solitons are independent of the dP/dμ. For the in-phase dipole gap soliton in semi-infinite gap and the out-of-phase dipole gap solitons in first gap,they are always unstable. For the out-of-phase dipole gap soliton in semi-infinite gap and the in-phase dipole gap solitons in first gap,they are stable far away from the band,while unstable near the band.

    Our work provides a qualitative and quantitative guideline towards the more detailed understanding of the possible dynamical features and of the interactions between the coherent structures of strongly correlated states of BEC in 2D optical lattices.

    猜你喜歡
    張娟王靜小貝
    貝小貝哭了
    Clinical study of warm needling moxibustion combined with entecavir in the treatment of compensated cirrhosis due to chronic hepatitis B
    夢寐以求的生日禮物
    The Management Methods And Thinking Of Personnel Files
    客聯(lián)(2021年9期)2021-11-07 19:21:33
    The Development of Contemporary Oil Painting Art
    青年生活(2019年16期)2019-10-21 01:46:49
    小貝趕集記
    兒童時代(2017年5期)2017-07-18 11:55:35
    Income Inequality in Developing Countries
    商情(2017年17期)2017-06-10 12:27:58
    Let it Go隨它吧
    小貝螺的大世界
    小貝當(dāng)刮臉
    成人精品一区二区免费| 国产精品九九99| 精品第一国产精品| 国产欧美日韩一区二区三| 国产成人免费观看mmmm| 国产精品九九99| 日韩制服丝袜自拍偷拍| 淫妇啪啪啪对白视频| 在线观看免费日韩欧美大片| 欧美在线黄色| 亚洲专区字幕在线| 久久天躁狠狠躁夜夜2o2o| 国产高清激情床上av| 亚洲一区中文字幕在线| 美女扒开内裤让男人捅视频| 狠狠狠狠99中文字幕| 日本五十路高清| 在线观看舔阴道视频| 18在线观看网站| 少妇被粗大的猛进出69影院| 午夜激情av网站| 国产激情久久老熟女| 国产一区二区三区综合在线观看| 少妇粗大呻吟视频| 制服人妻中文乱码| 中文字幕人妻丝袜制服| 搡老岳熟女国产| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产中文字幕在线视频| 免费看十八禁软件| 天堂√8在线中文| 搡老熟女国产l中国老女人| 精品国产一区二区久久| 日韩视频一区二区在线观看| 新久久久久国产一级毛片| 亚洲午夜理论影院| 国产精品久久电影中文字幕 | 涩涩av久久男人的天堂| 亚洲成人免费电影在线观看| 亚洲成人手机| 多毛熟女@视频| 国产精品美女特级片免费视频播放器 | 三上悠亚av全集在线观看| 99国产精品99久久久久| 亚洲精品国产色婷婷电影| 国产精品久久电影中文字幕 | 女人爽到高潮嗷嗷叫在线视频| 亚洲国产欧美网| 美女福利国产在线| 女人精品久久久久毛片| 午夜精品在线福利| 宅男免费午夜| 超碰97精品在线观看| www日本在线高清视频| 人人妻人人添人人爽欧美一区卜| 国产97色在线日韩免费| 午夜亚洲福利在线播放| 亚洲精品美女久久av网站| 999久久久精品免费观看国产| 一级毛片精品| 50天的宝宝边吃奶边哭怎么回事| 三级毛片av免费| 中文字幕人妻丝袜制服| 中文字幕人妻丝袜一区二区| 一a级毛片在线观看| 91成人精品电影| 国产欧美日韩一区二区三区在线| 侵犯人妻中文字幕一二三四区| 亚洲国产精品合色在线| 9191精品国产免费久久| 丰满迷人的少妇在线观看| 可以免费在线观看a视频的电影网站| 夜夜夜夜夜久久久久| 久久精品亚洲av国产电影网| 日韩欧美免费精品| 精品亚洲成国产av| 久热爱精品视频在线9| 在线国产一区二区在线| √禁漫天堂资源中文www| 嫁个100分男人电影在线观看| 天堂动漫精品| 啦啦啦免费观看视频1| 女性被躁到高潮视频| 最新美女视频免费是黄的| 99久久精品国产亚洲精品| 桃红色精品国产亚洲av| 日韩视频一区二区在线观看| 波多野结衣av一区二区av| 欧美日韩精品网址| 日日摸夜夜添夜夜添小说| 国产精品久久电影中文字幕 | 午夜免费成人在线视频| 激情在线观看视频在线高清 | 国产深夜福利视频在线观看| 美女国产高潮福利片在线看| 看黄色毛片网站| 国产精品欧美亚洲77777| 国产精品亚洲av一区麻豆| 高清黄色对白视频在线免费看| 俄罗斯特黄特色一大片| 日韩精品免费视频一区二区三区| 黑人巨大精品欧美一区二区蜜桃| a级毛片黄视频| 搡老熟女国产l中国老女人| 久久午夜综合久久蜜桃| 久99久视频精品免费| 精品国产乱码久久久久久男人| 黄片小视频在线播放| 亚洲精品久久成人aⅴ小说| 国产深夜福利视频在线观看| 建设人人有责人人尽责人人享有的| 在线国产一区二区在线| 9色porny在线观看| 精品欧美一区二区三区在线| 色综合欧美亚洲国产小说| 少妇裸体淫交视频免费看高清 | 欧美精品人与动牲交sv欧美| 99精品欧美一区二区三区四区| 一区二区日韩欧美中文字幕| 欧美日韩精品网址| 麻豆成人av在线观看| 欧美+亚洲+日韩+国产| 老司机深夜福利视频在线观看| 黄网站色视频无遮挡免费观看| 一边摸一边做爽爽视频免费| 中文字幕制服av| 日韩制服丝袜自拍偷拍| 日韩欧美国产一区二区入口| 精品久久久精品久久久| 亚洲 欧美一区二区三区| 两人在一起打扑克的视频| 1024视频免费在线观看| 男女免费视频国产| 一个人免费在线观看的高清视频| 国产午夜精品久久久久久| 欧美黑人欧美精品刺激| 国产主播在线观看一区二区| 欧美 日韩 精品 国产| 国产精品一区二区精品视频观看| 99久久99久久久精品蜜桃| 夜夜躁狠狠躁天天躁| 免费久久久久久久精品成人欧美视频| 老司机午夜十八禁免费视频| 18在线观看网站| 久久人妻熟女aⅴ| 亚洲精品粉嫩美女一区| 在线免费观看的www视频| 国产一区在线观看成人免费| 在线观看一区二区三区激情| 亚洲视频免费观看视频| 欧美日本中文国产一区发布| 人妻一区二区av| 亚洲五月色婷婷综合| 高清欧美精品videossex| 中文字幕制服av| 色婷婷久久久亚洲欧美| 9191精品国产免费久久| av电影中文网址| av免费在线观看网站| 妹子高潮喷水视频| 国产亚洲一区二区精品| 12—13女人毛片做爰片一| 亚洲 国产 在线| 一区福利在线观看| 国产精品香港三级国产av潘金莲| 久久天堂一区二区三区四区| 久久久久久久久久久久大奶| 久久香蕉激情| 免费av中文字幕在线| 久久精品亚洲精品国产色婷小说| 村上凉子中文字幕在线| 精品人妻熟女毛片av久久网站| 国产黄色免费在线视频| 色尼玛亚洲综合影院| 午夜激情av网站| 亚洲国产看品久久| 亚洲,欧美精品.| 捣出白浆h1v1| 女性生殖器流出的白浆| 欧美精品一区二区免费开放| 黄色片一级片一级黄色片| 成人国语在线视频| 少妇的丰满在线观看| 欧美久久黑人一区二区| 成年女人毛片免费观看观看9 | 亚洲av日韩在线播放| 窝窝影院91人妻| www.熟女人妻精品国产| 欧美一级毛片孕妇| 在线播放国产精品三级| 欧美性长视频在线观看| tube8黄色片| 午夜福利在线免费观看网站| 亚洲七黄色美女视频| 在线观看免费日韩欧美大片| 国产精品美女特级片免费视频播放器 | 国产精品免费视频内射| 欧美精品亚洲一区二区| 精品一区二区三区视频在线观看免费 | 欧美日韩亚洲综合一区二区三区_| 亚洲av熟女| 成人手机av| 国产成人一区二区三区免费视频网站| 亚洲成国产人片在线观看| 国产精品 欧美亚洲| 欧美大码av| 一区二区三区国产精品乱码| 99久久99久久久精品蜜桃| 国产成人av激情在线播放| 黄色丝袜av网址大全| 国产片内射在线| 黄色毛片三级朝国网站| 色婷婷久久久亚洲欧美| 香蕉丝袜av| 母亲3免费完整高清在线观看| 中文字幕人妻熟女乱码| 欧美另类亚洲清纯唯美| 欧美成人免费av一区二区三区 | 一本大道久久a久久精品| 久久久精品国产亚洲av高清涩受| 亚洲av片天天在线观看| 亚洲 国产 在线| 国产精品二区激情视频| 国产成人精品在线电影| 女人精品久久久久毛片| 日韩免费高清中文字幕av| 成年人午夜在线观看视频| 99热网站在线观看| 看免费av毛片| 亚洲av片天天在线观看| 女同久久另类99精品国产91| 不卡一级毛片| 久久精品熟女亚洲av麻豆精品| 热99久久久久精品小说推荐| 国产男女超爽视频在线观看| 最新的欧美精品一区二区| 色婷婷久久久亚洲欧美| 欧美精品人与动牲交sv欧美| 色婷婷av一区二区三区视频| 精品国产乱子伦一区二区三区| 精品国内亚洲2022精品成人 | 免费日韩欧美在线观看| 精品人妻1区二区| 美女高潮到喷水免费观看| 电影成人av| 亚洲人成伊人成综合网2020| 美女扒开内裤让男人捅视频| 国产单亲对白刺激| 久久久久久久久免费视频了| 欧美日韩亚洲综合一区二区三区_| 亚洲专区国产一区二区| 啦啦啦 在线观看视频| 久久香蕉国产精品| 啪啪无遮挡十八禁网站| 亚洲人成电影观看| 亚洲第一欧美日韩一区二区三区| 操美女的视频在线观看| xxxhd国产人妻xxx| 国产成+人综合+亚洲专区| 三上悠亚av全集在线观看| av视频免费观看在线观看| 亚洲黑人精品在线| av福利片在线| 老司机福利观看| 国产精品一区二区免费欧美| 757午夜福利合集在线观看| 久久久久久久国产电影| 1024香蕉在线观看| 精品乱码久久久久久99久播| 亚洲五月婷婷丁香| 色综合欧美亚洲国产小说| 99精品在免费线老司机午夜| 伦理电影免费视频| 亚洲va日本ⅴa欧美va伊人久久| 老熟女久久久| av视频免费观看在线观看| 国产xxxxx性猛交| 久久精品熟女亚洲av麻豆精品| 久久久久精品人妻al黑| 在线看a的网站| 女警被强在线播放| 黄色a级毛片大全视频| 岛国在线观看网站| 十八禁高潮呻吟视频| 日本精品一区二区三区蜜桃| 天天操日日干夜夜撸| 日本vs欧美在线观看视频| 精品少妇久久久久久888优播| 在线观看午夜福利视频| √禁漫天堂资源中文www| 黄色成人免费大全| 亚洲av电影在线进入| 欧美大码av| 国产一区有黄有色的免费视频| 波多野结衣一区麻豆| 欧美日韩av久久| 99精国产麻豆久久婷婷| 久久青草综合色| 午夜福利一区二区在线看| 99久久精品国产亚洲精品| 欧美日韩黄片免| 天天操日日干夜夜撸| 韩国av一区二区三区四区| 国产欧美日韩精品亚洲av| 欧洲精品卡2卡3卡4卡5卡区| 免费少妇av软件| 黄色片一级片一级黄色片| 亚洲欧美激情在线| 母亲3免费完整高清在线观看| 在线国产一区二区在线| 日韩欧美三级三区| 成人18禁在线播放| 在线观看一区二区三区激情| 免费不卡黄色视频| 欧美日韩亚洲高清精品| 1024视频免费在线观看| 欧美性长视频在线观看| 亚洲人成77777在线视频| 搡老乐熟女国产| 日韩成人在线观看一区二区三区| 免费观看a级毛片全部| 777久久人妻少妇嫩草av网站| 亚洲性夜色夜夜综合| 久久影院123| 精品视频人人做人人爽| 国产97色在线日韩免费| 女人被躁到高潮嗷嗷叫费观| 亚洲五月婷婷丁香| 亚洲专区国产一区二区| 黄频高清免费视频| 亚洲在线自拍视频| 久久天堂一区二区三区四区| 国产人伦9x9x在线观看| 久久精品国产亚洲av香蕉五月 | 成年人免费黄色播放视频| 午夜成年电影在线免费观看| 搡老岳熟女国产| 亚洲人成伊人成综合网2020| 男人舔女人的私密视频| 欧美日韩视频精品一区| 国产亚洲av高清不卡| 亚洲国产欧美一区二区综合| 亚洲午夜精品一区,二区,三区| 丝袜美腿诱惑在线| 一区二区三区国产精品乱码| 亚洲中文av在线| 夜夜躁狠狠躁天天躁| av欧美777| 久久久久视频综合| 波多野结衣av一区二区av| 国产成人av教育| 国产欧美日韩一区二区三区在线| 国产欧美日韩综合在线一区二区| 欧美 日韩 精品 国产| 男人的好看免费观看在线视频 | 久久国产精品男人的天堂亚洲| 国产精品美女特级片免费视频播放器 | 色94色欧美一区二区| 黄色怎么调成土黄色| 黄片小视频在线播放| www.熟女人妻精品国产| 亚洲精品久久成人aⅴ小说| 精品久久蜜臀av无| 欧美精品亚洲一区二区| 国产一区在线观看成人免费| 国产精品久久久久久精品古装| 新久久久久国产一级毛片| 欧美日韩乱码在线| 国产成人欧美在线观看 | 欧美日韩瑟瑟在线播放| 国产三级黄色录像| 老司机福利观看| 纯流量卡能插随身wifi吗| 高清黄色对白视频在线免费看| 高潮久久久久久久久久久不卡| 99re在线观看精品视频| 久久午夜亚洲精品久久| 一本综合久久免费| 亚洲成av片中文字幕在线观看| 国产单亲对白刺激| 国产高清视频在线播放一区| 老鸭窝网址在线观看| 在线观看免费高清a一片| 亚洲精品久久成人aⅴ小说| 美女 人体艺术 gogo| 久久国产精品影院| 日韩视频一区二区在线观看| 欧美精品av麻豆av| 国产野战对白在线观看| 欧美日韩视频精品一区| 成年版毛片免费区| 欧美日韩中文字幕国产精品一区二区三区 | 国产免费男女视频| 国产精品二区激情视频| 欧美一级毛片孕妇| 日韩免费av在线播放| 免费看十八禁软件| 亚洲一区二区三区不卡视频| 亚洲av成人av| 国产精品 欧美亚洲| 99国产极品粉嫩在线观看| 最新的欧美精品一区二区| 好看av亚洲va欧美ⅴa在| 午夜福利,免费看| 精品欧美一区二区三区在线| 三上悠亚av全集在线观看| 国产黄色免费在线视频| 久久精品亚洲精品国产色婷小说| 亚洲成av片中文字幕在线观看| 国产亚洲欧美精品永久| 动漫黄色视频在线观看| 交换朋友夫妻互换小说| 国产99久久九九免费精品| 亚洲国产看品久久| 一进一出抽搐动态| av免费在线观看网站| 19禁男女啪啪无遮挡网站| av一本久久久久| 亚洲美女黄片视频| 麻豆国产av国片精品| 亚洲国产精品合色在线| 久久久水蜜桃国产精品网| 亚洲精品国产色婷婷电影| 免费看a级黄色片| 麻豆国产av国片精品| 99国产精品一区二区三区| 啦啦啦 在线观看视频| 搡老乐熟女国产| 国产真人三级小视频在线观看| 成年人黄色毛片网站| 国产精品1区2区在线观看. | 日韩一卡2卡3卡4卡2021年| 欧美亚洲日本最大视频资源| 中亚洲国语对白在线视频| 亚洲成人免费电影在线观看| 免费观看a级毛片全部| 亚洲av成人不卡在线观看播放网| 18禁裸乳无遮挡动漫免费视频| 久久久国产一区二区| 啦啦啦免费观看视频1| 国产精品免费大片| 一区二区日韩欧美中文字幕| 精品国产乱子伦一区二区三区| 久久久国产精品麻豆| 中文字幕精品免费在线观看视频| 亚洲一码二码三码区别大吗| 激情视频va一区二区三区| 久久香蕉激情| 极品少妇高潮喷水抽搐| √禁漫天堂资源中文www| 欧美黑人欧美精品刺激| 啦啦啦在线免费观看视频4| 99国产综合亚洲精品| 欧美 亚洲 国产 日韩一| 欧美在线一区亚洲| 亚洲精品成人av观看孕妇| 91国产中文字幕| www.自偷自拍.com| 99久久综合精品五月天人人| 热99久久久久精品小说推荐| 啦啦啦视频在线资源免费观看| 久久精品人人爽人人爽视色| 成年人午夜在线观看视频| 又大又爽又粗| 国产精品久久久久久精品古装| 精品一区二区三区四区五区乱码| 波多野结衣av一区二区av| 男女高潮啪啪啪动态图| 捣出白浆h1v1| 免费在线观看视频国产中文字幕亚洲| 啪啪无遮挡十八禁网站| 12—13女人毛片做爰片一| 夜夜夜夜夜久久久久| av欧美777| 高潮久久久久久久久久久不卡| 天天添夜夜摸| 亚洲色图av天堂| 久久久久国内视频| 精品电影一区二区在线| 精品久久蜜臀av无| 亚洲av美国av| 国产成人精品久久二区二区免费| 国产有黄有色有爽视频| 下体分泌物呈黄色| 别揉我奶头~嗯~啊~动态视频| 日本五十路高清| 午夜成年电影在线免费观看| 国产午夜精品久久久久久| 王馨瑶露胸无遮挡在线观看| cao死你这个sao货| 飞空精品影院首页| 婷婷丁香在线五月| 宅男免费午夜| 欧美乱码精品一区二区三区| 久久久久国内视频| 亚洲伊人色综图| 午夜成年电影在线免费观看| 久久香蕉国产精品| 欧美不卡视频在线免费观看 | 精品国内亚洲2022精品成人 | 91精品国产国语对白视频| 欧美日韩乱码在线| 国产亚洲一区二区精品| 日韩一卡2卡3卡4卡2021年| 久热这里只有精品99| 国产在线精品亚洲第一网站| 日韩视频一区二区在线观看| 精品少妇一区二区三区视频日本电影| 国产欧美日韩综合在线一区二区| 国产真人三级小视频在线观看| 黄色a级毛片大全视频| 欧美日韩亚洲综合一区二区三区_| 黑人巨大精品欧美一区二区mp4| 亚洲人成电影免费在线| avwww免费| 国产精品电影一区二区三区 | 亚洲精品国产区一区二| 久久精品国产99精品国产亚洲性色 | 久久午夜综合久久蜜桃| 女性生殖器流出的白浆| 亚洲自偷自拍图片 自拍| 日韩大码丰满熟妇| 自线自在国产av| 亚洲综合色网址| 高清视频免费观看一区二区| 咕卡用的链子| 中国美女看黄片| 亚洲全国av大片| 欧美日韩亚洲国产一区二区在线观看 | 高清毛片免费观看视频网站 | 91麻豆精品激情在线观看国产 | 国产一区有黄有色的免费视频| 欧美久久黑人一区二区| 久久久国产精品麻豆| 国产精品98久久久久久宅男小说| 男人的好看免费观看在线视频 | 18禁美女被吸乳视频| 妹子高潮喷水视频| 久久久久久亚洲精品国产蜜桃av| 日韩欧美三级三区| 国产免费男女视频| 欧美乱妇无乱码| 亚洲成人免费av在线播放| 国产日韩一区二区三区精品不卡| 久久国产精品大桥未久av| 曰老女人黄片| 热99re8久久精品国产| 欧美黑人精品巨大| 国产1区2区3区精品| 电影成人av| av线在线观看网站| 国产精品一区二区在线不卡| 国产一区二区三区在线臀色熟女 | 91麻豆精品激情在线观看国产 | 丝瓜视频免费看黄片| svipshipincom国产片| 老司机影院毛片| 一级黄色大片毛片| 女人精品久久久久毛片| 久久中文看片网| 美国免费a级毛片| 香蕉国产在线看| 国产精品久久久av美女十八| 我的亚洲天堂| 天天操日日干夜夜撸| 成人国产一区最新在线观看| 日韩欧美在线二视频 | 手机成人av网站| 欧美日韩福利视频一区二区| 久久精品熟女亚洲av麻豆精品| 一区二区三区国产精品乱码| 丝袜人妻中文字幕| 国产日韩欧美亚洲二区| 国产aⅴ精品一区二区三区波| 黄色视频,在线免费观看| 精品一区二区三区视频在线观看免费 | 午夜精品久久久久久毛片777| 日韩 欧美 亚洲 中文字幕| 国产激情欧美一区二区| 久久精品91无色码中文字幕| 精品亚洲成国产av| 久久国产精品人妻蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 久久婷婷成人综合色麻豆| 一级片'在线观看视频| 色94色欧美一区二区| 日韩中文字幕欧美一区二区| 婷婷丁香在线五月| 90打野战视频偷拍视频| 精品国产亚洲在线| 99久久综合精品五月天人人| 国产精品.久久久| 亚洲三区欧美一区| 久久久久视频综合| 成人免费观看视频高清| 精品少妇一区二区三区视频日本电影| 亚洲成国产人片在线观看| 久久久久精品人妻al黑| 国产精品免费视频内射| 女人被狂操c到高潮| 99精品欧美一区二区三区四区| 下体分泌物呈黄色| 免费不卡黄色视频| 亚洲熟妇熟女久久| 九色亚洲精品在线播放| 国产成人欧美在线观看 | 成年人午夜在线观看视频| 一进一出抽搐动态| 欧美黄色片欧美黄色片| 欧美日韩av久久| 在线av久久热| 欧美日本中文国产一区发布| 国产精品98久久久久久宅男小说| 成人国语在线视频| 久久久久久久午夜电影 |