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

    A Study on the Correlation Between Poles and Cuts in ππ Scattering?

    2019-11-07 02:58:32LingYunDai戴凌云XianWeiKang康現(xiàn)偉TaoLuo羅濤andUlfMeiner
    Communications in Theoretical Physics 2019年11期

    Ling-Yun Dai (戴凌云), Xian-Wei Kang (康現(xiàn)偉), Tao Luo (羅濤), and Ulf-G.Mei?ner

    1School of Physics and Electronics, Hunan University, Changsha 410082, China

    2College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China

    3Key Laboratory of Nuclear Physics and Ion-beam Application(MOE)and Institute of Modern Physics,Fudan University,Shanghai 200443, China

    4Helmholtz Institut für Strahlen- und Kernphysik and Bethe Center for Theoretical Physics, Universit?t Bonn, D-53115 Bonn, Germany

    5Institute for Advanced Simulation, Institut für Kernphysik and Jülich Center for Hadron Physics, Forschungszentrum Jülich, D-52425 Jülich, Germany

    Abstract In this paper we propose a dispersive method to describe two-body scattering with unitarity imposed.This approach is applied to elastic ππ scattering.The amplitudes keep single-channel unitarity and describe the experimental data well, and the low-energy amplitudes are consistent with that of chiral perturbation theory.The pole locations of the σ, f0(980), ρ(770) and f2(1270) and their couplings to ππ are obtained.A virtual state appearing in the isospin-two S-wave is confirmed.The correlations between the left (and right) hand cut and the poles are discussed.Our results show that the poles are more sensitive to the right hand cut rather than the left hand cut.The proposed method could be used to study other two-body scattering processes.

    Key words:dispersion relations, partial-wave analysis, chiral lagrangian, meson production

    1 Introduction

    In a two-body scattering system, for example two hadrons, the general principals that we know are unitarity,analyticity,crossing,the discrete symmetries,etc.The resonances that appear as the intermediate states in such system are important.Among them the lightest scalar mesons, related toππscattering, have the same quantum numbers as the QCD vacuum and are rather interesting,for some early references, see Refs.[1–3].Theππscattering amplitude is also crucial to clarify the hadronic contribution to the anomalous magnetic moment of the muon, see e.g.Refs.[4–5].To study the resonances in a given scattering process, one needs dispersion relations to continue the amplitude from the reals-axis (the physical region)to the complex-splane(see e.g.,Refs.[6–9])where the pole locations and their couplings are extracted.Following this method,some work on the light scalars can be found in Refs.[10–15], where the accurate pole locations and residues of theσandκmesons are given.

    For the dispersive methods, a key problem is how to determine the left hand cut (l.h.c.) and the right hand cut (r.h.c.), with the unitarity kept at the same time.In Refs.[7–8]the l.h.c.is estimated by crossed-channel exchange of resonances, where chiral effective field theory(χEFT) is used to calculate the amplitude.And the contribution of r.h.c.is represented by an Omnés function,with unitarity kept.In the well-known Roy equations,crossing symmetry and analyticity are perfectly combined together as the l.h.c.is represented by the unitary cuts of the partial waves.The single channel unitarity is also well imposed by keeping the real part of the partial wave amplitudes the same as what is calculated by the phase shift directly, which could be obtained by fitting to the experimental data in some analyses.Until now, Roy and Roy-Steiner equations certainly give the most accurate description of the two-body scattering amplitude and the information of resonances appearing as the intermediate states, such asππ,πKscattering and the pole locations and residues of theρ,σ,f0(980) andκ, etc., see e.g.Refs.[13–15].In addition, Ref.[14]shows that the l.h.c.can not be ignored for the determination of the pole location of theσ.By removing the parabola term of the l.h.c.,theσpole location is changed by about 15% accordingly,while the unitarity is violated due to the removal of the l.h.c.And thus the method to get the poles on the second Riemann sheet, calculated from the zeros of the S-matrix,is not reliable any more, as the method is based on the continuation implemented by unitarity.Here,we focus on obtaining a quantitative relation between cuts and poles,with unitarity imposed and the l.h.c.and r.h.c.are correlated with each other.

    This paper is organized as follows:In Sec.2 we establish a dispersive method based on the phase.In the physical region we also represent the amplitudes by an Omn`es function of the phase above threshold.In Sec.3 we fit theππscattering amplitudes up to 1 GeV in a model-independent way, including theIJ=00,02,11,20 waves, whereIdenotes the total isospin andJthe angular momentum.The fit results are the same as those given by the Omnés function representation and comparable with those of chiral perturbation theory (χPT) in the low-energy region.The poles and couplings are also extracted.In Sec.4 we give the estimation of the relation between poles and cuts, including both the l.h.c.and the r.h.c.We end with a brief summary.

    2 Scattering Amplitude Formalism

    2.1 A Dispersive Representation

    The two-body partial wave scattering amplitude can be written as:

    withφ(s) the phase andf(s) a real function.By writing a dispersion relation for lnT(s), one has:

    Here,s0is chosen at a specific point where the amplitude is real, and “L” denotes the l.h.c.and “R” stands for the r.h.c.The amplitude turns into

    On the other hand, unitarity is a general principal required for the scattering amplitude.In the single channel case one has

    wheresis in the elastic region andρ(s)is the phase space factor.Substituting Eq.(3) into Eq.(4), we obtain a representation (in the elastic region) for a single channel scattering amplitude

    Also, the Omnés function of the phase[16]for the l.h.c.is correlated with that of the r.h.c.which is again valid in the elastic region.A simple way to get the two-body scattering amplitude proceeds in two steps:First, we follow Eq.(5) to fit the Omnés function of the r.h.c.to experimental data, and then use Eq.(6)and other constraints below the threshold to fit the Omnés function of the l.h.c.Note that Eq.(5)does not only work for the single channel case, but also for the coupled channel case in the physical region.

    2.2 On ππ Scattering

    In the equations above, the threshold factor is not included.Considering such factors, we need to change the amplitudes into:

    Here and in what follows, we takeππscatering as an example.Thus one haszJI=4Mπ2for the P-, D-, and higher partial waves,andzIJis the Adler zero for the S-waves.nJis one for S- and P-waves and two for D waves.We define a reduced amplitude

    and again we can write a dispersion relation for ln(s),so that we have

    Here,s0could be chosen from the range[0,4M2π].For the r.h.c., we cut off the integration somewhere in the high energy region, see discussions in the next sections.We have

    TheTIJ(s0) could be fixed byχPT or scattering lengths,or other low-energy constraints.For simplicity, we chooses0=0.Combining unitarity, embodied by Eq.(5), we have a correlation between Omnés functions of l.h.c.and r.h.c.in the elastic region

    This is similar to Eq.(6).Substituting Eq.(11) into Eq.(10), we still have Eq.(5).Since we know theππscattering amplitudes well in the region [4M2π, 2 GeV2]andχPT describes the amplitudes well in the low-energy region, we have to fit the l.h.c.to both Eq.(6) andχPT.

    3 Phenomenology

    3.1 ππ Scattering Amplitudes

    For theππscattering amplitude, we can parameterize the phase caused by the l.h.c.by a conformal mapping

    with

    Notice that Imω(s)behaves asarounds=0,which is consistent with that ofχPT,see Ref.[12]and references therein.

    As concerns the r.h.c., it is less known in the high energy region.However, these distant r.h.c.should have less important effects in the low-energy region, especially in the regions ≤1 GeV2.We choose three kinds of?IJR(s) to test the stability and uncertainty caused by the distant r.h.c.In Case A, the phases are cut off ats=2.25 GeV2.[7]In Case B, the phases are given by Ref.[7],up tos=22 GeV2.In Case C,the phases/Omn`es functions of the r.h.c.are given by Refs.[17–19]and references therein, up tos=22 GeV2.Here, the phases are fitted to the experimental data in Refs.[20–21]up t√o=2 GeV and constrained by unitarity up to=4 GeV.Notice that in Case A and B the phase of the isospin-one P-wave is given by CFDIV in Ref.[22],and we continue it to the higher energy region by means of the function

    with

    The function (and also its first derivative) is smooth at the pointsR.We setk=1,n=2,sR=1.42GeV2andφ11∞=160?, which is close toφ11R((1.4 GeV)2)=170?and ensures that the phase in the high energy region behaves smoothly.The upper limits of the integration of the r.h.c.of isospin-one P-wave are the same as the other partial waves.The three kinds of Omnes functions have different magnitudes.For example, in isospin zero S-wave the difference between Cases A/B and C is 17/6 percent ats=0.5 GeV2.And that is 2/17 percent for isospin one P-wave.

    The parameters of our fits for all the Cases are given in Table 1.ThecIJnare determined by the following procedure.In the elastic region, we choose one or two “mesh points”, equaling to the number of the coefficientscIJn.Combining Eqs.(11) and (12), we can build a matrix and solve forcIJn.We adjust the “mesh points” andsIJLto make the solutions to be consistent with that of Eq.(5)in the physical region.For instance, we choose one “mesh point” ass1=0.7 GeV2in Case A, and two points ass1=0.2 GeV2ands2=0.6 GeV2in Cases B and C.Other partial waves are treated similarly.Note that we do not choose more mesh points, as the correlation between each terms of dispersion relation will become stronger and we will have larger number ofcIJn, resulting in more twisted amplitudes.This strategy gives a good description of the amplitudes, with unitarity kept.See the fit results shown in Fig.1.Here all the partial waves refer to Case B,in which the phase is cut off ats=2.25 GeV2.The amplitudes from the other Cases are quite close to this one, except for the inelastic region and the distant l.h.c.(≤?0.4 GeV2).Our fit, both the real part (black solid line) and imaginary (black dotted line) part of the amplitudes shown in Fig.1, is indistinguishable from that given by the K-Matrix[7]or CFDIV.[22]Note that the amplitudes given by Eq.(5) are exactly the same as those of the K-Matrix or CFDIV fromππthreshold to the inelastic threshold.This implies that the unitarity is respected.To test it quantitatively, we define

    ?TIJ(sn)is the difference between our amplitude and that of Eq.(5).We use this equation to refine our solutions,that is, theand “mesh points” are selected out to makesmall.Here we choosesn=(0.1–0.9) GeV2for the S-waves,sn=(0.1–0.8) GeV2for the P-wave,andsn=(0.1-1.0) GeV2for the D-wave, with step of 0.1 GeV2.These points are located between theππand the inelastic thresholds.From here on all the steps are chosen to be 0.1 GeV2(or 0.1 GeV for).We find thatand=1.4%.The violation of unitarity is rather small at elastic region.Notice that if we remove out the Omnes function of the right hand phase, only the l.h.c.part(timesT(0)) remains and it is a smooth real function aboves >0.Since the contribution of the l.h.c.is determined rather well in the elastic region, it would be natrual to keep working well not faraway from the inelastic threshld.This is why the amplitudes could be described well up to 1.2 GeV for S- and P-waves.These are enough for us to discuss the poles and cuts ofσ,ρa(bǔ)ndf0(980).

    ForTIJ(0),χPT could be used to fix it.The analytical SU(3) 1-loopχPT amplitudes of each partial waves,are recalculated and given in Appendix A.The low-energy constants are given by Ref.[23].Those of SU(2) 2-loopχPT amplitudes are given by Refs.[24–25]and references therein.All the values ofTIJ(0) in Table 1 are chosen to be close to the prediction ofχPT or our earlier analyses,[17?18]as well as to minimize.In the isospin-zero S-wave, the magnitude of ourT00(0) is a bit larger than that ofχPT.This is consistent with what is known about this scattering length, where the one-loopχPT calculation gives a smaller result than what is obtained by dispersive methods, Roy equations or in experiment, see e.g.the review.[23]A better comparison would be given with the 2-loopχPT amplitudes.

    Fig.1 (Color online) Fit of the ππ scattering amplitudes for Case B.Notice that =sgn(s) The solid lines denote the real part of the amplitudes and the dotted, dashed, dash-dotted and dash-dot-dotted lines denote the imaginary part.The black lines are from our fit.The red lines are from a K-matrix fit[7] for the isospin-zero S-wave, and the violet lines are from CFDIV[22] for other waves.The borders of the cyan and green bands in the low-energy region are from SU(2) and SU(3) χPT, respectively.The CERN-Munich data are from Ref.[20], and the OPE and OPE-DP data are from Ref.[21].

    Table 1 The parameters for each fit.“-” means the absence of the corresponding quantity.For comparison,we also give the one-loop χPT results.

    In the isospin-zero D-wave, theT02(0) varies more in the different Cases.The reason is that some fine-tuning is needed as the inelastic r.h.c.is difficult to be implemented well.The amplitudes given by Eq.(5) are much different from that of CFDIV in the inelastic region where thef2(1270) appears.Notice further that the value ofT0D(0)is very small, one order smaller than that of the other waves.

    3.2 Pole Locations and Couplings

    With these amplitudes given by a dispersion relation,the information of the poles can be extracted.It is worth pointing out that the main goal of this section is to check the reliability of our representation, that is, whether our model is good on extracting the poles and residues in the(deep) complex plane or not.The polesRand its coupling/residuegfππon the second Riemann sheet are defined as

    Note that the continuation of theT(s) amplitude to the second Riemann sheet is based on unitarity,

    The poles and couplings/residues for Cases A, B, C are given in Table 2.

    Table 2 The pole locations and residues given by our fits.The notation “2S v.s.” denotes the virtual state in the isospin-two S-wave.

    All the poles and residues of the different Cases are close to each other, and also to those given by Roy[14]or Roy-like[22]equations,?Here we stress that we do not intend to extract more accurate pole locations than those of Roy or Roy-like equations.Note we o give an explicit description on the l.h.c., while it is hidden in the Roy or Roy-like equations.It has already been discussed in Ref.[26], within a unitarized χPT method.Here we use dispersion approach and re-confirm it, but we do not have the extra poles caused by unitarization.Recently, Lattice QCD gives negative scattering length as ?0.0412(08)(16),[27] ?0.04430 (25)(40).[28] ?0.04430 (2)Mπ? 1.[29] These values are consistent with that of χPT,[23] the Roy equations matched to χPT[30] and a dispersive analysis.[22]except for the pole location of thef2(1270).The reason is that thef2(1270) is located outside the elastic unitary cut ofT0D(s), while Eq.(11)only works in the elastic region.For this partial wave one needs a more dedicated method to study, including coupled-channel unitarity.For the poles of theσ, the differences between the different Cases is also a bit larger than those of other resonances such as theρ(770) and thef0(980).This is because theσis far away from the real axis.For the virtual state in the isospin-two S-wave, the poles and residues are a bit different from Cases A and B to Case C.This situation is comparable with that ofT2S(0),where in Cases A and BT2S(0) is 0.055 and in Case C it is 0.060, respectively.

    In addition, we also find that there exists a virtual state in the isospin-two S-wave very close tos=0.?Here we stress that we do not intend to extract more accurate pole locations than those of Roy or Roy-like equations.Note we o give an explicit description on the l.h.c., while it is hidden in the Roy or Roy-like equations.It has already been discussed in Ref.[26], within a unitarized χPT method.Here we use dispersion approach and re-confirm it, but we do not have the extra poles caused by unitarization.Recently, Lattice QCD gives negative scattering length as ?0.0412(08)(16),[27] ?0.04430 (25)(40).[28] ?0.04430 (2)Mπ? 1.[29] These values are consistent with that of χPT,[23] the Roy equations matched to χPT[30] and a dispersive analysis.[22]According to Eq.(18), the virtual state at the zero of the S-matrix below the threshold.This zero equals to the intersection point between two lines:T2S(s)and i/2ρ(s).As shown in Fig.2 the line of theT2S(s)and the line of i/2ρ(s)will always intersect with each other and the crossing point always lies in the energy region of [0,sa], wheresais the Adler zero.This is the virtual state.Since the scattering length is negative and the Adler zero (only one) is below threshold, one would expect that the amplitude ofT2S(s),froms=4M2πtos=0, will always cross the real axis ofsand arrive at the positive vertical axis.In all events, it will intersect with that of the i/2ρ(s).Thus the existence of the virtual state is confirmed.This inference is modelindependent, only the sign of the scattering length,§Here we stress that we do not intend to extract more accurate pole locations than those of Roy or Roy-like equations.Note we o give an explicit description on the l.h.c., while it is hidden in the Roy or Roy-like equations.It has already been discussed in Ref.[26], within a unitarized χPT method.Here we use dispersion approach and re-confirm it, but we do not have the extra poles caused by unitarization.Recently, Lattice QCD gives negative scattering length as ?0.0412(08)(16),[27] ?0.04430 (25)(40).[28] ?0.04430 (2)Mπ? 1.[29] These values are consistent with that of χPT,[23] the Roy equations matched to χPT[30] and a dispersive analysis.[22]the Adler zero, and analyticity are relevant.

    For a general discussion of the virtual state arising from a bare discrete state in the quantum mechanical scattering, we recommend readers to read[31?32]and references therein.We suggest that the isospin-two S-wave amplitude could be checked in the future measurement of Λ+c →Σ?π+π+.Its branching ratio[33]is large enough.

    Fig.2 (Color online)The lines of T2S(s)for the different Cases and i/2ρ(s).Note that all of them are real.The intersection point corresponds to the virtual state.

    The average values of the poles and residues of all the Cases define our central values.The deviations of the different Cases to the central values are used to estimate the uncertainties.The results are shown in Table 3.These are very similar from those of previous analyses.[7,14?15,34?35]Thef2(1270) has a much larger uncertainty compared to the other resonances, just as discussed before.The residues of all resonances have roughly similar magnitude at the region [0.25,0.55]GeV, except for that of the virtual state in the isospin-two S-wave,which is much weaker.But their phases are quite different.The phases ofρ(770)andf2(1270) are close to zero, while those of theσandf0(980) are close to?90?, and the virtual state one is close to 90?.This may imply thatρ(770) andf2(1270)are normal ˉqqstates but that theσandf0(980)have large molecular components.

    Table 3 The pole locations and residues by taking the averages between Cases A, B and C as discussed in the text.

    3.3 The Correlation Between Poles and Cuts

    It is interesting to find the correlation between the poles and cuts.We focus here on the isospin-zero S-wave and isospin-one P-wave, as thef2(1270) is far away from the l.h.c.and the virtual state is too close to the l.h.c.Also,the light scalars are more difficult to understand.All the fits of different Cases about these two partial waves are shown in Fig.3.In our approach only unitarity is used to constrain the amplitudes, but the low-energy amplitudes are consistent with those ofχPT.Only in Cases B and for the isospin-zero S-wave, theT0S(s) amplitude is inconsistent with that ofχPT at

    Fig.3 (Color online) Comparison of different solutions of the ππ scattering amplitudes.The solid lines are the real part of the amplitudes and other lines are the imaginary part.The green lines are from Case A,the cyan lines are from Case B, and the black lines are from Case C.The red lines are from K-Matrix[7] for isospin 0 S-wave,and the violet lines are from CFDIV[22] for isospin 1 P-wave.Note that the lines of K-matrix and/or CDFIV are overlapped with our fits in the elastic region, or even a bit further in the inelastic region.

    To see the variation of the l.h.c.in the different solutions,we apply Eq.(19)on ImTIJ(sn),with=?0.1 GeV to?0.6 GeV.Notice that at=0 all of the l.h.c.are zero and behave asthis partly ensures the l.h.c.to be consistent with that ofχPT in the low-energy region.We fix the average value of all Cases as the central value, and calculate the relative deviation for each point.At last we avarage these relative deviations to estimate the variation of the cuts.The variation of cuts and poles are defined as

    withspthe pole on the second Riemann Sheet.Finally,we collect the uncertainties in Table 4.And we define the correlation between poles and cuts as

    The simple meaning of the correlationCpoleis to answer the following question:When the cut is changed by 100%,how much would the pole location be changed?

    To test the correlation between poles and the r.h.c.,we simply setφtest(s)=1.04φ(s),and check the variation of poles and cuts, respectively.The relative uncertainty of the r.h.c.is also estimated by Eq.(19), withsn=(0.1–0.9) GeV2for the isospin-zero S-wave andsn=(0.1–0.8) GeV2for the isospin-one P-wave.The relative uncertainty of the poles and the correlation are calculated in the same way as that of the l.h.c., see Eqs.(19) and (20).

    Table 4 The correlation between l.h.c.and poles (represented by superscript‘L’),and between r.h.c.and poles(represented by superscript “R”).

    From Table 4, we find thatCRσis roughly two orders larger than that ofCLσ,thoughσis rather close to the l.h.c.Comparing to Ref.[14], which has roughly 15% contribution from l.h.c.,we have a rather smaller contribution from the l.h.c.,caused by the constraint of unitarity on the l.h.c.Also,is roughly three orders larger than that of, andis roughly one order larger than that of.These indicate that the correlation between the unitarity cut and the poles is much larger than that of the l.h.c.and poles.Note that in our case the l.h.c.is not arbitrary but correlated with the r.h.c.,constrained by unitarity and analyticity, see Eq.(11).For each Case,CLσis larger than.This is not surprising as theσis much closer to the l.h.c.Also,CRσis larger than.The reason is that theσis farther away from the real axis, the uncertainty of the pole is larger as the amplitude is continued from the physical region deeper into the complex-splane.It is interesting to see that in averageCLσis roughly two times larger thanAnd for the distance between these poles and l.h.c.(simply sets=0), |sσ| is one half of that of|sρ|, this tells us that the correlation between poles and l.h.c.is inversely proportional to their distance.In contrast,CRσis roughly one order larger than.For the distance between these poles and r.h.c.(simply sets=Respole),|Imsσ|is two times larger than|Imsρ|, this tells us that the correlation between poles and r.h.c.is proportional to their distance.These conclusions are still kept when comparing theσand thef0(980).The discussion of the correlation between poles and cuts reveals that our model, keeping analyticity and unitarity, give strong constrains on poles and l.h.c., even if the pole is a bit faraway from the real axis.

    4 Summary

    We proposed a dispersive method to calculate the twobody scattering amplitude.It is based on the Omn`es function of the phase, including that of the left hand cut and the right hand cut.The input of the r.h.c.is given by three kinds of parametrizations,and the l.h.c.is solved by Eq.(11), with unitarity and analyticity respected.The pion-pionIJ=00,02,11,20 waves are fitted within our method and the poles and locations are extracted.They are stable except for that of thef2(1270),which lies in the inelastic region.The r.h.c.has much larger contribution to the poles comparing to that of the l.h.c.This method could be useful for the studies of strong interactions in two-body scattering.For instance, it could be applied in theπηscattering, where the correlation between l.h.c.and r.h.c.could help us to constrain the phase shift.Also theππscattering amplitudes obtained here could be used for the future studies when one hasππfinal state interactions, see e.g.Refs.[36–40], and/or to multi-pions, see e.g.Refs.[41–42].

    Appendix A:Analytical Amplitudes of Partial Waves Within χPT

    The analytical 1-loop amplitudes of these partial waves within of SU(3)χPT are recalculated.For reader’s convenience, they are given below.We have the IJ=00 waves up toO(p4):

    Here the superscript ofTin the bracket means the chiral order, and subscripts represent for Isospin and spin, respectively.Note that for reader’s convenience we also give the analytical forms of the imaginary part (r.h.c.) of the amplitudes.TheI=2 S wave is

    TheI=1 P wave is

    And theI=0 D wave is

    It should be noted that in all these partial waves, 2L4+L5and 2L6+L8appear together.[43]TheA,B,Hfunctions are given as below

    withρ(t,m)=Notice that our amplitudes are calculated in the formalism ofwhile that of Refs.[44–47]is done in?1.The relation between our LECs (Li) and that of the latter one () is=Li+Γi/32π2.

    Acknowledgments

    We are grateful to Zhi-Yong Zhou for helpful discussions and for supplying us the files of the SU(2) 2-loopχPT amplitudes.

    亚洲国产欧美人成| 亚洲欧美成人精品一区二区| 色视频www国产| 婷婷色综合www| 菩萨蛮人人尽说江南好唐韦庄| www.av在线官网国产| 国产亚洲av片在线观看秒播厂| 亚洲精品成人av观看孕妇| 欧美成人a在线观看| 99久久精品热视频| 男女免费视频国产| 又爽又黄a免费视频| 亚洲欧美清纯卡通| 丰满乱子伦码专区| 亚洲熟女精品中文字幕| 亚洲欧美日韩东京热| 精品亚洲成a人片在线观看 | 久久99热6这里只有精品| 高清av免费在线| 国产一区亚洲一区在线观看| 最近最新中文字幕免费大全7| 国产精品偷伦视频观看了| 高清午夜精品一区二区三区| 精品亚洲成国产av| av不卡在线播放| 美女主播在线视频| 人人妻人人澡人人爽人人夜夜| 美女中出高潮动态图| 国产乱来视频区| 亚洲国产最新在线播放| 高清不卡的av网站| 青青草视频在线视频观看| 男女国产视频网站| 国产男女内射视频| 美女主播在线视频| 亚洲第一区二区三区不卡| 久久亚洲国产成人精品v| 中文字幕亚洲精品专区| 国产精品熟女久久久久浪| 国产精品三级大全| 亚洲欧美日韩无卡精品| 国产精品久久久久久久电影| 欧美老熟妇乱子伦牲交| 亚洲精品一区蜜桃| 妹子高潮喷水视频| 丝袜脚勾引网站| 欧美日韩国产mv在线观看视频 | 久久久精品免费免费高清| 久久99热6这里只有精品| 亚洲综合精品二区| 亚洲综合精品二区| 777米奇影视久久| 日本与韩国留学比较| 欧美日韩综合久久久久久| 亚洲av中文字字幕乱码综合| 久久精品久久久久久久性| 精品酒店卫生间| 国产精品99久久久久久久久| 亚洲国产精品国产精品| 韩国av在线不卡| 亚洲国产毛片av蜜桃av| 青春草亚洲视频在线观看| 国产免费又黄又爽又色| 91在线精品国自产拍蜜月| 国产一区有黄有色的免费视频| 欧美国产精品一级二级三级 | 欧美精品国产亚洲| 国产精品一区二区在线观看99| 国产成人freesex在线| av卡一久久| 国产爱豆传媒在线观看| 极品教师在线视频| 亚洲国产欧美在线一区| 伦精品一区二区三区| 国模一区二区三区四区视频| 日韩一本色道免费dvd| 色婷婷久久久亚洲欧美| 久久久久精品性色| 少妇熟女欧美另类| 丰满迷人的少妇在线观看| 久久久精品94久久精品| 亚州av有码| 我的女老师完整版在线观看| 午夜激情久久久久久久| 99九九线精品视频在线观看视频| 欧美xxxx性猛交bbbb| 一级爰片在线观看| 中文欧美无线码| 成人美女网站在线观看视频| 国产一区二区在线观看日韩| 狂野欧美激情性bbbbbb| 日本免费在线观看一区| 国产av国产精品国产| 亚洲,欧美,日韩| 激情五月婷婷亚洲| 久热这里只有精品99| av免费在线看不卡| 久热久热在线精品观看| 久久热精品热| 久久久久网色| 国产熟女欧美一区二区| av卡一久久| 伦理电影大哥的女人| 久久 成人 亚洲| 少妇高潮的动态图| 日本黄色日本黄色录像| 人妻 亚洲 视频| 97在线视频观看| 18+在线观看网站| 成人漫画全彩无遮挡| 国产精品99久久久久久久久| 国产av国产精品国产| 亚洲最大成人中文| 久久青草综合色| 国产黄片美女视频| 黑人猛操日本美女一级片| 精品久久久精品久久久| 晚上一个人看的免费电影| 久久婷婷青草| av在线蜜桃| 97超视频在线观看视频| 偷拍熟女少妇极品色| 国产精品一区www在线观看| 老师上课跳d突然被开到最大视频| 成人黄色视频免费在线看| 99热这里只有是精品50| 青青草视频在线视频观看| 欧美精品一区二区免费开放| 精品人妻一区二区三区麻豆| 我要看日韩黄色一级片| 亚洲电影在线观看av| 精品一区二区三区视频在线| 五月伊人婷婷丁香| 午夜免费观看性视频| 少妇猛男粗大的猛烈进出视频| 成人毛片a级毛片在线播放| 香蕉精品网在线| 日韩av免费高清视频| 久久女婷五月综合色啪小说| 久久精品夜色国产| 99热网站在线观看| 欧美人与善性xxx| 国产老妇伦熟女老妇高清| 欧美bdsm另类| 尾随美女入室| 国产淫片久久久久久久久| 国产精品秋霞免费鲁丝片| 尾随美女入室| 国产成人aa在线观看| 七月丁香在线播放| 最近最新中文字幕免费大全7| 亚洲精品久久午夜乱码| 下体分泌物呈黄色| 亚洲国产精品国产精品| 亚洲性久久影院| 大话2 男鬼变身卡| 久久精品国产自在天天线| 伦理电影免费视频| 欧美日韩国产mv在线观看视频 | 国产黄色免费在线视频| 国产精品人妻久久久影院| 啦啦啦视频在线资源免费观看| 亚洲欧美日韩另类电影网站 | 亚洲色图综合在线观看| 在线观看av片永久免费下载| 老女人水多毛片| 一二三四中文在线观看免费高清| 亚洲精品日韩在线中文字幕| 国产日韩欧美在线精品| 午夜福利视频精品| 亚洲欧美精品专区久久| 亚洲国产毛片av蜜桃av| 久久韩国三级中文字幕| av免费观看日本| 欧美高清成人免费视频www| 日韩免费高清中文字幕av| 少妇被粗大猛烈的视频| 中文资源天堂在线| 女性被躁到高潮视频| 久久精品久久久久久噜噜老黄| 亚洲av中文字字幕乱码综合| 日本欧美国产在线视频| 精品国产三级普通话版| 免费不卡的大黄色大毛片视频在线观看| 国产欧美亚洲国产| 亚洲成人一二三区av| 激情五月婷婷亚洲| 亚洲欧美精品自产自拍| 偷拍熟女少妇极品色| 免费少妇av软件| 日韩中字成人| 三级国产精品片| 永久网站在线| 午夜福利视频精品| 中文资源天堂在线| 日韩av免费高清视频| 亚洲精品一区蜜桃| av免费观看日本| 少妇熟女欧美另类| 1000部很黄的大片| 日韩欧美精品免费久久| 亚洲欧美日韩东京热| 又粗又硬又长又爽又黄的视频| 亚洲国产成人一精品久久久| 日韩强制内射视频| 男女边吃奶边做爰视频| 男女边吃奶边做爰视频| 免费av中文字幕在线| 免费av中文字幕在线| 国产亚洲最大av| 日韩强制内射视频| 青春草亚洲视频在线观看| 国产精品女同一区二区软件| 我要看日韩黄色一级片| 美女xxoo啪啪120秒动态图| 丰满迷人的少妇在线观看| 免费大片18禁| av国产精品久久久久影院| 国产高清有码在线观看视频| 免费观看无遮挡的男女| 免费观看无遮挡的男女| 亚洲成人手机| 性色avwww在线观看| 亚洲欧美清纯卡通| 国产老妇伦熟女老妇高清| 丝袜脚勾引网站| 少妇熟女欧美另类| 久久精品国产亚洲av天美| 久久久久久久国产电影| 国产爽快片一区二区三区| 精品亚洲乱码少妇综合久久| 精品亚洲成国产av| 一级毛片aaaaaa免费看小| 国产高潮美女av| 夜夜爽夜夜爽视频| 欧美一级a爱片免费观看看| 久久这里有精品视频免费| 午夜视频国产福利| 人妻一区二区av| 国产在线免费精品| .国产精品久久| 欧美xxⅹ黑人| 日韩av不卡免费在线播放| av在线观看视频网站免费| 国产一区二区三区av在线| 国内揄拍国产精品人妻在线| 精品一区二区免费观看| 99久久中文字幕三级久久日本| av又黄又爽大尺度在线免费看| 91精品国产国语对白视频| 丰满少妇做爰视频| 99久久精品国产国产毛片| 国产精品久久久久久久久免| 尤物成人国产欧美一区二区三区| 国产一区有黄有色的免费视频| 精品视频人人做人人爽| 80岁老熟妇乱子伦牲交| 国产成人午夜福利电影在线观看| 欧美激情国产日韩精品一区| 国产在线免费精品| 成人免费观看视频高清| 制服丝袜香蕉在线| 1000部很黄的大片| 久久99精品国语久久久| 视频区图区小说| 一本—道久久a久久精品蜜桃钙片| 亚洲内射少妇av| 在线观看美女被高潮喷水网站| av国产精品久久久久影院| 日韩欧美精品免费久久| 丝袜脚勾引网站| 成年免费大片在线观看| 亚洲av不卡在线观看| 日日啪夜夜爽| 日韩欧美一区视频在线观看 | 尾随美女入室| 天天躁日日操中文字幕| 熟女人妻精品中文字幕| 在线播放无遮挡| 91aial.com中文字幕在线观看| 少妇裸体淫交视频免费看高清| 国产深夜福利视频在线观看| 伦理电影大哥的女人| 亚洲av.av天堂| 王馨瑶露胸无遮挡在线观看| 欧美极品一区二区三区四区| 精品国产露脸久久av麻豆| 亚洲国产精品专区欧美| 精品久久国产蜜桃| 免费高清在线观看视频在线观看| 搡老乐熟女国产| 人妻制服诱惑在线中文字幕| 精品亚洲成a人片在线观看 | 伦精品一区二区三区| 亚洲最大成人中文| 一个人看的www免费观看视频| 免费av不卡在线播放| 国产精品熟女久久久久浪| 国产男人的电影天堂91| 成人国产av品久久久| 黄片无遮挡物在线观看| 国产午夜精品久久久久久一区二区三区| 国产精品福利在线免费观看| 欧美日韩视频精品一区| 视频区图区小说| 爱豆传媒免费全集在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 搡女人真爽免费视频火全软件| 久久久成人免费电影| 亚洲av男天堂| 97超视频在线观看视频| 国产成人精品婷婷| 国产亚洲最大av| av国产精品久久久久影院| 18禁在线无遮挡免费观看视频| 欧美日韩国产mv在线观看视频 | 亚洲国产日韩一区二区| 亚洲经典国产精华液单| 成人亚洲欧美一区二区av| 国产精品久久久久久久电影| 国产女主播在线喷水免费视频网站| 精品人妻偷拍中文字幕| 91aial.com中文字幕在线观看| 少妇丰满av| 亚洲av福利一区| 高清不卡的av网站| 少妇 在线观看| 香蕉精品网在线| h视频一区二区三区| 在线精品无人区一区二区三 | 一区二区三区四区激情视频| 在线免费十八禁| 蜜桃久久精品国产亚洲av| av专区在线播放| 国产男女超爽视频在线观看| 国产亚洲5aaaaa淫片| 日韩中文字幕视频在线看片 | 80岁老熟妇乱子伦牲交| 精品国产一区二区三区久久久樱花 | 韩国av在线不卡| 国产成人a区在线观看| 久久精品国产亚洲av涩爱| 寂寞人妻少妇视频99o| 在线观看一区二区三区| 激情五月婷婷亚洲| 日本一二三区视频观看| 国产精品蜜桃在线观看| av视频免费观看在线观看| 国产av码专区亚洲av| av免费观看日本| 伦理电影免费视频| 热99国产精品久久久久久7| 99热这里只有精品一区| 一级毛片 在线播放| 国产精品不卡视频一区二区| 全区人妻精品视频| 免费人妻精品一区二区三区视频| 啦啦啦在线观看免费高清www| 国产成人免费无遮挡视频| 国产91av在线免费观看| 男女啪啪激烈高潮av片| 狂野欧美激情性xxxx在线观看| 热99国产精品久久久久久7| 亚洲美女搞黄在线观看| 人人妻人人爽人人添夜夜欢视频 | 51国产日韩欧美| 免费少妇av软件| 精品国产露脸久久av麻豆| 成人国产麻豆网| 免费观看性生交大片5| 秋霞在线观看毛片| 日韩电影二区| 人妻 亚洲 视频| 国产av精品麻豆| 深夜a级毛片| 成人特级av手机在线观看| 亚洲精品国产成人久久av| 亚洲欧美清纯卡通| 亚洲av免费高清在线观看| 香蕉精品网在线| 91狼人影院| 亚洲av中文av极速乱| 精品久久久久久久久av| 国模一区二区三区四区视频| 狠狠精品人妻久久久久久综合| 国产伦理片在线播放av一区| 国产乱人视频| 99热网站在线观看| 美女福利国产在线 | 久久精品熟女亚洲av麻豆精品| 精品国产一区二区三区久久久樱花 | 九九在线视频观看精品| 少妇人妻 视频| 免费观看无遮挡的男女| 特大巨黑吊av在线直播| 人妻一区二区av| 免费人妻精品一区二区三区视频| 日韩强制内射视频| 最近的中文字幕免费完整| 在线观看人妻少妇| 久久久久久久久久久免费av| 亚洲熟女精品中文字幕| 美女xxoo啪啪120秒动态图| 少妇 在线观看| 人妻 亚洲 视频| 天天躁夜夜躁狠狠久久av| 五月玫瑰六月丁香| 国产精品久久久久成人av| 中文在线观看免费www的网站| 亚洲成人av在线免费| 男女下面进入的视频免费午夜| 日韩成人av中文字幕在线观看| 亚洲av不卡在线观看| 国产欧美亚洲国产| 日韩中字成人| 久久韩国三级中文字幕| 蜜桃久久精品国产亚洲av| 女人久久www免费人成看片| 美女主播在线视频| 精品少妇久久久久久888优播| 免费看av在线观看网站| 久久av网站| 国产精品无大码| 亚洲精品久久午夜乱码| 中文字幕亚洲精品专区| 久久久久国产网址| 日韩伦理黄色片| 国精品久久久久久国模美| 中文字幕制服av| 国产白丝娇喘喷水9色精品| 日本wwww免费看| 最近手机中文字幕大全| 免费少妇av软件| 日韩欧美 国产精品| 成人综合一区亚洲| 日韩欧美 国产精品| 激情五月婷婷亚洲| 国产视频内射| 欧美+日韩+精品| 欧美bdsm另类| 性色avwww在线观看| 成人美女网站在线观看视频| 日本与韩国留学比较| 欧美一区二区亚洲| 如何舔出高潮| 成人亚洲精品一区在线观看 | 美女xxoo啪啪120秒动态图| 超碰97精品在线观看| 色综合色国产| 高清日韩中文字幕在线| 日本一二三区视频观看| 国产欧美日韩精品一区二区| 午夜日本视频在线| 日韩成人av中文字幕在线观看| 亚洲国产成人一精品久久久| 成人亚洲欧美一区二区av| 777米奇影视久久| 少妇 在线观看| 伦精品一区二区三区| 视频区图区小说| 免费观看的影片在线观看| 成人午夜精彩视频在线观看| 国产黄频视频在线观看| 女人十人毛片免费观看3o分钟| 高清在线视频一区二区三区| 国产免费一级a男人的天堂| 伦理电影免费视频| 免费看av在线观看网站| 国产黄片美女视频| 欧美精品国产亚洲| 麻豆乱淫一区二区| 欧美激情极品国产一区二区三区 | 亚洲av欧美aⅴ国产| 伦理电影免费视频| 欧美 日韩 精品 国产| 成人亚洲欧美一区二区av| 日本猛色少妇xxxxx猛交久久| 中文字幕久久专区| 欧美3d第一页| 黄色日韩在线| 成人国产av品久久久| 美女脱内裤让男人舔精品视频| 亚洲自偷自拍三级| 亚洲精品第二区| 大香蕉久久网| 亚洲av福利一区| 性色av一级| 熟女电影av网| 伦理电影免费视频| 久久久久久久国产电影| 高清午夜精品一区二区三区| 国产精品三级大全| 亚洲国产成人一精品久久久| av.在线天堂| 亚洲精品国产成人久久av| 亚洲国产成人一精品久久久| 免费观看av网站的网址| 成人影院久久| 高清在线视频一区二区三区| 亚洲av欧美aⅴ国产| 黄色怎么调成土黄色| 高清在线视频一区二区三区| 日韩一本色道免费dvd| 国产精品久久久久久av不卡| 老师上课跳d突然被开到最大视频| 免费观看在线日韩| 亚洲国产精品一区三区| 新久久久久国产一级毛片| 国产爱豆传媒在线观看| 久久久久久久精品精品| 乱码一卡2卡4卡精品| 2022亚洲国产成人精品| 久久午夜福利片| 久久毛片免费看一区二区三区| 国产精品av视频在线免费观看| 国产 精品1| 成人免费观看视频高清| 亚洲国产欧美人成| 欧美成人午夜免费资源| 91午夜精品亚洲一区二区三区| 六月丁香七月| 国产成人午夜福利电影在线观看| 搡女人真爽免费视频火全软件| 国产成人aa在线观看| 老司机影院毛片| 青春草亚洲视频在线观看| 尤物成人国产欧美一区二区三区| 亚洲精品视频女| 国产在视频线精品| 国产成人一区二区在线| 亚洲自偷自拍三级| 国产v大片淫在线免费观看| 成年美女黄网站色视频大全免费 | 国产一区有黄有色的免费视频| 日本av免费视频播放| 人妻少妇偷人精品九色| 九草在线视频观看| 男女边吃奶边做爰视频| 久久精品国产鲁丝片午夜精品| 日日啪夜夜撸| 热99国产精品久久久久久7| 国产在线视频一区二区| 国产国拍精品亚洲av在线观看| 国产爽快片一区二区三区| 亚洲精品自拍成人| 一本久久精品| 日韩,欧美,国产一区二区三区| 国产成人精品婷婷| 欧美成人午夜免费资源| 99久国产av精品国产电影| 日日摸夜夜添夜夜添av毛片| 日本vs欧美在线观看视频 | 亚洲av免费高清在线观看| 人妻制服诱惑在线中文字幕| 男女国产视频网站| 丝袜喷水一区| 最新中文字幕久久久久| 性色av一级| 午夜福利高清视频| videossex国产| 内射极品少妇av片p| 色哟哟·www| 国产色爽女视频免费观看| 日韩伦理黄色片| 国产亚洲一区二区精品| 日韩欧美 国产精品| 性色avwww在线观看| 九草在线视频观看| 亚洲av电影在线观看一区二区三区| 国产欧美日韩精品一区二区| 日本vs欧美在线观看视频 | 乱系列少妇在线播放| 一级毛片黄色毛片免费观看视频| 久久6这里有精品| 国产毛片在线视频| 精品一区在线观看国产| 一级毛片aaaaaa免费看小| av免费在线看不卡| 蜜桃亚洲精品一区二区三区| 亚洲激情五月婷婷啪啪| 五月天丁香电影| 欧美三级亚洲精品| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩东京热| 免费大片18禁| 熟女av电影| 亚洲精品国产色婷婷电影| 精品视频人人做人人爽| 国产精品人妻久久久影院| 尤物成人国产欧美一区二区三区| 免费观看a级毛片全部| 久久99热6这里只有精品| 22中文网久久字幕| 热99国产精品久久久久久7| 国产高清不卡午夜福利| 久久 成人 亚洲| 久久鲁丝午夜福利片| 美女xxoo啪啪120秒动态图| 国产午夜精品一二区理论片| 国产国拍精品亚洲av在线观看| 久久久久久久亚洲中文字幕| 18禁裸乳无遮挡动漫免费视频| 中文字幕人妻熟人妻熟丝袜美| av国产免费在线观看| 热re99久久精品国产66热6| 午夜福利在线在线| 秋霞在线观看毛片| 在线观看av片永久免费下载| 国产 精品1| 亚洲av欧美aⅴ国产| av女优亚洲男人天堂| 国产女主播在线喷水免费视频网站| 99久久精品国产国产毛片| 国产有黄有色有爽视频| 97热精品久久久久久|