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: aggregation processes in which small objects merge irreversibly to form larger clusters are ubiquitous in nature @xcite . for example , aggregation underlies the evolution of planetary systems in astrophysics @xcite , cloud formation and dust accumulation in atmospheric sciences @xcite , as well as polymer and gel formation in chemical physics @xcite . aggregation also plays a central role in the theory of percolation @xcite , fractal formation @xcite , and network growth @xcite . often , aggregation is driven by a constant injection of mass and consequently , the total mass grows indefinitely with time @xcite . this is the case in chemical kinetics , and in particular polymerization where the aggregation rate depends on the number of clusters . here , due to aggregation , the total number of clusters typically decreases with time or saturates at a finite value . in such aggregation processes , the system may undergo gelation : a giant cluster develops and eventually , it contains all of the mass in the system @xcite .    in this study , we are interested in complementary aggregation processes that describe the growth of random structures such as random trees and random graphs , relevant in computer and information science @xcite . here , the system consists of an ensemble of clusters where both the number of clusters and the total mass grow with time . for example , in the internet , clusters are autonomous systems ; injection models the creation of new autonomous systems and aggregation describes merger of different autonomous systems @xcite .    we investigate aggregation processes with mass injection where the total merger rate does not depend on the total number of clusters . our main finding is that in such situations there is condensation rather than gelation . the system develops condensates that contain a finite fraction of the mass . these macroscopic condensates co - exist with microscopic clusters that contain the rest of the mass in the system .    in our formulation , mass is injected at a constant rate and clusters merge at a constant total rate . we address the three classic kernels @xmath2 for aggregation between clusters if size @xmath3 and @xmath4 , respectively : the constant rate @xmath5 , the sum rate @xmath14 , and the product rate @xmath15 .    for the constant aggregation rate , the cluster size distribution decays exponentially with the cluster size . for both the sum and the product rates , where the aggregation rate grows with the aggregate size , the system undergoes a phase transition as a function of the injection rate . when the injection rate is smaller than some critical value , the system is in a condensate phase . a finite fraction of the total mass constitutes a macroscopic condensate , but the remaining fraction of mass is in the form of microscopic clusters . the condensate and the finite clusters coexist . when the injection rate is larger than the critical rate , there are only finite clusters .    in the cluster phase , the distribution of cluster size generally decays as a power - law with the cluster size . however , different behaviors emerge in the condensate phase . for the product rate , the size distribution falls - off exponentially at large size . for the sum rate , however , the behavior is always power - law . interestingly , in the latter case , the decay exponent is not monotonic . it decreases monotonically with the injection rate in the condensate phase but increases monotonically in the cluster phase . the rest of this paper is organized as follows . we introduce the model , describe some of its basic features , and formulate the master equation approach in sec . the constant rate , the sum rate , and the product rate are analyzed in sections iii , iv , v , respectively . generally , our focus is the size of the condensate and the tail of the cluster size distribution . we conclude in sec . technical derivations of two results in sections iv and v are presented in appendices  [ log ] and [ giant ] . in our model , there are two independent and competing processes : mass injection and merger of clusters . in the first process , monodisperse elemental clusters are added to the system . this injection process occurs at a constant rate . in the second process , two clusters are merged . the mass of the resulting cluster is equal to the sum of the two original masses . the total merger rate is constant as well , and since both processes occur with constant rates , we may set one of them to unity without loss of generality . we therefore set the mass injection rate to @xmath0 and the merger rate to one . also , since the injection is monodisperse , we set this injection size as the mass unit . initially , the system contains no clusters . the total mass and the total number of clusters follow directly from these definitions . each injection event increases the number of clusters by one and similarly , each merger event decreases the number of clusters by one . thus , the average total number of clusters , @xmath16 , satisfies @xmath17 and consequently , there is simple linear growth @xmath18 we restrict our attention to situations where the number of clusters grows with time , @xmath19 . merger events conserve the total mass , and hence , the mass changes only through injection events . with each injection event , unit - size mass is added to the system , and thus , the average total mass , @xmath20 , obeys . consequently , the total mass also grows linearly with time , @xmath21    we investigate the cluster size distribution . let @xmath22 be the average number of clusters of size @xmath23 at time @xmath24 . this quantity satisfies the master equation @xmath25 the initial condition is @xmath26 . the first term accounts for mass injection and the last two terms account for merger . the kernel @xmath2 is defined as the aggregation rate between two clusters with size @xmath3 and @xmath4 . the total merger rate is constant , thereby implying the following constraint on the aggregation rate @xmath27 the total density is of course @xmath28 . summing the master equation ( [ ck - eq ] ) , and using the constraint ( [ norm ] ) , we confirm the linear growth of the total density ( [ nt ] ) . similarly , the total mass is @xmath29 . only the first term in the master equation affects the evolution of the total mass , and summing the rate equations , we recover ( [ mt ] ) . we analyze the three classic aggregation rates : the constant aggregation rate @xmath30 , the sum aggregation rate @xmath6 , and the product aggregation rate @xmath31 @xcite , each representing a distinct growing random structure . generally , random structures such as random trees and random graphs are made of nodes interconnected by links . mass injection represents addition of isolated nodes , and merger represents addition of a link between two nodes @xcite . in the constant rate process , two random structures are picked at random and an added link connects the two . in the product rate process , two nodes , picked at random , are connected by an added link . finally , the sum rate is a hybrid of the constant and the product cases : the link connects a randomly selected node and a randomly selected structure . we note that the constant aggregation rate models an ensemble of random growing trees ( no cycles are formed ) , while the sum and the product rates model an ensemble of random growing graphs . for the constant aggregation rate , all pairs of clusters merge at the same rate , irrespective of their size . this is the simplest and the most widely used aggregation process with examples including fractal aggregates @xcite , domain growth @xcite , and random trees @xcite . we consider the following size - independent aggregation rate @xmath32 this constant satisfies the normalization ( [ norm ] ) . it decreases as the number of clusters increases so that the overall merger rate does not change with time .    substituting this constant rate into the master equation ( [ ck - eq ] ) , the cluster size density satisfies @xmath33 the linear growth of the total density ( [ nt ] ) and the total mass ( [ mt ] ) suggest that the cluster size distribution also grows linearly with time . indeed , the density of the smallest clusters obeys @xmath34 so that this quantity , too , grows linearly with time , @xmath35 . thus , we write the cluster size density as a product of the overall density @xmath16 and the time - independent cluster size distribution @xmath8 , @xmath36 the cluster size distribution is normalized , @xmath37 . this form is consistent with the initial condition , and it satisfies ( [ ck - eq - const ] ) when the cluster size distribution obeys the recursion equation @xmath38    we utilize the generating function technique @xcite to solve this equation . let @xmath39 be the generating function @xmath40 normalization implies @xmath41 . multiplying eq . ( [ ck - eq - const ] ) by @xmath42 and summing over @xmath23 , the generating function satisfies the quadratic equation @xmath43 . the solution is therefore @xmath44.\ ] ] we can confirm that @xmath45 , in agreement with the above expression for @xmath46 . writing the generating function as a power - series , the cluster size distribution is obtained explicitly , @xmath47^k\ ] ] where @xmath48 is the gamma function . using the asymptotic property @xmath49 as @xmath50 , we find @xmath51 with the constants @xmath52 and @xmath53 . therefore , the cluster size distribution decays exponentially , but there is an algebraic correction . as the injection rate approaches the merger rate , @xmath54 , the cluster size distribution becomes algebraic @xmath55 . this limiting case coincides with the well - known behavior for time - independent constant aggregation rates @xcite . aggregation rates proportional to the sum of the cluster sizes are relevant in polymerization @xcite , coagulation under shear - flows @xcite , and random graphs @xcite . subject to the constraint ( [ norm ] ) , the sum aggregation rate is @xmath56    in this case , the master equation ( [ ck - eq ] ) becomes @xmath57    we again seek a solution of the form ( [ ck - def ] ) . the cluster size distribution remains normalized , @xmath37 , and it obeys the following recursion relation @xmath58 with the constant @xmath59 this recursion relation can be manually solved to find @xmath60 , @xmath61 , etc . this procedure can be formally related to integer partitions following the solution procedure of ref . such a solution is useful only when the cluster size distribution decays sharply . although it is difficult to obtain an explicit analytic solution as in ( [ ck - sol - const ] ) , it is still possible to obtain many of the interesting properties of cluster size distribution from asymptotic analysis of the generating function . the generating function ( [ fz - def ] ) obeys the nonlinear ordinary differential equation @xmath62 derivatives of the generating function at @xmath63 are related to the moments of the cluster size distribution . for example , the average cluster size follows from the first derivative , @xmath64 .    differentiating ( [ fz - eq - sum ] ) and then substituting @xmath63 yields a quadratic equation for the average cluster size @xmath65 naively , one expects that the average cluster size is the ratio between the total mass and the total number of clusters , @xmath66 . indeed , this is a solution of ( [ kav - eq - sum ] ) . however , there is another solution @xmath67 . this solution is not physical when @xmath13 since the product @xmath68 can not exceed the total mass in the system . each solution is relevant in the appropriate range of the parameter @xmath0 , so that the full solution is @xmath69 this assertion is supported by analysis below . the total mass contained by the clusters , @xmath70 , is of course proportional to the average cluster size , . therefore , when @xmath11 , finite clusters contain only a fraction of the mass when @xmath11 , but they contain all of the mass when @xmath13 .    this behavior can be reconciled with mass conservation only if there is a condensate of mass @xmath71 that contains the remaining fraction @xmath72 of the mass in the system as follows @xmath73 thus the system undergoes a phase transition . when @xmath74 , there is a condensate that contains a finite fraction of the mass . this condensate coexists with the finite clusters that contain the rest of the mass . the condensate contains nearly all of the mass in the limit @xmath54 : when injection is very slow , the condensate `` preys '' on newly added mass . as the transition point is approached , the condensate mass vanishes , @xmath75 as @xmath76 . when @xmath13 , the system contains only ordinary clusters . the tail of the cluster size distribution can be evaluated from the @xmath77 behavior of the generating function . in general , @xmath39 may contain both a regular component and a singular component , @xmath78 . the regular component is a power series in @xmath79 . let us assume a singular behavior with the leading behavior @xmath80 as @xmath77 , @xmath81 in the limit @xmath77 , this form satisfies the governing equation ( [ fz - eq - sum ] ) when @xmath82 the algebraic form of the generating function implies an algebraic form for the tail of the cluster size distribution , @xmath83 as @xmath84 . the exponent is found by substituting ( [ kav - sol - sum ] ) into the relation ( [ gamma - eq - sum ] ) . therefore , there are two regimes of behavior @xmath85 we note two remarkable features . first , the cluster size density is algebraic both in the condensate phase and in the cluster phase . second , the characteristic exponent is a non - monotonic function of the injection rate : it decreases monotonically with @xmath0 in the condensate phase and it increases monotonically in the cluster phase ( fig . [ fig - gamma - sum ] ) .     versus @xmath0 for the sum rate.,scaledwidth=40.0% ] the exponent @xmath86 is minimal , @xmath87 , at the phase transition point , @xmath88 . for @xmath89 , mass conservation would be violated because the sum @xmath90 is divergent . the restriction @xmath91 justifies our previous choice of the nontrivial solution @xmath92 in eq . ( [ kav - sol - sum ] ) . the behavior at the phase transition point requires a special treatment . we find that the cluster size density decays slightly faster than @xmath93 , namely @xmath94 the derivation of this result is detailed in appendix a. with this logarithmic correction , the sum @xmath90 or the total mass in finite clusters , remain finite .    in ordinary gelation ( or percolation ) , the cluster size distribution is algebraic _ only _ at the critical point ; away from criticality , it has exponential tails . in the present case , the cluster size distribution exhibits a strikingly different behavior  it is algebraic everywhere while at the critical point there is a logarithmic correction .    the size of the condensate can be obtained directly by focusing on its dynamics . let us assume that there is a giant cluster in the system with mass @xmath71 . its growth rate is @xmath95 and substituting the aggregation rate ( [ kij - sum ] ) , we arrive at @xmath96 in the second step , we made the approximation as the condensate is much larger than the rest of the clusters . when the condensate contains a finite fraction of the mass , @xmath97 , then from ( [ m*-eq - sum ] ) , we find @xmath98 . solving this equation , we recover the condensate mass ( [ m*-sol - sum ] ) . this approach can also be used to derive the size of the largest cluster in the cluster phase . from ( [ m*-sol - sum ] ) , the quantity @xmath71 is negligible compared with @xmath99 , and therefore , the governing equation ( [ m*-eq - sum ] ) reduces to @xmath100 so that @xmath101 with @xmath12 when @xmath13 . therefore , the size of the largest cluster grows algebraically with time . this result can be alternatively derived using the algebraic behavior ( [ ck - tail - sum ] ) and the extreme statistics criterion @xmath102 . the product aggregation rate models polymerization and gelation @xcite , as well as random graphs @xcite , and in our case it has the following form @xmath103    the explicit rate equation for the product aggregation rate is @xmath104 since the total mass appears in the denominator in the rate equation , we use a different normalization @xmath105 the average number of clusters of size @xmath23 still grows linearly with time . the transformation ( [ ck - def - prod ] ) reduces the master equation ( [ ck - eq - prod ] ) to the nonlinear recursion equation @xmath106 given of the structure of this equation , it is convenient to use a different definition of the generating function @xmath107 the generating equation satisfies the very same equation ( [ fz - eq - sum ] ) with @xmath108 although the governing equation is the same , the boundary condition is different . since the product aggregation rate grows with the cluster size , we expect that again there are two phases : a condensate phase and a cluster phase . the total mass contained in finite clusters is given by @xmath109 . therefore , in the condensate phase @xmath110 while in the cluster phase @xmath41 . this change in the boundary condition results in a drastically different behavior . we first discuss the cluster phase where @xmath41 , and consequently , the analysis is a straightforward generalization of the above . the first derivative at @xmath63 again satisfies @xmath111 ^ 2-f'(1)+1=0 $ ] . at the critical point @xmath112 and therefore @xmath113 . solving this quadratic equation , the first derivative is @xmath114 the first derivative is now the ratio between the second and the first moments of the cluster size distribution , @xmath115 . the tail behavior follows from the singular component of the generating function as in ( [ fz - def - sum ] ) @xmath116 this again implies the power - law decay ( [ ck - tail - sum ] ) for the cluster size distribution . we note that the unit shift in the exponent is due to the different definition of the generating function ( [ fz - def - prod ] ) . from the governing equation ( [ fz - eq - sum ] ) , the exponent @xmath86 satisfies @xmath117 , and the decay exponent is @xmath118 similar to ( [ gamma - sol - sum ] ) , the characteristic exponent grows linearly with the injection rate , @xmath119 , at large injection rates ( figure [ fig - gamma - prod ] ) . the minimum value , @xmath120 , is achieved at the phase transition point . a more careful analysis is again required at the phase transition point ; we find ( the analysis is essentially the same as in the case of the sum rate , appendix  [ log ] ) that @xmath121 as @xmath122 . versus @xmath0 for the product rate.,scaledwidth=40.0% ]    in the condensate phase , we are able to obtain the mass of the condensate only in the vicinity of the phase transition , @xmath123,\ ] ] as @xmath124 . the derivation of this result resembles that of refs . @xcite ; it is detailed in appendix  [ giant ] . in contrast with the sum rate , the phase transition is now very gentle : all derivatives of the condensate mass vanish at the transition point @xmath125 @xcite . in practice , it may be difficult to locate the phase transition point ( figure [ fig - m * ] ) . we comment that similar behavior was recently found in several models of growing networks @xcite . we performed monte carlo simulations to obtain the condensate mass shown in figure [ fig - m * ] . in the simulations there are two elemental steps : injection with probability @xmath126 and aggregation with probability @xmath127 . in an injection step , a cluster with unit mass is added into the system . in an aggregation step , two clusters , picked with probability proportional to the product of the two masses , are merged . the simulations results represent a single run with a total mass of @xmath128 .    as a signature of the phase transition , the quantity @xmath129 has a discontinuity at the phase transition point . in the condensate phase , this quantity is obtained , using the fact that @xmath110 , directly from the governing equation ( [ fz - eq - sum ] ) , @xmath130 . in the cluster phase it is given by ( [ fprime - sol - prod ] ) . therefore , the ratio between the second and the first moments has a jump at @xmath113 , @xmath131 a similar jump occurs in the berezinskii - kosterlitz - thouless phase transition @xcite .     as a function of @xmath0 obtained from monte carlo simulations of the aggregation process.,scaledwidth=45.0% ]    finally , we study the tail of the cluster size distribution . for this purpose , it is convenient to make the transformation @xmath132 so that the generating function ( [ fz - def - prod ] ) is redefined @xmath133 substituting this definition into the recursion equation ( [ ck - eq - prod ] ) , the generating function obeys @xmath134 the denominator suggests that @xmath135 has a singularity at @xmath136 with @xmath137 . from the definition ( [ fw - def ] ) , this implies the exponential decay @xmath138 with @xmath139 . since we can not solve for the generating function , we can not locate this singularity explicitly . nevertheless , one can still deduce the behavior near this singularity using asymptotic analysis . again , we assume that the generating function has a regular component and a singular component near @xmath140 , @xmath141 substituting this form into ( [ fw - eq - prod ] ) , and equating powers of @xmath142 on both sides of the equation we obtain . using the definition ( [ fw - def ] ) , the leading behavior of the singular component @xmath143 implies an algebraic correction to the leading exponential decay , @xmath144 as @xmath84 . we conclude that the tail of the cluster size distribution decays as follows @xmath145    the characteristic size can be related to the condensate mass in the vicinity of the phase transition as the characteristic size is expected to diverge in this limit , @xmath146 . consequently , the singularity is located close to the origin , @xmath147 . let us estimate the behavior of the generating function near the origin . from the definition ( [ fw - def ] ) , we have @xmath148 also , from the governing equation we have @xmath149 in the limit @xmath150 as in ( [ jump ] ) . therefore , @xmath151 as @xmath152 and from the condition @xmath137 , the location of the singularity is @xmath153 . therefore , the characteristic size grows as follows @xmath154 in the vicinity of the phase transition point , @xmath150 . the tail behavior coincides with the critical behavior at sufficiently small sizes , @xmath155 for @xmath156 , and exponential decay , @xmath138 , occurs beyond that scale , @xmath157 . the two behaviors should of course match at @xmath158 : this implies the proportionality constant in ( [ ck - tail - prod ] ) , @xmath159 . in conclusion , we studied aggregation with constant injection of mass . in this process , the total number of clusters grows with time . for aggregation rates growing as the sum or the product of the cluster sizes , there are two phases : a condensate phase and a cluster phase . in the condensate phase , a condensate containing a finite fraction of the mass coexists with finite clusters , while in the cluster phase there are only finite clusters .    for the sum rate , the mass of the condensate is a linear function of the injection rate . also , the cluster size distribution decays algebraically in both phases and interestingly , the decay exponent is a non - monotonic function of the injection rate . for the product rate the condensate mass is extremely small in the vicinity of the phase transition point and consequently , the phase transition is very gentle . in this case , the tail of the cluster size distribution is exponential in the condensate phase but algebraic in the cluster phase . we comment that there are two frameworks for describing aggregation processes : the flory approach that allows interaction between giant and finite clusters @xcite and the stockmayer approach that allows for interactions between finite clusters only @xcite . we used the more challenging former approach as it is the appropriate approach for modeling growing random structures @xcite . the various aggregation processes correspond to different random growing structures , but our study focused only on the size of these structures . we note that this theoretical framework can be generalized to also study structural properties such as paths and cycles @xcite . there are a number of possible extensions of this work . we focused on the three classic aggregation rates where the generating function obeys closed equations . this framework does not allow derivation of the necessary conditions for the emergence of a condensate as a function of the aggregation rate . based on the sensitive algebraic behavior in both of the phases , we speculate that the sum rate may be the marginal case for condensation . * acknowledgments . * we thank the isaac newton institute for mathematical sciences ( cambridge , england ) and the max planck institute for physics of complex systems ( dresden , germany ) , where this research was largely performed , for their hospitality . we acknowledge financial support from doe grant de - ac52 - 06na25396 and nsf grant che-0532969 .    99    m.  von smoluchowski , z. phys . * 17 * , 557 ( 1916 ) . f.  leyvraz , phys .  rep . * 383 * , 95 ( 2003 ) . s.  chandrasekhar , rev . * 15 * , 1 ( 1943 ) . h.  pruppacher and j.  klett , _ microphysics of clouds and precipitations _ ( kluwer , dordrecht , 1998 ) .    r.  l.  drake , in _ topics of current aerosol research _ , edited by g.  m.  hidy and j.  r.  brock ( pergamon , new york , 1972 ) . s.  k.  frielander , _ smoke , dust and haze : fundamentals of aerosol behavior _ ( wiley , new york , 1977 ) . j.  h.  seinfeld , _ atmospheric chemistry and physics of air pollution _ ( wiley , boston , 1986 ) . p.  j.  flory , j. amer . * 63 * , 3083 ( 1941 ) . w.  h.  stockmayer , j. chem . phys . * 11 * , 45 ( 1943 ) . j.  flory , _ principles of polymer chemistry _ ( cornell university press , ithaca , 1953 ) . d.  stauffer and a.  aharony , _ introduction to percolation theory _ ( taylor & francis , london , 1992 ) .    t.  a.  witten and l.  m.  sander , phys . rev . lett . * 47 * , 1400 ( 1981 ) . r.  albert and b.  a.  barabasi , rev . mod . phys . * 74 * , 47 ( 2002 ) . s.  n.  dorogovtsev and j.  f.  f.  mendes , adv . in phys . * 51 * , 1079 ( 2002 ) .    g.  b.  field and w.  c.  saslaw , astrophys . j. * 142 * , 568 ( 1965 )    w.  h.  white , j. colloid interface sci . * 87 * , 204 ( 1982 ) . j.  g.  crump and j.  h.  seinfeld , j. colloid interface sci . * 90 * , 469 ( 1982 ) . e.  m.  hendriks and m.  h.  ernst , j. colloid interface sci . * 97 * , 196 ( 1984 ) . d.  aldous and b.  pittel , random struct . * 17 * , 79 ( 2000 ) . z.  cheng , s.  redner , and f.  leyvraz , phys . * 62 * , 2321 ( 1989 ) . r. m. ziff , e. m. hendriks , and m. h. ernst , phys . lett . * 49 * , 593 ( 1982 ) . g.  j. van dongen , j. phys . a * 20 * , 1889 ( 1987 ) . a.  a.  lushnikov , phys . lett . * 93 * , 198302 ( 2004 ) . b.  ballobs , _ random graphs _ ( academic press , london , 1985 ) . d.  j.  aldous , bernoulli * 5 * , 3 ( 1999 ) . s.  janson , t.  uczak , and a.  rucinski , _ random graphs _ ( john wiley & sons , new york , 2000 ) . e.  ben - naim and p.  l.  krapivsky , phys . e * 71 * , 026129 ( 2005 ) . m.  fayed , p.  l.  krapivsky , j.  byers , m.  crovella , s.  redner , and d.  finkel , computer commun .  rev . * 33:2 * , 41 ( 2003 ) . p.  meakin , _ fractals , scaling and growth far from equilibrium _ ( cambridge university press , new york 1998 ) . a.  j.  bray , adv . phys . * 43 * , 357 ( 1994 ) . h.  m.  mahmoud , _ evolution of random search trees _ ( john wiley , new york , 1992 ) . e.  ben - naim , p.  l.  krapivsky , and s.  m.  majumdar , phys . rev . e * 64 * , 035101 ( 2001 ) . r.  l.  graham , d.  e.  knuth , and o.  patashnik , _ concrete mathematics : a foundation for computer science _ ( reading , mass . : addison - wesley , 1989 ) . r.  m.  ziff , j.  stat . * 23 * , 241 ( 1980 ) . e.  ben - naim and p.  l.  krapivsky , phys . e * 73 * , 03119 ( 2006 ) . p.  l.  krapivsky and b.  derrida , physica a * 340 * , 714 ( 2004 ) . s.  n.  dorogovtsev , j.  f.  f.  mendes , and a.  n.  samukhin , phys . e * 64 * , 066110 ( 2001 ) . v.  l.  berezinskii , sov . jetp * 32 * , 493 ( 1970 ) . j.  m.  kosterlitz and d.  j.  thouless , j. phys . c * 6 * , 1181 ( 1973 ) . d.  s.  callaway , j.  e.  hopcroft , j.  m.  kleinberg , m.  e.  j.  newman , and s.  h.  strogatz , phys . e * 64 * , 041902 ( 2001 ) . l.  krapivsky and s.  redner , computer networks * 39 * , 261 ( 2002 ) . j.  kim , p.  l.  krapivsky , b.  kahng , and s.  redner , phys . e * 66 * , 055102 ( 2002 ) . m.  bauer and d.  bernard , j. stat .  phys . * 111 * , 703 ( 2003 ) .    c.  coulomb and m.  bauer , eur . j. b. * 35 * , 377 ( 2003 ) . e.  ben - naim and p.  l.  krapivsky , j. phys . cond . mat . * 17 * , s4249 ( 2005 ) . at the phase transition point @xmath88 we have @xmath87 and therefore the leading singular term becomes regular . this suggests to use instead where @xmath160 vanishes slower than any power of @xmath161 as @xmath77 . thus at the phase transition point we employ the following expansion @xmath162 substituting ( [ exp - mod ] ) into ( [ fz - eq - sum ] ) yields the differential equation @xmath163 whose ( implicit ) solution is @xmath164 in the limit @xmath77 , the integration constant is negligible compared with the logarithmic term in the limit @xmath77 and consequently , @xmath165 . indeed , @xmath166 vanishes slower than any power of @xmath161 as @xmath77 . thus @xmath167 inverting this expansion leads to eq . ( [ ns - crit ] ) @xcite . the mass of the condensate follows from the behavior of the generating function at @xmath63 , @xmath168 . to analyze the behavior near this region , we make the transformations      with these transformations , the equation for the generating function ( [ fz - eq - sum ] ) is transformed into the following first - order nonlinear differential equation @xmath170 in writing this equation , we kept only the leading order terms . writing , and integrating , we have @xmath171 where @xmath172 . the integration constant can be evaluated by taking the @xmath173 limit . using , the first two terms in the above equation approach @xmath174 in the limit @xmath173 . using , the last term approaches @xmath175 . hence @xmath176    since we are interested in the behavior near the phase transition point , we take the limit @xmath177 . in this limit , we can replace @xmath178 by @xmath179 and also , the quantity @xmath180 by @xmath181 . additionally , we may replace @xmath182 by @xmath183 . with these substitutions , ( [ integ ] ) becomes @xmath184    next , we evaluate the left - hand side precisely at the phase transition point , @xmath185 . the critical behavior is detailed in appendix a. substituting @xmath186 and @xmath187 into ( [ uz ] ) , gives @xmath188 substituting this into ( [ g - int - near ] ) we obtain the condensate mass in the vicinity of the phase transition ( [ m*-prod ] ) .
Abstract:

