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: the study of the cosmic microwave background ( cmb ) radiation has had a long history . three aspects of the cmb might be considered : its existence , its spectrum , and its anisotropies . by firmly establishing that the universe expanded from an initially hot , dense state , the _ existence _ of the cmb underpins our entire cosmological framework . it has been recognized from the beginning as one of the pillars of the hot big bang cosmologies . the spectrum of the cmb is the most precise blackbody spectrum in nature , from which many inferences can be made . although this discovery is less than a decade old , its impact on models of the early universe been discussed extensively elsewhere , ( e.g.  @xcite ) . in this paper we would like to consider the _ anisotropies _ in the cmb , the small fluctuations imprinted on the sky by the progenitors of the large - scale structure seen in the distribution of galaxies today .    in the roughly seven years since the _ cobe _  dmr team announced the first detection of anisotropies in the cmb ( @xcite ) , more than a dozen groups have reported detections , covering the full range of frequencies and a wide range of angular scales ( see @xcite ) . due in large part to a dramatic increase in detector sensitivity , mapping the cmb anisotropy has become almost routine . our confidence in the results has grown as multiple observations by the same teams over a period of years , and then later by different experiments at different frequencies and sites , reproduced the same features on the sky and confirmed their black body nature .    over the same period much progress has been made in data analysis techniques and in the theoretical interpretation of cmb data . better physical understanding of the anisotropy generation has lead to faster algorithms for its computation ( e.g.  @xcite ) applicable to an impressively wide range of theories . the high precision calculations and accurate measurements of the anisotropy have spawned numerous ideas in data analysis , with a full likelihood analysis of mega - pixel cmb maps now within reach ( @xcite ) . however , since much of the progress has been incremental , it is not always obvious just how far we have really advanced . it therefore seems appropriate to take stock and ask the question :    _ what has the cmb ever done for us ? _ from a broad perspective , the main impact of cmb anisotropies has been to shrink substantially the range of cosmological models under active discussion . this is not always easy to see , since the _ number _ of models proposed at any time seems to be determined more by the number of theorists working in the field than by any constraints provided by the data . moreover , it sometimes seems that no class of model has been ruled out . however , looking back a decade in the literature makes it clear that this is not true .    even before _ cobe _ , the high level of isotropy of the cmb was perhaps the best possible evidence that the large - scale properties of the universe were well described by the friedman - robertson - walker metric . the assumption of homogeneity and isotropy , initially made for purely aesthetic reasons , turned out to be an extremely good approximation to the real universe . as the limits on anisotropy became stronger and stronger , the number of models based on anything but the frw metric became fewer and fewer . currently popular models assume that the matter in the frw universe is composed _ mostly _ of cold dark matter ( cdm ) , with smaller admixtures of baryons and perhaps massive neutrinos , plus curvature and/or vacuum components . for these cdm - inspired models , the cmb data have been instrumental in narrowing the range of possibilities , and most popular flavours of cdm now give remarkably similar predictions . models dominated by hot dark matter , already in trouble before _ cobe _ , are no longer discussed . two other classes of models , namely defects and isocurvature models , have not been ruled out definitively , but they are now very much on the defensive against the weight of data . explosion models ( @xcite ) , super - conducting cosmic string models ( @xcite ) , and late - time phase transition models ( @xcite ) have essentially vanished . figure  [ fig : current ] shows the current state of cmb measurements . included are all detections we are aware of that have been published or submitted for publication in 1998 . the results have been averaged in 12 bins , equally spaced in @xmath0 for clarity , and we have omitted the upper limits on smaller angular scales , most of which are off the right of the plot with our chosen @xmath1-axis range . this figure is meant to be indicative only . more statistically rigorous approaches exist for combining data sets ( e.g.  @xcite ) , and such methods should certainly be used for determining precise constraints on models . however , fig . [ fig : current ] gives approximately the correct visual impression for the combined constraining power of today s data . below we list two sets of statements that we believe are supported by the data : the first set contains ` fundamental truths ' about the universe ; and the second contains statements that will be fundamental truths if confirmed , but that for the present must be regarded more tentatively .    here is the ` a ' list :    *  gravitational instability in a dark matter dominated universe grew today s structure *  the universe ( re)combined *  there is an excess of temperature fluctuations at roughly the predicted angular scale *  the polarization of the cmb anisotropy is small    and the ` b ' list :    *  something like inflation produced adiabatic fluctuations *  the large - scale structure of space - time appears to be simple *  the gravity wave contribution to the anisotropy is not large *  there are constraints on non - standard physics at @xmath2    we now discuss these in turn , distinguishing between those demonstrated by _ cobe _  alone , and those demonstrated by the measurements at smaller angular scales that have been made since _ cobe_.    * 2.a.1 gravitational instability *    perhaps the most useful result of the _ cobe _  anisotropy data is the normalization of models of structure formation at large - angles , where the fluctuations in the matter and photons are expected to be in the linear regime . in today s favoured models of structure formation these large - angle anisotropies directly measure the amplitude of the gravitational potential on very large scales , allowing a theoretically clean and precise normalization of the matter power spectrum . this normalization , it turns out ( e.g.  @xcite ) , is in the right ball - park to explain the amplitude of galaxy clustering ( and with a little tuning of this or that parameter it is easy to get complete consistency ) . this is a vindication of our ideas that galaxies grew gradually under the action of gravitational instability .    before the _ cobe _ anisotropy was announced it was often claimed ( e.g.  @xcite ) that extra physics would be needed if the results turned out to yield yet more upper limits ; right up to the dmr announcement it was also commonly perceived that inflationary adiabatic models had difficulty having a high enough amplitude to form structure without violating cmb limits ( e.g.  @xcite ) . the fact that the anisotropies were measured at the levels predicted , in models with cold dark matter and adiabatic fluctuations , showed that there is no need to invoke extra magical processes to form structure by the present day . however , since the photons prevent baryonic matter from collapsing before recombination , we infer that the gravitational potentials had to be dominated by matter which was not prevented from collapsing by photon pressure , i.e.  matter that was not coupled to photons and was ` dark ' . the realization , from studies of the galaxy distribution in the local universe , that matter formed ` bottom up ' rather than ` top down ' constrains the velocity dispersion of the dominant dark matter component to be extremely small  the dark matter must be mostly _ cold_.    * 2.a.2 recombination *    here we are moving beyond simply an interpretation of the _ cobe _  data , and looking at the large number of detections of anisotropy at degree and sub - degree scales ( see fig .  1 ) . early reionization of the universe gives increased optical depth to thomson scattering from the present back to the epoch of reionization . the extreme case is a universe which did not ( re)combine at all and remained ionized for all time . multiple scattering erases existing anisotropies on scales smaller the horizon . thus reionization leads to damping of primordial anisotropies on small scales ( @xcite ) . the presence of fluctuations at @xmath3 is clear evidence that the universe was not reionized at a very early epoch . we can be confident that the universe recombined at @xmath4 , then remained largely neutral until some redshift @xmath5 , after which it was largely ionized ( as implied by the absence of gunn - peterson absorption in the spectra of high-@xmath6 quasars ) . the precise value of @xmath5 derived from fits to the data depends on the cosmological model , but is typically @xmath7 ( @xcite ) .    * 2.a.3 degree scale power *    we believe that fig . 1 shows a peak in power in the anisotropies at scales around a degree . the _ precise _ position of this peak , how high it might be , and whether it contains any substructure , are not so clear ( see e.g.  @xcite ) . however , it is striking that this feature is in the general location of the main acoustic peak predicted by currently favoured models , based on the angular size of the horizon at last scattering . it is worth stressing that this prediction was made more than a decade before the experiments were performed ( see for example @xcite ) . we expect the location of the peak to be determined definitively quite soon , by upcoming ground based and balloon experiments , interferometers and _ map _ , leading to very strong observational constraints on the angular diameter distance back to last scattering ( @xmath8 ) .    * 2.a.4 polarization *    it is a fundamental prediction of the gravitational instability paradigm that the cmb anisotropy is linearly polarized . in inflationary cdm - like models the level of polarization is a few percent of the anisotropy , and thus extremely small in absolute terms . there are already many limits on the polarization of cmb anisotropy ( see @xcite for a list ) , however they are all nearly an order of magnitude larger than the theoretical predictions . the fact that the cmb is not ` very ' polarized tells us important information about the conditions at the last scattering epoch . that the cmb is not very circularly polarized , for example , indicates that there were no large magnetic fields present at last scattering ( see also  2.b.4 ) , although we are only aware of very stringent upper limits at the smallest angular scales ( @xcite ) . * 2.b.1 inflation *    we put this item at the very top of our ` b'-list since we feel the weight of evidence is becoming very strong for something akin to inflation ( for a discussion of whether inflation is really a testable theory , see @xcite ) . to avoid semantic arguments , it is important at the outset to be clear about the meaning of ` inflation ' . here we refer to a period of accelerated expansion in the early universe . this is the only known mechanism for making an isotropic and homogeneous universe , and at the same time generates apparently acausal adiabatic fluctuations , i.e.  fluctuations in spatial curvature on scales larger than the hubble - length at a particular epoch . we do not intend ` inflation ' to carry the additional baggage of an inflaton field with a well - defined potential , connected with particle physics , etc . although ultimately we would all like to see the mechanism of inflation find a realization in a well motivated theory of fundamental physics . the amplitude and power spectrum of cmb anisotropies from degree - scales up to the largest scales probed by _ cobe _  seem to indicate that super - horizon size adiabatic fluctuations exist . our first hint comes from the normalization of the large - scale anisotropies relative to the matter ( see e.g.  discussion in @xcite ) . on dimensional grounds we expect that the amplitude of the temperature fluctuations be @xmath9 where @xmath10 is the large - scale gravitational potential . in adiabatic models a cancellation ( @xcite ) between intrinsic anisotropies and gravitational redshifts means that the coefficient is reduced to @xmath11 , i.e.  @xmath12 ( @xcite ) . in the simplest isocurvature models the coefficient is @xmath13 . since , as we mentioned before , our currently popular theories ` work ' , there is little room to absorb a factor of 6 in relative normalization . of course this alone is not proof of adiabatic fluctuations .    our next piece of observational evidence is the angular scale of the ` peak ' in power . the structure of the peaks ( locations , separations , relative heights ) is a strong discriminator between adiabatic and isocurvature models ( @xcite ) . in almost all isocurvature models the peak is shifted to smaller angular scales . since we observe excess power at about the right place for adiabatic fluctuations in a flat universe , there is little room for either spatial curvature or isocurvature fluctuations ( and the combination is particularly disfavoured ! ) . since the current evidence for a peak , in contrast to a rise , is modest we have put this in our ` b'-list . the observational situation is likely to change rapidly . in the future we can hope that detection of polarization on degree scales will finally pin down the fluctuation type beyond any argument ( @xcite ) , but this is a difficult measurement due to the low levels of signal . thus there is reasonable evidence for adiabatic fluctuations in a spatially flat universe . the latter has long been hailed as a ` prediction ' of inflation . the former is also tantamount to a ` proof ' of inflation , in the sense that the only _ causal _ means for generating nearly scale - invariant adiabatic fluctuations is a period when @xmath14 in the early universe ( see e.g.  @xcite ) . of course this condition is neither entirely necessary nor sufficient . on the sufficiency side , it is no doubt possible to imagine inflationary models which have fluctuations of an entirely different character , but it would seem pathological to deliberately avoid explaining density perturbations . and on the necessary side , one could in principle imagine some early universe physics which somehow mimics the effects of inflation by producing super - horizon adiabatic modes , and yet is not inflation . we would argue that this is a purely semantic distinction : if it looks like inflation and smells like inflation , then let s call it inflation while leaving open the possibility that current inflationary ideas may one day be shown to be part of some better paradigm . in the same vein it may also be argued that some planck - era physics somehow generates apparently acausal modes . again we would say that is either isomorphic with inflation , or simply an attempt to push the question of initial conditions into the realm of metaphysics .    * 2.b.2 space - time structure *    we have already mentioned that the extreme isotropy of the cmb is a strong indication that the frw metric is an excellent approximation to the large - scale properties of space - time . strong quantitative limits on the rotation and shear of space - time for specific bianchi models have been obtained from the _ cobe _  data ( @xcite ) . and limits on the geometry for general models can be placed at the @xmath15 level ( @xcite ) . cmb anisotropies probe the universe on the largest accessible scales , and so they also constrain things like the large - scale topology . there are quite stringent constraints in the simplest background models ( @xcite ) . however , in principle there may yet be observational consequences for compact topologies , in an open universe in particular ( @xcite ) . exactly how stringent the current constraints are , for general classes of cosmology on the largest scales , is still a matter of debate . nevertheless , we probably know at this point that the universe is nt _ very _ strange on gpc scales , quite an advance over our previous ignorance . * 2.b.3 gravity waves *    if whatever produces the initial density perturbations does nt discriminate on the basis of perturbation type we would expect that scalar , vector and tensor fluctuations would be produced at early times in roughly equal amounts . the vector modes , representing fluid vorticity , decay with time and so would not be present after a few expansion times . thus we would expect today to see only scalar ( density ) perturbations and tensor ( gravity wave ) perturbations . both of these types of perturbation would give rise to large - angle anisotropies , though only the former will seed large - scale structure . due to the aforementioned close consistency between the amplitude of the clustering on galaxy scales and the anisotropy seen by _ cobe _ there is a limit to how much the gravity wave signal can contribute to _ cobe_. roughly speaking , the tensor to scalar ratio @xmath16 ( see @xcite ) . if the tensor perturbations are not too different from scale - invariant this means that the possibility of seeing primordial gravity waves with detectors such as ligo or lisa is small ( @xcite ) .    as has been argued by lyth ( @xcite ) , the low - level of gravity waves is good news for our current ideas about realizing inflation in simple particle physics inspired models . in the most popular models today , the scalar modes are expected to dominate over the tensor modes by many orders of magnitude . the expectation is therefore that the tensor signal may not be measurable with any existing or planned experiments , or conversely that a positive detection of gravity waves would have profound implications for our ideas about inflation . however , for the time being , the constraints on the gravity wave contribution have not reached the level where we learn much about early universe physics  that will await future experiments . * 2.b.4 physics at @xmath2 *    it is possible to use the fact that the cmb anisotropies are largely as expected to limit the magnitude of any surprises at the last - scattering epoch . the arguments are much akin to those using the observed abundances of the light elements and big bang nucleosynthesis theory to limit ` exotic ' physics at early times . if something ` exotic ' would dramatically alter the theoretical predictions , it can be strongly constrained . a great many possible physical effects have been studied , but here we will list only a few things for which it is already possible to place observational bounds . strong limits exist on domains of anti - matter ( @xcite ) , particle decays near @xmath2 ( @xcite ) , primordial voids from an early phase transition ( @xcite ) and primordial magnetic fields ( @xcite ) , among other things . it was stated in the early 1960s , shortly before the discovery of the cmb , that there were only 2  facts in cosmology ( by peter scheuer , see @xcite ) . in a similar spirit , we have argued that there are perhaps 4 facts and 4 half facts currently known from cmb anisotropies . it has been recognized for some time that these anisotropies may answer some of our most fundamental questions about the universe . the current cmb data already indicate that gravitational instability , in a mostly cold dark matter dominated universe , amplified initially small adiabatic fluctuations into the large - scale structure that we see today . there is the potential to show what inflationary - like process happened in the early universe . and ultimately , the precise shape of the angular power spectrum holds the key to determining many of the fundamental cosmological parameters , either directly or in combination with other measurements . however , while it is interesting to track progress in this field and to speculate on what it all means , it seems clear that theorists have had long enough to manoeuvre that the present data no longer strongly constrain any popular cosmological model . with the coming of long duration balloon flights , the imminent launch of the _ map _  satellite , and the commissioning of three new cmb interferometers , we expect that to change . the boomerang team has already had a successful long duration balloon flight , and the analysis of that data set is eagerly awaited . similar flights will undoubtedly follow , along with other large data sets from new ground - based experiments . the race is on , since _ map _  is scheduled for late 2000 . a little later , sometime around 2006 , will see the launch of the planck surveyor . planck should supply us with essentially cosmic - variance limited information on all the angular scales relevant to primary anisotropies , over the full range of relevant frequencies . figure  2 is an estimate of how well the power spectrum might be constrained after _ map _  and after planck . with the proliferation of high precision data future ` a ' and ` b ' lists will be correspondingly longer and more detailed . our attempt at prognostication is represented in our list ` c ' :    *  cosmological parameters will be precisely determined *  polarization will be measured over a range of scales *  we will learn about early universe physics *  we will learn much about non - linear astrophysics    item  1 is in fact already happening , as discussed earlier . however , the current set of anisotropy data is not very constraining , since there is enough parameter freedom to fit models with quite different values of any individual parameter ( @xcite ) . this situation will undoubtedly improve in the future ( unless of course _ none _ of the current models fits the data , which is surely the most exciting prospect of all ! ) . certainly some degenerate parameter combinations will continue to exist in the model space ( particularly in models with the same ` angular diameter distance ' ) , but these degeneracies can be broken through combinations with other astrophysical data sets ( @xcite ) . if systematic errors can be kept under control , the combination of planck and data from redshift surveys will be particularly powerful at determining the cosmology . item  2 will be difficult , but we have no doubt that it will happen . _ map _  may yield some information , how much is difficult to estimate without more insight into the foregrounds . currently planned ground - based experiments may also give detections . and planck should provide polarization measurements over a reasonable range of scales . however , a full investigation of cmb polarization ( and certainly the ` curl ' or @xmath17-mode component produced by tensors ) may have to await an experiment even beyond planck . item  3 potentially involves information from both 1 and 2 . ultimately we will learn something about high energy physics through understanding the way in which fluctuations were laid down in the early universe , whether this involves discriminating any tensor component , measuring a changing spectral index , non - gaussian signatures , or something else . since the relevant energies are so far beyond what is achievable in particle accelerators , it is likely that cosmological phenomena will be the only way of constraining such models for quite some time . in addition to the ` initial conditions ' , the evolution of the fluctuations will provide us with information on the properties of the dark matter in the universe which may tie in directly to particle physics theories at the electroweak scale . item  4 includes a whole suite of potentially measurable effects , which can be thought of as processing the primary anisotropies . examples include gravitational lensing , non - linear potential growth , sunyaev - zeldovich effects , details of the reionization process , and extragalactic sources . there is a grey area between what is considered cosmic signal and what is considered a ` foreground ' . but whatever you call it , there is little doubt that data from the planck mission , for example , are likely to be mined for many years for the additional astrophysical information they contain . we expect rapid experimental progress in the next few years , and we trust that theoretical effort will be similarly feverish ( @xcite ) . as a result , there will no doubt be more physical processes uncovered which affect cmb anisotropies . at present the cdm - dominated inflationary paradigm looks like it s in pretty good shape . our ` c ' list may end up being quite inaccurate , and we can even imagine trouble for some entries in our ` b ' list . however , the spectral information from the cmb , together with the ` a ' list , provides a very solid foundation for the physics which generates the anisotropies . therefore we are confident that whatever proves to be the ultimate such list , a thorough investigation of cmb anisotropies will hold the key to learning about the background space - time and formation of large - scale structure in the universe .      99 barrow , j.d . , ferreira , p.g . , & silk , j , 1997 phys . rev . lett . , 78 , 3610 [ astro - ph/9701063 ] barrow , j.d liddle , a.r . , 1997 gen . , 29 , 1503 [ gr - qc/9705048 ] bartlett , j.g . , blanchard , a. , le dour , m. , douspis , m. , & barbosa , d. , 1998 in ` fundamental parameters in cosmology ' , moriond proceedings , in press [ astro - ph/9804330 ] bennett , c.l . , et al . , 1996 apj , 464 , l1 , astro - ph/9601067 bennett , c.l . , turner , m.s . , & white , m. , 1997 physics today 50 , 32 bond , j.r . , 1996 in ` microwave background anisotropies ' , proc . of xxxi recontres de moriond [ astro - ph/9610119 ] bond , j.r . , efstathiou , g. , & white , s.d.m . , 1992 mnras , 258 , 1p bond , j.r . , jaffe , a.h . , & knox , l.e . , 1998 submitted to apj [ astro - ph/9808264 ] borden , d. , ostriker , j.p . , & weinberg , d.h . , 1989 apj , 345 , 607 bunn , e.f . , ferreira , p. , & silk , j. , 1996 phys . rev . 77 , 2883 [ astro - ph/9605123 ] caldwell , r.r . , kamionkowski , m. , & wadley , l. , 1999 phys . , d59 , 027101 [ astro - ph/9807319 ] carr , b. , & ikeuchi , s. , 1985 mnras , 213 , 497 cornish , n.j . , spergel , d.n . , & starkman , g.d . , 1998 phys . d57 , 5982 [ astro - ph/9708225 ] de oliveira - costa , a. , smoot , g. , & starobinsky , a. , 1996 apj , 468 , 457 [ astro - ph/9510109 ] doroshkevich , a.g . , zeldovich , ya.b . , & sunyaev , r.a . , 1978 sov . astron . , 22 , 523 eisenstein , d.j . , hu , w. , & tegmark , m. , 1998 submitted to apj [ astro - ph/9807130 ] frieman , j. , hill , c.t . , & watkins , r. , 1992 phys . d46 , 1226 fuller , g. , & schramm , d.n . , 1992 phys . d45 , 2595 gooding , a.k . , park , c. , spergel , d.n . , turok , n. , & gott , r. , 1992 apj , 393 , 42 hancock , s. , rocha , g. , lasenby , a.n . , & guttirrez , c.m . , 1998 mnras , 294 , l1 [ astro - ph/9708254 ] hill , c.t . , schramm , d.n . , & fry , j. , 1989 comments . nucl . part . phys . , 19 , 25 hu , w. seljak , u. , white , m. , & zaldarriaga , m. , 1998 phys . rev . , d57 , 3290 [ astro - ph/9709066 ] hu , w. spergel , d.n . , & white , m. , 1997 phys . rev . , d55 , 3288 [ astro - ph/9605193 ] hu , w. , & white , m. , 1996 phys . lett . , 77 , 1687 [ astro - ph/9602020 ] hu , w. , & white , m. , 1997 apj , 479 , 568 [ astro - ph/9609079 ] hu , w. , & white , m. , 1997 newa , 2 , 323 [ astro - ph/9706147 ] hu , y. , turner , m.s . , & weinberg , e. , 1994 phys . d49 , 3830 [ astro - ph/9302002 ] ikeuchi , s. , 1981 pasj , 33 , 211 jaffe , a.h . , stebbins , a. , & frieman , j. , 1994 apj , 420 , 9 [ astro - ph/9301011 ] kinney , w.h . , kolb , e.w . , & turner , m.s . , 1997 phys . lett . , 79 , 2620 [ astro - ph/9704070 ] kogut , a. , hinshaw , g. , & banday , a.j . , 1997 phys . d55 , 1901 [ astro - ph/9701090 ] kolb , e.w . , & turner , m.s . , 1990 ` the early universe ' , addison - wesley , redwood city , p.@xmath18 krauss , l.m . , & white , m. , 1992 phys . rev . lett . , 69 , 869 levin , j. , barrow , j.d . , bunn , e.f . , & silk , j. , 1997 phys . 79 , 974 [ astro - ph/9702242 ] liddle , a.r . , 1994 phys . , d49 , 3805 [ astro - ph/9307036 ] liddle , a.r . , 1995 phys . , d51 , 5347 [ astro - ph/9410083 ] lineweaver , c.h . , 1998 apj , 505 , l69 [ astro - ph/9805326 ] longair , m.s . , 1993 qjras , 34 , 157 lyth , d.h . , 1997 phys . 78 , 1861 [ hep - ph/9606387 ] markevich , a.v . , & starobinsky , a.a . , 1996 astron . , 22 , 431 nordberg , h.p . , & smoot , g.f . , 1998 submitted to apj [ astro - ph/9805123 ] oh , s.p . , spergel , d.n . , & hinshaw g. , 1999 apj , in press [ astro - ph/9805339 ] ostriker , j.p . , & cowie , l. , 1981 apj , 243 , l127 ostriker , j.p . , & strassler , m.j . , 1989 apj , 338 , 579 ostriker , j.p . , & thompson , c. , 1987 apj , 323 , l97 ostriker , j.p . , thompson , c. , & witten , e. , 1986 phys . lett . , b280 , 231 partridge , r.b . , richards , e.a . , fomalont , e.b . , kellerman , k.i . , & windhorst , r.a . , 1997 apj , 483 , 38 pierpaoli , e. , & bonometto , s.a . , 1998 submitted to mnras [ astro - ph/9806037 ] press , w.h . , ryden , b. , & spergel , d. , 1990 phys . rev . lett . , 64 , 1084 python , m. , 1979 quoted in ` the life of brian ' , mandarin paperbacks , london ( 1993 ) sachs , r.k . , & wolfe , a.m. , 1967 apj , 147 , 73 sakai , n. , sugiyama , n. , & yokoyama , j. , 1997 preprint [ astro - ph/9712153 ] salopek , d.s . , 1995 in ` current topics in astrofundamental physics : the early universe ' , ed . n.  sanchez , a.  zichichi , kluwer , p.@xmath19 [ astro - ph/9412012 ] scott , d. , silk , j. , & white , m. , 1995 science , 268 , 829 scott , d. , & white , m. , 1996 comments on astrophysics , 18 , 289 [ astro - ph/9601170 ] seljak , u. , & zaldarriaga , m. , 1996 apj , 469 , 437 [ astro - ph/9603033 ] smoot , g.f . , 1991 in ` observational tests of cosmological inflation ' , ed .  t.  shanks , et al . , kluwer , p.@xmath20 smoot , g.f . , et al . , 1992 apj , 396 , l1 smoot , g.f . , & scott , d. , 1998 in caso c. , et al . j. , c3 , 1 , the review of particle physics . p.@xmath21 [ astro - ph/9711069 ] souradeep , t. , pogosyan , d. , & bond , j.r . , 1998 in ` fundamental parameters in cosmology ' , proc . of xxxiii recontres de moriond , in press [ astro - ph/9804042 ] stevens , d. , scott , d. , & silk , j. , 1993 phys . lett . , 71 , 20 subramanian , k. , & barrow , j.d . , 1998 phys . 81 , 3575 [ astro - ph/9803261 ] sugiyama , n. , silk , j. , & vittorio , n. , 1993 apj , 419 , l1 tegmark , m. , 1998 submitted to apj [ astro - ph/9809201 ] turner , m.s . , 1997 phys . d55 , 435 vishniac , e.t . , ostriker , j.p . , & bertschinger , e. , 1985 apj , 291 , 399 wandel , a. , 1985 apj , 294 , 385 wasserman , i. , 1986 phys . lett . , 57 , 2234 weinberg , d.h . , dekel , a. , & ostriker , j.p . , 1989 apj , 336 , 9 white , m. , 1998 apj , 506 , 495 [ astro - ph/9802295 ] white , m. , & bunn , e.f . , 1995 apj , 450 , 477 [ astro - ph/9503054 ] white , m. , & hu , w. , 1997 a&a , 321 , 8 [ astro - ph/9609105 ] zibin , j. , scott , d. , & white , m. , 1999 submitted to phys . d [ astro - ph/9901028 ]
Abstract:

Answer: the cosmic background explorer ( _ cobe _ ) satellite , and the differential microwave radiometer ( dmr ) experiment in particular , was extraordinarily successful . however , the dmr results were announced about 7 years ago , during which time a great deal more has been learned about anisotropies in the cosmic microwave background ( cmb ) . we assess the current state of knowledge , and discuss where we might be going . the cmb experiments currently being designed and built , including long - duration balloons , interferometers , and two space missions , promise to address several fundamental cosmological issues . we present our evaluation of what we already know , what we are beginning to learn now , and what the future may bring . _ all right . but apart from the sanitation , the medicine , education , wine , public order , irrigation , roads , the fresh water system , and public health  what have the romans ever done for us ? _


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: the large - order behavior of a perturbative expansion in gauge theories is inevitably dominated by the factorial growth of renormalon diagrams [ 1 - 4 ] . in the case of quantum chromodynamics ( qcd ) , the coefficients of perturbative expansions in the qcd coupling @xmath1 can increase dramatically even at low orders . this fact , together with the apparent freedom in the choice of renormalization scheme and renormalization scales , limits the predictive power of perturbative calculations , even in applications involving large momentum transfers , where @xmath1 is effectively small . a number of theoretical approaches have been developed to reorganize the perturbative expansions in an effort to improve the predictability of the perturbative qcd ( pqcd ) . for example , optimized scale and scheme choices have been proposed , such as the method of effective charges ( ech ) [ 5 ] , the principle of minimal sensitivity ( pms ) [ 6 ] , and the brodsky - lepage - mackenize ( blm ) scale - setting prescription [ 7 ] and its generalizations [ 8 - 20 ] . in ref.[4 ] , the resummation of the formally divergent renormalon series and the parametrization of related higher - twist power - suppressed contributions has been given . in general , a factorially divergent renormalon series arises when one integrates over the logarithmically running coupling @xmath2 in a loop diagram . such contributions do not occur in conformally invariant theories which have a constant coupling . of course , in the physical theory , the qcd coupling does run .    among the fundamental predictions of qcd are asymptotic scaling laws for large - angle exclusive processes [ 21 - 28 ] . qcd counting rules were formalized in refs.[22,23 ] . these reactions probe hadronic constituents at large relative momenta , or equivalently , the hadronic wave function at short distances . in particular , the meson wave function determines the leading higher - twist contribution to meson production at high @xmath0 . the hadronic wave functions in terms of quark and gluon degrees of freedom play an important role in the quantum chromodynamics predictions for hadronic processes . in the perturbative qcd theory , the hadronic distribution amplitudes and structure functions which enter exclusive and inclusive processes via the factorization theorems at high momentum transfers can be determined by the hadronic wave functions , and therefore they are the underlying links between hadronic phenomena in qcd at large ( nonperturbative ) and small distances ( perturbative ) . if the hadronic wave functions were accurately known , then we could calculate the hadronic distribution amplitude and structure functions for exclusive and inclusive processes in qcd . the study of hadron production at large transverse momenta in hadronic interactions is a valuable testing ground of the perturbative regime of qcd , providing information on both the parton distribution functions ( pdfs ) in the proton , and the parton - to - hadron fragmentation functions ( ffs)[29 ] . the frozen coupling constant approach can be applied for investigation , not only for exclusive processes , but also for the calculation of higher - twist contributions to some inclusive processes , for example as large -@xmath3 meson photoproduction [ 30 ] , two - jet+meson production in the electron - positron annihilation [ 31 ] . in the frozen coupling constant approach was used refs.[30,31 ] for calculation of integrals , such as @xmath4 according to ref.[7 ] , in pqcd calculations , the argument @xmath5 of the running coupling constant should be taken equal to the square of the momentum transfer of a hard gluon in a corresponding feynman diagram , in both the renormalization and factorization scale . but defined in this way , @xmath6 suffers from infrared singularities . for example in ref.[32 ] , @xmath5 equals to @xmath7 and @xmath8 , where @xmath9 , @xmath10 are the subprocess s mandelstam invariants . therefore , in the soft regions @xmath11 , @xmath12 , the integrals in ( 1.1 ) diverge and we need some regularization methods of @xmath13 in these regions for their calculation . in ref.[33 ] , the authors investigated the phenomenology of infrared renormalons in inclusive processes . the dispersive approach has been devised to extend properly modified perturbation theory calculations towards the low - energy region [ 34 ] . connections between power corrections for the three deep inelastic scattering sum rules have also been explored in ref.[35 ] . investigation of the infrared renormalon effects in various inclusive and exclusive processes is one of the most important and interesting problems in the perturbative qcd . it is known that infrared renormalons are responsible for factorial growth of coefficients in perturbative series for the physical quantities . but , these divergent series can be resummed by means of the borel transformation [ 1 ] and the principal value prescription [ 36 ] , and effects of infrared renormalons can be taken into account by a scale - setting procedure @xmath14 at the one - loop order results . technically , all - order resummation of infrared renormalons corresponds to the calculation of the one - loop feynman diagrams with the running coupling constant @xmath15 at the vertices or , alternatively , to calculation of the same diagrams with nonzero gluon mass . studies of infrared renormalon problems have also opened new prospects for evaluation of power - suppressed corrections to processes characteristics [ 37 ] . power corrections can also be obtained by means of the landau - pole free expression for the qcd coupling constant . the most simple and elaborated variant of the dispersive approach , the shirkov and solovtsov analytic perturbation theory , was formulated in ref.[38 ] .    by taking these points into account , it may be argued that the analysis of the next - to - leading order higher - twist effects on the dependence of the pion wave function in pion production at proton - proton collisions by the running coupling approach , are significant from both theoretical and experimental[39 ] points of view .    in this work we apply the running coupling approach[40 ] in order to compute the effects of the infrared renormalons on the meson production in proton - proton collisions in next - to - leading order . this approach was also employed previously[41,42 ] to calculate the inclusive meson production in proton - proton and photon - photon collisions . the running coupling method in next - to - leading order pion and kaon electromagnetic form factor was computed in [ 43 ] .    as we know , power - suppressed contributions to exclusive processes in qcd , which are commonly referred to as higher - twist corrections . the higher - twist approximation describes the multiple scattering of a parton as power corrections to the leading - twist cross section . we will show that higher - twist terms contribute substantially to the inclusive meson cross section at moderate transverse momenta . in addition , we shall demonstrate that higher - twist reactions necessarily dominate in the kinematic limit where the transverse momentum approaches the phase - spase boundary . a precise measurement of the inclusive charged pion production cross section at @xmath16 and @xmath17 is important for the proton - proton collisions program at the relativistic heavy ion collider ( rhic ) at the brookhaven national laboratory . another important aspect of this study is the choice of the meson model wave functions . in this respect , the contribution of the higher - twist feynman diagrams to a pion production cross section in proton - proton collisions is been computed by using various pion wave functions . also , higher - twist contributions which are calculated by the running coupling constant and frozen coupling constant approaches are been estimated and compared to each other . within this context , this paper is organized as follows : in sec.[ht ] , we provide formulas for the calculation of the contribution of the high twist diagrams . in sec . [ ir ] we present formulas and an analysis of the next - to - leading order higher - twist effects on the dependence of the pion wave function by the running coupling constant approach . in sec . [ lt ] , we provide formulas for the calculation of the contribution of the leading - twist diagrams . in sec . [ results ] , we give the numerical results for the cross section and discuss the dependence of the cross section on the pion wave functions . we present our conclusions in sec . the higher - twist feynman diagrams , which describe the subprocess @xmath18 for the pion production in the proton - proton collision are shown in fig.1 . in the higher - twist diagrams , the pion of a proton quark is directly observed . their @xmath19 power suppression is caused by a hard gluon exchange between pion constituents . the amplitude for this subprocess can be found by means of the brodsky - lepage formula [ 27 ] :    @xmath20    in eq.(2.1 ) , @xmath21 is the sum of the graphs contributing to the hard - scattering part of the subprocess . the hard - scattering part for the subprocess under consideration is @xmath22 , where a quark and antiquark form a pseudoscalar , color - singlet state @xmath23 . here @xmath24 is the pion wave function , i.e. , the probability amplitude for finding the valence @xmath25 fock state in the meson carry fractions @xmath26 and @xmath27 , @xmath28 . remarkably , this factorization is gauge invariant and only requires that the momentum transfers in @xmath21 be large compared to the intrinsic mass scales of qcd . since the distribution amplitude and the hard - scattering amplitude are defined without reference to the perturbation theory , the factorization is valid to leading order in @xmath29 , independent of the convergence of perturbative expansions . the hard - scattering amplitude @xmath21 can be calculated in perturbation theory and represented as a series in the qcd running coupling constant @xmath30 . the function @xmath31 is intrinsically nonperturbative , but its evolution can be calculated perturbatively . in our calculation , we have neglected the pion and the proton masses . turning to extracting the contributions of the higher - twist subprocesses , there are many kinds of leading - twist subprocesses in @xmath32 collisions as the background of the higher - twist subprocess @xmath33 , such as @xmath34 , @xmath35 , @xmath36 etc . the contributions from these leading - twist subprocesses strongly depend on some phenomenological factors , for example , quark and gluon distribution functions in the proton and fragmentation functions of various constituents , etc . most of these factors have not been well determined , neither theoretically nor experimentally . thus they cause very large uncertainty in the computation of the cross section of process @xmath37 . in general , the magnitude of this uncertainty is much larger than the sum of all the higher - twist contributions , so it is very difficult to extract the higher - twist contributions . the production of a hadron at large transverse momentum , @xmath0 , in a hadronic collisions is conventionally analyzed within the framework of perturbative qcd by convoluting the leading - twist @xmath38 hard subprocess cross sections with evolved structure and fragmentation functions . the most important discriminant of the twist of a perturbative qcd subprocess in a hard hadronic collision is the scaling of the inclusive invariant cross section [ 22 - 25,44 ] ,    @xmath39    at fixed @xmath40 and center - of - mass angle @xmath41 . as we know in the original parton model [ 45 ] the power fall - off is simply @xmath42 since the @xmath43 subprocess amplitude for point - like partons is scale invariant , and there is no dimensional parameter as in a conformal theory . the detected hadron @xmath44 can be produced directly in the hard subprocess reaction as in an exclusive reaction , then such direct higher - twist processes can give a significant contribution since there is no suppression from jet fragmentation at large momentum fraction carried by the hadron , @xmath45 , and the trigger hadron is produced without any waste of energy[44 ] . the mandelstam invariant variables for subprocesses @xmath18 are defined as @xmath46    in our calculation , we have also neglected the quark masses . we have aimed to calculate the pion production cross section and to fix the differences due to the use of various pion model functions . we have used five different wave functions : the asymptotic wave function ( asy ) , the chernyak - zhitnitsky wave function [ 28,46 ] , the cleo wave function [ 47 ] , the braun - filyanov pion wave functions [ 48 ] and the bakulev - mikhailov - stefanis pion wave function [ 49 ] . it should be noted that the wave functions of pions also are developed in refs.[50 - 52 ] by the dubna group :    @xmath47,\ ] ] @xmath48,\ ] ] @xmath49,\ ] ] @xmath50,\ ] ] @xmath51 @xmath52 @xmath53 where @xmath54 is the pion decay constant . here , we have denoted by @xmath55 , the longitudinal fractional momentum carried by the quark within the meson . then , @xmath56 and @xmath57 . the pion wave function is symmetric under the replacement @xmath58 . several important nonperturbative tools have been developed which allow specific predictions for the hadronic wave functions directly from theory and experiments . the qcd sum - rule technique and lattice gauge theory provide constraints on the moments of the hadronic distribution amplitude . however , the correct pion wave function is still an open problem in qcd . it is known that the pion wave function can be expanded over the eigenfunctions of the one - loop brodsky - lepage equation , _ i.e. _ , in terms of the gegenbauer polynomials @xmath59    @xmath60.\ ] ]    the evolution of the wave function on the factorization scale @xmath61 is governed by the functions @xmath62 , @xmath63^{\gamma_n/\beta_0},\ ] ] @xmath64    in eq.(2.6 ) , @xmath65 are anomalous dimensions defined by the expression ,    @xmath66.\ ] ]    the constants @xmath67 are input parameters that form the shape of the wave functions and which can be extracted from experimental data or obtained from the nonperturbative qcd computations at the normalization point @xmath68 . the qcd coupling constant @xmath13 at the one - loop approximation is given by the expression ,    @xmath69    here , @xmath70 is the fundamental qcd scale parameter , @xmath71 is the qcd beta function one - loop coefficient , respectively , @xmath72 the cross section for the higher - twist subprocess @xmath73 is given by the expression @xmath74 ^ 2}{{\hat s}^3}\left[\frac{1}{{\hat u}^2}+\frac{1}{{\hat t}^2}\right],\ ] ] where @xmath75+e_2\hat u\int_{0}^{1}dx\left[\frac{\alpha_{s}(\hat q_2 ^ 2)\phi_{\pi}(x,\hat q_2 ^ 2)}{1-x}\right].\ ] ] here @xmath76and @xmath77 , represent the momentum squared carried by the hard gluon in fig.1 , @xmath78 is the charge of @xmath79 and @xmath80 . the higher - twist contribution to the large-@xmath3 pion production cross section in the process @xmath81 is [ 53 ] : @xmath82 @xmath83 @xmath84 @xmath85 @xmath86 @xmath87 @xmath88 @xmath89 @xmath90    where @xmath91  is the transverse mass of pion , which is given by @xmath92 let us first consider the frozen coupling approach . in this approach , we take the four - momentum square @xmath93 of the hard gluon to be equal the pion s transverse momentum square @xmath94 . in this case , the qcd coupling constant @xmath95 in the integral ( 2.10 ) does not depend on the integration variable . after this substitution calculation of integral ( 2.10 ) becomes easy . hence , the effective cross section obtained after substitution of the integral ( 2.10 ) into the expression ( 2.9 ) is referred as the frozen coupling effective cross section . we will denote the higher - twist cross section obtained using the frozen coupling constant approximation by @xmath96 . we have extracted the following higher - twist subprocesses contributing to the two covariant cross sections in eq.(2.11 ) @xmath97 by charge conjugation invariance , we also have @xmath98 in this section , we shall calculate the integral ( 2.10 ) using the running coupling constant method and also discuss the problem of normalization of the higher - twist process cross section in the context of the same approach .    as is seen from ( 2.10 ) , in general , one has to take into account not only the dependence of @xmath99 on the scale @xmath100 , but also an evolution of @xmath101 with @xmath100 . the meson wave function evolves in accordance with a bethe - salpeter - type equation . therefore , it is worth noting that , the renormalization scale ( argument of @xmath95 ) should be equal to @xmath102 , @xmath103 , whereas the factorization scale [ @xmath61 in @xmath104 is taken independent from @xmath105 , we assume @xmath106 . such a approximation does not considerably change the numerical results , but the phenomenon considered in this article ( effect of infrared renormalons ) becomes transparent . the main problem in our investigation is the calculation of the integral in ( 2.10 ) by the running coupling constant approach . this integral in the framework of the running coupling approach takes the form    @xmath107    the @xmath108 has the infrared singularity at @xmath109 , if @xmath110 or @xmath111 , if @xmath112 and as a result integral @xmath113 diverges ( the pole associated with the denominator of the integrand is fictitious , because @xmath114 , and therefore , the singularity of the integrand at @xmath115 is caused only by @xmath116 ) . for the regularization of the integral , we express the running coupling at scaling variable @xmath116 with the aid of the renormalization group equation in terms of the fixed one @xmath13 . the renormalization group equation for the running coupling @xmath117 has the form [ 36 ]    @xmath118 ^ 2-\frac{\beta_{1}}{16}[\alpha(\lambda q^2)]^3\ ] ]    where @xmath119 the solution of eq.(3.2 ) , with the initial condition @xmath120 is [ 35 ]    @xmath121^{-1}.\ ] ]    this transcendental equation can be solved iteratively by keeping the leading @xmath122 and next - to - leading @xmath123 powers . for @xmath124 these terms are given by @xmath125}{[1+\ln\lambda / t]^2}\ ] ]    the first term in eq.(3.4 ) is the solution of the renormalization group eq.(3.2 ) with leading power accuracy , whereas the whole expression ( 3.4 ) is the solution of eq.(3.2 ) with next - to - leading power accuracy . after substituting eq.(3.4 ) into eq.(2.10 ) we get : @xmath126 @xmath127 @xmath128 @xmath129^{\gamma_{n}/\beta_{0}}c_{n}^{3/2}(2x-1)\right]}{(1-x ) ( 1+\ln\lambda / t_{1})}-\ ] ] @xmath130^{\gamma_{n}/\beta_{0}}c_{n}^{3/2}(2x-1)\right ] \ln(1+\ln\lambda / t_{1})}{(1-x ) ( 1+\ln\lambda / t_{1})^2}+\ ] ] @xmath131^{\gamma_{n}/\beta_{0}}c_{n}^{3/2}(2x-1)\right]}{(1-x ) ( 1+\ln\lambda / t_{2})}-\ ] ] @xmath132^{\gamma_{n}/\beta_{0}}c_{n}^{3/2}(2x-1)\right]\ln(1+\ln\lambda / t_{2})}{(1-x ) ( 1+\ln\lambda / t_{2})^2}\ ] ]    where @xmath133 and @xmath134 . it should be note that in [ 41 ] , we used @xmath135 . the integral ( 3.5 ) is common and , of course , still divergent , but now it is recast into a form , which is suitable for calculation . using the running coupling constant approach , this integral may be found as a perturbative series in @xmath13 @xmath136 the expression coefficients @xmath137 can be written as power series in the number of light quark flavors or , equivalently , as a series in power of @xmath71 : @xmath138 the coefficients @xmath139 of this series demonstrate factorial growth @xmath140 , which might indicate an infrared renormalon nature of divergences in the integral ( 3.5 ) and corresponding series ( 3.6 ) . the procedure for dealing with such ill - defined series is well known : one has to perform the borel transform of the series [ 15 ] @xmath141(u)=\sum_{n=0}^{\infty}\frac{d_n}{n ! } u^n,\ ] ] then invert @xmath142(u)$ ] to obtain the resummed expression ( the borel sum ) @xmath143 . after this we can find directly the resummed expression for @xmath144 . the change of the variable @xmath105 to @xmath145 , as @xmath146 . then @xmath147 @xmath148    the basic theoretical problem is how to define the borel sum @xmath142(u)$ ] of the integral in eq.(3.7 ) for the quantities we are interested in . in qcd this problem is usually solved using perturbative methods and calculating the corresponding multiloop feynman diagrams with a one - gluon line , dressed by the chains of fermion bubbles .    for the calculation the expression ( 3.7 ) we will apply the integral representations of @xmath149 and @xmath150 [ 54,55 ] . after this operation , formula ( 3.7 ) is simplified and we can extract the borel sum of the perturbative series ( 3.6 ) and the corresponding borel transform in dependence from the wave functions of the meson , respectively . also after such manipulations the obtained expression can be used for numerical computations . it is convenient to use the following integral representation for @xmath149 and @xmath151 @xmath152 and @xmath153 where @xmath154 is the euler - mascheroni constant . after inserting eq.(3.8 ) and ( 3.9 ) into ( 3.7 ) . then , we obtain @xmath155 @xmath156 @xmath157 @xmath158 in the case of @xmath159 for @xmath143 , we get @xmath160-\frac{4\sqrt{3 } \pi f_{\pi } e_{1 } \beta_{1 } \hat t}{\beta_{0}^3 } \cdot\ ] ] @xmath161(1-c-\ln u-\ln { t_1})u+\ ] ] @xmath162-\frac{4\sqrt{3 } \pi f_{\pi } e_{2}\beta_{1 } \hat u}{\beta_{0}^3}\cdot\ ] ] @xmath163(1-c-\ln u-\ln { t_2})u.\ ] ] in the case of the @xmath164 wave function , we find @xmath165 @xmath166^{50/81 } \left[\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}- \left.\frac{20}{4-u}\right]\right]-\ ] ] @xmath167 @xmath166^{50/81 } \left[\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}- \left.\frac{20}{4-u}\right]\right](1-c-\ln u-\ln { t_1})u+\ ] ] @xmath168 @xmath166^{50/81 } \left[\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}- \left.\frac{20}{4-u}\right]\right]-\ ] ] @xmath169 @xmath166^{50/81 } \left[\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}- \left.\frac{20}{4-u}\right]\right](1-c-\ln u-\ln { t_2})u\ ] ]    in the case of the @xmath170 wave function , we get @xmath171^{50/81}\cdot\ ] ] @xmath172 - 0.4125\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right]-\ ] ] @xmath174^{50/81}\cdot\ ] ] @xmath172 - 0.4125\left[\frac{\alpha_{s}(q^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right ] ( 1-c-\ln u-\ln { t_1})u+\ ] ] @xmath175^{50/81}\cdot\ ] ] @xmath172 - 0.4125\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right]-\ ] ] @xmath176^{50/81}\cdot\ ] ] @xmath172 - 0.4125\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right ] ( 1-c-\ln u-\ln { t_2})u\ ] ]    also , in the case of the @xmath177 wave function , we get @xmath178^{50/81}\cdot\ ] ] @xmath172 - 0.244\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right]-\ ] ] @xmath179^{50/81}\cdot\ ] ] @xmath172 - 0.244\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right](1-c-\ln u-\ln { t_1})u+\ ] ] @xmath180^{50/81}\cdot\ ] ] @xmath172 - 0.244\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right]-\ ] ] @xmath181^{50/81}\cdot\ ] ] @xmath172 - 0.244\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right ] ( 1-c-\ln u-\ln { t_2})u\ ] ]    equations(3.11)-(3.14 ) is nothing more than the borel sum of the perturbative series ( 3.6 ) , and the corresponding borel transform in the case @xmath159 is @xmath141(u)=\frac{1}{1-u}-\frac{1}{2-u},\ ] ] in the case @xmath164 is @xmath141(u)=\frac{1}{1-u}-\frac{1}{2-u}+0.84\left(\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right)^{50/81 } \left(\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}-\frac{20}{4-u}\right),\ ] ]    in the case @xmath170 is @xmath141(u)=\frac{1}{1-u}-\frac{1}{2-u}+0.405\left(\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right)^{50/81 } \left(\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}-\frac{20}{4-u}\right)-\ ] ] @xmath182 and in the case @xmath177 is @xmath141(u)=\frac{1}{1-u}-\frac{1}{2-u}+0.282\left(\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right)^{50/81 } \left(\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}-\frac{20}{4-u}\right)-\ ] ] @xmath183 the series ( 3.6 ) can be recovered by means of the following formula @xmath184(u)\mid_{u=0}\ ] ] the borel transform @xmath142(u)$ ] has poles on the real @xmath185 axis at @xmath186 which confirms our conclusion concerning the infrared renormalon nature of divergences in ( 3.6 ) . to remove them from eqs.(3.11)-(3.14 ) we applied the principal value prescription . hence , the effective cross section obtained after substitution of the expressions ( 3.10 - 3.14 ) into the expression ( 2.10 ) is referred as the running coupling effective cross section . we will denote the higher - twist cross section obtained using the running coupling constant approach by @xmath187 regarding the higher - twist corrections to the pion production cross section , a comparison of our results with leading - twist contributions is crucial . we take two leading - twist subprocesses for the pion production:(1 ) quark - antiquark annihilation @xmath188 , in which the @xmath189 meson is indirectly emitted from the gluon , @xmath190 and ( 2 ) quark - gluon fusion , @xmath191 , with subsequent fragmentation of the final quark into a meson , @xmath192 . the corresponding cross sections are obtained in @xmath193 @xmath194 for the leading - twist contribution , we find @xmath195 @xmath196 where @xmath197 @xmath198 and @xmath199 represents gluon and quark fragmentation functions into a meson containing gluon and quark of the same flavor . in the leading - twist subprocess , the @xmath189 meson is indirectly emitted from the gluon and quark with the fractional momentum @xmath45 . the @xmath200 function can be expressed in terms of the parton kinematic variables , and the @xmath45 integration can then be done . the final form for the cross section is @xmath201 @xmath202 @xmath203 @xmath204 in this section , we discuss the numerical results for next - to - leading order higher - twist effects with higher - twist contributions calculated in the context of the running coupling and frozen coupling approaches on the dependence of the chosen meson wave functions in the process @xmath205 . in the calculations , we use the asymptotic wave function @xmath206 , the chernyak - zhitnitsky @xmath207 , the cleo pion wave function [ 41 ] , the braun - filyanov pion wave functions [ 42 ] , and the bakulev - mikhailov - stefanis pion wave function [ 43 ] . for the higher - twist subprocess , we take @xmath208 and we have extracted the following four higher - twist subprocesses contributing to @xmath209 cross sections : @xmath210 @xmath211 , @xmath212 , @xmath213 contributing to cross sections . for the dominant leading - twist subprocess for the pion production , we take the quark - antiquark annihilation @xmath188 , in which the @xmath189 meson is indirectly emitted from the gluon and quark - gluon fusion , @xmath214 , with subsequent fragmentation of the final quark into a meson , @xmath215 . as an example for the quark distribution function inside the proton , the mrst2003c package [ 56 ] has been used . the higher twist subprocesses probe the meson wave functions over a large range of @xmath61 squared momentum transfer , carried by the gluon . therefore , in the diagram given in fig.1 we take @xmath216 , @xmath103 , which we have obtained directly from the higher - twist subprocesses diagrams . the same @xmath61 has been used as an argument of @xmath30 in the calculation of each diagram . the results of our numerical calculations are plotted in figs.2 - 15 . first of all , it is very interesting to compare the resummed higher - twist cross sections with the ones obtained in the framework of the frozen coupling approach . in figs.2 - 4 we show the dependence of next - to - leading order higher - twist cross sections @xmath217 calculated in the context of the running coupling constant approach and the ratios @xmath218/@xmath219 , @xmath217/@xmath220 as a function of the pion transverse momentum @xmath3 for different pion wave functions at @xmath221 . it is seen that the values of cross sections @xmath217 , and ratios for fixed @xmath222 and @xmath223 depend on the choice of the pion wave function . as seen from fig.2 the next - to - leading order higher - twist differential cross section is monotically decreasing with an increase in the transverse momentum of the pion . in figs.3 - 5 , we show the dependence of the ratios @xmath218/@xmath219 , @xmath217/@xmath220 , @xmath217/@xmath224 as a function of the pion transverse momentum @xmath3 for different pion wave functions . here @xmath217 , @xmath225 , @xmath224 are the higher - twist cross sections calculated in the context of the running coupling method , in the framework of the frozen coupling approach and is the leading - twist cross section , respectively . as seen from fig.5 , in the region @xmath226 next - to - leading order higher - twist cross section calculated in the context of the running coupling method is suppressed by about 1 - 2 orders of magnitude relative to the leading - twist cross section , but in the region @xmath227 is comparable with the cross section of leading - twist . in figs.6 - 8 , we have depicted next - to - leading order higher - twist cross sections@xmath217 , ratios @xmath228 , @xmath229 as a function of the rapidity @xmath222 of the pion at @xmath16 and @xmath230 . at @xmath16 and @xmath230 , the pion rapidity lies in the region @xmath231 . as seen from figs.6 - 8 , next - to - leading order higher - twist cross section and ratios have a different distinctive . in the region ( @xmath232 ) the next - to - leading order higher - twist cross section @xmath217 , is suppressed by about 1 - 2 order of magnitude relative to the leading - twist cross section . in the region ( @xmath233 ) , the ratio for all wave functions increase with an increase of the @xmath222 rapidity of the pion and has a maximum approximately at the point @xmath234 . besides that , the ratio decreases with an increase in the @xmath222 rapidity of the pion . as is seen from figs.7 - 8 , the ratio @xmath235 very sensitive to the choice of the meson wave functions . but , as seen from fig.8 , the ratio @xmath217/@xmath224 for all wave functions has a minimum approximately at the point @xmath234 . also , the distinction between @xmath236 with @xmath237 , @xmath238 , @xmath239 and @xmath240 have been calculated . for example , in the case of @xmath16 , @xmath221 , the distinction between @xmath236 with @xmath241(i = cleo , cz , bf , bms ) as a function of the pion transverse momentum @xmath3 is shown in table [ table1 ] . thus , the distinction between @xmath236 and @xmath242 is maximum at @xmath243 , with @xmath244 at @xmath245 but the distinction between @xmath236 with @xmath246 is minimum at @xmath245 , with @xmath244 at @xmath243 and increase with an increase in @xmath0 . such a behavior of @xmath235 may be explained by reducing all moments of the pion model wave functions to those of @xmath159 for high @xmath61 . also , we have calculated the distinction between @xmath236 with @xmath237 , @xmath238 , @xmath239 and @xmath240 as a function of the rapidity @xmath222 of the pion . for example , in the case of @xmath247 , @xmath248 the distinction between @xmath236 with @xmath241(i = cleo , cz , bf , bms ) as a function of the rapidity @xmath222 of the pion is presented in table [ table2 ] . we have also carried out comparative calculations in the center - of - mass energy @xmath17 . the results of our numerical calculations in the center - of - mass energies @xmath249 are plotted in figs.9 - 15 . analysis of our calculations at the center - of - mass energies @xmath16 and @xmath17 , show that with the increase in beam energy values of the cross sections , ratio @xmath250 , and contributions of next - to - leading order higher - twist to the cross section decrease by about 1 - 2 order . therefore the experimental investigation of higher - twist effects include renormalon effects conveniently in low energy . on the other hand , the higher - twist corrections and ratios @xmath235 a very sensitive to the choice of the pion wave function . analysis of our calculations at the center - of - mass energies @xmath16 and @xmath17 , show that values of the next - to - leading order cross sections decrease by about @xmath251 percent of magnitude relative to the leading - order cross sections . also , the distinction between @xmath236 with @xmath237 , @xmath238 , @xmath239 and @xmath240 have been calculated . for example , in the case of @xmath17 , @xmath221 , the distinction between @xmath236 with @xmath241(i = cleo , cz , bf , bms ) as a function of the pion transverse momentum @xmath3 is shown in table [ table3 ] . thus , the distinction between @xmath236 with @xmath241,(i= cleo , bf , ) is maximum at @xmath252 , with @xmath244 at @xmath253 , but the distinction between @xmath236 with @xmath238 , @xmath237 , @xmath239 is minimum at @xmath253 , with @xmath244 at @xmath254 and increase with an increase in @xmath0 . also , we have calculated the distinction between @xmath236 with @xmath237 , @xmath238 , @xmath239 and @xmath240as a function of the rapidity @xmath222 of the pion . for example , in the case of @xmath255 , @xmath256 the distinction between @xmath236 with @xmath241 as a function of the rapidity @xmath222 of the pion is presented in table [ table4 ] . the calculations show that the ratio @xmath241/@xmath236 , ( i = cleo , cz , bf , bms ) for all values of the transverse momentum @xmath0 of the pion identically equivalent to ratio @xmath257/@xmath258 .    in our calculations of the next - to - leading order higher - twist cross section of the process the dependence of the transverse momentum of meson appears in the range of @xmath259 . therefore , higher - twist cross section obtained in our paper should be observable at rhic in this work we have calculated the inclusive meson production via higher - twist mechanism and obtained the expressions for the subprocess @xmath260 cross section for mesons with symmetric wave functions . for calculation of the next - to - leading order cross section we have applied the running coupling constant method and revealed infrared renormalon poles in the cross section expression . infrared renormalon induced divergences have been regularized by the means of the principal value prescripton and the resummed expression ( the borel sum ) for the next - to - leading order higher - twist cross section has been found . in the resummed higher - twist cross section differs considerably from that found using the frozen coupling approach , some regions . also we have demonstrated that next - to - leading order higher - twist contributions to meson production cross section in the proton - proton collisions have important phenomenological consequences . our investigation enables us to conclude that the higher - twist pion production cross section in the proton - proton collisions depends on the form of the pion model wave functions and may be used for their study . analysis of our calculations shows that the magnitude of next - to - leading order cross sections calculated in the running coupling approach in some regions is larger than the leading - twist cross sections in 1 - 2 order . further investigations are needed in order to clarify the role of high twist effects in this process . we have demonstrated that the resummed result depends on the pion model wave functions used in calculations . the proton - proton collisions provide us with a new opportunity to probe a proton s internal structure . in particular , meson production in proton - proton collisions takes into account infrared renormalon effects : this opens a window toward new types of parton distributions which can not be measured by the deep inelastic lepton - proton scatterings . finally , we discuss the phenomenological consequences of possible higher - twist contributions to meson production in proton - proton collisions in next - to - leading order at rhic . future rhic measurements will provide further tests of the dynamics of large @xmath0 hadron production beyond leading twist . the work presented in this paper was completed while one of the authors , a.i.ahmadov , was visiting the hecap section of the abdus salam ictp , trieste , italy . he would like to express his gratitude to the members of the section especially to the head prof . s. randjbar - daemi for their hospitality . financial support by ictp is also gratefully acknowledged . hooft , in the why s of subnuclear physics , erice , 1977 , edited by a.  zichichi ( plenum , new york , 1979 ) , p.94 a. h.  mueller , nucl . phys . * b250 * , 327 ( 1985 ) ; phys . lett . * b308 * , 355 ( 1993 ) . v. i.  zakharov , nucl . b385 * , 452 ( 1992 ) . m.  beneke , phys . rep . * 317 * , 1 ( 1999 ) . g.  grunberg , phys . lett . * b95 * , 70 ( 1980 ) ; * b110 * , 501(e ) , ( 1982 ) ; phys . rev . * d29 * , 2315 ( 1984 ) . p. m.  stevenson , phys . rev . * d23 * , 2916 ( 1981 ) . s. j.  brodsky , g. p.  lepage and p. b.  mackenzie , phys . * d28 * , 228 ( 1983 ) ; g. p.  lepage and p. b.  mackenzie , ibid . * 48 * , 2250 ( 1993 ) . g.  grunberg and a. l.  kataev , phys . lett * b279 * , 352 ( 1992 ) . s. j.  brodsky and j.  rathsman , hep - ph/9906339 . h. j.  lu , and c. a. r. sa de  melo , phys . lett * b273 * , 260 ( 1991 ) ; * b285 * , 399(e ) ( 1992 ) ; h. j.  lu , phys . rev . * d45 * , 1217 ( 1992 ) . c. p.  lepage and p. b.  mackenzie phys . * d48 * , 2250 ( 1993 ) . s. j.  brodsky and h. j.  lu , phys . rev . * d51 * , 3652 ( 1995 ) . j.  rathsman , phys . rev . * d54 * , 3420 ( 1996 ) . m.  neubert , phys . rev . * d51 * , 5924 ( 1995 ) . m.  beneke and v. m.  braun , phys . lett . * b348 * , 513 ( 1995 ) ; p. ball , m.  beneke and v. m.  braun , nucl . phys . * b452 * , 424 ( 1995 ) ; m.  beneke , nucl . phys . * b405 * , 563 ( 1993 ) . c. n.  lovett - turner and c. j. maxwell , nucl . phys . b432 * , 147 ( 1994 ) . s. j.  brodsky , j.  ellis , e.  gardi , m.  karliner , and m. a.  samuel , phys . rev . * d56 * , 6980 ( 1997 ) . s. j.  brodsky , g. t.  gabadadze , a. l.  kataev , and h. j.  lu , phys . b372 * , 133 ( 1996 ) . s. j.  brodsky , m.  melles , and j.  rathsman , phys . rev . * d60 * , 096006 ( 1999 ) . s. j.  brodsky , m. s.  gill , m.  melles , and j.  rathsman , phys . rev . * d58 * , 116006 ( 1998 ) . s. v.  mikhailov , jhep 0706:009 ( 2007 ) . v. a.  matveev , r. m.  muradyan , and a. n.  tavkhelidze , lett . nuovo cimento 7 , 719 ( 1973 ) . s. j.  brodsky and g. r.  farrar , phys . lett . * 31 * , 1153 ( 1973 ) . j. f.  gunion , s. j.  brodsky , r.  blankenbecler , phys.rev . * d6 * , 2652 ( 1972 ) d. w. sivers , s. j. brodsky , r. blankenbecler , phys . rept . * 23 * , 1 ( 1976 ) . v. a.  matveev , l. a.  slepchenko , and a. n.  tavkhelidze , phys . lett . * b100 * , 75 ( 1981 ) . g. l.  lepage and s. j.  brodsky , phys . rev . * d22 * , 2157 ( 1980 ) . v. l.  chernyak and a. r.  zhitnitsky , phys . rep . * 112 * , 173 ( 1984 ) . franois  arleo and david  denterria , phys * d78 * : 094004 ( 2008 ) . j. a.  bagger and j. f.  gunion , phys * d25 * , 2287 ( 1982 ) . v. n.  baier and a. grozin , phys . lett . * b96 * , 181 ( 1980 ) ; s.  gupta , phys . rev . * d24 * , 1169 ( 1981 ) . a. i.  ahmadov , i.  boztosun , r. kh . muradov , a.  soylu and e. a.  dadashov , int . j. mod e15 * , 1209 ( 2006 ) . m.  maul , e.  stein , a.  schfer , l.  mankiewich , phys . lett.*b401*,100 ( 1997 ) . y. l.  dokshitzer , v. a.  khoze and s. i.  troyan , phys . rev . * d53 * , 89 ( 1996 ) . a. l.  kataev , mod . * a20 * , 2007 ( 2005 ) . h.  contopanagos and g.  sterman , nucl . phys . * b419 * , 77 ( 1994 ) . w.  greiner , s.  schramm and e.  stein , quantum chromodynamics , 2nd edn.(berlin , springer , 2002 ) , pp.551 . d. v.  shirkov and i. l.  solovtsov , phys . lett . * 79 * , 1209 ( 1997 ) . star collaboration ( j.adams et.al ) , phys b637 * , 161 ( 2006 ) . agaev , phys . lett . * b360 * , 117 ( 1995 ) ; * b369 * , 379(e ) ( 1996 ) . a. i.  ahmadov , coskun  aydin , sh . m.  nagiyev , yilmaz a.  hakan , and e. a.  dadashov , phys * d80 * , 016003 ( 2009 ) . a. i.  ahmadov , coskun  aydin , e. a.  dadashov and sh . m.  nagiyev , phys . rev . * d81 * , 054016 ( 2010 ) . agaev , modern physics letters * a13 * , 2637 ( 1998 ) . franois  arleo , stanley j.  brodsky , dae suang hwang , and anne m.sickles , phys * 105 * , 062002 ( 2010 ) . s. m.  berman , j. d.  bjorken and j. b.  kogut , phys . rev . * d18 * , 900 ( 1978 ) . v. l.  chernyak and a. r.  zhitnitsky , nucl . phys . * b201 * , 492 ( 1982 ) ; v. l.  chernyak , a. r.  zhitnitsky and i. r.  zhitnitsky , nucl . phys . * b204 * , 477 ( 1982 ) . a.  schmedding and o.  yakovlev , phys d62 * , 116002 ( 2000 ) . v. m.  braun and i. f.  filyanov , z. phys , * c44 * , 157 ( 1989 ) . a. p.  bakulev , s. v.  mikhailov and n. g.  stefanis , phys . lett , * b587 * , 91 ( 2004 ) . s. v.  mikhailov and a. v.  radyushkin , phys . rev . * d45 * , 1754 ( 1992 ) . bakulev , s. v.  mikhailov , z. phys . * c68 * , 451 ( 1995 ) . a. e.  dorokhov , jetp lett . * 77 * , 63 ( 2003 ) . j. f.  owens , rev . phys . * 59 * , 465 ( 1987 ) . j. zinn - justin , phys . rept . * 70 * , 109 ( 1981 ) . a.  erdelyi , higher transcendental functions ( mcgrow - hill book company , new york , 1953 ) , vol.2 . mrst2003c.f can be obtained from http://durpdg.dur.ac.uk/hepdata/pdf.html . see also , a.  d.  martin , r.  g.  roberts , w.  j.  stirling and r. s.thorne , hep - ph/0307262 ; r. s.thorne , hep - ph/0309343 .
Abstract:

Answer: in this article , we investigate the next - to - leading order contribution of the higher - twist feynman diagrams to the large-@xmath0 inclusive pion production cross section in proton - proton collisions and present the general formulae for the higher - twist differential cross sections in the case of the running coupling and frozen coupling approaches . we compared the resummed next - to - leading order higher - twist cross sections with the ones obtained in the framework of the frozen coupling approach and leading - twist cross section . the structure of infrared renormalon singularities of the higher twist subprocess cross section and it s resummed expression ( the borel sum ) are found . it is shown that the resummed result depends on the choice of the meson wave functions used in the calculations . we discuss the phenomenological consequences of possible higher - twist contributions to the meson production in proton - proton collisions in next - to - leading order at rhic .