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: perturbation theory remains the most accurate and reliable tool of theoretical analysis in elementary particle phenomenology . even the most promising approach based on the direct calculation of physical observables on the lattice still can not compete with perturbation theory supplemented with some non - perturbative additions as condensates or renormalons . the accuracy of experimental data for tests of the standard model permanently improves and demands the improvement of the accuracy of theoretical predictions  @xcite . within perturbation theory this demand implies the calculation of feynman diagrams with an increasing number of loops ( see e.g.  @xcite ) . at present , feynman diagrams of rather general form can be obtained analytically at most up to three loops . any of these calculations is still unique and involves sophisticated algorithms and an extensive usage of computer resources ( see e.g.  @xcite ) . even the notion of the analytical calculation is changing now because the analytical results are expressed through rather special transcendental numbers which eventually have to be calculated numerically . the striking example is the calculation of massive three - loop bubbles  @xcite where the analytical results for the diagrams are represented through the quantity which is given by a two - fold infinite sum ( see also  @xcite ) . not all of the multiloop feynman diagrams are equally complicated . there are some topologies of @xmath3-loop diagrams or special kinematic regimes for the external momenta which are much simpler to compute than the general @xmath3-loop case ( see e.g.  @xcite ) . the increasing complication of loop calculations within perturbation theory is more and more circumvented by using different kinds of expansion for the cases where the kinematics allows one to find a small parameter . heavy quark effective theory ( hqet ) and threshold expansions for heavy particle production are recent examples for such a situation . these expansions are mostly used for phenomenological applications where it even suffices to give approximate numerical values for the unexpanded multiloop integrals . in some cases one can obtain an analytical expansion that can be used for the simplification of the calculations or for the evaluation of the exact parameters of the effective lagrangians through matching conditions .    in this paper we discuss a method to calculate the near threshold expansion of the spectral density for the class of multiloop sunset - type feynman diagrams termed water melon diagrams  @xcite . these diagrams has recently drawn some attention as a laboratory for testing some advanced methods in loop calculations  @xcite . however , water melon diagrams also have numerous phenomenological applications . the most important application of these diagrams ( as well as spectacle diagrams as their first generalization ) is the calculation of the effective potential in quantum field theory both at finite temperature and at zero temperature  @xcite . among other applications one can name the sum rule analysis of baryons in qcd both in the massless approximation and with finite mass ( the leading order correlator being the standard sunset diagram ) @xcite or the treatment of baryons in the large @xmath4 limit of qcd  @xcite which requires the calculation of @xmath5-loop water melon diagrams in the leading order . an interesting application is the calculation of dibaryon properties in operator product expansion within the sum rule approach which is important for understanding many - quark bags in the nucleus  @xcite . the water melon diagrams emerge in chiral perturbation theory  @xcite . the corresponding integrals appear for the correlator of the effective operators related to the mixing of neutral mesons in flavour dynamics  @xcite and in the sum rule analysis of hybrid mesons  @xcite . water melon diagrams constitute an important part of the contributions generated by the recurrence relations of the integration - by - part techniques for three - loop diagrams  @xcite ( for the recent progress see e.g.ref . it is a rather hot topic and a dynamically developing area of loop calculations where new results are frequently reported ( see e.g.refs . @xcite ) .    in our opinion , the method presented in ref . @xcite completely solves the problem of computing this class of diagrams . the method is simple and reduces the calculation of a @xmath3-loop water melon diagram to a one - dimensional integral which includes only well - known special functions in the integrand ( bessel functions of different kinds ) for any values of the internal masses . the technique is universal . the most interesting part of our analysis of water melon diagrams is the construction of the spectral decomposition of water melon diagrams , i.e.  the determination of the discontinuity across the physical cut in the complex plane of the squared momentum . in this context a novel technique for the direct construction of the spectral density of water melon diagrams based on an integral transform in configuration space was presented in ref . @xcite . in the present paper we develop some explicit expansions for water melon diagrams and compare them with the exact results . the purpose of this paper is a three - fold one , namely :    * to demonstrate the ease with which threshold expansions to arbitrary space - time dimensions and a general number of internal lines with arbitrary masses can be generated within a configuration space technique . the entire construction reduces to algebraic operations while the only integral encountered is a simple integral of the type of euler s gamma function . * to generate explicit forms of threshold expansions and to analyze their convergence properties . the explicit forms of threshold expansions for the simple sunset can be compared with threshold expansions that are obtained in momentum space . * to consider a case when one mass is much smaller than the others . this is an important generalization of the threshold expansion and can be analytically done within the configuration space technique . the present analysis of the threshold expansion of water melon diagrams shows the effectiveness of the configuration space technique for this topology class of feynman diagrams . on the other hand the obtained explicit results can be useful for a variety of applications which include the evaluation of water melon type diagrams . the technique can readily be generalized to some other simple topologies of @xmath3-loop diagrams . the paper is organized as follows . in sec . 2 we provide the tools necessary for the following calculations . sec .  3 deals with the threshold expansion . here we outline our main strategy for this expansion in configuration space and in the subsections that follow we give examples for the sunset diagram and the water melon diagram with four or more internal lines . in sec . 4 we concentrate on the strongly asymmetric case where one of the masses is much smaller than the others . we introduce the procedure of resummation of the contributions of this smallest mass and show explicitly how we get to a closed expression for the spectral density even in this case . we give examples for a strongly asymmetric mass arrangement for the water melon diagrams with two ( degenerate case ) , three , and four or more internal lines , and compare the resummed results with the pure threshold expansion as well as with the exact spectral density . in sec . 5 we discuss how to recover the non - analytic part of the polarization function through the spectral density near threshold . in sec . 6 we give our conclusions . we start with a brief outline of the technique which we use in this paper  @xcite . the polarization function @xmath6 of a water melon diagram including @xmath3 internal lines with masses @xmath7 , @xmath8 in configuration space is given by the product @xmath9 the propagator @xmath10 of a massive line with mass @xmath11 in @xmath12-dimensional ( euclidean ) space - time is given by @xmath13 where we write @xmath14 . @xmath15 is a mcdonald function ( a modified bessel function of the third kind , see e.g.  @xcite ) . in the limit @xmath16 the propagator in eq . ( [ prop ] ) simplifies to @xmath17 where @xmath18 is euler s gamma function . ( [ prod ] ) contains all information about the water melon diagrams and in this sense is the final result for the class of diagrams under consideration . it is of course known and was used since long ago  @xcite . of some particular interest is the spectral decomposition of the polarization function @xmath6 which is connected to the particle content of a given model . the spectrum of particles is contained in the function @xmath19 related to the polarization function through the dispersion representation @xmath20 the dispersion representation of the polarization function in configuration space reveals the analytic structure of the polarization function @xmath6 . for applications , however , one may need the fourier transform of the polarization function @xmath6 given by @xmath21 ( we use the same notation for the function and its fourier transform because we think that this will cause no confusion ) . because of the lorentz invariance of the propagator the angular integration in eq . ( [ curr ] ) can be done explicitly in @xmath12-dimensional space - time ( for a generalization to tensor propagators see ref . the result reads @xmath22 where @xmath23 is the usual bessel function and @xmath24 is the rotationally invariant measure on the unit sphere in the @xmath12-dimensional ( euclidean ) space - time . note that the polarization function @xmath6 as well as its fourier transform @xmath25 are only functions of the absolute value @xmath26 and @xmath27 . the same is of course valid also for the other occurring functions . to simplify the notation we often write @xmath28 and @xmath29 for these absolute values . note that propagators of particles with non - zero spin in configuration space representation can be obtained from the scalar propagator by differentiation with respect to the space - time point @xmath30 . this does not change the functional @xmath30-structure and causes only minor modification of the basic technique ( for details see ref . @xcite and references therein ) . our final representation of the fourier transform of a water melon diagram is given by the one - dimensional integral @xmath31 which is a special kind of integral transformation with a bessel function as a kernel . this integral transformation is known as the hankel transform  @xcite . the representation given by eq . ( [ exactpi ] ) is quite universal regardless of whether tensor structures are added or particles with vanishing momenta are radiated from any of the internal lines . note in this context that bessel functions are objects well - studied during the last two centuries . they therefore can be added to the list of elementary functions .    because the dispersion representation of the polarization function in configuration space ( or the spectral density of the corresponding polarization operator ) has the form given in eq . ( [ convo1 ] ) , the analytic structure of the polarization function @xmath6 can be determined directly in configuration space without having to compute its fourier transform first . the transformation in eq . ( [ convo1 ] ) turns out to be a particular example of the hankel transform , namely the @xmath32-transform  @xcite . as pointed out in ref . @xcite , the inverse @xmath32-transform in this case is given by @xmath33 where @xmath34 is a modified bessel function of the first kind and the integration runs along a vertical contour in the complex plane to the right of the right - most singularity of @xmath35 . the inverse transform given by eq . ( [ exactrho ] ) completely solves the problem of determining the spectral density of the general class of water melon diagrams with any number of internal lines and different masses by reducing it to the computation of a one - dimensional integral . for @xmath3 internal lines with equal masses @xmath11 the spectral density reads explicitly @xmath36 in contrast to our technique , in the standard , or momentum , representation the polarization function @xmath25 is calculated from a @xmath37-loop diagram with @xmath37 @xmath12-dimensional integrations over the entangled loop momenta which makes the computation difficult when the number of internal lines becomes large . with our method described in detail in ref . @xcite the @xmath38-dependence of the spectral density can be calculated by a one - fold numerical integration according to eq . ( [ exactrho ] ) . the numerical integration in eq . ( [ exactrho ] ) can be done for arbitrary space - time dimensions and a general number of lines with arbitrary masses . in this sense this is the most efficient representation for the spectral density of the water melon diagram . however , we can also develop an explicit expansion near the threshold with any desired accuracy . it does not require any complicated integrations at all . the corresponding expansion of eq . ( [ exactrho ] ) can then be compared with series expansions near the production threshold obtained with the traditional momentum space technique . we stress that we are only interested in the spectral density because it is the main object for physical applications ( see e.g.  refs . @xcite ) . for practical reasons we start with eq . ( [ exactpi ] ) . the polarization function @xmath25 as written in eq . ( [ exactpi ] ) is , in general , uv divergent . the divergence can be subtracted by using the power series expansion of the weight function @xmath39 to an appropriate order which will be added and subtracted to this weight . this leads to a @xmath40-dependent power series of divergent subtraction terms plus an uv finite subtracted integral ( see e.g.  @xcite ) . but because the subtraction terms do not contribute to the spectral density , we can avoid this subtraction at all . in order that the formally written expressions make sense they are supposed to be dimensionally regularized . we use the simplified or unorthodox dimensional regularization method for water melon diagrams ( see ref . @xcite ) . the threshold region of a water melon diagram is determined by the condition @xmath41 where @xmath42 is the euclidean momentum and @xmath43 is the threshold value for the spectral density . we introduce the minkowskian momentum @xmath44 defined by @xmath45 which is an analytic continuation to the physical cut . operationally this analytic continuation can be performed by replacing @xmath46 . to analyze the region near the threshold we use the parameter @xmath47 which takes complex values . the parameter @xmath48 is more convenient in euclidean domain while the parameter @xmath49 is the actual energy counted from threshold which is used in phenomenological applications . the spectral density as a function of @xmath50 is written as @xmath51 in the following . the analytic continuation of the fourier transform in eq . ( [ exactpi ] ) to the minkowskian domain has the form @xmath52 for the threshold expansion we have to analyze the large @xmath30 behaviour of the integrand . it is this region that saturates the integral in the limit @xmath53 or , equivalently , @xmath54 . it is convenient to perform the analysis in a basis where the integrand has a simple large @xmath30 behaviour . the most important part of the integrand is the bessel function @xmath55 which , however , contains both rising and falling branches at large @xmath30 . it resembles the situation with elementary trigonometric functions @xmath56 and @xmath57 to which the bessel function @xmath23 is rather close ( in a certain sense ) . indeed , @xmath57 ( or @xmath56 ) is a linear combination of exponentials , namely @xmath58 and has also both rising and falling branches at large pure imaginary argument : the exponentials show simple asymptotic behaviour @xmath59 at @xmath60 . the analogous statement is true for @xmath23 which can be written as a sum of two hankel functions , @xmath61 where @xmath62 . the hankel functions @xmath63 show simple asymptotic behaviour at infinity , @xmath64 accordingly we split up @xmath65 into @xmath66 with @xmath67 the two parts @xmath68 of the polarization function @xmath65 have completely different behaviour near threshold which allows one to analyze them independently . this observation makes the subsequent analysis straightforward . we first consider the contribution of the part @xmath69 . the behaviour at large @xmath30 is given by the asymptotic form of the functions which we simply write up to the leading terms as @xmath70 the large @xmath30 range of the integral ( above a reasonably large cutoff parameter @xmath71 ) has the general form @xmath72 where @xmath73 the right hand side of eq . ( [ pilap ] ) is an analytic function in @xmath48 in the vicinity of @xmath74 . it exhibits no cut or other singularities near the threshold and therefore does not contribute to the spectral density . we turn now to the second part @xmath75 . in contrast to the previous case , the integrand of this part contains @xmath76 which behaves like a rising exponential function at large @xmath30 , @xmath77 therefore the integral is represented by @xmath78 the function @xmath79 is non - analytic near @xmath74 because for @xmath80 the integrand in eq . ( [ pilam ] ) grows in the large @xmath30 region and the integral diverges at the upper limit . therefore the function which is determined by this integral is singular at @xmath80 ( @xmath81 ) and requires an interpretation for these values of the argument @xmath48 . in the complex @xmath48 plane with a cut along the negative axis the function is analytic . this cut corresponds to the physical positive energy cut . the discontinuity across the cut gives rise to the non - vanishing spectral density of the diagram . let us first discuss the analytic part @xmath69 of the diagram . this part reduces to a regular water melon diagram . indeed , using the relation @xmath82 between bessel functions of different kinds one can replace the hankel function @xmath83 with the mcdonald function @xmath84 . since the propagator of a massive particle ( massive line in the diagram ) is given by the mcdonald function up to a power in @xmath30 , this substitution shows that the weight function behaves like a propagator of an additional line with the `` mass '' @xmath44 . the explicit expression is given by @xmath85 @xmath86 is the polarization function of a new effective diagram which is equal to the initial polarization function multiplied by a propagator with @xmath44 as mass parameter . we thus end up with a vacuum bubble of the water melon type with one additional line compared to the initial diagram ( see fig .  [ fig1 ] ) . these diagrams have no singular behaviour at the production threshold @xmath87 . as mentioned above , @xmath69 is analytic in @xmath48 near the origin @xmath74 and can therefore be omitted in the calculation of the spectral density . all derivatives of @xmath88 with respect to @xmath48 are represented as vacuum bubbles with one additional line carrying rising indices . such diagrams can be efficiently calculated within the recurrence relation technique developed in ref . @xcite . therefore one is left with the second part @xmath75 containing the non - analytical contributions in @xmath48 which leads to a non - vanishing spectral density . still the integral in eq . ( [ pm ] ) can not be done analytically . in order to obtain an expansion for the spectral density near the threshold in an analytical form we make use of the asymptotic series expansion for the function @xmath6 which crucially simplifies the integrands but still preserves the singular structure of the integral in terms of the variable @xmath48 . the asymptotic series expansion to order @xmath89 of the main part of each propagator , i.e.  of the mcdonald function , is given by  @xcite @xmath90,\quad    ( \lambda , n):=\frac{\gamma(\lambda+n-1/2)}{n!\gamma(\lambda - n-1/2)}\ ] ] ( @xmath91 $ ] ) . therefore the asymptotic expansion of the function @xmath6 consists of an exponential factor @xmath92 and an inverse power series in @xmath30 up to an order @xmath93 which is closely related to @xmath89 . it is this asymptotic expansion that determines the singularity structure of the integral . we write the whole integral in the form of the sum of two terms , @xmath94 the integrand of the first term @xmath95 behaves as @xmath96 at large @xmath30 while the integrand of the second term accumulates all lower powers of the large @xmath30 expansion . note that only the large @xmath30 behaviour is essential for the near threshold expansion of the spectral density . this fact has been taken into account in eqs . ( [ pilap ] ) and  ( [ pilam ] ) where we introduced a cutoff @xmath71 . however , from the practical point of view the calculation of the regularized integrals with an explicit cutoff is inconvenient . the final result of the calculation  the spectral density of the diagram  is independent of the cutoff , but the integration is technically complicated if the cutoff is introduced . however , in extending the integration over the whole region of the variable @xmath30 without using the cutoff one immediately encounters divergences at small @xmath30 because the asymptotic expansion is invalid in the region near the origin , so one is not allowed to continue the integration to this region . the standard way to cope with such a situation is to introduce dimensional regularization . it allows one to deal with divergent expressions at intermediate stages of the calculation and is technically simple because it does not introduce any cutoff and therefore does not modify the integration region drastically . note that dimensional regularization does not necessarily regularize all divergences in this case ( in contrast to the standard case of ultraviolet divergences ) but nevertheless suffices for our purposes . we use a parameter @xmath97 to regularize the divergences at small @xmath30 . this regularization prescription is an unorthodox version of the dimensional regularization ( see e.g.refs . @xcite ) .    the first part @xmath95 in eq . ( [ pidias ] ) containing a difference of the polarization function and its asymptotic expansion as the integrand gives no contributions to the spectral density up to a given order of the expansion in @xmath48 . this is because the subtracted asymptotic series to order @xmath89 cancels the inverse power behaviour of the integrand to this degree @xmath89 . the integrand decreases sufficiently fast for large values of @xmath30 and the integral converges even at @xmath74 . therefore this term is inessential when the expansion of the spectral density is evaluated up to some order . one can readily determine the order of the expansion near @xmath74 at which a contribution to the spectral density appears in using some further simplifications of the integrand of the term @xmath95 in eq . ( [ pidias ] ) . namely , we replace the hankel function under the integration sign by its asymptotic series expansion . the resulting exponential factor @xmath98 can then be expanded in the parameter @xmath47 and integrated together with the finite inverse power series in @xmath30 . one obtains a finite power series in this parameter @xmath48 which leads to a non - regular term of order @xmath99 ( for instance , @xmath100 or @xmath101 ) . therefore the part @xmath95 is regular and gives no contribution to the spectral density up to the order @xmath99 . for this reason we concentrate on the expansion of the second part @xmath102 and find that only this part contains the contribution to the spectral density up to the @xmath89 . therefore the expansion of the spectral density at small @xmath50 is determined only by the integral @xmath102 of eq . ( [ pidias ] ) . this integral is still rather complicated to compute but we can go a step further in its analytical evaluation . indeed , since the singular behaviour of @xmath102 is determined by the behaviour at large @xmath30 , we can replace the first factor , i.e.  the hankel function , in the large @xmath30 region by its asymptotic expansion up to some order @xmath89 . we use @xmath103\ ] ] ( cf . ( [ kas ] ) for the notation ) to obtain a representation @xmath104 the index `` _ das _ '' stands for `` double asymptotic '' and indicates that the integrand in eq . ( [ pidas ] ) consists of a product of two asymptotic expansions : one for the polarization function @xmath6 and another for the hankel function @xmath105 as weight ( or kernel ) . both asymptotic expansions are straightforward and can be obtained from standard handbooks on bessel functions . we therefore arrive at our final result : the integration necessary for evaluating the near threshold expansion of the water melon diagrams reduces to integrals of the type of euler s gamma function , i.e.integrals containing exponentials and powers . indeed , the result of the expansion in eq . ( [ pidas ] ) is an exponential function @xmath106 times a power series in @xmath107 , namely @xmath108 where @xmath109 has already been defined in eq . ( [ defa ] ) and the coefficients @xmath110 are simple functions of the momentum @xmath44 and the masses @xmath7 . the expression in eq . ( [ serx ] ) can be integrated analytically using @xmath111 the result is @xmath112 this expression is our final representation for the part of the polarization function of a water melon diagram necessary for the calculation of the spectral density near the production threshold . it is also one of the main results of our paper . next we discuss the general structure of the expression in eq . ( [ serint ] ) in detail . in the case where @xmath109 takes integer values , these coefficients result in @xmath113-divergences for small values of @xmath97 . the powers of @xmath48 in eq . ( [ serint ] ) have to be expanded to first order in @xmath97 and give @xmath114 because of @xmath115 @xmath116 in eq . ( [ serint ] ) contributes to the spectral density . for half - integer values of @xmath109 the power of @xmath48 itself has a cut even for @xmath117 . the discontinuity is then given by @xmath118 our method to construct a threshold expansion thus simply reduces to the analytical calculation of the integral in eq . ( [ pidas ] ) which can be done for arbitrary dimension and an arbitrary number of lines with different masses . in the next subsections we use our technique to work out some specific examples which demonstrate both the simplicity and efficiency of our method . the polarization function represented by the sunset diagram with three propagators with equal masses @xmath11 in @xmath119 space - time dimensions is given by @xmath120 the exact spectral density is given by the integral representation in eq . ( [ exactrho ] ) which for this particular case reads @xmath121 in order to obtain a threshold expansion of the spectral density in eq . ( [ exactrho4 ] ) we use eq . ( [ serint ] ) to calculate the expansion of the appropriate part of the polarization function . to illustrate the procedure we present the explicit shape of the integrand in eq . ( [ pidas ] ) which is given by an asymptotic expansion at large @xmath30 , @xmath122 from eq . ( [ ex ] ) we can easily read off the coefficients @xmath110 that enter the expansion in eq . ( [ serx ] ) . the spectral density is obtained by performing the term - by - term integration of the series in eq . ( [ ex ] ) and by evaluating the discontinuity across the cut along the positive energy axis @xmath81 . the result reads @xmath123 where the notation @xmath124 , @xmath125 is used . the simplicity of the derivation is striking . by no cost it can be generalized to any number of lines , arbitrary masses , and any space - time dimension . the standard equal mass sunset is chosen for the definiteness only . it also allows us to compare our results with results available in the literature . ( [ pidas430 ] ) reproduces the expansion coefficients @xmath126 obtained in ref . @xcite ( the fourth column in table  1 of ref . @xcite ) by a direct integration in momentum space within the technique of region separation  @xcite . the case of the equal mass standard sunset diagram is the simplest one . there exists an analytical expression for the spectral density of the sunset diagram with three equal mass propagators in @xmath119 space - time dimensions in terms of elliptic integrals  @xcite ( see also ref  @xcite ) . this expression can be used for a comparison with our exact result in eq . ( [ exactrho4 ] ) or with the expansion in eq . ( [ pidas430 ] ) . however , we only present the result for @xmath127 in order to keep the resulting expressions in a reasonably short form ( cf.ref . @xcite ) . in @xmath127 space - time dimensions the spectral density for a sunset diagram with equal masses @xmath11 can be readily obtained . we just use the exact expression for the spectral density in the convolution representation  @xcite and proceed towards @xmath128 equal masses . the convolution function for two spectral densities in @xmath127 dimensional space - time ( @xmath129 ) reads @xmath130 the two spectral densities one has to convolute are the spectral density of a correlator with two equal masses and the spectral density of a single massive line . while the latter is given by @xmath131 , the former can be obtained from eq . ( [ convolut ] ) by inserting @xmath132 , @xmath133 so the convolution leads to @xmath134 now we use the relation ( cf .  ref . @xcite ) @xmath135 @xmath136 with @xmath137 and the definition of the complete elliptic integral of the first kind @xmath138 ( remark the difference in the definition ) for @xmath139 , @xmath140 , @xmath141 , and @xmath142 to perform the integration in eq . ( [ ellint ] ) . we obtain @xmath143 and finally end up with @xmath144 therefore the spectral density in terms of the energy @xmath50 reads ( see e.g.ref . @xcite ) @xmath145 by expanding the elliptic integral in terms of the threshold parameter @xmath50 one reproduces the threshold expansion in eq . ( [ pidas430 ] ) . the result for @xmath119 space - time dimensions is expressible by the elliptic integrals with some rational functions as factors that makes the result a bit longer . note that the representation in eq . ( [ ell ] ) is understood to be an analytical expression for the spectral density . however , it is a matter of taste whether the representation through the elliptic integrals as in eq .  ( [ ell ] ) is considered simpler ( or in a more analytical form ) than the integral representation in eq . ( [ exactrho ] ) . the only objection against the latter which one can find in the literature is that the bessel functions are complicated ( see e.g.  ref . but after more than a century of intensive investigation they are well - known and no more complicated than the square root of the fourth order polynomial which is used in eq . ( [ elldef ] ) to define the elliptic integral . the water melon diagrams with four or more propagators can not be easily done by using the momentum space technique because it requires the multiloop integration of entangled momenta . within the configuration space technique the generalization to any number of lines ( or loops ) is immediate by no effort . consider first a three - loop case of water melon diagrams ( also called banana diagrams  @xcite or basketball diagrams  @xcite ) . the polarization function of the equal mass water melon diagram with four propagators in @xmath119 space - time is given by @xmath146 the exact spectral density of this diagram can be obtained from eq . ( [ exactrho4 ] ) while the near threshold expansion can be found using eq . ( [ serint ] ) . we construct the expansion of the spectral density near threshold explicitly and compare it with the exact result . the expansion of the integrand ( cf . ( [ pidas ] ) ) reads @xmath147 after the integration and the calculation of the discontinuity we obtain the expansion of the spectral density in the form @xmath148 where @xmath124 and @xmath149 is the threshold value . one sees the difference with the previous three - line case . in eq . ( [ pidas440 ] ) the cut represents the square root branch while in the three - line case it was a logarithmic cut . one can easily figure out the reason for this by looking at the asymptotic structure of the integrand . for even number of lines ( i.e.  odd number of loops ) it is a square root branch , while for an odd number of lines ( even number of loops ) it is a logarithmic branch . this is true in even space - time dimensions . in the general case the structure of the cut depends on the dimensionality of the space - time as well . the general formula reads @xmath150 for @xmath119 space - time dimension ( i.e.  @xmath151 ) we can verify the result of ref . @xcite ( cf . ( [ genthr ] ) ) , @xmath152    numerically eq .  ( [ pidas440 ] ) reads @xmath153 where we have written down the coefficients up to three decimal places . it is difficult to say anything definite about the convergence of this series . by construction it is an asymptotic series . however , we stress that the practical ( or explicit ) convergence can always be checked by comparing series expansions like the one shown in eq . ( [ sernum ] ) with the exact spectral density given in eq . ( [ exactrho4 ] ) by numerical integration . we conclude this part of the paper by the statement that the spectral density of the water melon diagram is most efficiently calculable within the configuration space technique . whether it is the exact result or the expansion , the configuration space technique can readily deliver the desired result . the exact formula in eq . ( [ exactrho4 ] ) as well as the threshold expansion obtained from it can be used to calculate the spectral density for an arbitrarily large number of internal lines . we stress that the case of different masses does not lead to any complications within the configuration space technique : the exact formula in eq . ( [ exactrho ] ) and/or the near threshold expansion work equally well for any arrangement of masses . we do not present plots for general cases of different masses because they are not very illustrative , showing only the common threshold . however , there is some interesting kinematic regime for different masses which is important for applications and which , to our best knowledge , have not been touched earlier . an analytical solution for the expansion of the spectral density in this regime is given in the next section . the threshold expansion for equal ( or close ) masses breaks down for @xmath155 . the example is shown in fig .  [ fig2 ] for the @xmath119 proper sunset . however , if the masses are not equal , the region of the break - down of the expansion is determined by the mass with the smallest numerical value . the simplest example where one can see this phenomenon is the analytical expression for the spectral density of the simple loop ( degenerate water melon diagram ) with two different masses @xmath156 and @xmath157 . in @xmath119 space - time dimensions ( see e.g.  ref . @xcite ) one has @xmath158 where @xmath159 . the threshold expansion is obtained by expanding the right hand side of eq . ( [ rho42 m ] ) in @xmath50 for small values of @xmath50 . if @xmath157 is much smaller than @xmath156 , the expansion breaks down at @xmath160 . the break - down of the series expansion can also be observed in more general cases . if one of the masses ( which we call @xmath161 ) is much smaller than the other masses , the threshold expansion is only valid in a very limited region @xmath162 .    to generalize the expansion and extend it to the region of @xmath163 one has to treat the smallest mass exactly . in this case one can use a method which we call the resummation of the smallest mass contributions . below we describe the resummation technique . we start with the representation @xmath164 which is the part of the polarization function contributing to the spectral density . the integrand in eq . ( [ pipas ] ) has the form @xmath165 where the asymptotic expansions are substituted for all the propagators except for the one with the small mass @xmath161 . this is indicated by the index `` _ pas _ '' in eq . ( [ pipas ] ) which stands for `` partial asymptotic '' . the main technical observation leading to the generalization of the expansion method is that @xmath166 is still analytically computable in a closed form . indeed , the genuine integral to compute has the form @xmath167 where @xmath168 and @xmath169 . the integral @xmath166 in eq . ( [ pipas ] ) is thus expressible in terms of hypergeometric functions  @xcite . for constructing the spectral density , being our main concern as mentioned before , one has to find the discontinuity of the right hand side of eq . ( [ intm0 ] ) . there are several ways to do this . we proceed by applying the discontinuity operation to the integrand of the integral representation of the hypergeometric function . the resulting integrals are calculated again in terms of hypergeometric functions . indeed , @xmath170 where @xmath171 . the final expression in eq . ( [ genuine ] ) completely solves the problem of the generalization of the near threshold expansion technique . for integer values of @xmath172 there are no singular gamma functions ( with negative integer argument ) . therefore we can lift up the regularization and set @xmath117 when using this expression . we thus have found a direct transition from the polarization function as expressed through the integral to the spectral density in terms of one hypergeometric function for each genuine integral . there is no need to use the recurrence relations available for hypergeometric functions .    in the following subsections we give explicit examples for @xmath119 and compare with the exact result in eq . ( [ exactrho4 ] ) and the pure expansion near the threshold . in the following the standard threshold expansion without resummation is called the pure threshold expansion . we start with a ( over)simplified example of the two - line diagram with masses @xmath11 and @xmath173 in four space - time dimensions . in this example the expansion of the spectral density and its generalized expansion can be readily compared analytically with the exact result in eq . ( [ rho42 m ] ) . this is the feature that justifies our discussion in this section . the results for the spectral density of this diagram are shown in fig . the solid curve displays the exact result obtained by using eq . ( [ exactrho4 ] ) ( which reproduces the analytical expression in eq . ( [ rho42 m ] ) ) . we compare this result with the two expansions . the pure expansion of the spectral density near threshold ( the second order asymptotic expansion should suffice to show the general features in a short and concise form ) is given by @xmath174 where @xmath175 . as mentioned above , this series breaks down for @xmath176 ( see eq . ( [ rho42 m ] ) ) . if we look at the dotted curves in fig . [ fig3 ] this becomes obvious . here we have plotted the series expansions up to the fourth order with the mass arrangement @xmath177 . the dashed lines represent the resummation of the smallest mass contributions . the analytical expression for the spectral density of the polarization function in eq .  ( [ pipas ] ) for the generalized asymptotic expansion based on eq . ( [ genuine ] ) is given by @xmath178 we have set the regularization parameter @xmath117 because the spectral density is finite . with @xmath117 the resulting expressions for the hypergeometric functions in eq . ( [ genuine ] ) simplify . the first term in the curly brackets of eq . ( [ pipas42 t ] ) is obviously equal to @xmath179 in this limit because the first parameter of the hypergeometric function vanishes for @xmath117 . however , we keep eq . ( [ pipas42 t ] ) in its given form to show the structure of the contributions . the generalized threshold expansion has the form @xmath180 where the functions @xmath181 represent effects of the resummation of the smallest mass and are not polynomials in the threshold parameter parameter @xmath50 . in the simple two - line case the hypergeometric functions reduce to elementary functions . for instance , @xmath182 higher order contributions are given by hypergeometric functions with larger numerical values of the parameters . they can be simplified by using gaussian recurrence relations for hypergeometric functions ( see e.g.ref . @xcite ) .    the convergence of the expansion in eq . ( [ pipas42 t ] ) breaks down only at @xmath183 . the resummation leads to an essential improvement of the convergence in comparison with the pure threshold expansion . in fig . [ fig4 ] we show the same curves divided by the leading order term . this representation is more convenient for the diagrams which we will discuss in following subsections . note that eq . ( [ resumsm ] ) does not lead to the exact function in eq . ( [ rho42 m ] ) because terms of order @xmath184 stemming from the difference part @xmath95 of the correlator are missing . it simply corrects the behaviour of the coefficient functions by the small mass contributions .      here we analyze the sunset diagram with two equal masses @xmath11 and a third mass @xmath173 ( @xmath177 ) . the exact result obtained by using eq . ( [ exactrho4 ] ) and normalized to the leading order term is shown in fig . [ fig5 ] as the solid curve . the pure expansion near the threshold reads @xmath185 it is shown by the dotted curves in fig . [ fig5 ] . in case of the resummation of the smallest mass contributions we obtain hypergeometric functions which do not obviously reduce to elementary functions in this case . the result for the spectral density within the asymptotic expansion up to the second order in eq . ( [ pipas ] ) is given by @xmath186 we see that the dashed curves in fig . [ fig5 ] that represent the result in eq . ( [ pipas43 t ] ) approximate the exact curve much better than the dotted curves . with this example we conclude our consideration of the strongly asymmetric case and at the same time show the way to the multi - line water melon diagrams which can be treated in an analogous manner . the result for the exact expression obtained by using eq . ( [ exactrho4 ] ) is shown in fig . [ fig6 ] as a solid line , normalized to the leading order term . the dotted lines represent the results for the pure expansion near threshold which is given by @xmath187 the asymptotic expansion to the second order in eq . ( [ pipas ] ) gives @xmath188 in fig . [ fig6 ] one can see how the expansion improves if the resummation of the smallest mass contributions ( displayed as dashed lines ) is performed . the result of this section is quite general and applicable to all cases of one small mass . for some particular arrangement of masses one can obtain even simpler expressions as discussed in the next section . in this section we obtain a result for the resummation of the smallest mass effects along a different route , namely , via the convolution of spectral densities . however , this method works in a narrower kinematic region than the method described in the previous section . in @xmath119 space - time dimensions , the convolution weight is given by @xmath189 the upper limit of the integration is determined by the requirement of positivity of the the square root argument . the zeros of the square root with respect to @xmath190 are given by @xmath191 , and the demand @xmath192 together with @xmath193 leads to @xmath194 or @xmath195 . the physical region is the latter one . with @xmath196 for the spectral density of the single small mass line we obtain @xmath197 where the low limit of integration is @xmath198 . we insert @xmath199 and @xmath200 and obtain @xmath201 where @xmath202 . for this function we use the threshold expansion in @xmath203 as expansion parameter . for small @xmath204 the threshold expansion inserted for @xmath205 is valid because @xmath204 implies @xmath206 . the described procedure can be extended to the case of a very light sub - block of the diagram , e.g.  a light fish diagram . in this case we have to replace @xmath207 by the spectral density of the light sub - diagram which is well - known . the analytic structure of water melon diagrams is completely fixed by the dispersion representation . therefore we have focussed on the computation of the spectral density as the basic quantity important both for applications and the theoretical investigation of the diagram . however , with an analytical expression for the spectral density @xmath19 at hand we can readily reconstruct the non - analytic piece of the polarization function in momentum space by using the dispersion relation @xmath208 we rewrite this equation in terms of threshold parameters according to @xmath209 and @xmath199 and obtain @xmath210 uv singularities can be removed by subtraction or by dimensional regularization . we again use the unorthodox dimensional regularization prescription . for a general form of the threshold expansion @xmath211 we have to calculate integrals of the form @xmath212 only the powers @xmath213 contribute to the singular part of the polarization function . expressions like the one presented in eq . ( [ restore ] ) then allow one to restore that part of the polarization function @xmath65 which has singularities near the threshold . we have discussed the configuration space technique for the calculation of @xmath3-line two - point diagrams termed water melon diagrams . this technique allows one to calculate the spectral density for arbitrary space - time dimensions and any number of internal lines with arbitrary mass values . within this technique one can use either an exact representation as one - fold integral or an expansion near the production threshold . we have developed an efficient method for constructing the near threshold expansion of water melon diagrams that uses only asymptotic expansions of bessel functions which are well - known and simple functions . we have considered a strongly asymmetric mass arrangement where one mass is much smaller than the others . in this case the `` pure '' threshold expansion which is done in terms of the threshold parameter @xmath50 breaks down at the energies in the vicinity of the smallest mass . in order to extend the approximation to higher energies we have introduced a generalized threshold expansion which exactly resums all contributions of the smallest mass . we have presented closed expressions for the generalized expansion in several particular cases and demonstrated the improvement of the convergence gained by the resummation of the smallest mass contributions . the particular kinematic regime of this case could be treated analytically because it reduced to the evaluation of the one - fold integral in terms of hypergeometric functions . the discontinuity across the physical cut for the generalized expansion has been found in terms of hypergeometric functions as well .    to conclude , we stress that the configuration space technique is a powerful and convenient tool for investigating different properties of water melon diagrams . the practical convenience of our method is demonstrated in fig . [ fig7 ] where we have plotted the spectral density for a four - line water melon diagram in @xmath214 , @xmath215 , and @xmath216 space - time dimensions .    * acknowledgments * + we thank andrei davydychev and jrgen krner for interesting discussions and careful reading the manuscript . a.a . pivovarov acknowledges a valuable discussion with v.a . smirnov on the region separation technique in threshold expansions . the work of a.a . pivovarov is supported in part by the volkswagen foundation under contract no . i/73611 and by the russian fund for basic research under contract 99 - 01 - 00091 . s.  groote gratefully acknowledges a grant given by the graduiertenkolleg `` eichtheorien  experimentelle tests und theoretische grundlagen '' .    99 particle data group , `` review of particle properties '' , eur . j. * c3 * ( 1998 ) 1 k.g . chetyrkin , j.h . khn and a.  kwiatkowski , phys . * 277 * ( 1996 ) 189 k.g . chetyrkin , j.h . khn and m.  steinhauser , nucl . * b505 * ( 1997 ) 213 d.j . broadhurst , eur . j.  * c8 * ( 1999 ) 311 o.m . ogreid , p.  osland , `` more series related to the euler series '' , + report no .  bergen-1999 - 04 , hep - th/9904206 r.  coquereaux , phys .  rev . * d23 * ( 1981 ) 2276 ; + v.i . zakharov , nucl . * b385 * ( 1992 ) 452 ; + d.j .  broadhurst and a.g . grozin , phys . rev .  * d52 * ( 1995 ) 4082 s.  groote , j.g . krner and a.a . pivovarov , nucl . * b542 * ( 1999 ) 515 ; eur . j.  * c11 * ( 1999 ) 279 ; phys . * 443 b * ( 1998 ) 269 f.a . berends , m.  buza , m.  bhm and r.  scharf , z.  phys . * c63 * ( 1994 ) 227 f.a . berends , a.i . davydychev and n.i . ussyukina , phys . lett .  * 426 b * ( 1998 ) 95 s.  coleman and e.  weinberg , phys . * d7 * ( 1973 ) 1888 ; + r.  jackiw , phys .  rev . * d9 * ( 1974 ) 1686 ; + r.  jackiw and s.templeton , phys . * d23 * ( 1981 ) 2291 a.k . rajantie , nucl . * b480 * ( 1996 ) 729 j.m . chung and b.k . chung , `` renormalization group improvement of the effective potential in massive @xmath217 theory : nnnlo logarithm resummation '' , report no . mit - ctp-2929 , hep - th/9911196 j.o . andersen , e.  braaten and m.  strickland , `` the massive thermal basketball diagram '' , hep - ph/0002048 j.m . chung and b.k . chung , phys . * d59 * ( 1999 ) 105014 d.j . gross , r.d . pisarski and l.g . yaffe , rev . phys .  * 53 * ( 1981 ) 43 ; + t.  appelquist and r.d . pisarski , phys . * d23 * ( 1981 ) 2305 b.l . ioffe , nucl . * b188 * ( 1981 ) 317 n.v . krasnikov , a.a .  pivovarov and a.n . tavkhelidze , + z.  phys . * c19 * ( 1983 ) 301 ; jetp lett . * 36 * ( 1982 ) 333 a.a . ovchinnikov , a.a . pivovarov and l.r . surguladze , + yad .  fiz .  * 48 * ( 1988 ) 562 ; int . j.  mod . * a6 * ( 1991 ) 2025 s.  groote , j.  g.  krner and a.  a.  pivovarov , phys . * d61 * ( 2000 ) 071501 g.  t hooft , nucl . phys .  * b72 * ( 1974 ) 461 e.  witten , nucl . phys .  * b160 * ( 1979 ) 57 s.a . larin _ et al . _ , yad .  fiz . * 44 * ( 1986 ) 1066 s.a . kulagin and a.a . pivovarov , yad . * 45 * ( 1987 ) 952 p.  post and k.  schilcher , phys . * 79 * ( 1997 ) 4088 j.  gasser and m.e . sainio , eur . j.  * c6 * ( 1999 ) 297 s.  narison and a.a . pivovarov , phys . * 327 b * ( 1994 ) 341 i.i . balitsky , d.i . dyakonov and a.v . yung , phys . * 112 b * ( 1982 ) 71 d.j . broadhurst , z.  phys . * c54 * ( 1992 ) 599 p.a . baikov , `` advanced method of solving recurrence relations for multiloop feynman integrals '' , talk presented at the @xmath218 international workshop on new computing techniques in physics research ( aihenp 99 ) , heraklion , greece , 1216 april , 1999 n.e . ligterink , `` solving multiloop feynman diagrams using light front coordinates '' , report no .  ect-99 - 016 , hep - ph/9911411 a.i . davydychev and v.a . smirnov , nucl . * b554 * ( 1999 ) 391 b.  kastening , h.  kleinert , `` efficient algorithm for perturbative calculation of multiloop feynman integrals '' , quant - ph/9909017 g.n . watson , `` theory of bessel functions '' , cambridge , 1944 e.  mendels , nuovo  cim .  * 45 a * ( 1978 ) 87 c.s . meijer , proc . amsterdam  akad . ( 1940 ) 599 ; 702 a.  erdelyi ( ed . ) , `` tables of integral transformations '' , + volume  2 , bateman manuscript project , 1954 a.h . hoang _ et al . _ , `` top - antitop pair production close to threshold : synopsis of recent nnlo results , '' report no . slac - pub-8369 , hep - ph/0001286 ; + a.  czarnecky and k.  melnikov , phys . * 80 * ( 1998 ) 2531 ; + m.  beneke , a.  signer and v.a . smirnov , phys . * 80 * ( 1998 ) 2535 ; + k.  melnikov and a.  yelkhovsky , nucl . * b528 * ( 1998 ) 59 ; + a.a . penin and a.a . pivovarov , hep - ph/9904278 ; nucl . * b549 * ( 1999 ) 217 ; nucl . * b550 * ( 1999 ) 375 ; phys . lett .  * 443 b * ( 1998 ) 264 ; + o.i . yakovlev , phys . lett .  * 457 b * ( 1999 ) 170 ; + a.h .  hoang and t.  teubner , * d60 * ( 1999 ) 114027 ; + t.  nagano , a.  ota and y.  sumino , phys . * d60 * ( 1999 ) 114014 n.n . bogoliubov and d.v . shirkov , `` quantum fields '' , benjamin , 1983 a.a . pivovarov , phys . * 236 b * ( 1990 ) 214 ; phys . lett .  * 263 b * ( 1991 ) 282 v.a . smirnov , phys . * b404 * ( 1997 ) 101 ; + m.  beneke and v.a . smirnov , nucl . * b522 * ( 1998 ) 321 ; + v.a . smirnov and e.r . rakhmetov , teor .  mat . 120 * ( 1999 ) 64 a.p . prudnikov , yu.a . brychkov and o.i . marichev , + `` integrals and series '' , vol .  2 , gordon and breach , new york , 1990 b.almgren , arkiv fr fysik , * 38 * ( 1967 ) 161 ; + s.bauberger , f.a . berends , m.  bhm and m.  buza , nucl . * b434 * ( 1995 ) 383 s.  groote , j.  g.  krner and a.  a.  pivovarov , phys . * d60 * ( 1999 ) 061701 a.a . pivovarov and v.f . tokarev , yad .  fiz . * 41 * ( 1985 ) 524 i.s . gradshteyn and i.m . ryzhik , + `` tables of integrals , series , and products '' , academic press , 1994 m.  abramowitz , i.a . stegun ( eds . ) , `` handbook of mathematical functions '' , dover publ . , new york , 9th printing , 1970
Abstract:

Answer: the near threshold expansion of generalized sunset - type ( water melon ) diagrams with arbitrary masses is constructed by using a configuration space technique . we present analytical expressions for the expansion of the spectral density near threshold and compare it with the exact expression obtained earlier using the method of the hankel transform . we formulate a generalized threshold expansion with partial resummation of the small mass corrections for the strongly asymmetric case where one particle in the intermediate state is much lighter than the others . mz - th/00 - 07 + hep - ph/0003115 + march 2000     * threshold expansion of feynman diagrams + within a configuration space technique * + s.  groote@xmath0 and a.a . pivovarov@xmath1 + @xmath0 institut fr physik , johannes - gutenberg - universitt , + staudinger weg 7 , d-55099 mainz , germany + @xmath2 institute for nuclear research of the + russian academy of sciences , moscow 117312