Answer: we investigate aggregation driven by mass injection . in this stochastic process , mass is added with constant rate @xmath0 and clusters merge at a constant total rate @xmath1 , so that both the total number of clusters and the total mass steadily grow with time . analytic results are presented for the three classic aggregation rates @xmath2 between clusters of size @xmath3 and @xmath4 . when @xmath5 , the cluster size distribution decays exponentially . when @xmath6 or @xmath7 , there are two phases : ( i ) a condensate phase with a condensate containing a finite fraction of the mass in the system as well as finite clusters , and ( ii ) a cluster phase with finite clusters only . for @xmath6 , the cluster size distribution , @xmath8 , has a power - law tail , @xmath9 in either phase . the exponent is a non - monotonic function of the injection rate : @xmath10 in the condensate phase , @xmath11 , and @xmath12 in the cluster phase , @xmath13 .


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: recent experimental @xcite and theoretical @xcite investigations have demonstrated that ( almost ) unipolar , high power electromagnetic field pulses are a useful new spectroscopic tool which is particularly well suited for investigating the dynamics of weakly bound rydberg electrons . in current experiments the duration @xmath0 of these half - cycle pulses ( hcps ) extends from the subpicosecond to the nanosecond regime . several features distinguish the interaction of rydberg electrons with this new type of electromagnetic radiation from the interaction with conventional laser light or microwave fields . first of all , and in contrast to the optical case , the hcp interacts with the rydberg electron at every point of its orbit . furthermore , due to its unipolar nature the influence of the hcp on the electron can often be described in relatively simple and intuitive terms , e.g. , with the help of the so - called sudden - impact approximation where the effect of the hcp is modeled as an instantaneous momentum change . finally , the energy transfer between an oscillating electromagnetic field and a rydberg electron is always governed by an approximate energy conservation whose energy uncertainty is typically small in comparison with the amount of energy transferred . as hcps are almost unipolar such an approximate energy conservation does not hold and arbitrary amounts of energies can be transferred to a charged particle from zero up to a maximum amount of the order of @xmath1 . these characteristic features give rise to interesting novel phenomena which have been explored in a number of recent investigations . so far these studies have concentrated on total and energy - resolved ionization probabilities @xcite . theoretical work in this context has focused on numerical approaches @xcite that often employ the sudden - impact approximation and on one - dimensional models @xcite which try to capture the most significant physical effects occuring along the direction of polarization of a linearly polarized hcp .    in this paper a multi - dimensional semiclassical description of excitation of a weakly bound rydberg electron by a half - cycle pulse is developed . it quantitatively connects the quantum evolution of the rydberg electron to its underlying classical dynamics . the presented theoretical treatment is nonperturbative and it is applicable to arbitrary spatial and temporal hcp shapes . in particular , the approach is not restricted to the sudden - impact approximation . within this theoretical framework measurable transition probabilities are represented as coherent sums of probability amplitudes which are associated with corresponding classical trajectories of the excited rydberg electron . this method is very accurate numerically especially for highly excited rydberg states with large principal quantum numbers and for hcps whose transferred energy is small in comparison with typical ground - state ionization energies . it is precisely this dynamical regime which usually causes severe numerical problems in fully quantum mechanical calculations due to the large spatial extension of rydberg states and the presence of the coulomb singularity close to the nucleus . on the basis of this theoretical treatment we demonstrate that angular distributions of ionization probabilities contain a wealth of information about the ionization dynamics of the rydberg electron which can not be obtained from energy - resolved ionization spectra alone . in the sudden - impact approximation , for example , in which the pulse duration of the exciting hcp is small in comparison with the classical kepler period of the rydberg electron energy- and angle - resolved ionization spectra are dominated by oscillatory structures and semiclassical catastrophes of the rainbow and glory type . these structures can be explained in a clear and intuitive way in terms of interferences between probability amplitudes of a few classical electronic trajectories . the study of the influence of finite pulse durations indicates that some parts of these oscillatory structures depend strongly on the duration of the exciting hcp . a main purpose of our subsequent discussion is thus the systematic exploration of these latter effects which can not be described appropriately by the sudden - impact approximation . first results of our studies have also been presented in ref . @xcite . this paper is organized as follows : in sec . ii a general ( multidimensional ) semiclassical description of excitation of a weakly bound rydberg electron by an hcp is developed . this general approach is valid for arbitrary spatial and temporal dependences of the hcp . starting from this general approach the simplifications arising in the sudden - impact approximation are discussed in detail . in sec . iii we first present a characterization of the various types of energy- and angle - resolved ionization spectra that can be obtained in the sudden - impact approximation . these spectra may be classified according to the appearing semiclassical scattering phenomena . subsequently , it is investigated how a finite duration of the exciting hcp influences the behavior of the spectra . the comparison with the corresponding results obtained in the sudden - impact approximation demonstrate which dynamical aspects are particularly sensitive to pulse duration effects . finally , in sec . iv a brief summary and conclusions are given . in the first part of this section a nonperturbative semiclassical framework is developed that describes the excitation of a weakly bound rydberg electron by an hcp and that is valid for arbitrary spatial and temporal pulse shapes . in the second part this treatment is specialized to the case of short hcps whose influence on a rydberg electron can be described within the sudden - impact approximation . hartree atomic units are used with @xmath2 ( @xmath3 and @xmath4 denote the electronic mass and charge ) . let us consider a typical process in which a weakly bound rydberg electron is excited by a half - cycle pulse . the dynamics of the electronic wavefunction @xmath5 is described by the time - dependent schrdinger equation @xmath6 with the hamiltonian @xmath7 ^ 2   + v({\bf{x } } ) . \label{schr}\ ] ] outside the core region the effective potential @xmath8 experienced by the rydberg electron is of the form @xmath9 whereas inside it is modified due to the presence of the other electrons and of the nucleus . the electromagnetic vector potential of the exciting hcp is denoted @xmath10 and fulfills the relations @xmath11 and @xmath12 . initially the rydberg electron is supposed to be prepared in an energy eigenstate @xmath13 with principal quantum number @xmath14 and angular quantum numbers @xmath15 and @xmath16 . well inside the classically allowed region and for low values of the angular momenta , i.e. @xmath17 , this state can be approximated semiclassically by @xcite @xmath18 \label{initial}\end{aligned}\ ] ] with the classical eikonal @xmath19 the local radial electronic momentum in the coulomb potential of the ionic core is given by @xmath20 with the initial energy @xmath21^{-1}.\ ] ] while the rapidly oscillating part of the wave function ( [ initial ] ) is determined by the eikonal @xmath22 its slowly varying amplitude is given by @xmath23 with the spherical harmonic @xmath24 . the spherical angles of @xmath25 are denoted @xmath26 and @xmath27 , respectively . equation ( [ initial ] ) is valid for electronic distances @xmath28 from the nucleus which are located well outside the core region . all electron correlation effects are localized inside this core region which has an extension of a few bohr radii . they can be taken into account within the framework of quantum defect theory @xcite . in the case of an inert ionic core which is considered here for the sake of simplicity these core effects can be described by a single quantum defect @xmath29 . this quantum defect is approximately energy independent sufficiently close to the ionization threshold . generalizations to more complicated situations in which electronic core excitations and channel couplings have to be taken into account are possible within the framework of multi - channel quantum defect theory @xcite but will not be considered in this work .    in eq . ( [ initial ] ) the initial state @xmath30 is represented by a sum of contributions each of which involves a slowly varying amplitude , i.e.  @xmath31 , and an exponential function with a large , imaginary - valued argument , i.e.  @xmath32 . within the semiclassical treatment one associates with each of these contributions a three - dimensional lagrangian manifold in phase space @xcite . the special form of the initial condition of eq . ( [ initial ] ) implies that this lagrangian manifold has two branches , namely @xmath33 as @xmath34 depends only on the radial electronic distance @xmath35 one finds @xmath36 as long as @xmath37 . semiclassically , the solution of the time - dependent schrdinger equation @xmath38 is determined by all solutions @xmath39 and @xmath40 of the classical equations of motion with hamiltonian @xmath41 ^ 2   + v({\bf{x } } ) \label{hamclass}\end{aligned}\ ] ] whose initial conditions are located in this lagrangian manifold @xmath42 . the probability amplitude @xmath43 is determined by all those classical trajectories @xmath44 , @xmath45 which reach point @xmath25 at time @xmath46 , i.e. @xcite @xmath47}\langle{\bf{x}}_0^{(j)}|\psi\rangle_{t=0}. \label{solution}\end{aligned}\ ] ] for each contributing trajectory @xmath48 this probability amplitude is determined by three characteristic classical properties , namely its classical action @xmath49 the determinant of its jacobi field @xcite @xmath50 and its morse index @xmath51 . thereby , @xmath52 denotes the classical lagrange function associated with the hamiltonian of eq . ( [ hamclass ] ) . the determinant of the jacobi field characterizes the stability properties of the contributing classical trajectory . the morse index @xmath53 is numerically equal to the numbers of zeroes of this determinant times their multiplicities @xcite .    from eq . ( [ solution ] ) one can evaluate the probability amplitude of measuring an ionized electron with final momentum @xmath54 . with the help of the stationary phase approximation one obtains the result @xmath55}\langle { \bf{x}}|\psi\rangle_t   \nonumber\\ & = & \sum_j p_j^{(cl ) }   e^{i[s_j({\bf{p}}^{(f ) } ) + w_j({\bf{x}}^{(j)}_0 ) - i\pi{\tilde \mu}_j/2]}. \label{semi}\end{aligned}\ ] ] this ionization amplitude is represented as a sum of all trajectories @xmath48 which start at positions @xmath56 with momenta @xmath57 and which assume the final momentum @xmath54 long after the interaction with the hcp , i.e. , at @xmath58 . in eq . ( [ semi ] ) the quantity @xmath59 denotes the classical contribution of trajectory @xmath48 to the angle- and energy - resolved ionization probability . the classical action @xmath60 is given by @xmath61 the phase contribution originating from the initial state is determined by @xmath62 \end{aligned}\ ] ] where one has to choose the plus or minus sign depending on whether the initial radial momentum of the rydberg electron is positive or negative . the maslov index @xmath63 is equal to the number of zeroes of the determinant @xmath64 times their multiplicities @xcite . in terms of the ionization amplitude of eq . ( [ semi ] ) the energy- and angle - resolved ionization probability is given by @xmath65 with @xmath66 and with @xmath67 and @xmath68 denoting the spherical angles of the final momentum @xmath54 . equations ( [ semi ] ) and ( [ diff ] ) are main results of this section . they yield a complete and numerically accurate ( compare with sec . iii a ) semiclassical description of the influence of an arbitrary hcp on a weakly bound rydberg electron that is prepared initially in an energy eigenstate . these equations are based on maslov s multidimensional generalization of the jwkb method as applied to the time - dependent hamiltonian of eq . ( [ schr ] ) @xcite . it should be mentioned that there exists also a wealth of alternative semiclassical methods to approach initial value problems . these methods have been pioneered by e.  heller @xcite and w.  h.  miller @xcite and have been generalized recently in various directions @xcite . these types of semiclassical approximations are mainly based on the dynamics of gaussian wave packets @xcite and so far they have been applied predominantly to problems with explicitely time - independent hamiltonians .    if the initial state of the rydberg electron is not an energy eigenstate of the form of eq . ( [ initial ] ) almost all of the considerations of this section still apply . it is only the classical lagrangian manifold associated with the initial state that has to be modified appropriately . in the case of a spatially localized electronic rydberg wave packet , for example , which can be described by a quantum state of the form @xmath69 with a large eikonal @xmath70 and a slowly varying amplitude @xmath71 the associated lagrangian manifold is given by @xmath72 contrary to the initial state considered in eq . ( [ initial ] ) this lagrangian manifold has only one branch so that the structure of the interferences appearing in eq . ( [ semi ] ) is expected to be changed significantly . finally , we would like to mention a useful scaling property of the classical hamiltonian ( [ hamclass ] ) . if spatial variations of the pulses are negligible the classical trajectories @xmath73 and @xmath74 induced by two hcps @xmath75 and @xmath76 are related by the scaling relation @xmath77 these scaling properties are exploited in sec . iii a where they simplify the classification of the types of spectra obtainable in the sudden - ionization approximation . if the pulse duration @xmath0 of the exciting hcp is small in comparison with the classical kepler period @xmath78 of the initially prepared rydberg electron , the determination of @xmath5 can be simplified considerably . neglecting for the sake of clarity any spatial dependence the vector potential of the applied hcp can be approximated by @xmath79 with the characteristic momentum @xmath80 thereby , @xmath81 denotes the hcp field strength . thus the solution @xmath5 of the time - dependent schrdinger equation with hamiltonian ( [ schr ] ) is approximately given by @xmath82 . within this sudden - impact approximation the rydberg electron experiences a sudden change of its momentum by the amount @xmath83 and subsequently evolves under the influence of the atomic hamiltonian @xmath84 without being further affected by the hcp . in this case ionization probabilities can also be evaluated quantum mechanically by means of a partial wave analysis .    in the case of a hydrogen atom , for example , which is prepared initially in the state @xmath85 the angle- and energy - resolved ionization probability is given by the partial wave expansion @xmath86 thereby @xmath87 and @xmath88 are the legendre polynomial and the spherical bessel function of order @xmath89 , respectively . the regular , energy normalized radial coulomb wave function is denoted @xmath90 and @xmath91 is the coulomb phase shift . the radial wave function of the initially prepared rydberg state @xmath85 is denoted @xmath92 . in eq . ( [ partial ] ) it has been assumed that the hcp is linearly polarized and that the spherical angles @xmath93 of the final momentum @xmath54 are measured with respect to the direction of polarization . alternatively , @xmath5 can also be determined semiclassically with the approach developed in sec . ii a. for this purpose one starts from the electronic state which is modified by the sudden momentum change due to the hcp . immediately after the application of the hcp this state is given by @xmath94 of eq . ( [ initial ] ) , the classical lagrangian manifold which is associated with this state is given by @xmath95 in analogy to eq . ( [ solution ] ) the probability amplitude @xmath43 is determined by all solutions of the classical equations of motion with hamiltonian @xmath96 whose initial conditions are located in this lagrangian manifold . as this hamiltonian is explicitly time independent the energy of these trajectories is conserved during their time evolution . this implies that outside the core region where @xmath97 the initial positions @xmath98 on this lagrangian manifold which yield ionizing trajectories with final energy @xmath99 are determined by energy conservation , i.e. @xmath100 ^ 2 - \frac{1}{|{\bf x}_0 | } & = & \epsilon^{(f)}.\end{aligned}\ ] ] inserting eq . ( [ eikonal ] ) and assuming a linearly polarized hcp with @xmath101 , @xmath102 , yields @xmath103 with @xmath104 denoting the polar angle of the initial position @xmath105 . equation ( [ energy ] ) determines the initial positions @xmath105 of the classical trajectories which ionize with asymptotic energy @xmath99 . according to eq . ( [ energy ] ) there is one value of @xmath106 for each final energy @xmath99 and initial angle @xmath104 . at the critical final energy @xmath107 the local radial momentum of eq . ( [ momentum ] ) vanishes and the primitive semiclassical approximation for the initial state @xmath13 is no longer applicable . in this special case a proper semiclassical description of the initial state can be obtained with the help of uniform semiclassical methods @xcite . this peculiar property suggests that one can at least distinguish two different energy regimes in the excitation of a rydberg electron by a hcp , namely a region of low final energies for which @xmath108 and a high energy regime in which @xmath109 . for a fixed value of the final ionization energy @xmath99 the initial positions of the classical trajectories which ionize into a particular final angle @xmath67 have to be determined by the solutions of the classical equations of motion with hamiltonian ( [ hacl ] ) . in the special case of a hydrogen atom these relations can be evaluated analytically @xcite . a comprehensive discussion of energy- and angle - resolved ionization probabilities which is based on the semiclassical results of this section is presented in sec .  iii . in this section energy- and angle - resolved ionization probabilities are discussed on the basis of the semiclassical approach developed in sec . ii . in the first part characteristic predictions of the sudden - impact approximation are discussed in the semiclassical limit . the second part discusses modifications of these ionization spectra that arise from a finite duration of the exciting hcp . this part gives a general overview of the characteristic properties of angle- and energy - resolved ionization spectra in the sudden - impact approximation for the case @xmath110 . as it will be shown below these spectra can be described very accurately with the help of semiclassical methods . thus the characteristic qualitative aspects of their behavior may be understood by studying the underlying classical dynamics . the classical trajectories pertaining to the parameter sets @xmath111 and @xmath112 , @xmath113 , are related by the scaling transformation of eq . ( [ scaling ] ) . therefore it is sufficient to examine the classical dynamics as a function of the parameters @xmath114 and @xmath115 , for example . the result of a corresponding study is summarized in fig . [ fig1 ] . in this diagram various regions in the @xmath116-plane are identified that give rise to different types of behavior of the ionization spectra . region ( i ) pertains to parameters with @xmath117 . according to eq . ( [ energy ] ) , in this case classical trajectories start with negative ( positive ) radial momenta @xmath118 from the positive ( negative ) @xmath119-half plane . in order to illustrate the typical behavior of the resulting spectra we consider as an example the parameter values @xmath120 , @xmath121 . figure [ fig2](a ) depicts the initial positions of the classical trajectories in the @xmath122-plane and their respective final angles for the initial principal quantum number @xmath123 ( which implies @xmath124 mev ) . to show the distribution of the scattering angles clearly , in fig . [ fig2](b ) the  deflection function " is displayed , i.e. , the dependence of the final angle @xmath125 on the polar angle @xmath126 of the corresponding initial position . from figs . [ fig2](a ) and ( b ) it is apparent that three different initial positions pertain to each final emission angle @xmath125 . the relative significance of these trajectory classes for the ionization signal can be inferred from fig . [ fig2](c ) where their classical weights @xmath127 ( i.e. , the jacobi determinants ) are depicted as a function of the final angle @xmath67 . for small and very large final angles the contributions of all three classes are relevant whereas in the intermediate region the weight of class ( ii ) becomes negligibly small . this decrease is due to the fact that class ( ii ) includes trajectories with initial positions @xmath105 in the vicinity of @xmath128 and the behavior of such trajectories depends very sensitively on the initial conditions . as shown in fig . [ fig2](d ) the existence of different contributing trajectories with approximately equal weight leads to pronounced interference effects in the ionization spectrum . in particular , for small angles one observes rapid oscillations which gradually diminish as @xmath67 is increased . this diminution is directly related to the peculiar @xmath67-dependence of the classical weight of trajectory class ( ii ) that was described above . an important feature of figs . [ fig2](a ) and ( b ) is the existence of trajectories with @xmath129 whose initial positions do not lie on the @xmath119 axis . consequently , due to the rotational symmetry of the system around this axis , there are infinitely many classical trajectories that contribute to the ionization amplitude at these final angles so that the classical probabilities @xmath130 diverge there . this phenomenon is known as semiclassical glory scattering @xcite . for a proper semiclassical description the primitive ionization amplitude of eq . ( [ semi ] ) has to be regularized with the help of uniformization methods @xcite . thereby , one replaces the contributions of the trajectory classes @xmath131 and @xmath132 at both sides of the divergence by the uniform expression @xmath133 } \sqrt { \frac{\pi(f_a - f_b)}{4 } } \ { ( g_a + g_b)\times\nonumber\\ & & j_0[(f_a - f_b)/2 ] + i ( g_a - g_b ) j_1[(f_a - f_b)/2 ] \ } \label{uni1}\end{aligned}\ ] ] with @xmath134 and @xmath135 ( @xmath136 ) . the trajectories are labeled in such a way that @xmath137 so that @xmath138 . the bessel functions of order zero and one are denoted @xmath139 and @xmath140 . the contribution of the third trajectory can be described by the primitive semiclassical expression given in eq . ( [ semi ] ) . in fig . [ fig2](d ) the quantum mechanical ionization spectrum ( full curve ) is compared to the uniform ( dashed ) and the primitive ( dotted ) semiclassical approximation . all spectra in this section are calculated for an initial hydrogenic energy eigenstate with @xmath141 . the uniform curve describes the quantum mechanical result very well over the whole angular range whereas the primitive approximation diverges in small regions around @xmath142 and @xmath143 . the characteristic qualitative features of the example described above , i.e. , the existence of three contributing trajectories for all angles @xmath144 , their relative weights and the appearance of forward and backward glory effects are characteristic for all spectra represented by region ( i ) of fig . [ fig1 ] . the second region of particular importance in fig . [ fig1 ] is denoted ( ii ) . this region is bounded from below and above by the curves @xmath145 and @xmath146 with @xmath147 and @xmath148 ) and ( [ curve2 ] ) can be approximated as @xmath149 and @xmath150 , respectively . the spectra in region ( ii ) are characterized by the appearance of a semiclassical rainbow scattering phenomenon and a backward glory whereas the forward glory present in the spectra of region ( i ) has disappeared . again this behavior is most clearly illustrated by studying a particular example . to this end , the parameter values @xmath151 , @xmath152 are used ( with @xmath123 in fig.[fig3 ] which implies @xmath153 mev ) . in contrast to the previous case ( and , more generally , as soon as @xmath154 ) trajectories originating from the positive ( negative ) @xmath119-half plane now start with positive ( negative ) radial momentum . figure [ fig3](a ) depicts the location of initial positions in the @xmath122-plane and their corresponding final angles whereas fig.[fig3](b ) shows again the deflection function . apparently , for final angles larger than the  rainbow angle " @xmath155 the contributions of three different trajectory classes have to be taken into account . at @xmath155 , however , two of these classes coalesce which gives rise to the semiclassical rainbow scattering phenomenon @xcite . again the classical weights of the corresponding classical trajectories diverge at this angle [ cf . [ fig3](c)].to determine the transition amplitude the primitive approximation has to be replaced by a uniformized expression . using the notation of expression ( [ uni1 ] ) and defining @xmath156^{2/3}$ ] this uniform transition amplitude takes the form @xmath157 } \ { ( g_1 + g_2)\times\nonumber\\ & & \xi^{1/4}{\rm ai}(-\xi ) + i ( g_1 - g_2 ) \xi^{-1/4}{\rm ai}'(-\xi ) \}. \label{uni2}\end{aligned}\ ] ] in ( [ uni2 ] ) @xmath158 and @xmath159 denote the airy function and its derivative . quantum mechanically , the rainbow scattering phenomenon does not cease abruptly at the rainbow angle but extends smoothly into the classically forbidden region . this may be taken into account semiclassically by extrapolating the quantities @xmath160 , @xmath161 , and @xmath162 of expression ( [ uni2 ] ) to values @xmath163 . furthermore , figs . [ fig3](a)-(c ) show that the forward glory has disappeared but that the backward glory is still present . the classical weight associated with the different types of trajectories can be inferred from fig . [ fig3](c ) . figure [ fig3](d ) shows that the most prominent difference between the spectra of cases ( i ) and ( ii ) is that the latter start for small angles with a smooth non - oscillating part . oscillations set in only around @xmath155 . if , for a fixed value of @xmath164 , the final energy @xmath115 is increased the rainbow angle @xmath155 grows towards the value of @xmath143 and the interference effects in the spectrum are more and more diminished . if the line given by eq . ( [ curve2 ] ) is crossed , rainbow and backward glory scattering vanish completely . in this region ( iv ) the spectrum is determined by the contributions of a single trajectory class and interferences are no longer visible . however , fig . [ fig1 ] also indicates that in this dynamical region the angle - integrated ionization probability @xmath165 is very small so that this type of behavior is only of little practical relevance , in general . regions ( i ) and ( ii ) are separated by a small area ( iii ) in which the classical dynamics is more complicated . in this transition region besides the backward glory both forward glory and rainbow scattering characteristic of regions ( i ) and ( ii ) are present . in addition , however , a second classical rainbow angle appears at a very small value of @xmath67 . a proper semiclassical description of this behavior would require uniformization based on a higher - order catastrophe @xcite . figures [ fig2](a ) and [ fig3](a ) indicate that there is a strong correlation between the final direction into which the electron is emitted and its position before application of the pulse . this correlation might be useful for probing the spatial distribution of localized wave packets . in order to assess the perspectives of such an approach we plot in fig . [ fig6 ] several ionization spectra obtained in the sudden - impact approximation . they pertain to a radial rydberg wave packet detected at different stages of its time evolution . for the calculation , the wave packet was assumed to be generated by a weak and short laser pulse as described in ref.@xcite and the electronic angular momentum was taken as @xmath141 to show the effects in question most clearly . the spectral envelope @xmath166 of the laser pulse was chosen as @xmath167 $ ] with the mean excited energy @xmath168 a.u . corresponding to a mean quantum number @xmath169 and an energy spread determined by @xmath170 . figure [ fig6](a ) shows the spectra ( on a non - logarithmic scale ) for @xmath153 mev and various delay times between the exciting laser pulse and the hcp with @xmath171 a.u . all delay times are chosen such that the wave packet has not yet reached its outer turning point and is still moving outward . the inset depicts some corresponding radial distributions @xmath172 for the wave packet . the connection between these distributions and the form of the spectra can be made with the help of the right part of fig . [ fig3](a ) where @xmath173 . as long as @xmath172 is peaked at a distance @xmath174 less than ca . 4000 a.u . the spectrum attains its maximum at an angle @xmath175 whose value agrees with what can be inferred directly from @xmath174 and fig.[fig3](a ) . for larger values of @xmath174 the maximum is reached at @xmath176 but decreases as the wave packet moves beyond 4000 a.u .    for the case of an inward moving wave packet which is examined in fig.[fig6 ] the spectra are essentially peaked around the rainbow angle @xmath155 independent of the delay time . this observation agrees again with the classical picture obtained from fig . [ fig3 ] as the associated weight of the electron trajectories diverges at this angle . in order to achieve a connection between the maxima in the spectrum and the electronic probability distribution similar to fig . [ fig6](a ) one should investigate spectra for final energies @xmath177 [ cf . [ fig2](a ) ] . in summary , [ fig6 ] demonstrate that the spatial distribution of the wave packet is indeed reflected in the ionization spectrum . in order to obtain a complete picture of the distribution one will have to analyze a set of spectra with different parameters . a particular advantage of the use of hcps , compared , e.g. , to optical methods , is that they can interact with the wave packet at every point of its orbit .      in this subsection it is discussed to which extent the results of sec.iii a can be applied to the study of angle- and energy - resolved ionization spectra from hcps of finite duration . it is natural to examine this problem within a semiclassical approach as the theoretical framework discussed in sec . ii can be applied equally well to both instantaneous and time - dependent pulses and its quantitative accuracy has been demonstrated in sec . iii a. thus the main task consists in studying the classical electronic dynamics . in a fully quantum - mechanical treatment the partial wave expansion of sec.ii b would no longer be applicable . one would have to turn to a numerical integration of the time - dependent schrdinger equation which is significantly more complicated due to the coulomb singularity at the nucleus and due to the large spatial extension of highly excited rydberg states .    in the present context our main question is to which extent the essential results of the sudden - impact approximation can be recovered in experiments with finite - duration hcps . it is therefore reasonable to concentrate the study on short pulses with durations @xmath0 up to the order of , say , one tenth of the rydberg electron s orbit time @xmath178 . typically this corresponds to pulse lengths around 1 ps which can easily be produced in the laboratory . before discussing in detail an example that illustrates characteristic effects originating from finite - duration hcps let us briefly summarize the main results of our investigations :    \(i ) the behavior of trajectories starting close to the nucleus is modified significantly even by very short pulses whereas the dynamics of orbits originating at a large distance from the core is changed only very gradually . this is of course due to the different initial velocities of the trajectories : far way from the nucleus the electron moves slowly and the pulse can still be considered as almost instantaneous . from figs.[fig2 ] and [ fig3 ] it follows that the part of the spectrum with small final angles @xmath67 , where the ionization probability is concentrated and where forward glory and rainbow scattering can be observed , is related to trajectories with large initial distances from the nucleus . thus these essential regions of the spectrum can be expected to be relatively insensitive to effects of finite pulse durations . on the other hand , the behavior of the spectrum at large final angles is determined by trajectories starting close to the core . in this region of the spectrum , where the ionization probability is relatively small , pulse duration effects are most prominent . in particular , in our examples we even find backward glory scattering to vanish . \(ii ) the deviations of the spectrum from the sudden - impact approximation are mainly determined by the pulse duration whereas details of the pulse shape seem to be of less importance . we compared the classical dynamics and spectra resulting from two types of hcps the shapes of which were almost rectangular and a half - period sine wave , respectively . these choices should represent typical variations of realistic hcp shapes . the differences in the resulting spectra for these two pulse forms ( calculated for identical pulse durations and integrated field strengths ) were small in comparison to the deviations from the sudden - impact approximation ( see figs .  [ fig4](c ) and [ fig5 ] ) . \(iii ) in the sudden - impact approximation there are classical trajectories for _ all _ values of final energies @xmath179 . with increasing pulse duration , however , we find that the range of accessible final energies becomes more and more reduced in accordance with the energy - time uncertainty principle mentioned in the introduction . in order to discuss effects of finite pulse duration in more detail we consider as an example the ionization spectrum for an hcp with an almost rectangular pulse shape given by @xmath180 $ ] . the parameters @xmath181 and @xmath182 are determined so that the integrated field strength and the pulse duration fulfill @xmath151 and @xmath183 , respectively . again , for the calculation an initial hydrogenic energy eigenstate with @xmath141 and with principal quantum number @xmath123 was chosen . figure [ fig4 ] shows an analysis of the corresponding classical dynamics and of the resulting ionization spectrum for a final energy of @xmath152 ( compare to fig . [ fig3 ] ) . figure [ fig4](a ) details the changes in the classical dynamics . three major modifications can be observed :    \(i ) the final angles @xmath67 of trajectories with initially outgoing radial momentum do not cover the full range @xmath184 as in the sudden - impact approximation . instead , they are restricted to values @xmath185 with @xmath186 in the present case . here , this does not lead to major changes in the spectrum . this is because in the angular region where the deviations of this trajectory class from the sudden - impact case become significant the dominant contributions already originate from trajectory classes ( i ) and ( ii ) [ see figs.[fig3](b ) and [ fig4](b ) ] . however , due to the absence of trajectory class ( iii ) for @xmath187 the small short - period modulations of the spectrum of fig . [ fig3](c ) , which are visible for @xmath67 between @xmath188 and @xmath189 , are no longer present in fig . [ fig4](c ) . \(ii ) a new class of trajectories appears which has no counterpart in the sudden - impact approximation . these orbits have initially incoming radial momenta and originate from the half - plane @xmath173 . their final angles vary between @xmath190 and @xmath191 . the maximum final angle is not attained at @xmath192 but at a finite distance from the nucleus so that there is a rainbow effect connected with this class of trajectories [ cf . [ fig4](b ) ] . however , due to the low classical weight of these trajectories they show up in the spectrum most prominently in the form of small modulations in the vicinity of @xmath176 . \(iii ) the initial positions of trajectories which start with incoming radial momenta from the half - plane @xmath193 are changed significantly . in addition , the final angle @xmath194 is reached only for initial positions with @xmath195 , so that backward glory scattering has disappeared . the changes in the ionization spectrum at large angles which are due to these effects represent the most significant influence of finite pulse durations .    in fig . [ fig4](c ) the resulting ionization spectrum ( full curve ) is compared to a spectrum which has been calculated for an hcp of identical integrated field strength and pulse duration but with shape @xmath196 ( dashed curve ) . the spectrum for the sudden - impact approximation is also shown ( dotted curve ) . as mentioned above , for these pulse durations the influence of the hcp shape is only secondary even for large final angles . however , as the pulse length is increased to values larger than ca . 0.1 @xmath197 these influences eventually become more significant . figure [ fig5 ] shows the results of an investigation similar to fig . [ fig4](c ) for a final energy of @xmath198 ( @xmath199 mev ) which corresponds to region ( i ) of the sudden - impact approximation ( cf .  fig . [ fig1 ] ) . as in the previous case , it is found that the main features of the spectrum of the sudden - impact approximation are left almost unchanged . the classical dynamics change similarly to the previous example . in particular , the new trajectory class still exists , but in the rapidly oscillating part of the spectrum its contribution is almost undiscernible . a general semiclassical treatment of excitation of weakly bound rydberg electrons by intense hcps has been presented . thus a quantitative connection between observable transition probabilities and the underlying classical dynamics of the excited rydberg electron has been established . this approach is numerically accurate for high principal quantum numbers of the rydberg electron and in cases in which the energy transfer from the exciting hcp to the rydberg electron is small in comparison with typical ionization energies of low - lying states . the underlying classical dynamics of the rydberg electron yields a clear and natural explanation for the oscillatory structures which govern the resulting energy- and angle - resolved ionization spectra . these structures arise from interference between probability amplitudes associated with different ionizing trajectories . thereby two types of semiclassical catastrophes appear , namely glory scattering for electron emission in the forward and backward direction of the polarization of the hcp and a rainbow phenomenon . parts of these oscillatory structures depend strongly on the ratio between the pulse duration of the exciting hcp and the classical kepler period of the initially prepared rydberg state . the amplitude and phase information which is contained in the semiclassical probability amplitudes might be useful for the reconstruction of quantum states of rydberg electrons . a first indication of the usefulness of such an approach was given by the investigation of the spectra pertaining to radial wave packets . g.  a.  acknowledges support by the deutsche forschungsgemeinschaft within the spp  zeitabhngige phnomene und methoden . " o.  z.  is supported by the u.s . office of naval research under contract no . 14 - 91-j1205 and by the u.s .  army research office . r.  r.  jones , d.  you , and p.  h.  bucksbaum , phys . 70 * , 1236 ( 1993 ) ; g.  m.  lankhuijzen and l.  d.  noordam , phys .  rev .  lett . * 74 * , 35 ( 1994 ) ; n.  e.  tielking and r.  r.  jones , phys . a. * 52 * , 1371 ( 1995 ) ; r.  b.  vrijen , g.  m.  lankhuijzen , and l.  d.  noordam , phys . lett . * 79 * , 617 ( 1997 ) . r.  r.  jones , phys . * 76 * , 3927 ( 1996 ) . c.  raman _ et al . _ , phys .  rev .  lett . * 76 * , 2436 ( 1996 ) . c.  o.  reinhold , h.  shaw , and j.   burgdrfer , j.  phys . b. * 27 * , l469 ( 1994 ) ; k.  j.  la gattuta and p.  b.  lerner , phys . a. * 49 * , r1547 ( 1994 ) ; a.   bugacov _ et al . _ , phys . a. * 51 * , 1490 ( 1995 ) ; c.  o.   reinhold _ et al . _ , j.  phys . b. * 28 * , l457 ( 1995 ) ; c.  o.  reinhold _ et al . _ , phys . a. * 54 * , r33 ( 1996 ) ; m.  t.   frey _ et al . _ , phys . a. * 55 * , r865 ( 1997 ) ; r.   gebarowski , j.  phys . b. * 30 * , 2143 ( 1997 ) ; s.  yoshida _ et al . _ , phys . a. * 58 * , 2229 ( 1998 ) . m.  mallalieu and shih - i chu , chem . * 258 * , 37 ( 1996 ) . c.  d.  schwieters and j.  b.  delos , phys . a * 51 * , 1023 ( 1995 ) ; phys .  rev . a * 51 * , 1030 ( 1995 ) . g.  alber and o.   zobay , to be published in phys . rev .  a. h.  a.  bethe and e.  salpeter , _ quantum mechanics of one- and two- electron atoms _ ( plenum , new york , 1977 ) . m.  j.  seaton , rep . 46 * , 167 ( 1983 ) . u.  fano and a.  r.  p.  rau , _ atomic collision and spectra _ ( academic , new york , 1986 ) . maslov and m.  v.  fedoriuk , _ semiclassical approximation in quantum mechanics _ ( reidel , boston , 1981 ) . j.  b.  delos , adv . phys .  * 65 * , 161 ( 1986 ) . l.  s.  schulman , _ techniques and applications of path integration _ ( wiley , new york , 1981 ) . e.  j.  heller , j.  chem . * 75 * , 2923 ( 1981 ) . w.  h.  miller , adv . phys .  * 25 * , 69 ( 1974 ) . m.  f.  herman and e.  kluk , chem . phys . , 27 ( 1984 ) . k.  g.  kay , j.  chem .  phys . * 100 * , 4377 ( 1994 ) . g.  campolieti and p.  brumer , phys . a * 50 * , 997 ( 1994 ) . f.  grossmann , comments at . phys . , 141 ( 1999 ) . m.  v.  berry and k.  e.  mount , rep . .  phys . * 35 * , 315 ( 1972 ) . l.  d.  landau and e.  m.  lifshitz , _ mechanics _ ( pergamon , oxford , new york , 1976 ) . relations ( [ curve1 ] ) and ( [ curve2 ] ) can be derived by noting that @xmath200 changes sign when the boundary between regions ( ii ) and ( iii ) [ or ( iv ) , respectively ] is crossed . g.  alber , h.  ritsch , and p.  zoller , phys . a * 34 * , 1058 ( 1986 ) .
Abstract:

Answer: the interaction of a weakly bound rydberg electron with an electromagnetic half - cycle pulse ( hcp ) is described with the help of a multidimensional semiclassical treatment . this approach relates the quantum evolution of the electron to its underlying classical dynamics . the method is nonperturbative and is valid for arbitrary spatial and temporal shapes of the applied hcp . on the basis of this approach angle- and energy - resolved spectra resulting from the ionization of rydberg atoms by hcps are analyzed . the different types of spectra obtainable in the sudden - impact approximation are characterized in terms of the appearing semiclassical scattering phenomena . typical modifications of the spectra originating from finite pulse effects are discussed .