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

    PV-reduction of sunset topology with auxiliary vector

    2022-10-22 08:14:50BoFengandTingfeiLi
    Communications in Theoretical Physics 2022年9期

    Bo Feng and Tingfei Li

    1 Beijing Computational Science Research Center,Beijing 100084,China

    2 Zhejiang Institute of Modern Physics,Zhejiang University,Hangzhou,310027,China

    3 Center of Mathematical Science,Zhejiang University,Hangzhou,310027,China

    4 Peng Huanwu Center for Fundamental Theory,Hefei,Anhui,230026,China

    Abstract The Passarino—Veltman(PV)reduction method has proven to be very useful for the computation of general one-loop integrals.However,not much progress has been made when it is applied to higher loops.Recently,we have improved the PV-reduction method by introducing an auxiliary vector.In this paper,we apply our new method to the simplest two-loop integrals,i.e.,the sunset topology.We show how to use differential operators to establish algebraic recursion relations for reduction coefficients.Our algorithm can be easily applied to the reduction of integrals with arbitrary high-rank tensor structures.We demonstrate the efficiency of our algorithm by computing the reduction with the total tensor rank up to four.

    Keywords: PV-reduction,sunset topology,loop

    1.Introduction

    Calculation of Feynman integrals at the multi-loop level is of great importance for both perturbation quantum field theory and particle experiments.The general strategy is to expand arbitrary integrals to finite simpler integrals (called Master Integrals)and then do the integration of master integrals only.Thus,tensor reduction of Feynman integrals is one of the critical steps for various realistic calculations in the Standard Model.There is a large body of research on one-loop reduction [1—25].For reducing multi-loop integrals,various ideas and methods have been developed.There are some works on two-loop tensor reduction for some special integrals[26—31],while an algorithm for the reduction of general Feynman integrals was proposed in [32—37].Among these methods,the Integration-By-Parts (IBP) method [36,37]is widely used and has been implemented by some powerful programs such as FIRE,LITERED,KIRA,etc [38—45].Although its popularity,the number of IBP identities grows very fast with the increasing number of mass scalars and higher rank of tensor structure,it will become a bottleneck when dealing with complicated physical processes.Even at the one-loop level,reducing a general tensor massive pentagon with the IBP method is difficult.Thus it is always welcome to find a more efficient reduction method.

    PV-reduction [3]for one-loop integrals has played a significant role in general one-loop computations in history.However,when trying to generalize the method to higher loops,it becomes hard,and there are not many works in the literature.In [26,46],the tensor two-loop integrals with only two external legs are reduced to scalar integrals with the PVreduction method.However,an extra massless propagator appears.Integrals with new topologies have been introduced as master integrals to address the problem.This results against the intuitive picture: Reduction is achieved by removing one or more propagators from the original topology.Thus the reduction and master integrals should be constrained to the original topology and its sub-topologies.

    Recently,we have generalized the original PV-reduction method for one-loop integrals by introducing an auxiliary vector R in[47,48].Our improved PV-reduction method can be easily carried out and compute reduction coefficients analytically by employing algebraic recursion relations.With these analytical expressions,one can do many things.For example,by taking limits properly,we can solve the reduction with propagators having general higher powers [49].When external momenta are not general,i.e.,the Gram determinant becomes zero,master integrals in the basis will not be independent anymore.We can systematically study their degeneration [50].

    Encouraged by the results in [47,48],we want to see if such an improved PV-reduction method can be transplanted to higher loops.When moving to two loops,some complexities arise.The first is that we need to determine the master integrals.For one-loop,the master integrals are trivially known.For higher loops,the choice of master integrals becomes an art.With proper choice,one can reduce the computation greatly[51].The second is that when intuitively generalizing our method to two loops,we will meet the irreducible scalar products,which cannot be reduced.In fact,these two complexities are related to each other.To avoid unnecessary complexity,we will focus on the simplest nontrivial two-loop integrals,i.e.,the sunset topology in this paper.With this example,we will show how to generalize our new PV-reduction method to two loops.

    This paper is structured as follows.In section 2,we briefly review our previous work on one-loop tensor reduction and discuss the main idea for the reduction of sunset integrals.In section 3,we derive recursion relations for sunset integrals using differential operators and the proper choice of the master integrals in our algorithm.In section 4,we demonstrate our algorithm successfully to sunset integrals with the total rank from one to four.In section 5,we give some discussions and the plan for further studying.Technical details are collected in the appendix.

    2.Brief discussion of tensor reduction with auxiliary vector

    In this section,we will give a brief discussion of the general idea of using the auxiliary vector R to do tensor reduction.We first review the method presented in[47,48]and then move to the two-loop sunset integrals.For the tensor reduction of a general one-loop integral (to avoid confusion,we denote Pjas propagators for one-loop integrals while Djfor two-loop integrals)

    one can recover its tensor structure by multiplying each index with an auxiliary vectorRi,μi.Furthermore,we can combine these RitoR=xi Rito simplify the expression (2.1) to

    The good point of using (2.2) instead of (2.1) is that we have avoided the complicated Lorentz tensor structure in the reduction process,while they can easily be retained using the expansion of R and taking corresponding coefficients of xi's.Since we are using dimensional regularization,we like to keep the D as a general parameter.For D=d-2∈,the master basis is the scalar integrals with n≤d+1 propagators,i.e.,the reduction results can be written as

    where we denoteas the integrals got by removing propagatorsPi1,Pi2,…,Piafrom In.For simplicity,we denote

    By the explicit permutation symmetry in (2.1),we only need to calculatewhile other reduction coefficients can be got by proper replacements and momentum shifting.With the introduction of auxiliary vector R in (2.2),one can expand the reduction coefficients according to their tensor structure

    By acting two types of differential operatorson (2.2),we get the recursion relations for the expansion coeffciients

    Table 1.Number of expansion coefficients and independent equations for tensor bubble.

    is the n×n rescaled Gram matrix and T=diag(a1+1,a2+1,…,an+1)is a diagonal matrix.The vectors in the equations (2.6) and (2.7) are defined as

    and

    whereindicates to drop the index ai.With the known boundary conditions,one can obtain all reduction coefficients by applying the recursions (2.6) and (2.7) iteratively.

    One important point of the above reduction method is that we need to use bothDiandT operators to completely fix unknown coefficients in (2.5).One simple explanation is that in (2.5) coefficients there are (n+1) indices (ai,i=1,…,n and rank m),so naively(n+1)relations are needed.A more accurate explanation is a little tricky.Let us take the bubble,i.e.,n=1 as an example.As shown in table 1,for rank m,there areunknown coefficients while the number of independentD-type equations isThus,only using theD-type operators we can fix all coefficients with odd rank m=(2k+1) by these known terms with lower rank m≤2k.However,for m=2k+2,Nc-ND=1,thus just usingD-type relations is not enough and we need to adopt theT operator to provide one extra independent relation.In fact,T provides more than just one relation,but other relations are not independent to these coming fromD and can be taken as the consistent check of the reduction method.

    In this paper,we want to generalize our reduction method from one-loop to two-loop integrals.Again,with a simplified tensor structure,we define the general tensor integrals of sunset topology as5For simplicity,we denote the scalar integral Ia1,a 2, a3 ≡ .

    where the propagators are

    In this paper,we consider only the reduction with all ai=1.For ai≥2,one can use the same strategy presented in[49].For simplicity,we denote

    Similar to what we do for one-loop case,we can expand reduction coefficients according to their tensor structure6We will discuss the choice of master integrals J in section 3.3.

    where we have written C,J to emphasize they are vector,while other kinematic variables are

    In(2.13),the expansion coefficientsare actually vectors with components corresponding to different master integrals Ji.The summation over indices{i1,i2,j}≥0 admits to the restriction 0≤(rl-il-j)/2 ∈N,l=1,2.We make the indices free by setting all=0for inappropriate indices throughout the paper,so we will drop the prime in other summations.Similar to the strategy for one-loop case,we solve expansion coefficients from lower rank levels to higher rank levels iteratively,where the rank level is defined as the total rank (r1+r2).

    When we do the tensor reduction,we will reach the sub-topologies where one of the propagators in(2.11)is removed.For these cases,the integrals are reduced to the product of two one-loop tadpoles.Their tensor reduction has been solved in appendix B.To simplify notation,we denote

    and their reductions are written as

    In (2.13),there are five indices for the expansion coefficients,i.e.,i1,i2,j,r1,r2.With the experience from one-loop integrals,every index needs a recursion relation,so in total five kinds of differential operators are needed to give sufficient recursion relations.From the form of (2.2),we can construct threeT-type operators

    By applying these operators on the tensor integrals (2.2),we produce combinations?12,?22,?1·?2in the numerator,respectively.Using the following algebraic decomposition

    we can write the result as a sum of terms with one-lower or two-lower rank levels and terms for lower topologies,thus we can establish three recursion relations for these operators with expansion (2.13).

    Having discussed the operatorsTi,it is natural to think about following twoD-type operatorsKμ·andKμ·.However,when acting on (2.2),it gives the combinations ?1·K and ?2·K.These two factors do not have simple algebraic decompositions like(2.18).Thus it is not so easy to find corresponding recursion relations.One key input in our paper is that the recursion relation ofD-type operators can be established if we consider the reduction of the sub-one-loop integrals first.As we will show,such a recursion relation is highly nontrivial and more discussions will be given later.

    3.Recursion relations for tensor integrals of sunset topology

    In this section,we first derive the recursion relations for expansion coefficients of threeT-type operators.Then we establish another two recursion relations ofD-type operators by reducing the left/right sub-one-loop first.As we will see,in the process,one needs to use a highly nontrivial relation of the reduction coefficients of tensor bubbles.Combining the recursion relations fromT-type andD-type operators,we get five relations.When using them to solve expansion coefficients,we find that not all coefficients can be fixed.It indicates that some integrals with nontrivial numerators should be recognized as master integrals.

    After introducing two auxiliary vectors R1,R2,we can define differential operators(where we use the shorthand ?A=)

    Using the two expressions and external Lorentz vectors,we can get the following Lorentz invariant combinations,i.e.,threeT-type operatorsT00,T00′,T0′0′defined in (2.17) and sixD-type operators

    It is easy to use the definition (3.1) to find their action on reduction coefficients,for example,

    Employing them,we can write down recursion relations for expansion coefficients.

    3.1.Recursion relations from T-type operators

    Here,we derive the recursion relations for threeT-type operators.Let us start with the action ofT00.It is easy to check

    by using the algebraic relation (2.18).Plugging the expansion (2.13) into both sides,we find the LHS of the equation is

    and the RHS is

    where the definition ofhas been given in(2.16).Comparing both sides,we obtain the recursion relation for operator T00with r1≥2

    Similarly,one can derive the recursion relations for operatorT0′0′with r2≥2

    Finally,using the algebraic relation

    with

    we derive recursion relations of operatorT00′for r1≥1,r2≥1 as

    where we have used the shorthand

    Employing the expansion (2.13),we find the LHS of (3.11) is

    and the RHS is

    Comparing both sides of equation (3.11),we have

    By now,we have obtained threeT-type relations(3.7),(3.8)and(3.15).As discussed in the introduction,these relations are not sufficient to solve all expansion coefficients iteratively.In table 2,we list the number of expansion coefficients (each vector→is counted as one coefficient)and the independent equations given by threeT-type recursions.Comparing two tables,we find the number of reminding unknown terms is universal:Nα-NT=1.So for each rank level r ≡(r1+r2)there reminds(r+1) unknown expansion coefficients.We need to find more relations to determine them.

    3.2.Recursions relation from D-type operators

    Having used theT-type operators,from the experience of one-loop integrals,it is obvious that we need to consider theD-type operators.Naively actingD10andD10′defined in(3.2)globally to(2.2),we get K·?iin numerator,which are irreducible scalar products and cannot be reduced further algebraically.How to get out of the deadlock?Let us look back to the two-loop integrals.Instead of doing the loop integration together,one can think it as two times integration:first,we integrate the ?1,then we carry out the integration of ?2.Using this idea,we can first do the tensor reduction of sub-one-loop integrals.Let us see what we can get.Let us start with the sub-one-loop integrals.From appendix A,the r-rank tensor reduction of one-loop bubbles

    can be written as a summation of (r-1)-rank and (r-2)-rank bubbles as well as the contribution of tadpoles,whereP0≡?2-M02,P1≡ (?-K)2-M12.The striking point is that the Gram determinant K2only appears once in the denominator,i.e.,

    where we have defined

    For example

    where the pole of s11only appears in the overall factor.Now we insert(3.17)to the two loop integration by regarding propagator D1,D3as the propagators P0,P1,respectively,i.e.,we reduce ?1first

    where we have defined

    and

    thus in (3.19) we have

    The appearancein the denominator in(3.19)causes some trouble since it depends on ?2,while the original sunset integrals do not have such a denominator.In the paper[46],they keep this factor and regard the two-loop integrals within the denominator as a new master integral.However,according to the reduction idea,such a solution is somewhat surprising since we expect that the reduction is achieved by removing original propagators.Thus,we should reach the original topology and its sub-topologies at last.

    Since we want to keep the original reduction picture,we will try to get rid ofin (3.19).The idea is simple: one can multiply both sides with(D+r1-3)before integrating ?2and get

    Writing

    the LHS of the equation becomes (where we have suppressed the integral sign ∫1,2for simplicity)

    where in the last term we replace the new tensor structure 2K·?2by the action of differential operatorD10′on the basic form I [r1,r2+1](see(3.21)).One important point in(3.26)is that we have theD-type action on I [r1,r2+1],which is what we are looking for.

    Now we consider the RHS of (3.25).There are three terms:For the frist term,using the expression ofin(3.22)and the rewriting in(3.26),we can write it as

    Table 2.Number of expansion coefficients (left) and independentT-type equations (right).

    From the last equation in(3.28)one can see that the rank level of all terms will be less than(r1+r2+1),except the last term with rank level(r1+r2+1).For example,the termhas rank r1+r2.Similar analysis can be done for the other two terms in the RHS of(3.25),it is easy to see that they can be written as a combination of terms with rank level lower than (r1+r2+1).

    Collecting all terms together,the RHS of (3.25) is

    Rearranging (3.25),we arrive

    The recursion relation(3.30)is the wantedD-type recursion relation.The LHS of the equation can be written as a summation of expansion coefficients with rank (r1,r2+1),(r1-1,r2+2) of rank level (r1+r2+1) while the RHS can be written as a sum of expansion coefficients with rank level less than (r1+r2+1) and lower topologies.Thus,it is a recursion relation for rank level,instead of the explicit rank configuration (r1,r2).The reason is that under theD-type action,different rank configurations are mixed as at the LHS of (3.30).

    By our recursion assumption,the RHS is considered to be known,and we can just write it as

    The expression of (3.31) for explicit r1,r2can be found in appendix C.Putting the expansion (2.13) into the LHS and after

    Table 3.Number of expansion coefficients and independent equations for several rank levels.

    some algebra,finally we arrive at the algebraic recursion relation for unknown coefficients→

    Similarly,if we reduce ?2first,we have

    which is the dual expression of (3.30).Similarly,the lower-topology termin the equation can be got fromby regarding P0,P1as D2,D3,respectively,and multiplying withAgain,after plugging expansion(2.13)into the LHS of (3.33) and writing the RHS as

    we will get another algebraic recursion equation for expansion coefficients which is dual to (3.32).

    By now,we have derived twoD-type relations (3.30) and (3.33).Unlike theT-type relations which give a recursion relation for a particular rank configuration,D-type7will mix them: oneD-type relation (3.30) mixes ranks (r1,r2+1) and(r1-1,r2+2)while the otherD-type relation(3.33)mixes ranks(r1+1,r2)and(r1+2,r2-1).When using both types,one must do reduction level by level.As shown in table 3,one can solve all expansion coefficients for r1+r2>2 by combining T-type andD-type relations while there are still four expansion coefficients to be determined,which indicates we need more relations for r1+r2≤2 level.As will be discussed later,this observation gives us a hint of the choice of master integrals.

    3.3.Master integrals choice

    Up to now,we have avoided giving the explicit choice of master basis in (2.13),since all recursion relations ofT-type and D-type are independent of the choice.The only constraint is that they cannot contain auxiliary vectors.However,as pointed out in table 3,the reduction for some lower-rank tensor integrals is not completely clear.

    Let us consider them one by one.For rank level zero,it is the scalar sunset I1,1,1with no auxiliary vectors,so there are no differential operators that can act on it with nonzero results.Then we must regard it as a master integral,i.e.,J1≡I1,1,1.Next,we consider integrals with rank level one,for example,the integration

    According to its tensor structure,we should have

    The logic of the PV-reduction method is to solve B by multiplying K on both sides

    Now a key point appears.Unlike the one-loop case,the contraction ?1·K cannot be written as the combinations of Di.Thus,we cannot solve B as a function of external momentum K,masses,and space-time dimension D only.This observation is well known,and it is related to the concept of the irreducible scalar product(ISP).For a given L-loop integrals with E+1 external legs,there areindependent scalar products involving at least one loop momentum.If there are N propagators Di,N scalar products can be written as the linear combination of Di's,which leadsirreducible scalar products.For one-loop integrals,we have L=1 and N=E+1,so there is no ISP,and every tensor integral can be decomposed to the scalar basis.But for two-loop integrals,there are nontrivial ISPs.Thus there are multiple master integrals for a given topology.For our sunset topology,we have L=2,E=1,N=3 and NISP=2,which are given by (?1·K) and(?2·K).With the above analysis and the irreducibility of the LHS of(3.37),it is natural to take the LHS to be a master integral

    Similarly,we get another master integral

    In table 3,there are two undetermined expansion coefficients for rank level ne.After taking J2,J3above to be master integrals,they can be solved easily as shown in the next section.

    Then we consider the rank level two integrals,for example,

    According to its tensor structure,we have

    The only differential operator that can reduce the integral isT00′,which gives just one equation while there are two unknown expansion coefficients.As pointed out in table 3,there is only one undetermined expansion coefficient for rank level two.When combining the information of ISP,it is natural to take the integral obtained by acting withD10D10′,i.e.,

    to be another master integral.In fact,another two choices,i.e.,

    should be equivalently fine.

    Having chosen master integrals for sunset topology,we discuss the master integrals for lower topologies.These lower topologies are obtained by removing one of the sunset’s propagators,and the resulted topology is the product of two one-loop tadpole integrations.For one-loop integrals,we know the basis is the scalar integrals,thus it is easy to see that we have the following three master integrals

    One technical point is that the standard form of J5,J6should be written with the proper momentum shifting.For example,

    Our discussion on ISP and the choice of basis hints that the number N of master integrals is related to the number Nispof ISP when the masses and momenta are general.Another point is that when using the FIRE to do the reduction,it will take a different master basis.In appendix D we will present the non-degenerate transformation matrix between these two basis,thus our choice of basis is legitimate.

    Before ending this subsection,let us discuss briefly the number of master integrals when the kinematics and masses are not general,for example K2=0 or M1=M2.For K2=0,but M1,M2,M3are different,there are only four master integrals J1,J5,J6,J7and J2,J3,J4will be decomposed to a linear combination of them.We can get the proper reduction coefficients by taking the limit carefully,as discussed in[50].If M1=M2≠M(fèi)3,but K2≠0,we will have J5=J6and J2=J3,i.e.,there are only five master integrals.If the reduction coefficients are not singular under the limit M1→M2,we just add reduction coefficients together,i.e.,c2J2+c3J3→(c2+c3)J2.If the reduction coefficients are singular under the limit,we need to expand J3by other master integrals with expansion coefficients as the Taylor series of (M1-M2) as presented in [50].

    4.Examples

    Having laid out the general reduction frame in the previous sections,we present examples to demonstrate our algorithm.In this section,we will give the exact reduction process for the rank levels from one to four.The full expansion of coefficients is collected in appendix E and an attached Mathematica file.

    4.1.Rank level one

    Having chosen seven master integrals Ji,we start the reduction for tensor sunset with rank level one.There are two integrals

    Due to the total rank being less than 2,bothT-type andD-type recursions(3.30)and(3.33)cannot be used.However,we can apply operatorsD10andD10′on the two integrals,respectively.For example,the action ofD10on the integralgives

    where we have defined a constant vector

    The reason one can obtain the relation is that we choose J2as a master integral.

    For this case,the expansion (2.13) is

    Plugging it into the LHS of (4.2),we have the relation

    Comparing both sides,we can easily solve

    Similarly,we have

    4.2.Rank level two

    There are three integrals.By the expansion (2.13),we have

    Let us start from the rank (1,1).The onlyT-type recursion relation is given byT00′.Using (3.15),we have

    Then we consider the action ofD10D10′

    where={0,0,0,1,0,0,0}.Using the expansion (4.8),we find the LHS is

    Plugging it into (4.10),we have

    Combining (4.9) and (4.12) we can solve out

    Using the results of lower rank levels and lower topologies,we get

    From the above analysis,one can see the reason for taking master integral J4as discussed in the previous section.

    Next,we consider the rank (2,0).The onlyT-type relation is given byT00.Using (3.7) we have

    TheD-type relation (3.33) mix the rank (2,0) and the (1,1) and give

    Combining (4.16) and (4.15) and using (4.14),we can solve

    By proper replacement,we have

    4.3.Rank level three

    There are four rank configurations and the expansion (2.13) is

    We list 14 linear equations ofT-type andD-type relations for these ten expansion coefficients.There are only ten of them are linearly independent and the other four extra can be used as the self-consistence check.We listT-type equations first,

    ·T-type relations

    — rank (3,0)

    — rank (2,1)

    — rank (1,2)

    — rank (0,3)

    These six independent equations restrict the number of unknown expansion coefficients to 4.For every rank configuration,the number ofT equations is always one less than the number of expansion coefficients,which agrees with our discussion.

    Now we consider theD relations for rank level three,which are represented by four arcs in figure 1,where two red arcs for reducing ?1first and the other two orange arcs for reducing ?2first.We list them below:

    ·D-type relations by reducing ?1first

    — Mix rank (0,3) and (1,2)

    The reason for having two equations is that even with the removal of one R2byD operator,we are still left two different tensor structures,i.e.,s0′0′and.

    — Mix rank (1,2) and (2,1)

    ·D-type relations by reducing ?2first

    — Mix rank (3,0) and (2,1)

    — Mix rank (1,2) and (2,1)

    One can pick ten independent equations from the relations above and solve all expansion coefficients.Denoting

    we have linear equations=bwith matrix W

    and the known vector b which corresponds to lower rank level and lower topologies.We solve these linear equations and present the final results in appendix E.

    4.4.Rank level four

    There are five rank configurations and the expansion (2.13) is

    One can get 34 linear equations fromT-type andD-type relations for these 20 expansion coefficients.We can pick 20 independent equations among them.There are 18T-type equations

    For 16D-type equations,we have

    · reducing ?1first

    · reducing ?2first

    Picking 20 independent equations,we can solve all expansion coefficients of rank level four.Since its expression is long,we present them in an attached Mathematica file.All coefficients have been checked with the analytic output of FIRE[38,39,41].

    Before ending this section,let us summarize the algorithm to reduce tensor integrals of sunset topology to the master basis.The number of expansion coefficients for a particular rank level r is given by

    whereNc(m)=+1is the number of expansion coefficients for a one-loop rank-m tensor bubble.The recursion procedure can be nicely represented in figure 1.

    5.Conclusion and outlook

    In this paper,we have taken the first step in generalizing our improved PV-reduction method with auxiliary vector R to twoloop integrals,i.e.,we have carried out the tensor reduction for the simplest sunset topology.For two-loop integrals,theT-type relations can be established straightforwardly,while theD-type relations are rather nontrivial.We need to consider the reduction of its sub-one-loop integrals first.Combining the two types of relations,one can solve all expansion coefficients with the proper choice of master basis.

    One of our motivations for the paper is to find an alternative method for higher loop tensor reduction other than the IBP method.Comparing these two methods,one can find some similarities.All IBP relations are generated by following six primary relations·Kjwith i=1,2 and Kj=(?1,?2,K).Similarly,in our algorithm,there are also six operatorswith i=1,2 and=(R1,R2,K).The one-to-one mapping between these two types of manipulations seems to hint at a deep connection between these two methods although we are not clear at this moment.As pointed out in the paper,in our algorithm we have avoided the appearance of higher power propagators,thus the number of equations one needs to solve is much fewer than that using the standard IBP method.However,in[52]an idea has been suggested to avoid the higher power of propagators in the IBP method.It will be interesting to see if such an idea can be used in our algorithm.

    There are several interesting questions coming from this work.The first one is that,in section 3,we have discussed the choice of master integrals to be evaluating the ISP.We hope that for arbitrary two-loop or higher loop integrals with general masses and momenta K,one can determine the master basis from this perspective.

    The second question is the following.As emphasized several times,the nontrivialD-type relations depend crucially on the reduction relation (3.17) of one-loop integrals.We have tested the reduction of triangles and found similar reduction relation also holds for several examples.We find such a reduction relation holds for general one-loop integrals with the proof being presented in [50].

    The third problem is obvious.In this paper,we focus on the sunset topology.There are other topologies for two-loop integrals.One needs to check if the algorithm present in this paper can be generalized to these cases.With results from higher topologies,we can reduce sunset with propagators having higher power using the idea given in[49].We hope to give a positive answer to these interesting questions in our future work.

    Acknowledgments

    This work is supported by Chinese NSF funding under Grant Nos.11935013,11947301,and 12047502 (Peng Huanwu Center).

    Appendix A.The reduction of one-loop tadpoles and bubbles

    In this part,we will collect the reduction of one-loop tensor tadpoles and bubbles,which have been extensively used in this paper.

    For the one-loop rank r tensor tadpole

    the reduction result is given by

    where

    Figure 1.Algorithm for tensor reduction of sunset topology,where every point represents a particular tensor configuration (r1,r2).After using theT-type relations,only one unknown coefficient is left for each point.To determine the last unknown piece,we need to useD-type relations to relate different tensor configurations with the same rank level.We have used the red lines to represent relations coming from D-type relations for ?1,and orange lines,for ?2 respectively(here we drawD-type relations for r1+r2≤5).Each line provides a nontrivial relation among unknown coefficients.One can see that for rank level r≥3,the number of lines is equal or larger than the number of unknown coefficients,thus they will give enough equations.

    For the one-loop rank r tensor bubble

    the reduction result is given by

    where

    If we collect these coefficients into a listusing the notations

    we will have the following results:

    · For the rank r=1:

    · For the rank r=2:

    For our application,we have separated terms according to the power of pole s11.One crucial point is that we can writeC(22)as the combination of lower-rank coefficients and the contributions from lower topologies (tadpoles)

    · For the rank r=3:The three reduction coefficients are

    By matching the terms with different powers of pole s11,we can expandC(23)as

    One technical point of expansion (A.12) is that there are other different expansions,for example

    The difference between (A.13) and (A.12) is whether the coefficient of the bubble in the last term is zero or not.For simplicity of application,we prefer the form of(A.12).The choice of bubble coefficient to be zero for the remaining term has fixed the expansion uniquely.

    · For the rank r=4:After some algebra,we have

    · For the rank r=5:

    From these examples,one can observe the general expansion for any rank r >0 to be

    Let us start fromD-type recursion relations,i.e.,the differential equation generated byD1for the bubble.Acting with K?Ron both sides of (A.12) and using 2K·?=D0-D1+f1.The LHS is

    The termcan be reduced to scalar tadpole without any poles from Gram determent s11,and for simplicity,we just denote the contribution of tadpole topology as O(tad).The RHS becomes

    To require (A.18) to be equal to (A.17),one needs the following equation

    which is just (A.16) with r=r0-1.

    Then we consider theT-type recursion.Acting withT on the LHS of (A.16),we get

    while acting on the RHS we get

    One can prove (A.20) is equal to (A.21) by employing the recursion relation (A.16) for rank r0-2.

    One crucial point of the above proof is that,in principle,one can establish recursion relations offor general r.Thus we can solve it and get the complete relation (A.16) with explicit expression for.This result will give another recursive way to analytically compute reduction coefficients,which may be more efficient than the one given in [48].Furthermore,knowing,one may write down the analytic expression ofin(3.32)and its dual form.Then we can get the algebraic recursion relation for all coefficients for any rank level,just like what has been done in[47,48].It is obvious that with higher and higher rank levels,all available programs,like FIRE,Kira,Reduce [40—43]will become harder and harder to handle analytically.However,the method proposed in this paper can still work with fewer efforts.

    Appendix B.Reduction of lower topologies

    Here we give the reduction results for integrals got by removing the ith propagator in the tensor sunsets,i.e.,.We denote

    from which we can easily read outin(B.1).Forbecause the subtlety discussed in(3.45),we need to shift ?1in the numerator.The computation details are:

    Similarly,we have

    Appendix C.Lower rank/topology terms in D-type relation

    To write downD-type recursion relation for given r1,r2,we need calculatedefined in(3.31).Here we list some results.For the rank (1,r2) and (r1,1),the general forms are

    One can also write down expressions for(2,r2),(3,r2)etc.The expression will become longer and longer.Since in examples the rank level is up to four,we just list the terms we used in the examples

    Appendix D.Transfer matrix between tensor MIs and scalar MIs

    In section 3.3,we have discussed the choice of master basis.The basis chosen by FIRE6 is

    and the basis used in this paper is

    The transform matrix J=TF is

    One can check that with general masses and K,it is non-degenerate.

    Appendix E.Results

    Here we collect the reduction results from rank levels one to three for reference.It is too long to write in text form for rank level four,so we collect them in an attached Mathematica file.

    · rank r1+r2=1

    · rank r1+r2=2

    · rank r1+r2=3

    There are ten expansion coefficients for rank-3 level.

    精品福利观看| 久久久成人免费电影| 国产又黄又爽又无遮挡在线| 亚洲成人久久爱视频| 国产精品嫩草影院av在线观看 | 免费搜索国产男女视频| 欧美三级亚洲精品| 国产白丝娇喘喷水9色精品| 国产伦精品一区二区三区四那| 国产精品久久久久久久久免| 91精品国产九色| 人妻久久中文字幕网| 在线免费观看不下载黄p国产 | 国产探花在线观看一区二区| 国产69精品久久久久777片| 三级毛片av免费| 性插视频无遮挡在线免费观看| 日韩高清综合在线| 高清毛片免费观看视频网站| 啪啪无遮挡十八禁网站| 国产伦人伦偷精品视频| 亚洲黑人精品在线| 国产一区二区在线av高清观看| 黄片wwwwww| 神马国产精品三级电影在线观看| 国产白丝娇喘喷水9色精品| 国产伦一二天堂av在线观看| 免费在线观看日本一区| 麻豆一二三区av精品| 亚洲性久久影院| 欧美日韩国产亚洲二区| 亚洲最大成人中文| 97碰自拍视频| 少妇高潮的动态图| 少妇高潮的动态图| h日本视频在线播放| av.在线天堂| 日本免费一区二区三区高清不卡| 在线观看av片永久免费下载| 亚洲 国产 在线| 在线观看av片永久免费下载| 亚洲成人中文字幕在线播放| 国产91精品成人一区二区三区| 国产精品一区二区三区四区久久| 床上黄色一级片| 男女边吃奶边做爰视频| 久久久久久久午夜电影| 美女xxoo啪啪120秒动态图| 日韩欧美免费精品| 97热精品久久久久久| 97热精品久久久久久| 亚洲国产色片| 亚洲欧美激情综合另类| 韩国av在线不卡| 午夜激情福利司机影院| 蜜桃亚洲精品一区二区三区| 亚洲成a人片在线一区二区| 午夜激情欧美在线| 中文字幕精品亚洲无线码一区| 国产亚洲av嫩草精品影院| 国产精品久久视频播放| 亚洲四区av| 精品人妻熟女av久视频| 亚洲va在线va天堂va国产| 欧美黑人巨大hd| 色播亚洲综合网| 免费大片18禁| 18禁在线播放成人免费| 观看美女的网站| 少妇的逼好多水| 欧美日本亚洲视频在线播放| 久久热精品热| 国语自产精品视频在线第100页| 黄色丝袜av网址大全| 一本久久中文字幕| 免费看美女性在线毛片视频| 88av欧美| 国产精品国产三级国产av玫瑰| 精品国产三级普通话版| 久久精品国产自在天天线| 少妇人妻一区二区三区视频| 村上凉子中文字幕在线| 国产三级在线视频| av天堂在线播放| 日韩精品青青久久久久久| 国产精品无大码| 3wmmmm亚洲av在线观看| 日韩欧美一区二区三区在线观看| 国产一区二区三区av在线 | 嫩草影院精品99| 99在线人妻在线中文字幕| 亚洲av成人av| 国产精品人妻久久久影院| 国产国拍精品亚洲av在线观看| 久久久成人免费电影| 他把我摸到了高潮在线观看| 91麻豆精品激情在线观看国产| 香蕉av资源在线| 色吧在线观看| 女生性感内裤真人,穿戴方法视频| 有码 亚洲区| 琪琪午夜伦伦电影理论片6080| 九九热线精品视视频播放| 俄罗斯特黄特色一大片| 久久久久久久午夜电影| 九九热线精品视视频播放| 女生性感内裤真人,穿戴方法视频| 丰满乱子伦码专区| 在线观看美女被高潮喷水网站| 亚洲国产欧美人成| 亚洲av日韩精品久久久久久密| 亚洲av成人av| 我要看日韩黄色一级片| 嫩草影院入口| 久久久国产成人免费| 久久草成人影院| 高清在线国产一区| 可以在线观看毛片的网站| 国产免费男女视频| 色尼玛亚洲综合影院| 91麻豆av在线| 亚洲第一电影网av| 淫妇啪啪啪对白视频| 99热这里只有是精品50| 免费一级毛片在线播放高清视频| 亚洲性久久影院| 又紧又爽又黄一区二区| 亚洲专区中文字幕在线| 国产 一区精品| 少妇的逼好多水| www.色视频.com| 国产三级中文精品| 日韩欧美精品v在线| 热99在线观看视频| 国产私拍福利视频在线观看| 免费av毛片视频| 九色成人免费人妻av| 国产精品一区二区三区四区久久| av天堂中文字幕网| 1000部很黄的大片| 精品人妻1区二区| 国产一区二区在线观看日韩| 国产成年人精品一区二区| 久久精品综合一区二区三区| 少妇裸体淫交视频免费看高清| a在线观看视频网站| x7x7x7水蜜桃| 人妻久久中文字幕网| 国产精品一区二区免费欧美| 搡老熟女国产l中国老女人| 国产精品98久久久久久宅男小说| 精品人妻偷拍中文字幕| 一级a爱片免费观看的视频| 美女黄网站色视频| a级毛片免费高清观看在线播放| 啪啪无遮挡十八禁网站| 久久午夜亚洲精品久久| 人妻少妇偷人精品九色| 中文字幕精品亚洲无线码一区| 中出人妻视频一区二区| 日本 欧美在线| 中文字幕免费在线视频6| 舔av片在线| 国产精华一区二区三区| 国产精品98久久久久久宅男小说| 午夜日韩欧美国产| 欧美日本视频| av.在线天堂| avwww免费| 久久久久久国产a免费观看| 99九九线精品视频在线观看视频| 亚洲黑人精品在线| 亚洲人与动物交配视频| 日韩欧美国产在线观看| 国产探花极品一区二区| 亚洲精品亚洲一区二区| 九九久久精品国产亚洲av麻豆| 成人特级av手机在线观看| 最近最新免费中文字幕在线| 麻豆精品久久久久久蜜桃| 国内精品久久久久精免费| 亚洲最大成人中文| 亚洲中文字幕一区二区三区有码在线看| 日日啪夜夜撸| 他把我摸到了高潮在线观看| 变态另类丝袜制服| 亚洲精品影视一区二区三区av| 免费看日本二区| 真人做人爱边吃奶动态| 亚洲成人免费电影在线观看| 深夜a级毛片| 舔av片在线| 国产精品一区www在线观看 | 国产精品久久久久久亚洲av鲁大| 国产91精品成人一区二区三区| 日韩av在线大香蕉| 国产精品一区www在线观看 | 久久婷婷人人爽人人干人人爱| 国产一区二区三区在线臀色熟女| 久久久国产成人免费| 午夜久久久久精精品| 亚洲午夜理论影院| 国产老妇女一区| 一级a爱片免费观看的视频| 国产精品爽爽va在线观看网站| 久久午夜亚洲精品久久| 亚洲avbb在线观看| 亚洲中文日韩欧美视频| 亚洲人成网站在线播放欧美日韩| 亚洲精品成人久久久久久| 国产精品日韩av在线免费观看| 天天躁日日操中文字幕| 成年女人毛片免费观看观看9| 99久久精品一区二区三区| 又紧又爽又黄一区二区| 欧美日韩精品成人综合77777| a级毛片免费高清观看在线播放| 高清毛片免费观看视频网站| 婷婷丁香在线五月| 成人综合一区亚洲| 日韩,欧美,国产一区二区三区 | 久久久午夜欧美精品| a级毛片免费高清观看在线播放| 日本 欧美在线| 久久精品国产亚洲网站| 最近最新中文字幕大全电影3| 一个人看视频在线观看www免费| 一区二区三区高清视频在线| 香蕉av资源在线| 丝袜美腿在线中文| 国产麻豆成人av免费视频| 国产 一区 欧美 日韩| 欧美日韩瑟瑟在线播放| 日韩一本色道免费dvd| 五月玫瑰六月丁香| 欧美黑人巨大hd| 日韩强制内射视频| 国产精品嫩草影院av在线观看 | av专区在线播放| 欧美日韩黄片免| 日本五十路高清| 国产午夜精品久久久久久一区二区三区 | 色精品久久人妻99蜜桃| 麻豆精品久久久久久蜜桃| 麻豆一二三区av精品| 亚洲一区二区三区色噜噜| 男女做爰动态图高潮gif福利片| 久久久久久久精品吃奶| 夜夜夜夜夜久久久久| 免费av观看视频| 久久热精品热| 琪琪午夜伦伦电影理论片6080| 亚洲成人久久爱视频| 最近最新中文字幕大全电影3| 国产精品乱码一区二三区的特点| 精品人妻视频免费看| 高清日韩中文字幕在线| 亚洲国产精品成人综合色| 亚洲成人久久爱视频| 欧美一级a爱片免费观看看| 全区人妻精品视频| 久久99热这里只有精品18| 美女xxoo啪啪120秒动态图| 啦啦啦啦在线视频资源| 麻豆一二三区av精品| 精品一区二区三区视频在线观看免费| 国产蜜桃级精品一区二区三区| 丝袜美腿在线中文| 中出人妻视频一区二区| 国产探花极品一区二区| 村上凉子中文字幕在线| 成人特级av手机在线观看| 可以在线观看的亚洲视频| 国产探花在线观看一区二区| 日本黄色视频三级网站网址| 男女做爰动态图高潮gif福利片| 69av精品久久久久久| 国产v大片淫在线免费观看| av在线亚洲专区| 亚洲真实伦在线观看| 舔av片在线| 亚洲一区高清亚洲精品| 尾随美女入室| 日韩欧美国产一区二区入口| 男人舔奶头视频| 亚洲av美国av| 日韩欧美精品免费久久| 久久精品国产鲁丝片午夜精品 | 免费观看人在逋| 男人的好看免费观看在线视频| 久9热在线精品视频| 久久精品国产亚洲网站| 又爽又黄无遮挡网站| 免费黄网站久久成人精品| 级片在线观看| 国产精品电影一区二区三区| netflix在线观看网站| 亚洲,欧美,日韩| 看免费成人av毛片| 夜夜看夜夜爽夜夜摸| 99久久久亚洲精品蜜臀av| 久久欧美精品欧美久久欧美| 一级a爱片免费观看的视频| 成人国产麻豆网| 美女高潮喷水抽搐中文字幕| 黄色视频,在线免费观看| 性插视频无遮挡在线免费观看| 在现免费观看毛片| 国产日本99.免费观看| 91精品国产九色| 最近最新中文字幕大全电影3| 亚洲天堂国产精品一区在线| 久久国产精品人妻蜜桃| 少妇的逼水好多| av福利片在线观看| 中文字幕av成人在线电影| 日韩欧美精品v在线| 两个人视频免费观看高清| 欧美另类亚洲清纯唯美| 欧美激情久久久久久爽电影| 欧美bdsm另类| 久久久久国内视频| 国产色爽女视频免费观看| 国产精品一区二区性色av| 亚洲内射少妇av| 无人区码免费观看不卡| 免费一级毛片在线播放高清视频| 免费在线观看影片大全网站| 久久精品综合一区二区三区| 99久久精品国产国产毛片| 成人一区二区视频在线观看| 悠悠久久av| 欧美日韩中文字幕国产精品一区二区三区| 久9热在线精品视频| 窝窝影院91人妻| 2021天堂中文幕一二区在线观| 久久久久久久久久黄片| 亚洲国产色片| 免费观看精品视频网站| 麻豆久久精品国产亚洲av| 国产一区二区亚洲精品在线观看| 精品人妻熟女av久视频| 成人精品一区二区免费| 自拍偷自拍亚洲精品老妇| 一区二区三区四区激情视频 | 91午夜精品亚洲一区二区三区 | 亚洲在线自拍视频| 夜夜看夜夜爽夜夜摸| avwww免费| 日本一二三区视频观看| 亚洲在线自拍视频| 91久久精品电影网| 亚洲性夜色夜夜综合| 直男gayav资源| 亚洲成人久久性| 超碰av人人做人人爽久久| 亚洲人成网站在线播| 天天一区二区日本电影三级| 91麻豆精品激情在线观看国产| 欧美色视频一区免费| 国产亚洲精品综合一区在线观看| 日韩中字成人| 一区二区三区免费毛片| 免费人成在线观看视频色| 搡老岳熟女国产| 日韩欧美精品免费久久| 波多野结衣巨乳人妻| 97人妻精品一区二区三区麻豆| 色精品久久人妻99蜜桃| 亚洲欧美日韩高清在线视频| 久久热精品热| 国产黄色小视频在线观看| 国产综合懂色| 人妻夜夜爽99麻豆av| 欧美一区二区亚洲| 国产精品久久电影中文字幕| 精品一区二区三区视频在线观看免费| 国产一区二区在线av高清观看| 欧美丝袜亚洲另类 | 俄罗斯特黄特色一大片| 中文字幕免费在线视频6| 亚洲成人久久性| 狠狠狠狠99中文字幕| 国产高清三级在线| 婷婷精品国产亚洲av在线| 亚洲黑人精品在线| 免费在线观看日本一区| 国产三级在线视频| 中亚洲国语对白在线视频| 窝窝影院91人妻| 91av网一区二区| 99riav亚洲国产免费| 搡老熟女国产l中国老女人| 丰满的人妻完整版| 国产精品一区www在线观看 | 亚洲性夜色夜夜综合| 内地一区二区视频在线| 日韩欧美 国产精品| 在线观看美女被高潮喷水网站| 观看免费一级毛片| av天堂在线播放| 日日干狠狠操夜夜爽| 啦啦啦韩国在线观看视频| 亚洲真实伦在线观看| 亚洲av日韩精品久久久久久密| 国产视频内射| 久久香蕉精品热| 在线观看美女被高潮喷水网站| 日韩欧美国产在线观看| 午夜免费成人在线视频| 一个人观看的视频www高清免费观看| 级片在线观看| 人人妻,人人澡人人爽秒播| 91在线精品国自产拍蜜月| 国产成人a区在线观看| 亚洲精品在线观看二区| videossex国产| 午夜视频国产福利| www.色视频.com| 久久久久性生活片| 亚洲人成网站在线播放欧美日韩| 极品教师在线免费播放| 黄色丝袜av网址大全| 日韩欧美国产一区二区入口| 午夜福利在线观看吧| 亚洲午夜理论影院| 久久久久久久久久黄片| 免费看日本二区| 22中文网久久字幕| 亚洲七黄色美女视频| 中文字幕人妻熟人妻熟丝袜美| 不卡一级毛片| 免费观看人在逋| 99久久成人亚洲精品观看| 深夜a级毛片| 在线免费十八禁| 国语自产精品视频在线第100页| 99在线人妻在线中文字幕| 亚洲欧美精品综合久久99| 12—13女人毛片做爰片一| 中国美女看黄片| 老司机深夜福利视频在线观看| АⅤ资源中文在线天堂| 国产单亲对白刺激| 亚洲国产高清在线一区二区三| 国产高清有码在线观看视频| 国产亚洲精品av在线| 国产精品一及| 亚洲精品456在线播放app | 免费搜索国产男女视频| 亚洲成av人片在线播放无| 天堂√8在线中文| 日本在线视频免费播放| 人人妻人人澡欧美一区二区| 成年人黄色毛片网站| 免费无遮挡裸体视频| 乱人视频在线观看| 乱码一卡2卡4卡精品| 午夜福利欧美成人| 国产成人影院久久av| 国内精品久久久久精免费| www日本黄色视频网| x7x7x7水蜜桃| 免费观看在线日韩| 欧美zozozo另类| 亚洲午夜理论影院| 亚洲精品色激情综合| 亚洲精品乱码久久久v下载方式| 亚洲美女黄片视频| 午夜爱爱视频在线播放| 亚洲电影在线观看av| 中文字幕熟女人妻在线| 人妻丰满熟妇av一区二区三区| 欧美bdsm另类| 美女 人体艺术 gogo| 精品一区二区三区视频在线观看免费| 69av精品久久久久久| 赤兔流量卡办理| 我的女老师完整版在线观看| 91麻豆av在线| 久久久午夜欧美精品| 精品一区二区三区视频在线| av福利片在线观看| 欧美日韩乱码在线| 午夜视频国产福利| 国产高清视频在线观看网站| 亚洲图色成人| 国产一区二区三区av在线 | 国产精品自产拍在线观看55亚洲| 日本一本二区三区精品| 久久精品久久久久久噜噜老黄 | 极品教师在线视频| 欧美又色又爽又黄视频| 午夜精品久久久久久毛片777| 九九在线视频观看精品| 日本 欧美在线| 伊人久久精品亚洲午夜| 综合色av麻豆| 淫妇啪啪啪对白视频| 久久99热这里只有精品18| 婷婷丁香在线五月| 午夜福利18| 亚洲成人久久爱视频| 不卡一级毛片| 亚洲七黄色美女视频| 国产午夜精品论理片| 亚洲精品粉嫩美女一区| 淫妇啪啪啪对白视频| 草草在线视频免费看| 91麻豆av在线| 无遮挡黄片免费观看| 麻豆精品久久久久久蜜桃| 免费在线观看影片大全网站| 日韩欧美精品免费久久| 成年版毛片免费区| 精品免费久久久久久久清纯| 日韩在线高清观看一区二区三区 | 成人性生交大片免费视频hd| 国产极品精品免费视频能看的| 日韩av在线大香蕉| 天天一区二区日本电影三级| 偷拍熟女少妇极品色| 亚洲精华国产精华液的使用体验 | 两性午夜刺激爽爽歪歪视频在线观看| 91久久精品国产一区二区三区| 亚洲人成网站高清观看| 99久久精品热视频| 亚洲人成网站在线播| 五月伊人婷婷丁香| 国产日本99.免费观看| 男女啪啪激烈高潮av片| 最近在线观看免费完整版| 精品一区二区三区人妻视频| 最后的刺客免费高清国语| 色综合色国产| 日韩av在线大香蕉| 日韩精品中文字幕看吧| 一夜夜www| 国产午夜精品论理片| 我要看日韩黄色一级片| 成人欧美大片| 亚洲精品粉嫩美女一区| 尤物成人国产欧美一区二区三区| 免费av毛片视频| 中文字幕免费在线视频6| 美女大奶头视频| 成人二区视频| 日本与韩国留学比较| x7x7x7水蜜桃| 久久欧美精品欧美久久欧美| 高清毛片免费观看视频网站| 欧美高清成人免费视频www| 男女边吃奶边做爰视频| 日韩欧美国产在线观看| av天堂在线播放| 热99re8久久精品国产| 99国产极品粉嫩在线观看| 一级黄片播放器| 国产男靠女视频免费网站| 在线免费十八禁| 97碰自拍视频| 久久久久性生活片| 午夜免费成人在线视频| 欧美色欧美亚洲另类二区| 欧美又色又爽又黄视频| 两人在一起打扑克的视频| 村上凉子中文字幕在线| 欧美中文日本在线观看视频| 日本色播在线视频| 亚洲av不卡在线观看| 日韩精品有码人妻一区| 久久久久久久久久久丰满 | 在线观看美女被高潮喷水网站| 91狼人影院| 一夜夜www| 色综合色国产| 国产精品久久视频播放| 在线观看免费视频日本深夜| 国产亚洲精品久久久com| 精品一区二区三区视频在线观看免费| 亚洲最大成人手机在线| 国产精品综合久久久久久久免费| 能在线免费观看的黄片| 午夜福利在线在线| 欧美中文日本在线观看视频| 3wmmmm亚洲av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 免费电影在线观看免费观看| 成人欧美大片| 嫩草影院精品99| 亚洲欧美日韩东京热| 性插视频无遮挡在线免费观看| 日韩中文字幕欧美一区二区| 真人一进一出gif抽搐免费| 中文字幕av成人在线电影| 亚洲成a人片在线一区二区| 日本色播在线视频| 亚洲久久久久久中文字幕| 搡女人真爽免费视频火全软件 | 桃红色精品国产亚洲av| 日韩欧美 国产精品| 51国产日韩欧美| 亚洲美女黄片视频| 欧美成人免费av一区二区三区| 九九久久精品国产亚洲av麻豆| 人妻久久中文字幕网| 久久久久久久久大av| 精品久久久久久久久久免费视频| 十八禁网站免费在线| 少妇的逼水好多| 日本免费a在线| 国产高清不卡午夜福利| 一本一本综合久久| 12—13女人毛片做爰片一| 国产美女午夜福利|