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

    Color-kinematics duality for Sudakov form factor in non-supersymmetric pure Yang–Mills theory

    2022-08-02 02:47:50ZeyuLiGangYangandJinxuanZhang
    Communications in Theoretical Physics 2022年6期

    Zeyu Li,Gang Yangand Jinxuan Zhang

    1 CAS Key Laboratory of Theoretical Physics,Institute of Theoretical Physics,Chinese Academy of Sciences,Beijing 100190,China

    2 School of Physical Sciences,University of Chinese Academy of Sciences,No.19A Yuquan Road,Beijing 100049,China

    Abstract We study the duality between color and kinematics for the Sudakov form factors of tr(F2) in non-supersymmetric pure Yang-Mills theory.We construct the integrands that manifest the color-kinematics duality up to two loops.The resulting numerators are given in terms of Lorentz products of momenta and polarization vectors,which have the same powers of loop momenta as that from the Feynman rules.The integrands are checked by d-dimensional unitarity cuts and are valid in any dimension.We find that massless-bubble and tadpole topologies are needed at two loops to realize the color-kinematics duality.Interestingly,the two-loop solution contains a large number of free parameters suggesting the duality may hold at higher loop orders.

    Keywords:scattering amplitudes,form factors,color-kinematics duality,unitarity method

    1.Introduction

    Significant progress has been made in the study of scattering amplitudes in the past thirty years,see e.g.[1–10]for review.These studies have not only important phenomenological applications but also have uncovered various new structures that are not obvious at all from the traditional Feynman diagram method.In this paper,we focus on one of such structures,the so-called color-kinematics (CK) duality discovered by Bern,Carrasco,and Johansson [11,12].The duality conjectures that there exists a trivalent Feynman-like diagrammatic representation of amplitudes in which the kinematic numerators satisfy the same algebraic relations as the color factors associated with the same graphs.Importantly,this duality indicates a deep connection between the gauge and gravity theories:via double copy[12,13],one can obtain gravitational amplitudes directly from gauge amplitudes once the latter are organized to respect the CK duality,which is also closely related to the KLT relations [14]and the CHY formula [15,16].See [17]for an extensive review of the duality and its applications.

    Although the CK duality has been understood at tree level using e.g.string theory monodromy relations[18,19]or gauge theory recursion relations [20],it remains a conjecture at loop level.In supersymmetric theories,the CK duality has been found to exist at high loop orders.For example,in the maximally supersymmetric N=4 super-Yang–Mills (SYM)theory,the four-gluon amplitude has been found to preserve the duality up to four loops [12,21]and for five-point amplitudes up to three loops[22].For the form factors of the stress-tensor multiplet inN=4 SYM,the CK-dual representation has been obtained up to five loops for the two-point Sudakov form factor [23,24]and up to four loops for the three-point form factor [25–27].High-loop constructions for half-maximally supersymmetric theories were also studied in[28–33].In the case of non-supersymmetric theories,some examples at one-and two-loop orders were found [34–42].

    Despite this much progress,it is far from clear to which extent the duality holds at general loop orders.For example,for the four-point amplitude inN=4 SYM,although a simple four-loop CK-dual solution has been obtained for a while [21],it has proven difficult to construct the five-loop integrand which could manifest full CK duality [43,44].In the non-supersymmetric gauge theories,the CK-dual representations of loop integrands are generally much harder to obtain.For example,for the four-gluon amplitude in pure YM theory,such a representation is so far only obtained for the case with identical four-dimensional external helicities at two loops [35].For the two-loop five-gluon amplitudes in pure YM theory with identical helicities,it was found that the numerators with twelve powers of loop momenta have to be used to realize the duality,which is far beyond the seven powers expected from Feynman diagrams [37].

    Figure 1.Trivalent graphs of the four-point tree amplitude.

    In this paper,we explore this duality in the non-supersymmetric pure Yang–Mills theory by considering the Sudakov form factor of the operator tr(F2) in the theory,which is defined as

    where piare on-shell massless momenta of external gluons and q=p1+p2is the off-shell momenta associated with the operator.Sudakov form factor plays a central role in the study of IR divergences of gauge theories [45–48].In pure YM theory,the two-loop Sudakov form factor of tr(F2)was firstly obtained to all orders in the dimensional regularization regulator ∈=(4-d)/2 using Feynman diagrams in [49](see also [50]).In this work,we will reproduce this result using the CK duality and unitarity method,and we would like to stress that our main concern here is the structure of the form factor integrand.

    We will start with constructing an ansatz of integrand by using CK duality,and then we solve the ansatz by applying the unitarity-cut method[51–53].Once the result satisfies a spanning set of cut constraints,it is a physically correct one.In this way,we obtain the integrand solutions up to two loops which not only manifest all dual Jacobi relations but also satisfy all possible cuts.Compared to the results in N=4 SYM [23],there are several important new features for the pure YM Sudakov form factor.First,the d-dimensional unitarity cuts are necessary to obtain the full results,and the resulting integrands also depend on the dimensional regulator ∈=(4-d)/2.Second,massless-bubble and tadpole topologies are needed in order to satisfy CK duality at two-loop order.In contrast,the two-loopN=4 Sudakov form factor is much simpler,which contains only a planar-ladder and a cross-ladder topology.Our pure YM results are given in terms of Lorentz products of momenta and polarization vectors,and they are valid in arbitrary dimensions.The CK-dual numerators have the same powers of loop momentum as expected from Feynman diagrams.Moreover,we find that the two-loop solution space still contains a large number of free parameters.These imply that the CK duality may promisingly hold for higher loop form factors in pure YM theory.

    The rest of the paper is organized as follows.In section 2,we construct the ansatz of the form factor integrands by imposing that they manifest the CK duality.In section 3,we apply unitarity cuts to solve for the coefficients in the ansatz and also study the constraints related to the tadpoles and massless bubbles.In section 4,we perform integral reductions and show that all parameters cancel after integral reductions.We also discuss the difference between the integrated results of the d-dimensional and 4-dimensional integrands.A summary and outlook are given in section 5.Several appendices provide some further details.In appendix A,we give the full set of two-loop dual Jacobi relations.The explicit expressions of the two-loop master numerators are given in appendix B.In appendix C,we provide explicit expressions for the integration-by-part(IBP)master integrals and the master coefficients.

    2.CK duality and ansatz construction

    In this section,we apply CK duality to construct the integrand ansatz for the Sudakov form factor of tr(F2)in pure YM theory.After giving a brief review of the duality,we will demonstrate how to generate the ansatz in the one-loop case step by step.Next,we perform a similar construction at two loops.

    2.1.Review of CK duality

    In SU(Nc) gauge theory,the structure constant is defined by

    in which Taare SU(Nc) generators and are normalized by tr(TaT b)=δab.By definition structure constants (2.1)satisfy the Jacobi relation:

    As depicted in section 3 of [11],one can represent the fullcolor four-point tree amplitude as

    which correspond to s,t,u-channels shown in figure 1.The color factors satisfy the Jacobi relation

    Figure 2.Loop diagrams related by Jacobi relation.

    It turns out that the numerators also satisfy a similar linear relation as:

    This correspondence is called CK duality.The relation (2.5)for the numerators is named as dual Jacobi relation or CK relation.

    One can impose similar CK duality at loop level as follows.For a certain loop diagram with color-dressed vertices,if we replace a four-point sub-diagram of it with the other two kinds of four-point diagrams in figure 1,we will obtain two new loop diagrams.This process is shown in figure 2.Note that except for the four-point sub-digram,the rest of the three diagrams is the same,so their color factors take the form:

    Obviously they satisfy the Jacobi relation:

    For these diagrams,we can postulate a dual numerator relation as

    As mentioned in the introduction,it is still a conjecture whether such a requirement can be fulfilled at general loop level.Our strategy is to assume such duality can be achieved and use this structure to construct an ansatz for the loop integrand.Then we check the ansatz by physical unitarity-cut constraints.Once the CK-dual ansatz has a solution consistent in all possible unitarity cuts,this then shows that the CK duality can be indeed realized at this order.Practically,this conjecture (as long as it works) can also greatly simplify the integrand construction.

    We will focus on the Sudakov form factor of tr(F2) in pure YM theory,as given in (1.1).At tree level,the form factor is given as

    with color factor

    and the color-stripped tree form factor is

    An l-loop full-color Sudakov form factor takes the general form as

    Figure 3.One-loop cubic graphs.

    with the meaning of each term explained as follows.The summation over σ2means taking the permutation of external legs into account.The summation over Γimeans to sum over all possible trivalent graphs.Si,coming from overcounting of contributions of the ith graph,is the symmetry factor of it.Ciis the color factor of the ith graph,given as products of structure constants dressed on every vertex.1/Di,adenotes the ath propagator of the ith graph.Finally,Niare the kinematic numerators,which are the main goal of the construction.The reader can find more details of the general strategy in e.g.[21,54,55].

    2.2.One-loop numerator ansatz

    In this subsection,we will construct the ansatz for the oneloop Sudakov form factor.We start with generating all oneloop cubic graphs.There are two topologies to consider,as shown in figure 3.

    We write the one-loop form factor as

    As in figure 3,we use Pito denote the propagator momenta,which will also be used as color indices for each edge.We choose the convention of color factors such that the three color indices of eachfollow a clockwise direction for each vertex shown in figure 3.This convention will be used throughout the paper.Explicitly,the color factors of the graphs in figure 3 are given as

    We would like to point out that C2is zero since the product ofandvanishes.However,as we will discuss below,it is important to keep it in the form of (2.14) for getting the correct sign factors for the dual Jacobi relations.The two kinematic numerators are defined as

    Figure 4.Map from CK operation results to convention.

    which are the main functions to be solved.

    2.2.1.Generating CK relations.In this part,we use the oneloop form factor to describe in detail how to generate dual Jacobi relations that will provide linear relations for the numerators.

    We start with picking up a propagator in the trivalent graph which will play the role of s-channel propagator in figure 2.In order to apply the dual Jacobi relation,such a propagator should not be connected to the operator leg because the operator vertex is not dressed with a structure constant.For the first topology in figure 3,the only propagator which can be considered is denoted by the red color in figure 4.The Jacobi relation leads to the t-channel and u-channel graphs in figure 4.At one-loop this will be the only dual Jacobi relation one needs to consider.

    Color Jacobi relation—We first consider the color Jacobi relation which will be important to determine the sign factors in the dual Jacobi relation for the kinematic numerators.We set Csas C1in(2.14)and define the other two color factors as

    which satisfy

    To extract the correct sign factors for the dual Jacobi relation,we can map the color indices in Ct/uto our convention in the form of (2.14).As shown in figure 4,the t-and u-channel topologies are isomorphic to the bubble and triangle topologies respectively.We define the map of color indices as φt/ushown in figure 4 which are

    As a result,we will obtain

    Note that the equality after Ct/uholds because φt/uonly rename color indices.Then we can directly compare mapped Ct/uwith Ci.The two plus signs are the results of the comparison.We have

    The same signs will be used for the kinematic numerators.We emphasize again that although Ct=C2=0 is correct in value,we need to use its form as in (2.19) to get the wanted sign.See also [26]for further examples and discussion on this point.

    Dual Jacobi relation—Next we derive the dual Jacobi relation for the numerator functions (2.15).Mapping to the momenta in the t-and u-channel diagrams in figure 4,we have the dual Jacobi relation:

    where the two plus signs in the first equation are in accord with the signs in the color relation (2.20).

    Choose master numerator—We can represent N2by a linear combination of N1as

    We will choose N1as the master numerator.Any groups of numerators which can express all numerators through CK relations can be called a set of CK master numerators.Corresponding cubic graphs of these numerators are called CK master topology.In this one-loop case,there is a unique choice of master numerator as N1.However,in general the choice of master numerators is not unique,such as in the twoloop case discussed below.Usually,we tend to choose simple topologies (such as planar topologies and symmetric topologies) and also minimize the number of master topologies,so that the ansatz construction will be simplified.

    Figure 5.One-loop master graph and the labeling of its propagators.

    2.2.2.Constructing ansatz for the master numerator.We need to further construct an ansatz for the master numerator N1.The numerator is given as a polynomial in terms of following Lorentz product variables:

    where the momenta Pifor propagators are shown in figure 5.The complete basis is

    We impose the power counting condition following Feynman diagrams.The one-loop diagram contains two trivalent-vertices and an operator vertex.Each trivalent-vertex contributes mass dimension one to the numerator,and the operator vertex contributes mass dimension two.For the complete Feynman rule of the form factor,see e.g.[56].So numerator needs to be mass dimension four in total.In addition,the numerator must depend on each polarization vector linearly.After collecting all permitted monomials,the ansatz of N1can be written as

    There are 26 monomials in total.Notice that for completeness,we include terms containing D1,2,such asa7D1K13K23anda4D1D2E,which will reduce to massless bubbles or tadpole,and they are zero after integration.As we will see,the coefficients of these kinds of terms can not be determined by the unitarity cut,and one can also exclude such terms in the ansatz from the beginning.However,this is not possible in the two-loop case,which will be discussed in section 3.3.1.

    Figure 6.Feynman diagrams which have d dependence in the numerators.

    There is another important comment on the coefficients ai.From the Feynman diagram calculation,two Feynman graphs as shown in figure 6 can produce terms that are linear in the space-time dimension parameter d in the numerator,due to the contraction of metric=δμμ=d.Therefore,we expect that the coefficients will in general be degree-one polynomials of d as

    where aj,0and aj,1are pure rational numbers.

    2.2.3.Symmetrization of the ansatz.Here we would like to impose a further constraint on the ansatz such that each numerator should respect the symmetry property of the corresponding diagram.More concretely,symmetry conditions come from graph self-isomorphism and we require that CiNigoes back to itself after the symmetry transformations.

    Considering the first graph in figure 3,along the horizontal axis,there is a symmetry transformation:

    For the second graph in figure 3,it has two symmetries,by exchanging bubble propagators or exchanging external legs,and they correspond to the following transformations respectively:

    These symmetry transformations impose the following three constraints:

    As discussed below(2.19),we can not simply set C2to 0 but should use the form with color indices shown explicitly.In terms of the numerator functions,they are

    Figure 7.Two-loop cubic graphs without massless bubble and tadpole.

    Figure 8.Two-loop cubic graphs with massless bubble and tadpole.

    These conditions will solve for 9 coefficients in N1,leaving 17 coefficients unsolved.

    2.3.Two-loop numerator ansatz

    In this subsection,we go on to construct CK-dual ansatz for the two-loop Sudakov form factor.Since the procedure is similar to the one-loop case,we will be brief in the discussion and mainly focus on the new features at two loops.

    To begin with,we generate all possible cubic graphs.There are 17 different graphs which are shown in figure 7 and figure 8.In particular,eight of them contain massless bubble or tadpole sub-graphs,which are all collected in figure 8.The two-loop form factor can be given the following ansatz form:

    As mentioned,the convention of color factors is chosen such that the color indices of eachfollow a clockwise direction for the corresponding vertex in figure 7 and figure 8.The kinematic numerators are defined as functions Ni[l1,l2,p1,p2],where l1,2are labeled in figures.

    2.3.1.Generating CK relations.Following the strategy in section 2.2.1,we have generated all dual Jacobi relations in appendix A.As mentioned before,the choice of the master topologies is not unique.We will choose the numerators of the first two graphs in figure 7 as master numerators,which are both planar.

    Using a set of 15 dual Jacobi relations,one can obtain all other 15 numerators from the two master numerators.An explicit set of such relations can be given as

    Figure 9.CK operation corresponding to the third CK relation.

    Figure 10.Two-loop propagator basis.

    where Nirepresents Ni[l1,l2,p1,p2]for brevity.The first seven relations can generate other seven numerators in figure 7,while the other eight relations can be used to generate all numerators in figure 8.As an example,we demonstrate the relation (2.34) as diagrams in figure 9.

    2.3.2.Constructing ansatz for master numerators.Next,we construct the ansatz for master numerators N1and N2.The momenta for propagator basis are labeled as Piin figure 10.Define the two-loop Lorentz products as

    The numerators are given as polynomials in terms of following Lorentz product basis,which contains 15 elements:

    We use power counting of the Feynman diagram:there are 4 three-vertices and an operator vertex,we require the numerators to have mass dimension 6.The ansatz for each master numerator turns out to have 444 different monomials in the most general form:

    where Miare the monomials similar in the one-loop ansatz (2.25).

    As in the one-loop case,the numerators from Feynman diagrams can also depend on d.An inspection of Feynman diagrams shows that only two Feynman diagrams,as shown in figure 11,can produce terms proportional to d2,coming from metric contractionsThus we expect the coefficients to be in general degree-two polynomials in d as

    where ak,αand bk,αare pure rational numbers.

    d2-dependence—One can note that the two Feynman diagrams in figure 11 have special loop structures that are related to trivalent topologies (1),(6),(7) and (9) in figure 7.Thus one may expect that only the numerators of these four topologies {N1,N6,N7,N9} have d2dependence.Using dual Jacobi relations,one can find that there is a chain of relations(2.35),(2.36),and(2.38),such that the numerators N6,N7,N9can be generated by only using N1(which is one of the master numerators).This raises the question that if it is possible only to let N1contain d2terms and set d2terms of the second master numerator N2equal zero.If possible,this may be used to simplify the ansatz and reduce the number of free parameters.However,as we will see in section 3.3.1,this is in general not possible.The unitarity-cut constraints together with CK duality require us to include d2terms in N2.

    Figure 11.Feynman diagrams which have a contribution of d2.

    Figure 12.Symmetry of tadpole.

    Reduced tadpoles and massless bubbles—In the ansatz consisting of full 444 terms in (2.49),we include monomials that can reduce the maximal topologies to sub-topologies that contain tadpoles and massless bubbles.Such terms correspond to scaleless integrals and will not contribute to the final form factor result.Moreover,they are not detectable by unitarity cuts.Thus we may wonder if one can exclude these terms from the ansatz for the master numerators.If possible,this can also simplify the ansatz for master numerators from the beginning.However,as we will see in section 3.3.1,this is in general not possible.The unitarity-cut constraints together with CK duality require us to include such terms in the master numerators.

    2.3.3.Symmetrization of the ansatz.Finally,like the oneloop case discussed in section 2.2.3,we impose symmetry constraints on the numerators.Since the tadpole topologies appear in the two-loop case for the first time,we discuss their feature in more detail.Each tadpole subgraph has the flip symmetry as shown in figure 12.Under the flip symmetry,the momentum flowing in the tadpole will change its direction and the color factor will produce an additional minus sign.Thus,using the topology (10) in figure 8 as an example,one has the constraint on the numerator

    from which N10must be made up only by terms with odd powers of tadpole-loop momentum.We point out that this symmetry property is actually automatically satisfied if the tadpole numerators are generated through dual Jacobi relations,see section 3.3.2 for further discussion.

    It turns out to be that 508 parameters of N1and N2can be solved after using symmetry constraints,and there are 380 parameters left.We point out that after this step,the numerators also satisfy all dual Jacobi relations as listed in appendix A.

    3.Solving CK ansatz

    Having obtained the CK ansatz,now we apply the unitarity cuts to constrain the ansatz so that it provides the physical result.We stress that one needs to apply d-dimensional cuts to get the complete result since we are studying the pure YM theory.The feature that the ansatz numerators depend on the dimension parameter d also implies that d-dimensional cuts are necessary.

    3.1.Review of d-dimensional cut

    The central idea of the unitarity method [51–53]is that by setting internal propagators to be on-shell as

    the amplitude or form factor will be factorized as a product of lower-order amplitudes or form factors,such as

    The physical result must be consistent in all possible cut channels.Note that since there is no sub-leading color contribution for Sudakov form factor at one and two loops,it is enough to consider planar cuts.(See also [26]for discussion on the non-planar cuts.)

    To apply d-dimensional planar cuts,the building blocks are color-stripped tree amplitudes or form factors that can be calculated by Feynman diagrams.All the expressions appear as Lorentz products which are valid in d dimensions.When we multiply them together,we need to sum over all the helicities,and in d-dimensional cut,this operation corresponds to the contraction of the polarization vectors for the internal gluons by using the following rule:

    Figure 13.Double cut for one-loop two-point form factor.

    where the ξμis a light-like reference momenta and the result after summation should be independent of the choice of it.By matching the resulting tree product expression and the cut of the ansatz integrand,one can fix the coefficients in the ansatz.

    We checked that after the helicity sum,the ξ-dependent terms in the expression vanish.As another check,when we set d=4,the tree products match the result obtained by spinor helicity formalism.Compared to the 4-dimensional cut using spinor helicity formalism,d-dimensional cuts can capture all terms valid for general dimensions.

    3.2.One-loop solution

    For the one-loop case,there is only one cut we need to consider,which is the double cut in figure 13.The one-loop Sudakov form factor will be factorized as the products of the two-point tree form factor and the four-point tree amplitude:

    where the integral measure is the two-particle phase-space measure which is defined as

    As discussed in the previous section,after the symmetry constraints there are 17 monomial coefficients left in the ansatz.When we match the CK ansatz under the cut and the tree product as in (3.4),we can further fix 7 coefficients and 10 of them remain unsolved.These coefficients can not be fixed by symmetry and unitarity.We find that in this solution space,we can let the numerator of the second topology be zero,whose color factor is zero and will not contribute to the final result.When we apply this extra constraint,there will be only 6 free coefficients.

    We present the final numerator of the triangle topology by separating it into two parts as

    where N1,0contains no free parameters and all the monomials with six unfixed coefficients are collected in N1,1.They are explicitly given as

    In (3.7),except for the first two terms,all other terms are proportional to ∈=(4-d)/2 and they can not be fixed in the 4-dimensional cut.For the terms depending on free coefficients collected in(3.8),we can easily see that all of them are proportional to D1or D2,which will reduce topologies with massless bubble or tadpole.They can not be detected by unitarity cuts and thus they remain as free parameters.Since they are scaleless integrals,they vanish after integration.

    3.3.Two-loop solution

    In this subsection,we will constrain the two-loop CK ansatz.We will first apply unitarity cuts.As we will see,there is a large solution space.Next we will discuss possible further constraints to reduce the solution space from tadpoles and massless bubbles.Finally,we discuss some features of the solution space.

    3.3.1.Unitarity constraints.At two loops,it is necessary to consider several different cuts.A complete spanning set of planar cuts are shown in figure 14.We mention that cut (b)and cut(c)give the same constraint because of the cyclic and reflection symmetry of the 4-point tree amplitude.

    As discussed in section 2.3,after requiring all the numerators to have the same symmetry as their topologies,we are able to constrain a lot of parameters and there will be 380 coefficients left.After applying all the cuts,there are 235 coefficients that remain unsolved.All of the remaining parameters should vanish after the integral reduction,which we have checked to be so.Before considering further constraints,let us address the two questions we raised during the construction of ansatz in section 2.3.2.

    On d2-dependence—First we consider the question if one can set the d2terms in N2to be zero.In the unitarity calculation we find that only in the cut of figure 14(b),the tree product will provide non-zero d2terms,which is consistent with the Feynman diagram analysis.Topologies that contribute to cut-(b) are (1),(6),(7) and (9) in figure 7.However,in the solution space after unitarity cuts,we find that one can not let all d2terms be zero in N2.This is an interesting feature that the CK-dual solution has a different ddependence structure from that of Feynman diagrams since all Feynman diagrams which can give contributions to the topology of N2can not produce d2terms The origin of this different structure is that here we require the integrand to not only pass unitarity cuts but also satisfy the requirement of CK duality.

    Figure 14.A spanning set of cuts for two-loop Sudakov form factor.

    Figure 15.Via CK operation term in diagram (a) reducing to a tadpole may generate a term in diagram (b) reducing neither to tadpole nor massless bubble.

    On reduced tadpoles and massless bubbles—In the oneloop case,all terms with remaining free parameters can be reduced to massless bubbles or tadpoles,as shown in(3.8)in section 3.2.One can exclude these terms when we construct the ansatz for the master numerator.We may wonder whether there is a similar feature at two loops,and this is the question that is also raised in section 2.3.2.To check this,we collect all terms that will reduce to tadpoles or massless bubbles in the solution space of master numerators,and then we ask if they can be set to zero in the solution space.We find that for the first planar-ladder master topology,one can indeed make its numerator N1contain no tadpoles and massless bubbles in the solution space,but for the second master we are not able to do so.This implies that one can not drop all the terms corresponding to reduced tadpoles or massless bubbles from the starting ansatz.One way to understand this is that a tadpole or massless bubble term in one topology may become a term having non-trivial contribution in other topologies by CK operation.We show an example in figure 15.Although such terms have no physical effect,they are needed in the numerators because of the requirement of CK duality.

    3.3.2.Constraints from tadpoles and massless bubbles.In this subsection,we consider the possible further constraints on the solution space by asking if one can set the numerators of trivalent topologies containing massless bubbles or tadpoles to be zero.It is important to first emphasize that here the tadpoles or massless bubbles should not be confused with the reduced tadpoles or massless bubbles discussed in the previous subsection.In this subsection,we consider the trivalent topologies with a maximal number of propagators,and the tadpole or massless-bubble topologies all refer to the diagrams in figure 8.A complete set of CK relations will generate these topologies.As we will discuss below,we find that some of them can be set to zero while others can not.They are shown in figure 16 and figure 17 respectively.

    Constraints from tadpoles—We consider tadpole topologies first.We find that all the tadpoles on the external legs,i.e.(a)and(c)in figure 16,their numerators N14and N12can be set to zero.These tadpoles are related to massless bubbles via dual Jacobi relations.For tadpoles related to massive bubbles,however,we can not let all of them be zero but only N15for figure 16(b).Requiring the above three tadpole numerators to be zero,we can solve for 91 parameters,leaving 144 parameters.

    We provide an interpretation of why we may not set all tadpoles to be zero as follows.Through CK relations,the numerator of a tadpole topology is related to the numerator of topology containing a sub-bubble,such as shown in figure 18.This means that the vanishing of a tadpole numerator will require the corresponding bubble diagram that generates it should only have even powers of the loop momentum flowing inside the bubble.

    This requirement can be too strong for a massive bubble since it also has to satisfy constraints directly from unitarity cuts.On the other hand,a massless bubble on the external leg is not constrained by unitarity cuts directly(although through CK relations it can be related to other topologies and is constrained by unitarity cuts indirectly).Thus the constraints are weaker in the latter case.

    Constraints from massless bubbles—Next we consider topologies that contain a bubble on external legs.Although the massless bubble integral is scaleless and becomes zero after integration in dimensional regularization,there is still a subtlety at the integrand level.For such trivalent topologies,there is an intermediate on-shell propagator which is divergent and makes the integrand apparently ill-defined.

    Figure 16.Tadpoles or massless bubbles which can be set to zero.

    Figure 17.Tadpoles or massless bubbles which can not be set to zero.

    Figure 18.The generation of tadpoles from sub-bubble graphs using dual Jacobi relations.The red lines represent the propagators on which we apply CK operation.

    As pointed out in [35],one can impose constraints such that the numerator has the same power-counting property as that of Feynman diagrams.The numerator of a massless bubble obtained from Feynman rules satisfies the property that it contains two powers of the l or pi,where l is the loop momentum flowing inside the massless bubble and piis the corresponding external momentum.After integration,they will always be proportional to pi2and cancel the on-shell propagator in the denominator,such as

    which vanishes for ∈<0.But the numerator constructed by CK duality in general may contain numerators not satisfying this property.For example,terms like (ε1·ε2) (p1·p2)3would be ill-defined even after integration since the divergent propagator 1/pi2still exist.

    In our case,there are two bubble-on-external-leg topologies,figure 16(d) and figure 17(d).For the first case,the whole numerator N17can be set to zero,which will solve for 40 parameters,leaving 104 parameters.For the second topologies,one can not set the full numerator N13to be zero.We impose the condition that the terms which are not proportional tol22,(p1-l2)2and ε1·l2should be zero.Terms that proportional to ε1·l2are allowed because after integration,l2will be replaced by p1and these terms will vanish because of the transverse condition ε1·p1=0.We find these conditions can indeed be satisfied for this topology.This can further solve for 22 parameters,leaving 82 free parameters in the solution space.

    At this point,we further ask if we can let some other numerators be zero.We find that N9and N6can be set to zero in the above solution space.This will further solve for 8 and 9 parameters respectively,leaving 65 free parameters.To summarize,in the solution space with 235 parameters,one can set N6,N9,N12,N14,N15and N17to be zero,and no other numerators can be set to zero.

    3.4.Origins of the remaining parameters

    After applying unitarity constraints and other conditions as discussed in previous subsections,we find the CK-dual solution space still contains a large number of free parameters.All these parameters cancel in the final physical results.Below we discuss the origin of these parameters.

    The first origin of the parameters is related to the fact that the ansatz we construct includes terms that can be reduced to tadpoles and massless bubbles.Since the unitarity cut can not detect such terms,many of such parameters will remain undetermined.In the one-loop case,we have seen in section 3.2 that all the remaining parameters in (3.8) are related to this type of origin.However,in the two-loop case,as discussed in section 3.3.1,one can not exclude all terms that are reduced to tadpoles and massless bubbles,because of the requirement of CK duality.

    Figure 19.Different trivalent topologies may reduce to the same sub-topology.

    The second origin is more non-trivial and it comes from the freedom of redistributing terms in different topologies.In the CKdual ansatz,all the integrals correspond to trivalent topologies which have a maximal number of propagators.If a term in the numerator is proportional to one of the propagators,this term can be reduced to a sub-topology.It is possible that different trivalent topologies may reduce to the same sub-topology and such an example is shown in figure 19.In other words,the same contribution of a certain sub-topology may be expressed in terms of different maximal topologies,and this fact leads to some degrees of freedom when expressing the integrand.

    We have checked that most of the free parameters are associated with the second origin.Actually,only two parameters(out of 235) are purely due to the first trivial origin,which reduces only to tadpole or massless bubble sub-topologies.Below we briefly comment on the physical meaning of the free parameters.Since all parameters cancel in the final result,they should be understood as ‘gauge’ parameters in the sense that different choices of the parameters give different but physically equivalent representations of the same form factor.Such kind of freedom generally exists for loop integrands of form factors or amplitudes.The real non-trivial point here is that we also require the integrand solution to satisfy the CK duality.In particular,different solutions in the solution space preserve all dual Jacobi relations.The free parameters in this sense can be understood as CK-preserving deformation parameters.We mention that a similar solution space with a large number of parameters was also observed for the three-point half-BPS form factors inN=4 SYM up to four-loops[25–27].As discussed in[25,26],the CKpreserving deformation is also related to the generalized gauge transformation associated with the operator insertion.The existence of such a deformation is due to the color-singlet nature of the operator and thus is a special feature for form factors.We refer the interested reader to [25,26]for more details on this point.

    The two master numerators in the solution space of 235 free parameters are explicitly given in appendix B where the free parameters are set to zero for simplicity.The full numerators with all free parameters are provided in the auxiliary files.

    4.Integral reduction and integration

    In this section we will show that all the parameters disappear after integral reduction.We also discuss the integrated results and focus on the differences between the results of the ddimensional integrand and 4-dimensional integrand.Since the one-loop case is relatively simple,we will mainly focus on the two-loop case.

    To simplify the integrand,we first evaluate the color factors.This can be done by first expanding the structure constants in terms of color trace products using(2.1)and then applying the contraction rule:

    For the two-loop case,after contracting internal color indices,all the color factors will either be proportional toor equal to zero.We find that only the first 5 topologies in figure 7 have non-zero color factors,so all other topologies will not contribute to the full form factor in gauge theory and hence we will omit them in the following discussion.3The topologies with zero color factors can have a non-trivial contribution in the study of double copy for gravitational quantities.We will not consider this in the present work.We define the five non-trivial color factors asciand ciare given as

    For completeness,we also give their symmetry factors coming from the isomorphism of graphs:

    The full integrand of the two-loop form factor can be written as:

    4.1.Integral reduction

    In this subsection,we do the integral reduction and check that the remaining parameters all cancel in the final results.Notice that our integrand contains terms like εi·lj.A common strategy would be to first reduce such terms by the PV reduction[57]and then perform IBP reduction[58,59].Here we will use an alternative way based on gauge-invariant basis projection.

    Since the final result must be gauge-invariant,one can expand the final result by a set of gauge-invariant basis.For a general discussion about gauge-invariant basis,one may refer to section 2 of[60].For a two-point form factor,there is only one gauge-invariant basis,which is defined as:

    and this is obviously equivalent to the tree-level Sudakov form factor in (2.9).The integrand can be expanded as:

    where f(pj,lk) contains only Lorentz product momentum variables and we can perform IBP directly.To get the expression of f(pj,lk) we multiply B0to both sides of (4.6)and sum over the polarization vectors by the rule in (3.3),which gives

    Here we point out a technical subtlety:since the integrand is not manifestly gauge-invariant,the reference momentum ξ introduced through the helicity sum (3.3) will not vanish in f(pj,lk).We take ξ as another external momentum and perform the IBP with a set of propagator basis including ξ (here we use the LiteRed package [61]).We find that the ξ dependence indeed disappears after IBP,which also provides a cross-check for the result.

    After collecting the coefficients of master integrals,we find the parameters cancel and the coefficients we obtain are consistent with the known result[49].For the convenience of discussing the integrated results,we also use the following form

    where

    and the expansions of I(l=1,2)in terms of master integrals are given in appendix C.

    4.2.Integrated results

    The integrand numerators we obtain are polynomials of spacetime dimension parameter d=4-2∈.The term depending on ∈will not be fixable by four-dimensional cuts.For one-loop amplitudes,such terms will contribute to rational terms after integration.One may wonder what is the contribution of the ∈-dependent terms for the Sudakov form factor.In this subsection we address this problem by considering the integrated results.

    Considering first the one-loop result,we separate the result into two terms

    Next for the two-loop result,we separate the integrand into three terms

    Unlike the one-loop case,now we can see that terms proportional to ∈in the integrand will provide contribution starting fromdivergence.This shows clearly that ddimensional cuts are necessary to get the correct form factor result.

    5.Summary and discussion

    In this paper,we study the CK duality for one and two-loop Sudakov form factors of tr(F2) in the non-supersymmetric pure YM theory.We find that massless bubble and tadpole type topologies are necessarily involved.Although they are zero after integration,they are required to preserve the CK duality.The CK-dual numerators have the same powers of loop momenta as expected from Feynman diagrams.Interestingly,despite the non-trivial constraints from both unitarity cuts and dual Jacobi relations,the two-loop solution space still contains a large number of free parameters.Our results imply that the CK duality may also hold in the pure YM theory at three loops or for higher-point form factors,which are certainly interesting to be explored further.

    Since the CK duality plays an important role in constructing gravitational amplitudes via double copy,it would be interesting to explore the double-copy of the form factors.For the Sudakov form factor studied in this paper,one may naively apply double copy.However,to have a consistent gravitational quantity,it is crucial to ensure that the double-copy quantity not only preserves the diffeomorphism invariance but is also consistent with all unitarity cuts in the gravitational theory.Some progress has been made recently for the double copy of tree-level form factors in[62]where it was found that certain spurious poles become physical poles after the double copy.The new tree-level double-copy results also provide important building blocks for the further study of unitarity cuts at the loop level.It would be highly interesting to study further the double copy for loop form factors which we leave for future works.

    Acknowledgments

    It is a pleasure to thank Guanda Lin and Siyuan Zhang for their discussions.This work is supported in part by the National Natural Science Foundation of China (Grants No.12175291,11 935 013,11 822 508,12 047 503),and by the Key Research Program of the Chinese Academy of Sciences,Grant NO.XDPB15.We also thank the support of the HPC Cluster of ITP-CAS.

    Appendix A.Complete two-loop dual Jacobi relations

    Here are all CK relations we have between 17 cubic graphs in figure 7 and figure 8.(A.1)–(A.7) can generate other seven numerators detectable by cuts.(A.8)–(A.15) are relations to generate numerators of graphs with a massless bubble or tadpole.(A.16)–(A.23) are other CK relations.They will not provide further constraints for the 15 numerators obtained from (A.1) to (A.15).

    Appendix B.Two-loop master numerators

    In this appendix,we provide the two-loop master numerators after using the unitarity constraints in section 3.3.1,which contain 235 unfixed parameters.Here we set all free parameters to zero for simplicity,and the full numerators with parameters are provided in the auxiliary files.

    Figure 20.Graphs corresponding to master integrals.

    Appendix C.Master integral

    The definition of I(l)is:

    After IBP,we can expand I(l)by master integrals:

    C.1.One-loop

    Master integral for one-loop is

    The coefficient of this master is

    C.2.Two-loop

    Master integrals for two-loop are below.Notice that the previous three integrals have exact expressions in terms of Gamma functions.Integral I4can only be expressed in hypergeometric functions and we do not list it.For details,you can refer to section 2 of [49].

    The coefficients of these masters for the two-loop form factor are

    ORCID iDs

    亚洲国产毛片av蜜桃av| 两人在一起打扑克的视频| 老司机午夜福利在线观看视频| 日本三级黄在线观看| 亚洲精品国产一区二区精华液| 国产精品秋霞免费鲁丝片| 51午夜福利影视在线观看| 无限看片的www在线观看| 91大片在线观看| 精品欧美一区二区三区在线| 久久亚洲真实| 日韩三级视频一区二区三区| 国产一卡二卡三卡精品| 免费人成视频x8x8入口观看| 免费在线观看亚洲国产| 怎么达到女性高潮| 动漫黄色视频在线观看| 亚洲片人在线观看| 99国产综合亚洲精品| 黄片大片在线免费观看| 久久久久国产一级毛片高清牌| 精品久久久久久成人av| 成年版毛片免费区| 色在线成人网| 午夜福利免费观看在线| 在线看a的网站| 国产欧美日韩一区二区精品| 精品久久蜜臀av无| 亚洲色图av天堂| 国产精品综合久久久久久久免费 | 精品熟女少妇八av免费久了| 欧美黄色淫秽网站| 欧美大码av| 午夜日韩欧美国产| 美女国产高潮福利片在线看| 亚洲熟女毛片儿| 天天影视国产精品| 在线天堂中文资源库| 日韩欧美免费精品| 男人的好看免费观看在线视频 | 午夜福利,免费看| 婷婷六月久久综合丁香| 黄色视频,在线免费观看| 久久亚洲真实| 首页视频小说图片口味搜索| 最新美女视频免费是黄的| 日日摸夜夜添夜夜添小说| 中文字幕精品免费在线观看视频| 夜夜看夜夜爽夜夜摸 | 日本一区二区免费在线视频| 正在播放国产对白刺激| 好看av亚洲va欧美ⅴa在| 欧美精品啪啪一区二区三区| 国产精品综合久久久久久久免费 | 亚洲欧美一区二区三区久久| 欧美+亚洲+日韩+国产| 一进一出好大好爽视频| 99国产精品一区二区蜜桃av| 国产97色在线日韩免费| 成年版毛片免费区| 亚洲第一av免费看| 国产精品av久久久久免费| 后天国语完整版免费观看| 91麻豆av在线| 亚洲成人久久性| 亚洲成av片中文字幕在线观看| 免费女性裸体啪啪无遮挡网站| 午夜福利在线观看吧| 国产亚洲欧美98| 黄片播放在线免费| 不卡av一区二区三区| av电影中文网址| 国产一区二区在线av高清观看| 搡老岳熟女国产| 欧美av亚洲av综合av国产av| 免费看a级黄色片| 国产99白浆流出| 热re99久久精品国产66热6| 激情视频va一区二区三区| 国产亚洲精品第一综合不卡| 无遮挡黄片免费观看| 亚洲av成人av| 老汉色av国产亚洲站长工具| 人人妻人人爽人人添夜夜欢视频| 午夜免费观看网址| 亚洲激情在线av| 国产欧美日韩精品亚洲av| 日韩三级视频一区二区三区| 国产免费现黄频在线看| 成人黄色视频免费在线看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品国产高清国产av| 黑人巨大精品欧美一区二区蜜桃| 久久香蕉国产精品| 久久久久国内视频| 高潮久久久久久久久久久不卡| 亚洲午夜精品一区,二区,三区| 国产野战对白在线观看| 午夜福利在线免费观看网站| 熟女少妇亚洲综合色aaa.| 亚洲国产欧美网| 身体一侧抽搐| 久久精品aⅴ一区二区三区四区| 1024视频免费在线观看| svipshipincom国产片| 欧美日韩av久久| av天堂在线播放| 啦啦啦 在线观看视频| 午夜免费激情av| 丰满迷人的少妇在线观看| 午夜激情av网站| 午夜影院日韩av| 一夜夜www| 久久这里只有精品19| 看免费av毛片| 亚洲三区欧美一区| 啦啦啦 在线观看视频| 黄片大片在线免费观看| 正在播放国产对白刺激| 大型黄色视频在线免费观看| 久久精品人人爽人人爽视色| bbb黄色大片| 免费av中文字幕在线| 老熟妇乱子伦视频在线观看| 国产野战对白在线观看| 日韩高清综合在线| 国内久久婷婷六月综合欲色啪| 一本大道久久a久久精品| 午夜视频精品福利| 日韩欧美免费精品| 大码成人一级视频| 手机成人av网站| av中文乱码字幕在线| 国产精品九九99| 日韩大码丰满熟妇| 9色porny在线观看| 国产亚洲精品第一综合不卡| 国产成人精品在线电影| 欧美成人免费av一区二区三区| 水蜜桃什么品种好| 免费av毛片视频| 日本免费a在线| 欧美日韩亚洲国产一区二区在线观看| 热99国产精品久久久久久7| 亚洲少妇的诱惑av| 午夜亚洲福利在线播放| 99久久综合精品五月天人人| 一本大道久久a久久精品| av中文乱码字幕在线| 午夜福利免费观看在线| 十分钟在线观看高清视频www| 老熟妇乱子伦视频在线观看| 国产欧美日韩一区二区三| 日韩欧美免费精品| 久久久久久亚洲精品国产蜜桃av| 老熟妇乱子伦视频在线观看| 亚洲视频免费观看视频| 中亚洲国语对白在线视频| 人人妻人人添人人爽欧美一区卜| 精品电影一区二区在线| 一区福利在线观看| 啪啪无遮挡十八禁网站| 欧美一级毛片孕妇| 免费在线观看视频国产中文字幕亚洲| 男女下面进入的视频免费午夜 | 高潮久久久久久久久久久不卡| 国产精品亚洲一级av第二区| 欧美午夜高清在线| 老司机午夜十八禁免费视频| 欧美人与性动交α欧美软件| 欧美黑人精品巨大| 国产熟女午夜一区二区三区| 99国产精品99久久久久| 国产欧美日韩一区二区三| 国产黄色免费在线视频| 丰满人妻熟妇乱又伦精品不卡| 色综合站精品国产| 国产一区二区三区综合在线观看| 男人的好看免费观看在线视频 | 国产一区二区三区视频了| 精品乱码久久久久久99久播| 国产单亲对白刺激| 国产免费男女视频| 国产av一区二区精品久久| 69精品国产乱码久久久| 亚洲av片天天在线观看| av天堂在线播放| 亚洲精华国产精华精| 国内久久婷婷六月综合欲色啪| 久热爱精品视频在线9| 男人舔女人下体高潮全视频| 精品福利永久在线观看| 国产成人啪精品午夜网站| 9热在线视频观看99| 亚洲av成人av| 欧美黄色淫秽网站| 一进一出抽搐gif免费好疼 | 久久精品国产99精品国产亚洲性色 | 国产亚洲精品第一综合不卡| 亚洲精华国产精华精| 日本三级黄在线观看| 国产乱人伦免费视频| 日韩精品免费视频一区二区三区| 亚洲一区二区三区欧美精品| 咕卡用的链子| 亚洲精品久久成人aⅴ小说| av网站免费在线观看视频| 国产单亲对白刺激| 亚洲专区中文字幕在线| 男女之事视频高清在线观看| 窝窝影院91人妻| 亚洲 国产 在线| 国产精品av久久久久免费| 国产高清videossex| 精品国产超薄肉色丝袜足j| 日本黄色日本黄色录像| 免费在线观看黄色视频的| 久久国产亚洲av麻豆专区| 日韩欧美国产一区二区入口| 这个男人来自地球电影免费观看| 国产激情欧美一区二区| 搡老岳熟女国产| 日日干狠狠操夜夜爽| 校园春色视频在线观看| 亚洲中文字幕日韩| 淫秽高清视频在线观看| 我的亚洲天堂| 亚洲 国产 在线| 18美女黄网站色大片免费观看| 精品第一国产精品| 日韩欧美国产一区二区入口| 天堂中文最新版在线下载| 久久热在线av| 亚洲专区中文字幕在线| 国产有黄有色有爽视频| 午夜免费鲁丝| 纯流量卡能插随身wifi吗| 午夜福利在线免费观看网站| 亚洲男人的天堂狠狠| 精品一区二区三区视频在线观看免费 | 欧美日本亚洲视频在线播放| 视频在线观看一区二区三区| 淫妇啪啪啪对白视频| 成年版毛片免费区| 亚洲狠狠婷婷综合久久图片| 国产成人啪精品午夜网站| 99久久人妻综合| 亚洲欧美一区二区三区黑人| 91大片在线观看| 亚洲av成人不卡在线观看播放网| 免费搜索国产男女视频| 日日夜夜操网爽| 国产又色又爽无遮挡免费看| 国产av又大| 一本大道久久a久久精品| 他把我摸到了高潮在线观看| 美女 人体艺术 gogo| 麻豆一二三区av精品| 精品国产美女av久久久久小说| 久久性视频一级片| 欧美日韩国产mv在线观看视频| 亚洲成人免费av在线播放| 久久午夜亚洲精品久久| 午夜福利影视在线免费观看| 日韩欧美一区视频在线观看| 国产单亲对白刺激| 精品国产超薄肉色丝袜足j| 亚洲久久久国产精品| 亚洲久久久国产精品| 国产av又大| 9191精品国产免费久久| 亚洲av日韩精品久久久久久密| 久久久国产成人精品二区 | 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区| 乱人伦中国视频| 欧美激情高清一区二区三区| xxxhd国产人妻xxx| 国产成人系列免费观看| 国产精品 国内视频| 咕卡用的链子| 99热只有精品国产| 天堂影院成人在线观看| 成人国产一区最新在线观看| 国产色视频综合| 国产人伦9x9x在线观看| 久久国产精品男人的天堂亚洲| 大香蕉久久成人网| 香蕉久久夜色| videosex国产| 国产高清国产精品国产三级| 男女午夜视频在线观看| 视频在线观看一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 在线观看一区二区三区| 国产极品粉嫩免费观看在线| 99国产精品99久久久久| 亚洲精品粉嫩美女一区| 亚洲国产中文字幕在线视频| 最好的美女福利视频网| 国产精品偷伦视频观看了| 亚洲国产欧美网| 天堂√8在线中文| 中文字幕色久视频| 老汉色∧v一级毛片| 亚洲精华国产精华精| 中文字幕av电影在线播放| 国产精品久久久av美女十八| 欧美成人性av电影在线观看| 美国免费a级毛片| 一进一出好大好爽视频| 国产激情久久老熟女| 亚洲国产欧美日韩在线播放| 国产伦人伦偷精品视频| 国产99久久九九免费精品| 99国产极品粉嫩在线观看| 久99久视频精品免费| 色婷婷久久久亚洲欧美| 久久精品影院6| 成人18禁高潮啪啪吃奶动态图| 欧美一区二区精品小视频在线| 男人的好看免费观看在线视频 | 日本 av在线| 男人操女人黄网站| 老熟妇乱子伦视频在线观看| 久久人妻福利社区极品人妻图片| 国产精品秋霞免费鲁丝片| 午夜福利免费观看在线| 久久精品aⅴ一区二区三区四区| 亚洲av成人av| 亚洲aⅴ乱码一区二区在线播放 | 男女做爰动态图高潮gif福利片 | 欧美日韩一级在线毛片| 美女福利国产在线| 国产高清视频在线播放一区| 岛国视频午夜一区免费看| 99riav亚洲国产免费| 伦理电影免费视频| 国产精品99久久99久久久不卡| av视频免费观看在线观看| 最近最新中文字幕大全免费视频| 午夜老司机福利片| 国产亚洲精品综合一区在线观看 | 亚洲精华国产精华精| 怎么达到女性高潮| 99国产精品免费福利视频| 少妇的丰满在线观看| 国产视频一区二区在线看| 黑丝袜美女国产一区| 国产高清国产精品国产三级| 无遮挡黄片免费观看| 18禁黄网站禁片午夜丰满| 欧美 亚洲 国产 日韩一| 日韩视频一区二区在线观看| 精品国产乱子伦一区二区三区| 久久人人97超碰香蕉20202| 操出白浆在线播放| 黄色成人免费大全| 99精国产麻豆久久婷婷| 搡老熟女国产l中国老女人| 免费av毛片视频| 久久久精品国产亚洲av高清涩受| 日韩免费av在线播放| 免费人成视频x8x8入口观看| 欧美乱码精品一区二区三区| 亚洲九九香蕉| 日日摸夜夜添夜夜添小说| 久热这里只有精品99| 十分钟在线观看高清视频www| 国产99白浆流出| 国产精品成人在线| 两性午夜刺激爽爽歪歪视频在线观看 | 啦啦啦免费观看视频1| 久99久视频精品免费| 黄色成人免费大全| 51午夜福利影视在线观看| 免费日韩欧美在线观看| 国产一区二区在线av高清观看| 国产精品 欧美亚洲| av视频免费观看在线观看| 色婷婷av一区二区三区视频| 天堂俺去俺来也www色官网| 国产精品偷伦视频观看了| 午夜福利影视在线免费观看| 男女之事视频高清在线观看| 老司机福利观看| 一级a爱视频在线免费观看| 国产熟女xx| 看片在线看免费视频| 亚洲中文字幕日韩| 午夜福利欧美成人| 亚洲中文av在线| 99久久国产精品久久久| 欧美日韩国产mv在线观看视频| 精品国产一区二区三区四区第35| 国产精华一区二区三区| av国产精品久久久久影院| 9热在线视频观看99| 久久人人精品亚洲av| 一个人免费在线观看的高清视频| 精品熟女少妇八av免费久了| 色播在线永久视频| 在线视频色国产色| 最近最新中文字幕大全免费视频| 国产精品 欧美亚洲| 最好的美女福利视频网| 日本免费a在线| 免费av中文字幕在线| 女同久久另类99精品国产91| 午夜福利,免费看| 99香蕉大伊视频| 精品一区二区三区四区五区乱码| 黄片播放在线免费| 国产精品免费视频内射| 国产成人影院久久av| 日韩精品中文字幕看吧| 国产免费现黄频在线看| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品美女久久av网站| 欧美中文综合在线视频| 免费少妇av软件| 午夜精品久久久久久毛片777| 亚洲专区中文字幕在线| 69av精品久久久久久| av免费在线观看网站| 黄频高清免费视频| 亚洲 欧美一区二区三区| 亚洲精品粉嫩美女一区| 黄色a级毛片大全视频| 国产午夜精品久久久久久| 制服人妻中文乱码| 亚洲第一av免费看| 国产高清videossex| 一a级毛片在线观看| 免费在线观看亚洲国产| 亚洲第一av免费看| 久久精品国产99精品国产亚洲性色 | 免费av毛片视频| 国产精品久久电影中文字幕| 国产精品免费一区二区三区在线| 国产精品永久免费网站| 人人妻,人人澡人人爽秒播| 久久久久久久精品吃奶| 欧美乱码精品一区二区三区| 欧美在线一区亚洲| 亚洲国产精品合色在线| 亚洲人成77777在线视频| 国产精品乱码一区二三区的特点 | 欧美成人性av电影在线观看| 亚洲av成人一区二区三| 人成视频在线观看免费观看| 国产免费男女视频| 亚洲av日韩精品久久久久久密| 国产精品野战在线观看 | 精品卡一卡二卡四卡免费| 老汉色av国产亚洲站长工具| 人人澡人人妻人| 成人亚洲精品一区在线观看| 99国产精品99久久久久| 国产免费现黄频在线看| 91在线观看av| 日韩av在线大香蕉| 美女国产高潮福利片在线看| 国产亚洲精品久久久久5区| 精品福利观看| 男女下面插进去视频免费观看| 9191精品国产免费久久| 久久久久久久精品吃奶| 精品福利观看| 麻豆国产av国片精品| 亚洲成人免费电影在线观看| 麻豆国产av国片精品| 欧美日韩黄片免| 视频在线观看一区二区三区| 欧美一区二区精品小视频在线| 99在线视频只有这里精品首页| 久久精品91蜜桃| 久久这里只有精品19| 大型黄色视频在线免费观看| www国产在线视频色| 午夜影院日韩av| 久久这里只有精品19| 亚洲国产看品久久| 91麻豆av在线| 亚洲国产毛片av蜜桃av| 最新美女视频免费是黄的| 午夜免费成人在线视频| 精品一区二区三区四区五区乱码| 欧美中文综合在线视频| 欧美日韩黄片免| 丝袜人妻中文字幕| av免费在线观看网站| 人成视频在线观看免费观看| 亚洲专区字幕在线| 国产有黄有色有爽视频| 咕卡用的链子| 天天添夜夜摸| 黑人操中国人逼视频| 欧美日韩中文字幕国产精品一区二区三区 | ponron亚洲| 18禁裸乳无遮挡免费网站照片 | 欧美精品亚洲一区二区| 99精品欧美一区二区三区四区| √禁漫天堂资源中文www| 大型黄色视频在线免费观看| 嫁个100分男人电影在线观看| 好男人电影高清在线观看| 18禁黄网站禁片午夜丰满| 久久狼人影院| 高清av免费在线| 男女做爰动态图高潮gif福利片 | 亚洲视频免费观看视频| 久久 成人 亚洲| 亚洲中文av在线| 欧美久久黑人一区二区| 国产一区二区三区在线臀色熟女 | 国产一卡二卡三卡精品| 99久久99久久久精品蜜桃| 国产一区二区在线av高清观看| 女性被躁到高潮视频| 亚洲一码二码三码区别大吗| 欧美激情久久久久久爽电影 | 成人手机av| 不卡一级毛片| 女人被狂操c到高潮| 成人三级黄色视频| 久久久国产一区二区| 亚洲熟妇熟女久久| 日日爽夜夜爽网站| 国产精品久久久久久人妻精品电影| 女性被躁到高潮视频| 天堂中文最新版在线下载| 超碰成人久久| 一区福利在线观看| 黑人欧美特级aaaaaa片| 在线免费观看的www视频| 日日夜夜操网爽| 亚洲精品粉嫩美女一区| 国产乱人伦免费视频| 久久精品91无色码中文字幕| 人人妻,人人澡人人爽秒播| e午夜精品久久久久久久| 欧美激情 高清一区二区三区| 国产成人av激情在线播放| 国产一卡二卡三卡精品| 高潮久久久久久久久久久不卡| 国产午夜精品久久久久久| 丝袜美足系列| 久久中文看片网| 女人被狂操c到高潮| 法律面前人人平等表现在哪些方面| 在线观看舔阴道视频| 亚洲男人天堂网一区| 高清av免费在线| 久久国产乱子伦精品免费另类| 精品熟女少妇八av免费久了| 在线观看免费视频日本深夜| а√天堂www在线а√下载| 香蕉久久夜色| 一级,二级,三级黄色视频| 欧美国产精品va在线观看不卡| 亚洲,欧美精品.| 操美女的视频在线观看| 日本黄色视频三级网站网址| 中文字幕另类日韩欧美亚洲嫩草| 人人澡人人妻人| 日韩欧美一区视频在线观看| 亚洲人成伊人成综合网2020| 色在线成人网| 日本一区二区免费在线视频| 夜夜看夜夜爽夜夜摸 | 国产午夜精品久久久久久| 黑人猛操日本美女一级片| 亚洲,欧美精品.| 91精品三级在线观看| 操美女的视频在线观看| 欧美在线黄色| 欧美成人免费av一区二区三区| 亚洲伊人色综图| 欧美日韩视频精品一区| a级毛片黄视频| 国产精品爽爽va在线观看网站 | 母亲3免费完整高清在线观看| 亚洲精品成人av观看孕妇| 大香蕉久久成人网| 国产一区在线观看成人免费| 久久国产亚洲av麻豆专区| 中文字幕另类日韩欧美亚洲嫩草| 国产精品国产高清国产av| 黄色片一级片一级黄色片| 国产成人系列免费观看| 日日爽夜夜爽网站| 午夜亚洲福利在线播放| 色老头精品视频在线观看| 精品一区二区三区四区五区乱码| 亚洲自拍偷在线| 一级毛片精品| 免费观看人在逋| 日韩中文字幕欧美一区二区| 久99久视频精品免费| 成人av一区二区三区在线看| 国产av一区在线观看免费| 中文字幕另类日韩欧美亚洲嫩草| 国产精品1区2区在线观看.| 女人精品久久久久毛片| 两个人免费观看高清视频| 国产亚洲精品一区二区www| 国内毛片毛片毛片毛片毛片| 在线十欧美十亚洲十日本专区| 亚洲一码二码三码区别大吗| 午夜激情av网站| 亚洲一区高清亚洲精品| 国产91精品成人一区二区三区| 国产激情欧美一区二区| 在线观看午夜福利视频|