Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: perturbative methods are indispensable in order to establish consistent theories of particle interactions , and to predict quantitatively their experimental manifestations . the anticipation of new phenomena in modern experiments and theoretical extensions of the standard model , requires cross - sections for complicated processes . this has driven a remarkable progress in the development of new computational methods . at the one - loop level , the calculation of cross - sections with five external legs is gradually becoming a routine activity ( e.g.  @xcite ) . at two - loops , there have been recent successful computations of amplitudes with four external legs and up to three kinematic scales ( e.g.  @xcite ) . at three - loops and beyond , amplitudes with up to one parametric variable have also been computed ( e.g.  @xcite ) . many new processes with higher final - state multiplicity , number of loops , and kinematic scales , have been identified to be important at the tev energy frontier . the aim of our paper is to provide a new method which can be used to compute loop amplitudes for such , more complicated , processes . loop integrations are cumbersome due to the presence of infrared singularities . loop integrals with many kinematic scales have , in addition , a complicated analytic structure with respect to their kinematic parameters . the tensor structure in gauge theories is also an issue , since it proliferates the number of terms . a general method for computing arbitrary loop integrals should extract their infrared and ultraviolet singularities , and treat simply kinematic discontinuities and threshold singularities . for practical applications , it should also be able to handle tensor integrals efficiently . there is no method which addresses satisfactorily all these issues ; known techniques can compute a limited number of amplitudes where some simplifications occur in special cases .    following a traditional approach to calculate loop amplitudes , one reduces the number of terms to a few master integrals and computes the latter analytically . one - loop integrals can be reduced using the classical method of passarino and veltman  @xcite . for generic multi - loop computations one can derive reduction identities from integration by parts  @xcite and the invariance of scalar integrals under lorentz transformations  @xcite . there is a large variety of approaches for evaluating the master integrals analytically . for example , one can compute integrals with a simple singularity structure and a small number of kinematic scales by integrating directly their feynman parameter representation . for more complicated cases one can use advanced techniques , such as the method of differential equations  @xcite . this approach can fail , however , if we apply it to complicated processes . the reduction algebra is hard , and the expressions of amplitudes in terms of master integrals may have spurious singularities which hamper their numerical evaluation . the extraction of @xmath1 poles in the master integrals , the evaluation of the coefficients of the @xmath1-expansion in terms of known analytic functions , and the analytic continuation of the latter to physically interesting kinematic regions are also involved . it is , thus , very well motivated to improve or replace the `` traditional scheme '' and to develop new automated methods . one - loop amplitudes can be entirely determined in four dimensions in terms of basic functions such as logarithms and polylogarithms that appear in the one - loop scalar box . new methods introduce sophisticated algorithms for the reduction to the basic functions ; by either numerical or analytical techniques , they control the appearance of spurious singular terms and minimize the size of the intermediate expressions  @xcite . recently , cross - sections for @xmath2 processes at nlo were computed  @xcite using such techniques . a method , which is inspired by techniques for real radiation at nlo , renders one - loop graphs numerically integrable  @xcite with the introduction of universal subtraction terms . a different approach uses unitarity , dualities , and analyticity properties for the determination of one - loop amplitudes  @xcite .    beyond one - loop , there has been a significant progress in the automation of the reduction methods to master integrals  @xcite . the infrared structure of multi - loop integrals is substantially more complex than at one - loop ; the development of methods for their numerical evaluation is more difficult . nevertheless , there is significant progress in this direction  @xcite . a powerful numerical method for multi - loop calculations is the method of sector decomposition  @xcite ; it simplifies recursively the singularities of feynman parameterizations and allows a straightforward expansion in @xmath1 . the method of sector decomposition has been very succesfully employed in the computation of several multi - leg integrals at one , two and three loops @xcite . this method has been introduced , recently , for the purely numerical evaluation of multi - loop amplitudes  @xcite . however , it is perplexing how to apply it for loop - integrals in non - euclidean regions .    in 1999 , smirnov  @xcite and , soon later , tausk  @xcite , introduced a new method for the evaluation of loop integrals . in their pioneering papers , smirnov and tausk @xcite computed analytically the first infrared divergent double box integrals . in  @xcite new two - loop integrals for @xmath3 massless processes were computed . in  @xcite the method was applied to double - box integrals with one additional mass - scale . the method was spectacularly applied in the computation of three - loop amplitudes  @xcite in @xmath4 supersymmetric yang - mills theory , shedding light to novel cross - order perturbative relations of the theory  @xcite . the smirnov - tausk method is based on a few simple ideas . starting from the feynman parameterization of a loop integral we can derive a new representation in terms of a multiple complex contour integral . such , mellin - barnes ( mb ) , parameterizations have enabled compicated loop calculations by using powerful methods for complex integration  @xcite . smirnov and tausk exploited a novel property of these representations . infrared divergences localize on simple poles inside the complex integration volume . we can isolate the divergent pieces of the integral at @xmath5 , by using the cauchy theorem . after subtracting the divergent residues , we can perform a taylor expansion in @xmath1 , and sum up the remaining infinite series of residues . finally , we can work to derive analytic expressions for the infinite sums in the coefficients of the expansion in terms of logarithms , generalized polylogarithms  @xcite , and more complicated functions . the smirnov - tausk method is very powerful ; however , it is laborious and intricate . the isolation of the divergent residues in multiple mellin - barnes integrals is convoluted . in addition , it is difficult to identify infinite sums in terms of polylogarithm functions with known analyticity properties . as a consequence , the analytic continuation in physical kinematic regions is also involved . due to these complications , the method has been applied to a few master integrals with a small number of kinematic scales . in this paper , we generalize the method to a broader spectrum of applications .    as a first task , we automate the procedure for the isolation of the divergent residues at @xmath0 . smirnov  @xcite and tausk  @xcite use different techniques for finding these residues . the approach of smirnov is very intuitive , but daedal . we have found that the technique described by tausk in ref . @xcite resembles closely to a programmable algorithm . we have used it as a guide and we have written computer programs which subtract the @xmath6 poles in arbitrary mellin - barnes integrals .    in our method , we avoid entirely the painstaking tasks of finding analytic expressions for infinite sums in terms of polylogarithms , and performing the analytic continuation in the arguments of polylogarithms . mellin - barnes representations are valid in kinematic regions where loop integrals may be complex - valued . we have found that , in a broad spectrum of applications , it is simple to calculate the representations numerically . the only analytic continuation that is ever required is that of logarithms with a single kinematic scale as an argument . an important goal of our method is to calculate loop amplitudes in realistic gauge theories . we have found that tensor integrals and , furthermore , diagrams which belong to the same topology can be calculated collectively . as we will show , the integrand of a representation for a scalar integral will be a product of gamma functions and powers of kinematic invariants , while that of a generic tensor will be the same integrand as in the scalar integral multiplied by a polynomial in the integration variables . the evaluation of polynomials is fast in a numerical program ; the computational cost for evaluating tensor integrals or loop diagrams is not significantly larger than evaluating scalar integrals . the only practical issue that we need to address , is the book - keeping of the various terms that contribute to the polynomial ; we present an efficient solution of this problem here .    in this paper , we apply our method to a number of examples . we , first , test our method in scalar and tensor integrals of the one - loop massless hexagon topology . the purpose of this computation is to introduce our method for tensor integrals and to demonstrate that we can tackle problems which are relevant in computations of physical amplitudes . we present here results for tensors through rank six , in both the euclidean and the physical region for @xmath7 processes . the numerical programs that we have constructed are suitable for the evaluation of the qcd amplitudes in four - jet production at hadron colliders . in the second set of examples , we compute scalar two and three - loop integrals which are known analytically in all kinematic regions : the massless planar  @xcite and cross  @xcite double - box , the massless double - box with one off - shell leg  @xcite , and the massless planar triple - box  @xcite . they serve to cross - check our algorithms and to demonstrate that we can easily reproduce state - of - the - art computations . our numerical results are in excellent agreement with the analytic expressions . we also present a number of new results that would require significant efforts for their computation with traditional approaches . we present , for the first time , double - box integrals with up to four kinematic parameters , and triple - box integrals with up to three kinematic parameters computed in all physical regions .    in section [ sec : mbrepresentations ] , we explain our technique for deriving mellin - barnes representations for loop integrals . in section [ sec : analytic ] , we describe our routines for performing an @xmath1-expansion of mellin - barnes representations for scalar integrals . we extend these results to loop integrals with tensor numerators in section [ sec : tensors ] and , in section [ sec : numerical ] , we present our methods for the numerical integration . section [ sec : results ] is devoted to present our results for several integrals at the one , two and three loop level . finally we present our conclusions . we start with a brief discussion on the derivation of mellin - barnes representations for loop integrals from their feynman parameterization . the construction of parameterizations is not unique , and various representations of the same integral may have quite different features . for example , they could have a different integral dimensionality , or they could be better or worse suited for numerical evaluation . a valid parameterization , however , can always be found . a pedagogical introduction to the topic of mellin - barnes representations is presented in ref . @xcite .    as a concrete example we consider the one - loop box integral with two adjacent external legs off - shell and massless propagators , @xmath8 with @xmath9 and @xmath10 , @xmath11 , @xmath12 , @xmath13 , @xmath14 . the powers of the propagators @xmath15 and the dimension @xmath16 are kept arbitrary and we will derive a mellin - barnes parameterization of the box - integral in this general case . in this section , we will seek values for @xmath17 where the representation is well defined and the integral is finite . in the next section , we will explain the technique for the analytic continuation to the values of the parameters , such as @xmath18 , where the integral develops divergences . our starting point is the feynman parameterization of the box integral : @xmath19 with @xmath20 . the main tool for obtaining the mellin - barnes representation of the above integral is the formula : @xmath21 where the following conditions are satisfied :    1 . @xmath22 are complex numbers with @xmath23 , 2 . the contour of integration is a straight line parallel to the imaginary axis separating the poles of @xmath24 and @xmath25 , i.e. @xmath26 . we can use the identity of eq . [ eq : mbrep ] to simplify the entangled denominator of eq . [ eq : box_2m_fp ] , with the cost of introducing new integrations in the complex plane . before we do so , we would like to point out an attractive feature of eq . [ eq : mbrep ] . the above two conditions guarantee that the integrand in the right hand side of eq .  [ eq : mbrep ] vanishes at infinity . we can then close the contour of integration at infinity and calculate the integral using the cauchy theorem . if we close the contour to the side of the positive real - axis we obtain @xmath27 this is the taylor expansion of the left hand side of eq . [ eq : mbrep ] in the region @xmath28 . equivalently , by closing the contour of integration to the side of the negative real axis , we obtain a taylor expansion in the complementary region @xmath29 . this is a particularly useful property : representations of feynman integrals which are derived by using the mellin - barnes decomposition are valid in all kinematic regions . in addition , @xmath22 can be complex ; therefore , we can account for the infinitesimal imaginary part @xmath30 that is assigned to invariant masses and kinematic parameters . 3 may be recursively applied to denominators with more than two terms , yielding : @xmath31 we are now ready to apply the decomposition of eq . [ eq : mbrep_multi ] to the denominator of eq . [ eq : box_2m_fp ] , introducing three contour integrations : @xmath32 where the contours must satisfy condition 2 above , i.e. they must be chosen in such a way that the arguments of the gamma functions are all positive . it is now straightforward to integrate out the feynman parameters , using @xmath33 to exchange the order of integrations , we must assume that the exponents of the feynman parameters satisfy @xmath34 with these conditions the gamma functions on the right hand side of eq.[eq : int_feynman ] have arguments with positive real part as well . the final result is : @xmath35 where @xmath36 and we have adopted the shorthand notation @xmath37 in general , mellin - barnes representations for a loop - integral are well defined for number of dimensions and powers of propagators which render the real parts of the arguments in all gamma functions positive . for example , it is impossible to satisfy positivity in all gamma function arguments for @xmath38 . this is in accordance with the expectation that the scalar box integral is infrared divergent in four dimensions . but , for example , the representation of eq . [ eq : box_2m_mb ] is well defined if we choose the contour @xmath39 and the parameters @xmath40 . in the next section we will detail the analytic continuation from values of parameters for which the mellin - barnes representation is well - defined to values for which the integral develops divergences . that is , a procedure which makes the divergences appear explicitly in the form of poles in the parameter space . from eq . [ eq : box_2m_mb ] we observe that it is simple to implement the analytic continuation of kinematic parameters in a numerical evaluation , since they always appear in a factorized form , e.g. @xmath41 we need the analytic continuation of simple logarithms ; these are evaluated trivially in all kinematic regions : @xmath42 it is clear from the above example that we can derive a mellin - barnes representation for any integral starting from its feynman parameterization . the form of the parameterization depends on the choice of feynman parameters , the order in which the @xmath43 terms appear in the right hand side of eq . [ eq : mbrep_multi ] , and the implementation of the delta function constraint for the feynman parameters . this arbitrariness is more pronounced beyond one - loop .    in section [ sec : tensors ] we will derive representations for tensor one - loop integrals which simplify their numerical evaluation . we aim to find similarly simple representations for multi - loop tensor integrals as well . for this purpose , it is convenient to use the one - loop representations as building blocks , by employing a re - insertion technique  @xcite . we explain this technique in our second example . we now derive a mellin - barnes representation for the massless planar double - box integral with two adjacent legs off - shell . the integral is defined : @xmath44 with @xmath45 and @xmath10 , @xmath11 , @xmath12 , @xmath13 , @xmath14 . for simplicity , we have set the powers of propagators to one . we perform , first , the @xmath46-loop integral and derive the representation for this integral reading it off from the one - loop result of eq . [ eq : box_2m_mb ] , @xmath47 the off - shell legs of the @xmath46-loop box correspond to the @xmath48 propagators while the @xmath49 invariant mass gives rise to a new propagator for the second integration : @xmath50 the propagators @xmath51 form a new one - loop box with two adjacent legs off - shell , and we can read off , once again , the representation for this second integration using eq . [ eq : box_2m_mb ] . the final result for the two - loop box with two adjacent legs off - shell is : @xmath52 we have produced a representation for a two - loop integral by embedding representations for one - loop integrals into other representations . it is obvious that we can use the re - insertion method for writing representations for arbitrary multi - loop integrals . representations for integrals with given kinematic scales can be used to derive , easily , results for simpler integrals where some of the scales are taken to zero . let us , for illustration , derive the representation for the double - box with one - leg off - shell and the on - shell double - box . we must first take the limit @xmath53 in eq . [ eq : dbox_2m_mb ] . the term @xmath54 is vanishing in this limit unless @xmath55 at the same time . if we use the cauchy theorem , we find that the wanted limit is given by taking the residue of the integrand at @xmath56 . therefore , the double - box with one - leg off - shell is : @xmath57 similarly , the on - shell double - box is given by : @xmath58 note that , in order to have a finite limit for a vanishing kinematic scale , the mellin - barnes representation should have terms in the integrand which behave as @xmath59 . these emerge naturally in non - trivial representations due to eq . [ eq : mbrep ] . in the previous section we have noted that the mellin - barnes representation of a loop integral is valid if appropriate poles of the gamma functions lay separated on the right and left of the integration contours . this condition guarantees the equivalence between the mellin - barnes integral and the original denominator in the loop integral . it often occurs that this condition can not be satisfied for values of the dimension parameter and the powers of the propagators which are relevant for a realistic application . let us consider , as an example , the one - loop box representation of eq . [ eq : box_2m_mb ] in the case of powers of propagators set to unity and the dimension in @xmath16 , @xmath60 this representation can only be valid for values of @xmath1 different from zero . for example , if we choose the contour @xmath61 we find that @xmath1 should be in the interval @xmath62 . we can use the cauchy theorem to obtain a representation in terms of a sum of contour integrals , which is valid at @xmath63 and the @xmath6 poles appear explicitly . the key point is that if the value of @xmath1 is chosen outside the allowed region , for example @xmath63 , some of the poles of the gamma functions will be on the wrong side of the contours ( see fig .  [ fig : pole_pos ] ) . to recover the original result , we must compare the position of the poles relative to the contours for values of @xmath1 inside the allowed region and the point outside . depending whether the poles crossed the contours from left or right , we should correct the representation by adding or subtracting the residue of the integrand on those poles . this simple observation can be easily cast into general purpose algorithms which extract the poles in @xmath1 of arbitrary mellin - barnes representations .    , and some of the poles cross to the other side of the contour . to recover a valid representation for the loop - integral , these poles sould be isolated using the cauchy theorem . , height=151 ]    let us choose a value of @xmath1 inside the allowed region and observe the position of the poles with respect to the contours for some @xmath1-dependent gamma functions in our example . the poles of @xmath64 and @xmath65 , at @xmath66 and @xmath67 , are situated to the right of the @xmath68 and @xmath69 contours , respectively , for all non - negative integers @xmath70 . let us now take @xmath0 . we observe that the first poles , for @xmath71 , moved to new positions , @xmath72 and @xmath73 , which are on the left side of the contours . the poles for @xmath74 remain on the right side . therefore , the representation is not valid at @xmath75 , since the contours separate poles which originate from the same gamma function . to obtain a valid expression for the integral at @xmath76 , we need to use the cauchy theorem to isolate the crossing poles ( @xmath71 ) . smirnov , in ref . @xcite , provides a number of pedagogical examples where this is done . the approach of smirnov is to deform the multiple contours , using the cauchy theorem , so that the `` offending '' poles never cross the contour of integration . tausk , in ref . @xcite , describes in detail this task as well , using as an explicit example the mellin - barnes representation for the two - loop cross box integral . the technique of tausk is different from the one of smirnov ; the main idea is to account for the poles that end up on the `` wrong '' side of the contour in a progressive way , as they cross the contours when changing continuously the value of @xmath1 . in this paper , we have followed the approach of tausk . we would like to recommend that the interested reader studies the pedagogical example of ref . @xcite . here we discuss our implementation of the algorithm which is depicted in the flow diagram of figure  [ fig : flow1 ] . consider the mellin - barnes representation @xmath77 of a loop integral . the integrand of @xmath77 depends initially on gamma functions . the application of the analytic continuation algorithm requires the iterative evaluation of residues . these may also give rise to psi functions , which have poles in the same positions the gamma functions do , and are treated identically for the purposes of analytic continuation . the contours @xmath78 of the representation are straight lines parallel to the imaginary axis , and are chosen so that the real parts of the arguments of all gamma functions are positive . the same condition determines an allowed region for @xmath1 , which we will assume that includes the point @xmath79 . for the sake of clarity , we will consider the case @xmath80 , but the method proceeds analogously when @xmath81 . the first step is to find the maximum value @xmath82 , with @xmath83 , for which the representation holds . at this value , one or some of the positions of the poles run into a contour . if the set of these poles , @xmath84 , is empty , i.e. no poles cross the contours when moving @xmath1 from @xmath85 to 0 , then the integral can be safely expanded around @xmath63 . so , let us focus on the more interesting case in which @xmath84 is non empty . here we have to distinguish between @xmath86 and @xmath87 .    if @xmath86 , at least one pole lies on the contours for the terminating value of @xmath1 . then , the representation has an unregulated singularity , and it is not feasible to perform an @xmath1 expansion and evaluate the contour integrals numerically . to remedy the situation we must change the contour of integration ; we use the cauchy theorem and displace one or more of the contours parallel to the imaginary axis . then , we attempt afresh the analytic continuation ; the poles should now hit the contour at @xmath87 . if this does not happen , due to an unfortunate choice of the contour we displaced , we simply iterate the procedure until we succeed .    when @xmath88 , we have to determine the residues in the poles that cross the contour . if there is only one pole in @xmath84 the situation is simple . let us look back at the specific example of the one - loop box . by applying the previous steps of the algorithm , we find that the first pole to cross a contour is at @xmath89 . we then choose one of the multiple integrals , e.g. the @xmath69 integral , and subtract the residue of the integrand on this pole , at @xmath90 , to the representation : @xmath91 the residue term in the above expression is a new multiple mellin - barnes integral with one integration variable less than in @xmath77 . in our example , the pole crosses the @xmath69 contour from the right side and we subtract its residue from the original representation . if a pole crosses the contour from the left , we should add its residue . at this point , we have achieved the analytic continuation in @xmath1 from @xmath85 to @xmath92 . we now repeat the algorithm setting as initial value @xmath93 , and perform the analytic continuation to the next value of @xmath1 where a new crossing occurs . we proceed iteratively , until we find all poles which cross the contours of integration in the terms of eq . [ eq : ancont ] and its descendants . the algorithm is more complicated if many poles cross simultaneously the contour at @xmath82 . in our one - loop box example , if we had chosen a symmetric contour @xmath94 then we would have two poles crossing at the same time : @xmath95 . one way to deal with this situation is to displace the contours in such a way that the degeneracy is removed , i.e. introduce asymmetric contours from the beginning . another alternative is to take combined residues in all the poles . to correctly account for all the terms in this case , we can imagine that the contours are only slightly displaced from their current position . then the poles will hit them one at the time and we can take residues successively . finally , we can restore the original position of the contours without affecting the position of the poles . there is , however , a caveat in taking multiple residues . let us consider the case of the one - loop box , where the poles @xmath96 and @xmath97 hit the contour at the same time . we take residues in sequence , considering first @xmath98 . we make the `` unfortunate '' choice to take the residue on the @xmath99 integration variable , and not on @xmath68 . the analytic continuation yields : @xmath100 then we consider the crossing of the @xmath101 pole on each of the two terms on the right hand side of the above expression . we find no problems in taking the residue of the first term . however , the position of the @xmath101 pole for the second term is now written : @xmath102 the pole @xmath101 is now independent of @xmath1 and gets nailed to the contour . this situation is similar to the one we encountered for @xmath86 , and we can deal with it in exactly the same way by displacing the contours . we apply the algorithm recursively , to all the terms which are generated by adding the residues of crossing poles to the original parameterization . at the end of this procedure , we obtain a sum of contour integrals plus some terms that do not contain any remaining integral due to taking residues on all integration variables . the method assures that the sum of all these contributions equals the original loop integral when @xmath103 . so we can expand in a power series of @xmath1 .    in this section , we have presented the algorithm assuming that only the dimension parameter @xmath1 is required to regulate the mellin - barnes parameterization . however , the method proceeds identically for integrals where multiple analytic regulators are required . this could be the case in non - covariant gauges or in integrals with irreducible numerators  @xcite or additional positive powers of propagators  @xcite . the actual implementation of the algorithm in a set of routines in both maple and mathematica is very efficient , and allows to perform the analytic continuation in the regulator parameters in a few seconds , for all the mellin - barnes representations that we have studied in this paper . mellin - barnes representations , with the exception of the massless two - loop diagonal box topology  @xcite , have been traditionally employed for the evaluation of master integrals . in this section , we introduce an efficient generalization of the method to tensor integrals ; our method does not require any reduction techniques to master integrals . common reduction methods produce extremely large algebraic expressions which may also have spurious singular denominators . here we derive efficient representations for one - loop tensor integrals , which are amenable to numerical integration . we can derive analogous representations for multi - loop integrals by using the re - insertion method of section  [ sec : mbrepresentations ] . we consider a generic one - loop tensor integral of rank @xmath104 with @xmath70 external legs : @xmath105 \left [ ( k+q_2)^2-m_2 ^ 2 \right ] \cdots \left [ ( k+q_n)^2-m_n^2 \right ] } \,\ ] ] where , @xmath106 and @xmath107 are external incoming momenta . we first introduce feynman parameters for the denominator . we obtain , @xmath108^n}\,,\ ] ] where @xmath109 and @xmath110\ ] ] now we perform the usual shift in the loop momentum , @xmath111 , and we obtain , @xmath112 } } { ( k^2-\delta)^n}\ , , \ ] ] where , we denote @xmath113 } = \sum_{\{j_1,\ldots , j_m\ } } k^{\mu_{j_1}}\cdots k^{\mu_{j_r } }   p^{\mu_{j_{r+1}}}\cdots p^{\mu_{j_m } } ,   \\   & & \hspace{2 cm } \{j_1,\ldots , j_m\ }   \in \mbox{permutations}(1,\ldots , m ) \nonumber .   \end{aligned}\ ] ] we find the standard feynman representation for the generic tensor one - loop @xmath70-point function , by integrating the loop momentum : @xmath114 } \delta^{\frac{r}{2 } } , \ , \end{aligned}\ ] ] where @xmath115 for @xmath116 odd , and @xmath117}$ ] for @xmath116 even . we observe that the sum , in the second line of the equation , is a polynomial in the feynman parameters , with tensor coefficients : @xmath118 } \delta^{\frac{r}{2 } } = \sum_{\{\nu_i > 0 \ } } c(\{\nu_i\ } ; \{p_i^{\mu_j}\ } ) x_1^{\nu_1 - 1}\cdots x_n^{\nu_n-1}\ ] ] for the scalar integral with unit powers of propagators , this sum is @xmath119 . we will use eq . [ eq : tensor - feynf ] to derive mellin - barnes representations for the tensor loop - integral . we would like to exploit the fact that the feynman representation of the scalar and tensor integrals have the same denominator : @xmath120 in section  [ sec : mbrepresentations ] , we have seen that we can obtain a mellin - barnes representation for a loop integral by decomposing the denominator of the feynman parameterization with eq . [ eq : mbrep_multi ] . we should anticipate that the mellin - barnes representations of tensors and scalars are , therefore , very similar . we can verify this intuition , if we consider one generic term in the polynomial of eq . [ eq : fp_polynomial ] , and follow the steps of section  [ sec : mbrepresentations ] . decomposing the denominator with eq . [ eq : mbrep_multi ] , we obtain @xmath121 where @xmath122 is the number of kinematic scales @xmath123 in the integral . the factors @xmath124 are universal for all terms in eq . [ eq : fp_polynomial ] , and originate from the dependence of @xmath125 on the feynman parameters . we now proceed to integrate out the feynman parameters using eq . [ eq : int_feynman ] . this integration yields gamma functions @xmath126 rewriting all gamma functions as above in eq . [ eq : fp_polynomial ] , we obtain , for the tensor integral , the gamma functions of the scalar integral multiplied by simple factors with mellin - barnes integration variables . the result for the one - loop tensor integrals depends on the topology and the rank of the tensors in a rather complicated manner . however , we find that the representation is of the following form @xmath127 the first line of the above representation contains all the terms which already appear in the analogous representation of the scalar integral with unit powers of propagators . the function @xmath128 is a _ polynomial in the mellin - barnes variables _ , with tensor coefficients in terms of the external momenta . this is very important for our strategy to evaluate tensor integrals and feynman diagrams . the polynomial in the numerator is smooth , and it does not affect the analytic continuation in @xmath1 . we can therefore perform the continuation collectively for all tensor integrals and diagrams of the same topology . we organize the evaluation of tensor integrals and diagrams of a topology as follows . we first find the mellin - barnes parameterization of the scalar integral of the topology , and multiply the integrand with a template function @xmath129 which we assume depends on all mellin - barnes variables . then , we apply the analytic continuation algorithm of section  [ sec : analytic ] keeping @xmath130 general . we expand in @xmath1 , and create numerical programs for the evaluation of the expansion for any smooth function @xmath130 . this part of the evaluation needs to be performed only once for a given topology . then we must identify the polynomial for the evaluation of the specific tensors or diagrams that we are interested in . we have written form  @xcite and mathematica programs which perform the steps that we described in this section , and derive the explicit functional form of the polynomials in eq . [ eq : tensor_mbrep ] . we then write numerical routines for their evaluation , and link them to the general routines for the @xmath1 expansion of the integrals of the topology . this approach is efficient for the application of our technique to tensor integrals . it turns out , however , that the expressions for the polynomials are quite long for high rank tensors . we have observed that we can reduce the size of the expressions substantially with a simple modification . the terms @xmath131 in eq . [ eq : tensor - feynf ] are lengthy , and when we integrate out the feynman parameters they give rise to large expressions in @xmath128 of eq . [ eq : tensor_mbrep ] . we rewrite eq . [ eq : tensor - feynf ] in a different way , @xmath132 } \,\end{aligned}\ ] ] then , we introduce a mellin - barnes decomposition for each denominator : @xmath133 and integrate out the feynman parameters . at the end of this procedure , we obtain a sum of @xmath134 + 1 $ ] mellin - barnes integrals for a tensor of rank @xmath104 . all integrals correspond to the mellin - barnes representation of the scalar integral with shifted dimension @xmath135 . each of them has a different polynomial @xmath136 in the numerator . however , these are substantially shorter than the polynomial @xmath128 in eq . [ eq : tensor_mbrep ] . it is worth to note that , with our method , we never introduce spurious singularities . the evaluation of the higher rank tensors and diagrams in gauge theories is very efficient , since we never require the manipulation of lengthy expressions . the polynomials @xmath136 can be large , however , they require a minimum amount of handling before generating numerical codes for their evaluation .    here we have presented the derivation of efficient representations for one - loop integrals . the method may be applied to multi - loop integrals using the one - loop tensor integral representations as building blocks . the analogous functions of @xmath136 in multi - loop tensors contain gamma functions in the denominator . however , they are smooth and maintain the properties which are required for a collective evaluation of all tensors of a topology . explicit applications with multi - loop tensor integrals will be examined in future work . after analytic continuation , the contour integrals can be safely expanded in power series of @xmath1 . the coefficients of the expansion are also contour integrals which , in general , contain in their integrand gamma functions , gamma function derivatives , simple logarithms , and powers of the kinematic parameters . in earlier works , these integrals were computed analytically . the standard approach is to evaluate them by means of the residue theorem , i.e. summing an infinite number of residues and expressing the power series as logarithms , polylogarithms and harmonic polylogarithms . this method is cumbersome for the evaluation of integrals with more than one mellin - barnes variable . in all cases presented in the literature , integrals with more than one dimension are reduced to one - dimensional integrals with a variety of clever , however not general , tricks . typical procedures include direct integration in simple situations , the application of barnes lemmas  @xcite , and the expansion in laurent series around an integer value of one of the variables . in a few cases , a careful inspection of the multi - dimensional integrals reveals exact cancellations between different components which yield significant simplifications . many integrals which are relevant for physical applications , such as the massless double box with two external legs off - shell and one - loop tensor integrals with six external legs , after the analytic continuation in @xmath1 , are expressed in terms of hundreds of mellin - barnes components ; some of them have a very high dimensionality . the analytic methods are not suitable for such problems . in our method , we choose to evaluate the mb integrals numerically by direct integration over the contours . the procedure we follow is , conceptually , straightforward . in practice , the implementation of the numerical evaluation requires a significant work in automating the book - keeping of the various terms . tensor integrals require significantly more complex book - keeping due to large expressions that multiply the gamma functions in the integrand of their mellin - barnes representation . our computer programs are organized in the following way . we first use maple , form  @xcite , and mathematica routines to derive mellin - barnes representations for the integrals that we need to compute . we then perform the analytic continuation of section  [ sec : analytic ] , and we expand the integrands in powers of @xmath1 using maple and mathematica . as the next step , we translate the integrands into fortran routines which evaluate them as complex quantities . there are many options to calculate numerically the sum of the mellin - barnes integrals which emerge from the routines for the @xmath1 expansion . for example , we can combine all contributions in a single common integrand , so that large numerical cancellations take place before integration . however , the finding of the peaks and the adaptation of the integration routines become less efficient . following a diametric approach , we could integrate all components separately and compute their sum at the end . then , the adaptation of the numerical integration is optimal , however , the numerical precision is sensitive to cancellations . in practice , we follow a hybrid approach and group together integrals with equivalent contours . we have achieved a reasonable efficiency and speed for all the computations that we present in this paper . since we aim to keep our integration method general , we do not attempt any special simplifications by , for example , applying barnes lemmas or exploiting symmetries of particular contours . for the numerical integration , we maintain the contours as straight lines parallel to the imaginary axis . we map them onto a real interval with a simple transformation , @xmath137 in all cases that we have considered , the integrands vanish rapidly when the integration variables take values away from the real line @xmath138 . the mellin - barnes integrals do not present problematic numerical instabilities , and converge rather fast . we perform the multidimensional integrations using the cuhre and vegas routines of the cuba library  @xcite . to gain some speed , we integrate the real and imaginary parts of the integrals concurrently , using the same grid . the grid is adapted to the peaks of the real part , however , the quality of the integration for the imaginary part is only mildly affected . on the other hand , there is another subtle point concerning the convergence of the imaginary parts . as they always start an order later in the expansion in @xmath1 , the analytic expressions for them would be simpler than the corresponding real ones . when integrating numerically this relative simplicity is reflected by a faster convergence when compared to the real pieces at the same order . the numerical results depend on the values for the kinematic scales of the integrals . with our approach , we use the same expressions for the evaluation of the integrals in all kinematic regions . as we explained earlier , the only analytic continuations in kinematic variables that we must perform , is for simple functions of the type @xmath139 and @xmath140 . these are implemented trivially , in a branched format , in numerical programs . in this paper , we present multi - loop integrals for fixed kinematic parameters . however , we have arranged our routines so that a numerical integration over the phase - space can be performed simultaneously with the mellin - barnes integrations . a problematic case for our numerical algorithms is the evaluation of loop integrals with massive internal propagators in non - euclidean kinematic regions . typically , mellin - barnes represenations for such integrals exhibit a very slow damping at @xmath141 in physical regions . the mapping of the integration region to a finite interval that we have considered in eq . [ eq : mapping ] is not adequate , and contour deformations together with more sophisticated algorithms for the evaluation of oscillatory integrals are required . we defer the study of such situations for a future project . we have achieved a full automatization of the code generation for the numerical evaluation . this has allowed us to apply our method to diverse problems . our approach can be improved and refined further in the future by considering , for example , a clever analysis of the terms before integration in order to avoid numerical cancellations or more sophisticated mappings which can smoothen the peak structure of the integrands . however , as we will see in the rest of the paper , our first , naive , implementation is already sufficient to obtain accurate results for complicated integrals in a reasonable time .    in what follows we present results for one , two , and three loop integrals . some of the integrals are known analytically , and we use them to verify our numerical programs . we also present new loop integrals , which would be extremely tedious to calculate with traditional methods . in this section we present results for one , two , and three - loop integrals that we obtain with our method . for each integral , we present the coefficients of its expansion in @xmath142 through the finite term : @xmath143 where @xmath116 is the depth of the leading pole in @xmath1 . we estimate the errors associated with our results for @xmath144 , by adding , in quadrature , the integration errors of all contributions from the residue decomposition . the first integral we consider is the massless one - loop hexagon , depicted in figure [ fig : hexagon ] .        for massless external legs , we derive an eight - dimensional mellin - barnes representation for this integral . due to the high dimensionality , the analytic continuation in @xmath1 involves the crossing of the contours of integration from many poles , and the mellin - barnes expression after the continuation consists of about two hundred terms . this example demonstrates the advantage of the automatic algorithm we described in section [ sec : analytic ] , since the book - keeping is done automatically and our routines perform the @xmath1 expansion in fractions of a minute .    after the analytic continuation and the expansion in series of @xmath1 , only integrals with low dimensionality contribute to the poles and finite pieces of the series . the integration over feynman parameters produces a denominator factor @xmath145 , which contributes factors of @xmath1 in the series expansion . the eight - dimensional component , in which no residues have been taken , drops out from the finite part . in addition , most of the residues in the analytic continuation do not give rise to poles in @xmath1 ; a few components only develop a singularity that compensates the @xmath1 in the overall factor . for the scalar hexagon , we find integrals with up to three dimensions which contribute to the poles and the finite term . this mechanism , which reduces the dimensionality of the representation in the leading coefficients of the @xmath1 expansion , is not specific to the hexagon topology , but is typical of mellin - barnes representations of scalar multi - leg integrals .    in table [ tb : hex - scalar.physical ] we present explicit results for the scalar hexagon in phase space points corresponding to the physical region for the @xmath146 process . to check the efficiency of our integration routines in all the phase - space , we have sampled hundreds of points . here , we only show representative results in the phase - space points of table [ tb : hexagon - phasepoints ] . these points were generated starting from explicit values for the momenta in four dimensions , so they have a vanishing gram determinant . .explicit values of the invariants for some of the phase space points used for the evaluation of the hexagon topology . the points correspond to the physical region of the @xmath146 scattering of massless particles . [ cols="^,^,^,^,^,^,^,^,^,^",options="header " , ]         we also present our results in the kinematical region corresponding to the process @xmath147 shown in figure [ fig:2b2mass ] , relevant , for instance , for the process @xmath148 . in figure [ fig:2off2bem ] we plot the values of the coefficients , for this configuration , up to the finite terms as a function of @xmath49 in the case of two equal masses , @xmath149 , where the energy scale is fixed by @xmath150 . figure [ fig:2off2bm10 ] corresponds to the same kinematical region but for @xmath150 , @xmath151 and @xmath152 . the error bars again lie within the points in the plot and are , for the finite pieces , typically less that 1% after a run of a couple of minutes per point on a 2.8ghz cpu . process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath149.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]     process where the two particles in the final state are massive . we plot the values of the coefficients as function of @xmath49 for @xmath150 , @xmath151 and @xmath152.,width=226 ]      at the three loop level , we start by considering the on - shell triple box ( figure [ fig:3bonshell ] ) , impressively computed in an analytic way in @xcite using the mb technique .    using the re - insertion method described in section ii , we get a 7-fold mb representation . the analytic continuation leaves us with contributions with up to 5 contour integrals for the constant pieces in the @xmath1 expansion , and poles up to @xmath153 . although the numerical integration is certainly more involved in this case , due to the depth of the poles in @xmath1 which generates very complicated expressions for the finite terms , it is relatively easy to achieve a precision of @xmath154 , as shown in figure [ fig : on3bem ] . as in the two loop case we show results for the physical region corresponding to @xmath155 , @xmath156 with @xmath157 and @xmath158 . we compare our results with the analytic computation of @xcite . we used the mathematica package hpl @xcite for the evaluation of the one - dimensional harmonic polylogarithms .         process . on the upper and lower left panel we plot the real and imaginary parts , respectively , of @xmath159 as a function of @xmath49 for fixed value of @xmath150 . on the right panels we show the corresponding ratios of the numerical calculation to the analytic result of @xcite for the same kinematics , the bands are given by the error in the numerical integrations.,width=302 ]     process . on the upper and lower left panel we plot the real and imaginary parts , respectively , of @xmath159 as a function of @xmath49 for fixed value of @xmath150 . on the right panels we show the corresponding ratios of the numerical calculation to the analytic result of @xcite for the same kinematics , the bands are given by the error in the numerical integrations.,width=302 ]     process . on the upper and lower left panel we plot the real and imaginary parts , respectively , of @xmath159 as a function of @xmath49 for fixed value of @xmath150 . on the right panels we show the corresponding ratios of the numerical calculation to the analytic result of @xcite for the same kinematics , the bands are given by the error in the numerical integrations.,width=302 ]     process . on the upper and lower left panel we plot the real and imaginary parts , respectively , of @xmath159 as a function of @xmath49 for fixed value of @xmath150 . on the right panels we show the corresponding ratios of the numerical calculation to the analytic result of @xcite for the same kinematics , the bands are given by the error in the numerical integrations.,width=302 ]      as we did for the double box topology , we consider now a further step in complication for the triple box , adding a mass in the external state . this is , as far as we know , the first evaluation of a three - loop box with three mass scales . we obtained a mb representation with 8 variables . as usual , the dimensionality of the problem is reduced after the analytic continuation and we found that , up to order @xmath160 , only up to six - fold integrals contribute . we present results for this novel integral in the physical region of the decay of a heavy particle , @xmath161 ( figure [ fig:3bdecay ] ) as we did in the two loops case . in figure [ fig:1off3bm ] we show the the corresponding coefficients of the expansion in @xmath1 as a function of the invariant @xmath162 for fixed values of the other two variables : @xmath163 and @xmath164 . with the notable exception of the left - most point in the real part of the constant term , the error bars are all contained in the plotted points . they are typically of the order of @xmath165 or better for the finite parts after 50 minutes , in average , per point on a 2.8ghz cpu .    , width=151 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     ( figure [ fig:3bdecay ] ) . we show the results for fixed values of @xmath163 and @xmath164 as a function of the remaining invariant @xmath162.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ]     process with the massive leg on the final state . we plot the real and imaginary parts as a function of the invariant @xmath158 for fixed value of @xmath166 and @xmath167.,width=226 ] in this paper , we have introduced a new method for the numerical evaluation of arbitrary loop integrals . we first derive mellin - barnes representations by either decomposing directly the feynman representation , or by using one - loop representations as building blocks for representations of multi - loop integrals . for a given representation , we find values of the space - time dimension and the powers of the propagators which yield a well defined representation . then we use automated programs to perform the analytic continuation to values of the parameters for which the integral may develop divergences . finally , we expand in @xmath1 and compute the integral coefficients of the expansion numerically . the method is based on the earlier work of smirnov  @xcite and tausk  @xcite . the novel features in our publication are : ( i ) the automation of the procedure for the @xmath1 expansion , ( ii ) the numerical evaluation of the coefficients of the expansion , avoiding cumbersome summations of multiple infinite series and the analytic continuation in the arguments of polylogarithms , ( iii ) and the efficient generalization of the method to tensor integrals . we believe that these new features render our method suitable for the practical evaluation of multi - loop amplitudes in gauge theories .    we have performed a number of explicit calculations to verify our algorithms and demonstrate the power of our method . we first recalculated complicated loop integrals which are only recently known in the literature . namely , we evaluated the planar and cross on - shell double boxes , the planar double - box with one leg off - shell , and the triple planar box with on - shell legs . in all cases , we have found an excellent agreement between our numerical results and the known analytic results . the analytic results require complicated analytic continuations of polylogarithms to non - euclidean kinematic regions . with our method , we were able to compute the integrals in all kinematic regions effortlessly .    in this paper , we have presented results for loop - integrals which were previously unknown . at two - loops , we have computed the double - box integral with two adjacent legs off - shell for independent external mass - parameters . this is one of the most complicated two - loop box integrals that enters the evaluation of two - loop amplitudes for heavy boson pair production at colliders . at three - loops , we have computed the planar triple - box with one off - shell leg ; this integral emerges in @xmath168 , or the production of a single heavy boson in association with a jet at hadron colliders at nnnlo in qcd . to the best of our knowledge , it is the first time that a two - loop box with four kinematic scales or a three - loop box with three scales are ever computed in all physical regions . many processes with six external legs are particularly important at the lhc . at present , there has been no nlo calculation of a cross - section for a hadron collider process with six external states , due to the lack of efficient methods for evaluating loop amplitudes . in this paper , we applied our method to the evaluation of tensors through rank six for the hexagon topology . the purpose of this application was two - fold ; first , we wanted to demonstrate that we were able to extend the method to tensor integrals economically , and second , to setup programs that are efficient for the evaluation of multi - scale one - loop amplitudes ( e.g. six parton qcd amplitudes ) . we have found that the evaluation of the tensor integrals is not significantly more complicated than the scalar integral . the programs that we developed can be immediately used for the evaluation of the six parton one - loop qcd amplitudes . our technical goals for the current study were to build the necessary programs , implement an efficient book - keeping platform , and to examine the scaling of computing intensity in applications with diverse features . we anticipate that our method can be improved and generalized even further , when we consider issues that were only briefly examined in this first study . for example , in preliminary investigations , we have found that it is possible to exploit further the cauchy theorem . if we specify a kinematic region , it turns out that some of the numerical integrals can be approximated very accurately by the sum of a finite number of residues . we expect that we can improve substantially our efficiency if we replace some of the integrations which involve kinematic scales with appropriate finite sums of such residues . in addition , we are investigating hybrid approaches , combining our method with ideas in ref . @xcite . the applications of our method are numerous . we believe that it will be particularly useful for precision calculations of observables in collider physics . we plan to apply our techniques to key nlo and nnlo calculations for the lhc and to answer more formal questions on the perturbative behavior of gauge theories in the near future . we are indebted to vladimir smirnov , bas tausk , simon schwarz and jeppe andersen for extensive discussions on the topic . we would like to thank aude gehrmann - de ridder , thomas gehrmann , gudrun heinrich , zoltan kunszt , zoltan nagy , and bas tausk for their suggestions and for reading the manuscript . we are grateful to thomas gehrmann for providing us results and numerical routines from refs . @xcite , and daniel matre for discussions about his program to evaluate harmonic polylogarithms  @xcite . this work was supported in part by the swiss national science foundation ( snf ) under contract number 200020 - 109162 and by the forschungskredit der universitt zrich . j.  campbell , r.  k.  ellis , f.  maltoni and s.  willenbrock , arxiv : hep - ph/0510362 . l.  w.  garland , t.  gehrmann , e.  w.  n.  glover , a.  koukoutsakis and e.  remiddi , nucl . b * 642 * , 227 ( 2002 ) [ arxiv : hep - ph/0206067 ] . g.  passarino and m.  j.  g.  veltman , nucl . b * 160 * , 151 ( 1979 ) . f.  v.  tkachov , phys . b * 100 * , 65 ( 1981 ) ; k.  g.  chetyrkin and f.  v.  tkachov , nucl . b * 192 * , 159 ( 1981 ) . t.  gehrmann and e.  remiddi , nucl . b * 580 * , 485 ( 2000 ) [ arxiv : hep - ph/9912329 ] . z.  bern , l.  j.  dixon and d.  a.  kosower , nucl . b * 412 * , 751 ( 1994 ) [ arxiv : hep - ph/9306240 ] . a.  v.  kotikov , phys . b * 254 * , 158 ( 1991 ) . w.  t.  giele and e.  w.  n.  glover , jhep * 0404 * , 029 ( 2004 ) [ arxiv : hep - ph/0402152 ] ; r.  k.  ellis , w.  t.  giele and g.  zanderighi , arxiv : hep - ph/0508308 ; r.  k.  ellis , w.  t.  giele and g.  zanderighi , phys . d * 72 * , 054018 ( 2005 ) [ arxiv : hep - ph/0506196 ] . a.  denner and s.  dittmaier , arxiv : hep - ph/0509141 ; a.  denner and s.  dittmaier , nucl . b * 658 * , 175 ( 2003 ) [ arxiv : hep - ph/0212259 ] ;    t.  binoth , j.  p.  guillet , g.  heinrich , e.  pilon and c.  schubert , arxiv : hep - ph/0504267 ; t.  binoth , g.  heinrich and n.  kauer , nucl .  phys . b * 654 * , 277 ( 2003 ) [ arxiv : hep - ph/0210023 ] ; t.  binoth , j.  p.  guillet , g.  heinrich and c.  schubert , nucl . b * 615 * , 385 ( 2001 ) [ arxiv : hep - ph/0106243 ] ; t.  binoth , j.  p.  guillet and g.  heinrich , nucl . b * 572 * , 361 ( 2000 ) [ arxiv : hep - ph/9911342 ] . a.  van hameren , j.  vollinga and s.  weinzierl , eur . j.  c * 41 * , 361 ( 2005 ) [ arxiv : hep - ph/0502165 ] . f.  del aguila and r.  pittau , jhep * 0407 * , 017 ( 2004 ) [ arxiv : hep - ph/0404120 ] . a.  denner , s.  dittmaier , m.  roth and l.  h.  wieders , arxiv : hep - ph/0505042 ; a.  denner , s.  dittmaier , m.  roth and l.  h.  wieders , phys . b * 612 * , 223 ( 2005 ) [ arxiv : hep - ph/0502063 ] . z.  nagy and d.  e.  soper , jhep * 0309 * , 055 ( 2003 ) [ arxiv : hep - ph/0308127 ] . z.  bern , l.  j.  dixon and d.  a.  kosower , arxiv : hep - ph/0507005 . s.  laporta , int . j.  mod . a * 15 * , 5087 ( 2000 ) [ arxiv : hep - ph/0102033 ] . c.  anastasiou and a.  lazopoulos , jhep * 0407 * , 046 ( 2004 ) [ arxiv : hep - ph/0404258 ] . a.  v.  smirnov and v.  a.  smirnov , arxiv : hep - lat/0509187 ; o.  v.  tarasov , nucl . instrum . a * 534 * , 293 ( 2004 ) [ arxiv : hep - ph/0403253 ] ; o.  v.  tarasov , acta phys . b * 29 * , 2655 ( 1998 ) [ arxiv : hep - ph/9812250 ] ; p.  a.  baikov , arxiv : hep - ph/0507053 ; p.  a.  baikov , nucl .  phys .  proc . * 116 * ( 2003 ) 378 . v.  a.  smirnov and m.  steinhauser , nucl . b * 672 * , 199 ( 2003 ) [ arxiv : hep - ph/0307088 ] . s.  actis , a.  ferroglia , g.  passarino , m.  passera and s.  uccirati , nucl . b * 703 * , 3 ( 2004 ) [ arxiv : hep - ph/0402132 ] ; a.  ferroglia , m.  passera , g.  passarino and s.  uccirati , nucl . b * 680 * , 199 ( 2004 ) [ arxiv : hep - ph/0311186 ] ; a.  ferroglia , m.  passera , g.  passarino and s.  uccirati , nucl . b * 650 * , 162 ( 2003 ) [ arxiv : hep - ph/0209219 ] ; g.  passarino and s.  uccirati , nucl . b * 629 * , 97 ( 2002 ) [ arxiv : hep - ph/0112004 ] . s.  p.  martin and d.  g.  robertson , arxiv : hep - ph/0501132 . m.  caffo , h.  czyz and e.  remiddi , nucl . suppl .   * 116 * , 422 ( 2003 ) [ arxiv : hep - ph/0211178 ] ; s.  laporta , phys . b * 504 * , 188 ( 2001 ) [ arxiv : hep - ph/0102032 ] . t.  binoth and g.  heinrich , nucl . b * 585 * , 741 ( 2000 ) [ arxiv : hep - ph/0004013 ] ; t.  binoth and g.  heinrich , nucl . b * 680 * , 375 ( 2004 ) [ arxiv : hep - ph/0305234 ] . k.  hepp , commun . * 2 * , 301 ( 1966 ) . m.  roth and a.  denner , `` high - energy approximation of one - loop feynman integrals , '' nucl . b * 479 * , 495 ( 1996 ) [ arxiv : hep - ph/9605420 ] . t.  binoth , g.  heinrich and n.  kauer , nucl . b * 654 * , 277 ( 2003 ) [ arxiv : hep - ph/0210023 ] . c.  anastasiou , k.  melnikov and f.  petriello , arxiv : hep - ph/0505069 . v.  a.  smirnov , phys . b * 460 * , 397 ( 1999 ) [ arxiv : hep - ph/9905323 ] . j.  b.  tausk , phys . b * 469 * , 225 ( 1999 ) [ arxiv : hep - ph/9909506 ] . v.  a.  smirnov and o.  l.  veretin , nucl . b * 566 * , 469 ( 2000 ) [ arxiv : hep - ph/9907385 ] . c.  anastasiou , j.  b.  tausk and m.  e.  tejeda - yeomans , nucl . suppl .   * 89 * , 262 ( 2000 ) [ arxiv : hep - ph/0005328 ] . v.  a.  smirnov , phys . b * 491 * , 130 ( 2000 ) [ arxiv : hep - ph/0007032 ] ; v.  a.  smirnov , phys . b * 500 * , 330 ( 2001 ) [ arxiv : hep - ph/0011056 ] ; v.  a.  smirnov , phys . b * 524 * , 129 ( 2002 ) [ arxiv : hep - ph/0111160 ] ; g.  heinrich and v.  a.  smirnov , phys . b * 598 * , 55 ( 2004 ) [ arxiv : hep - ph/0406053 ] . v.  a.  smirnov , phys . b * 567 * , 193 ( 2003 ) [ arxiv : hep - ph/0305142 ] . z.  bern , l.  j.  dixon and v.  a.  smirnov , phys . d * 72 * , 085001 ( 2005 ) [ arxiv : hep - th/0505205 ] . c.  anastasiou , z.  bern , l.  j.  dixon and d.  a.  kosower , phys . * 91 * , 251602 ( 2003 ) [ arxiv : hep - th/0309040 ] . e.  remiddi and j.  a.  m.  vermaseren , int . j.  mod . a * 15 * , 725 ( 2000 ) [ arxiv : hep - ph/9905237 ] . ussyukina , teor.mat.fiz . 22 ( 1975 ) 300 ; e.  e.  boos and a.  i.  davydychev , theor . * 89 * , 1052 ( 1991 ) [ teor .  mat . fiz .   * 89 * , 56 ( 1991 ) ] ; n.  i.  ussyukina and a.  i.  davydychev , phys . b * 305 * , 136 ( 1993 ) ; n.  i.  ussyukina and a.  i.  davydychev , phys . b * 298 * , 363 ( 1993 ) . t.  gehrmann and e.  remiddi , nucl . b * 601 * , 248 ( 2001 ) [ arxiv : hep - ph/0008287 ] . t.  gehrmann and e.  remiddi , comput . commun .   * 141 * , 296 ( 2001 ) [ arxiv : hep - ph/0107173 ] . t.  gehrmann and e.  remiddi , comput . phys .  commun . * 144 * , 200 ( 2002 ) [ arxiv : hep - ph/0111255 ] . t.  gehrmann and e.  remiddi , nucl . b * 640 * , 379 ( 2002 ) [ arxiv : hep - ph/0207020 ] .
Abstract:

Answer: we present a new method for the numerical evaluation of arbitrary loop integrals in dimensional regularization . we first derive mellin - barnes integral representations and apply an algorithmic technique , based on the cauchy theorem , to extract the divergent parts in the @xmath0 limit . we then perform an @xmath1-expansion and evaluate the integral coefficients of the expansion numerically . the method yields stable results in physical kinematic regions avoiding intricate analytic continuations . it can also be applied to evaluate both scalar and tensor integrals without employing reduction methods . we demonstrate our method with specific examples of infrared divergent integrals with many kinematic scales , such as two - loop and three - loop box integrals and tensor integrals of rank six for the one - loop hexagon topology .