• <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.

    成人国产麻豆网| 熟女人妻精品中文字幕| 欧美黑人欧美精品刺激| 99久久精品一区二区三区| 99九九线精品视频在线观看视频| 亚洲精品日韩av片在线观看| 日本-黄色视频高清免费观看| 国产亚洲精品久久久久久毛片| 国产免费一级a男人的天堂| 国产探花在线观看一区二区| www.www免费av| 欧美一区二区国产精品久久精品| 色综合亚洲欧美另类图片| 亚洲国产精品久久男人天堂| 久久精品久久久久久噜噜老黄 | 特大巨黑吊av在线直播| 男女视频在线观看网站免费| 97超视频在线观看视频| 国产乱人视频| 国内精品一区二区在线观看| 老师上课跳d突然被开到最大视频| 久9热在线精品视频| 成人特级av手机在线观看| 村上凉子中文字幕在线| 精品福利观看| 免费观看精品视频网站| 成人鲁丝片一二三区免费| 国产精品久久久久久av不卡| 亚洲国产欧洲综合997久久,| 久久欧美精品欧美久久欧美| 欧美一级a爱片免费观看看| 狠狠狠狠99中文字幕| 永久网站在线| 九九在线视频观看精品| 特大巨黑吊av在线直播| 美女黄网站色视频| 午夜免费男女啪啪视频观看 | 九九在线视频观看精品| 国产色爽女视频免费观看| 免费在线观看成人毛片| 欧美成人性av电影在线观看| 免费av毛片视频| av天堂在线播放| 在线观看免费视频日本深夜| 看黄色毛片网站| 中文字幕免费在线视频6| 能在线免费观看的黄片| 久久天躁狠狠躁夜夜2o2o| 精品久久久久久,| 成人鲁丝片一二三区免费| 99热这里只有精品一区| 日韩一区二区视频免费看| 成人一区二区视频在线观看| 一本久久中文字幕| 欧美色视频一区免费| 亚洲不卡免费看| 国产激情偷乱视频一区二区| 精品不卡国产一区二区三区| 波多野结衣高清无吗| 日韩中字成人| 波多野结衣高清作品| 在线观看av片永久免费下载| 国产一级毛片七仙女欲春2| 久久久久久伊人网av| 久久久久久久久久黄片| 久久久久久伊人网av| 国产精品女同一区二区软件 | 最近中文字幕高清免费大全6 | 久久久久久久久中文| 日韩精品中文字幕看吧| 国产免费一级a男人的天堂| 国产三级在线视频| 麻豆久久精品国产亚洲av| 精品久久久久久,| 神马国产精品三级电影在线观看| 欧美激情国产日韩精品一区| 久久精品国产99精品国产亚洲性色| 国内精品美女久久久久久| 中文资源天堂在线| 国产人妻一区二区三区在| 日韩欧美精品免费久久| 18+在线观看网站| 干丝袜人妻中文字幕| 性色avwww在线观看| 有码 亚洲区| 极品教师在线免费播放| 亚洲性夜色夜夜综合| 又粗又爽又猛毛片免费看| 精品不卡国产一区二区三区| 一个人看视频在线观看www免费| 成人综合一区亚洲| 蜜桃久久精品国产亚洲av| 性欧美人与动物交配| 黄色视频,在线免费观看| 日韩中字成人| 亚洲美女视频黄频| 亚洲av成人av| 黄色女人牲交| 午夜福利在线在线| 女同久久另类99精品国产91| 一个人看视频在线观看www免费| 亚洲色图av天堂| 亚洲在线观看片| 成年人黄色毛片网站| 啪啪无遮挡十八禁网站| 日韩强制内射视频| 日本精品一区二区三区蜜桃| 亚洲美女视频黄频| 欧美区成人在线视频| 又黄又爽又刺激的免费视频.| 国产 一区 欧美 日韩| av专区在线播放| 亚洲人成网站在线播| 99在线视频只有这里精品首页| 黄色视频,在线免费观看| 深爱激情五月婷婷| 久久亚洲真实| 国产成人福利小说| 久久精品国产亚洲网站| 丰满乱子伦码专区| 亚洲熟妇中文字幕五十中出| 在线免费观看的www视频| 亚洲最大成人中文| 中国美白少妇内射xxxbb| 精品不卡国产一区二区三区| av在线观看视频网站免费| 少妇被粗大猛烈的视频| 国产av麻豆久久久久久久| 女人被狂操c到高潮| 午夜爱爱视频在线播放| 亚洲人成网站在线播| 老女人水多毛片| 国产成人一区二区在线| 国产成人一区二区在线| 亚洲国产精品久久男人天堂| 99热精品在线国产| 久久九九热精品免费| 毛片一级片免费看久久久久 | АⅤ资源中文在线天堂| a级毛片免费高清观看在线播放| 国产欧美日韩精品亚洲av| 长腿黑丝高跟| 欧美黑人巨大hd| 99国产精品一区二区蜜桃av| 欧美xxxx黑人xx丫x性爽| 亚洲欧美日韩高清专用| 日本撒尿小便嘘嘘汇集6| 99热这里只有精品一区| 国产高清三级在线| 国产成人a区在线观看| 免费人成在线观看视频色| 又黄又爽又免费观看的视频| 久久精品人妻少妇| 欧美一区二区精品小视频在线| 亚洲真实伦在线观看| 精品国产三级普通话版| 国产91精品成人一区二区三区| 国产成年人精品一区二区| 国产精品日韩av在线免费观看| 搡老妇女老女人老熟妇| 69人妻影院| 桃红色精品国产亚洲av| 两人在一起打扑克的视频| 热99re8久久精品国产| 久久精品国产99精品国产亚洲性色| 男女那种视频在线观看| 午夜久久久久精精品| 国内精品久久久久精免费| 欧美激情久久久久久爽电影| 国产视频一区二区在线看| 精品不卡国产一区二区三区| 国产成人一区二区在线| 国产成人一区二区在线| 免费av不卡在线播放| av在线观看视频网站免费| 搡老岳熟女国产| 亚洲精品色激情综合| 亚洲最大成人av| 色精品久久人妻99蜜桃| 男人和女人高潮做爰伦理| videossex国产| 在现免费观看毛片| 日本撒尿小便嘘嘘汇集6| 国产综合懂色| 日韩,欧美,国产一区二区三区 | 一卡2卡三卡四卡精品乱码亚洲| 两人在一起打扑克的视频| 欧美一级a爱片免费观看看| netflix在线观看网站| 波多野结衣巨乳人妻| 尾随美女入室| 精品一区二区三区视频在线观看免费| 午夜免费成人在线视频| 99热网站在线观看| 美女被艹到高潮喷水动态| 又黄又爽又刺激的免费视频.| 淫妇啪啪啪对白视频| 99热这里只有是精品50| 日韩一区二区视频免费看| 国产高清三级在线| 大又大粗又爽又黄少妇毛片口| 中国美白少妇内射xxxbb| 日韩欧美国产一区二区入口| 十八禁网站免费在线| 欧美日韩黄片免| 中文资源天堂在线| 成人一区二区视频在线观看| 久久久久国产精品人妻aⅴ院| 白带黄色成豆腐渣| 婷婷精品国产亚洲av| 亚洲电影在线观看av| 一区二区三区高清视频在线| 婷婷精品国产亚洲av在线| 精华霜和精华液先用哪个| 亚洲美女视频黄频| 精品日产1卡2卡| 国模一区二区三区四区视频| 久久九九热精品免费| x7x7x7水蜜桃| 欧美高清性xxxxhd video| 人妻丰满熟妇av一区二区三区| 在线播放无遮挡| 亚洲国产精品合色在线| 成人特级黄色片久久久久久久| av在线观看视频网站免费| 欧美日韩瑟瑟在线播放| 精品久久久久久久人妻蜜臀av| 精品福利观看| 精品人妻视频免费看| 天堂动漫精品| 午夜a级毛片| 国产成年人精品一区二区| 午夜老司机福利剧场| 啦啦啦啦在线视频资源| 国产欧美日韩精品一区二区| 国产视频一区二区在线看| 搡老岳熟女国产| 听说在线观看完整版免费高清| 男人和女人高潮做爰伦理| 老师上课跳d突然被开到最大视频| 黄色一级大片看看| 日本色播在线视频| 亚洲精品国产成人久久av| 久9热在线精品视频| 国产 一区 欧美 日韩| 精品国内亚洲2022精品成人| 成人国产综合亚洲| 自拍偷自拍亚洲精品老妇| 成人特级av手机在线观看| 少妇人妻精品综合一区二区 | 欧美精品啪啪一区二区三区| 在线观看av片永久免费下载| 深夜精品福利| 丰满乱子伦码专区| 简卡轻食公司| 看片在线看免费视频| 一级黄色大片毛片| eeuss影院久久| 国产伦人伦偷精品视频| 色视频www国产| 国产色婷婷99| av.在线天堂| 日韩精品有码人妻一区| 搡老妇女老女人老熟妇| 亚洲一区高清亚洲精品| 国产成人a区在线观看| 欧美日韩黄片免| 亚洲av免费在线观看| 我要看日韩黄色一级片| 1024手机看黄色片| 亚洲久久久久久中文字幕| 国产精品不卡视频一区二区| 少妇高潮的动态图| 人妻制服诱惑在线中文字幕| 国产男靠女视频免费网站| 成年女人永久免费观看视频| 99久久精品国产国产毛片| 久久草成人影院| 国产精品野战在线观看| 亚洲,欧美,日韩| 欧美激情久久久久久爽电影| 日韩欧美精品v在线| av视频在线观看入口| 国内久久婷婷六月综合欲色啪| 99久久精品一区二区三区| 国产女主播在线喷水免费视频网站 | 色精品久久人妻99蜜桃| 99国产精品一区二区蜜桃av| 午夜免费激情av| 国产高清三级在线| 亚洲电影在线观看av| 国产精品免费一区二区三区在线| 听说在线观看完整版免费高清| 在线免费十八禁| 九九爱精品视频在线观看| 亚洲色图av天堂| 成人三级黄色视频| 搡老妇女老女人老熟妇| 一卡2卡三卡四卡精品乱码亚洲| 国内精品一区二区在线观看| 99久久久亚洲精品蜜臀av| 精品久久久久久久久av| 免费看美女性在线毛片视频| aaaaa片日本免费| 日韩国内少妇激情av| 男人舔奶头视频| 亚洲无线在线观看| 极品教师在线视频| 亚洲午夜理论影院| 一级a爱片免费观看的视频| 天天躁日日操中文字幕| 午夜免费激情av| 高清在线国产一区| 大又大粗又爽又黄少妇毛片口| 在线天堂最新版资源| 亚洲av美国av| 免费搜索国产男女视频| 欧美另类亚洲清纯唯美| 赤兔流量卡办理| 亚洲va日本ⅴa欧美va伊人久久| 午夜免费成人在线视频| 99九九线精品视频在线观看视频| 国产日本99.免费观看| 免费一级毛片在线播放高清视频| 身体一侧抽搐| 亚洲av美国av| 国产成人影院久久av| 在线播放国产精品三级| 亚洲va日本ⅴa欧美va伊人久久| 一a级毛片在线观看| 内射极品少妇av片p| 亚洲av二区三区四区| 日本精品一区二区三区蜜桃| 99久久精品国产国产毛片| 亚洲在线自拍视频| 黄色配什么色好看| 自拍偷自拍亚洲精品老妇| 五月玫瑰六月丁香| 国产真实乱freesex| 赤兔流量卡办理| 亚洲国产精品sss在线观看| 国产男人的电影天堂91| 真人一进一出gif抽搐免费| 精品99又大又爽又粗少妇毛片 | 亚洲最大成人av| 精品99又大又爽又粗少妇毛片 | 久久亚洲真实| 欧美高清性xxxxhd video| 免费av不卡在线播放| 久久中文看片网| 精品人妻熟女av久视频| 国产亚洲精品综合一区在线观看| 国产高清视频在线播放一区| 露出奶头的视频| 女人被狂操c到高潮| 国产白丝娇喘喷水9色精品| 一a级毛片在线观看| 亚洲aⅴ乱码一区二区在线播放| 日韩高清综合在线| 免费黄网站久久成人精品| 麻豆一二三区av精品| 自拍偷自拍亚洲精品老妇| 嫩草影院入口| 天堂网av新在线| 人人妻,人人澡人人爽秒播| 精品免费久久久久久久清纯| 亚洲国产精品久久男人天堂| 免费在线观看影片大全网站| 听说在线观看完整版免费高清| 国产69精品久久久久777片| 日韩一区二区视频免费看| 麻豆久久精品国产亚洲av| 97人妻精品一区二区三区麻豆| 色精品久久人妻99蜜桃| 免费人成在线观看视频色| 男人舔女人下体高潮全视频| 国产亚洲av嫩草精品影院| 免费人成在线观看视频色| 国产精品无大码| 欧美成人a在线观看| 老女人水多毛片| 亚洲精品亚洲一区二区| 欧美性猛交黑人性爽| 99久国产av精品| 999久久久精品免费观看国产| 九九热线精品视视频播放| 亚洲精品乱码久久久v下载方式| 全区人妻精品视频| 91久久精品国产一区二区三区| 成年免费大片在线观看| 日韩在线高清观看一区二区三区 | 亚洲美女视频黄频| 俺也久久电影网| 久久精品国产亚洲av天美| 看十八女毛片水多多多| 99热只有精品国产| 在线看三级毛片| 免费无遮挡裸体视频| 变态另类丝袜制服| 久久九九热精品免费| 在线播放无遮挡| 欧美+亚洲+日韩+国产| 极品教师在线免费播放| 日韩亚洲欧美综合| 欧美最黄视频在线播放免费| 婷婷精品国产亚洲av在线| 丝袜美腿在线中文| 欧美高清性xxxxhd video| 欧美一区二区精品小视频在线| 欧美性猛交╳xxx乱大交人| 精品人妻熟女av久视频| 91狼人影院| 深夜精品福利| 国产69精品久久久久777片| 日本成人三级电影网站| 国产男人的电影天堂91| 成人二区视频| 91久久精品电影网| 欧美性感艳星| 中文字幕精品亚洲无线码一区| 亚洲国产日韩欧美精品在线观看| 国产精品福利在线免费观看| 久久亚洲真实| 国产 一区 欧美 日韩| 99国产极品粉嫩在线观看| 日本 欧美在线| 欧美一级a爱片免费观看看| 欧美日本亚洲视频在线播放| 久久九九热精品免费| 搡老熟女国产l中国老女人| 99久久精品热视频| 麻豆精品久久久久久蜜桃| 在线观看免费视频日本深夜| 在线观看美女被高潮喷水网站| 免费观看精品视频网站| 亚洲精品成人久久久久久| 男人和女人高潮做爰伦理| 国产精品一区二区三区四区免费观看 | 麻豆一二三区av精品| 少妇丰满av| 国产美女午夜福利| 国产高清不卡午夜福利| 国产国拍精品亚洲av在线观看| 成人性生交大片免费视频hd| 免费看光身美女| 十八禁网站免费在线| 直男gayav资源| 亚洲精品日韩av片在线观看| 美女 人体艺术 gogo| 欧美日韩黄片免| 嫩草影院新地址| 日日摸夜夜添夜夜添小说| 草草在线视频免费看| 国内精品久久久久久久电影| 国产免费一级a男人的天堂| 亚洲五月天丁香| 在线国产一区二区在线| 亚洲av.av天堂| 成年女人毛片免费观看观看9| 999久久久精品免费观看国产| 精品久久久久久,| a级毛片a级免费在线| 精品99又大又爽又粗少妇毛片 | 国产欧美日韩一区二区精品| 在线a可以看的网站| 噜噜噜噜噜久久久久久91| 观看免费一级毛片| 搡老岳熟女国产| av女优亚洲男人天堂| 国产亚洲精品久久久久久毛片| 嫩草影院新地址| 久久久久久久久中文| 十八禁网站免费在线| 又黄又爽又刺激的免费视频.| 搡女人真爽免费视频火全软件 | 日韩强制内射视频| 老女人水多毛片| 日韩欧美一区二区三区在线观看| 91精品国产九色| 精品人妻1区二区| 男人和女人高潮做爰伦理| 三级男女做爰猛烈吃奶摸视频| а√天堂www在线а√下载| 看黄色毛片网站| 亚洲内射少妇av| 国产又黄又爽又无遮挡在线| 麻豆久久精品国产亚洲av| 校园人妻丝袜中文字幕| 国产成人a区在线观看| 亚洲美女视频黄频| 91麻豆av在线| 村上凉子中文字幕在线| 午夜激情福利司机影院| 亚洲一区二区三区色噜噜| 午夜福利在线观看吧| 天天一区二区日本电影三级| 久久久久久大精品| av黄色大香蕉| 男人舔女人下体高潮全视频| 欧美丝袜亚洲另类 | 99九九线精品视频在线观看视频| 啦啦啦观看免费观看视频高清| 日本 av在线| 麻豆精品久久久久久蜜桃| 欧美一区二区精品小视频在线| 日本欧美国产在线视频| 在线天堂最新版资源| videossex国产| 欧美绝顶高潮抽搐喷水| 黄色视频,在线免费观看| 人妻丰满熟妇av一区二区三区| 国产男人的电影天堂91| 国内精品久久久久精免费| 91麻豆精品激情在线观看国产| 老司机福利观看| 亚洲美女黄片视频| 国产麻豆成人av免费视频| 国产精品亚洲一级av第二区| 男女视频在线观看网站免费| 日韩,欧美,国产一区二区三区 | 校园人妻丝袜中文字幕| 91狼人影院| 欧美性猛交╳xxx乱大交人| 在线观看舔阴道视频| 偷拍熟女少妇极品色| 91狼人影院| 国内毛片毛片毛片毛片毛片| 韩国av一区二区三区四区| 我要搜黄色片| 亚洲欧美日韩东京热| 国产精品女同一区二区软件 | 国国产精品蜜臀av免费| 日日夜夜操网爽| 欧美日韩中文字幕国产精品一区二区三区| 国产精品电影一区二区三区| 免费高清视频大片| 麻豆av噜噜一区二区三区| 身体一侧抽搐| 亚洲最大成人手机在线| 一级黄色大片毛片| 最近最新免费中文字幕在线| 伊人久久精品亚洲午夜| 淫秽高清视频在线观看| 中文字幕久久专区| АⅤ资源中文在线天堂| 精品久久久久久久人妻蜜臀av| 色综合站精品国产| 少妇的逼水好多| 久久久久久大精品| 天堂av国产一区二区熟女人妻| 五月玫瑰六月丁香| 中文字幕久久专区| 琪琪午夜伦伦电影理论片6080| 午夜免费男女啪啪视频观看 | 国产精品精品国产色婷婷| 中文资源天堂在线| 制服丝袜大香蕉在线| 国产精品乱码一区二三区的特点| 欧美激情在线99| 小说图片视频综合网站| 久久久久九九精品影院| 少妇高潮的动态图| 亚洲欧美日韩卡通动漫| 看免费成人av毛片| 99久久无色码亚洲精品果冻| 又紧又爽又黄一区二区| 久久九九热精品免费| 国产高清视频在线播放一区| 国产成人av教育| 亚洲乱码一区二区免费版| 亚洲av不卡在线观看| 久久婷婷人人爽人人干人人爱| 亚洲天堂国产精品一区在线| 99热只有精品国产| 国产精品野战在线观看| 最近视频中文字幕2019在线8| 国产亚洲精品综合一区在线观看| 国产精品自产拍在线观看55亚洲| 亚洲国产日韩欧美精品在线观看| 超碰av人人做人人爽久久| 国产精品人妻久久久久久| 琪琪午夜伦伦电影理论片6080| 日韩一区二区视频免费看| 综合色av麻豆| 亚洲国产日韩欧美精品在线观看| 观看美女的网站| 波野结衣二区三区在线| 亚洲四区av| 一区二区三区高清视频在线| 午夜激情福利司机影院| 成人三级黄色视频| 亚洲第一电影网av| 91麻豆精品激情在线观看国产| 成人综合一区亚洲| 丰满乱子伦码专区| 日日摸夜夜添夜夜添av毛片 | av在线亚洲专区| 国产精品久久久久久av不卡| 日本黄色视频三级网站网址| 国产精品嫩草影院av在线观看 | 国国产精品蜜臀av免费| 国产精品亚洲美女久久久| 色播亚洲综合网| 精品一区二区三区视频在线观看免费| 午夜激情欧美在线| 亚洲精品成人久久久久久| 级片在线观看| 国产一区二区在线av高清观看| 搡老熟女国产l中国老女人| 少妇猛男粗大的猛烈进出视频 | 国产精品三级大全| 在线播放无遮挡| 亚洲黑人精品在线| 日本黄大片高清| 两人在一起打扑克的视频|