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: understanding the behavior of a particle diffusing in the presence of traps is important as it captures the essence of many physical processes , including diffusion controlled reactions where one of the reactants is immobile , trapping of excitons , etc . @xcite mapping of the diffusion problem with traps to the scalar elasticity problem also affords an insight into lattice vibration of systems such as binary alloys with contrasting elastic constants . from a theoretical point of view , this problem is analogous to the ideal chain in inhomogeneous media . the ideal chain problem is interesting as it is the gaussian limit of the self - avoiding walks and yet possesses a universality class of its own in inhomogeneous media , distinct from the usual random walk . these considerations resulted in many analytical and numerical attempts to gain understanding of the problem @xcite . typical quantities used to characterize diffusion in the presence of traps include @xmath0 , the probability that the diffusing particle returns to the starting point after time @xmath1 . @xmath0 is related to the number of distinct sites visited by the diffusing particle and is the laplace transform of the vibrational density of states of the corresponding scalar elasticity problem , which can often be measured experimentally by methods such as raman and neutron scattering @xcite .    for the case of diffusion in the presence of traps distributed randomly with no correlations , it has been proved rigorously by donsker and varadhan @xcite that the decay of @xmath0 with time is slower than exponential , indicating that , in the long time limit , the properties of the diffusing particle is dominated by the presence of large trap free regions which have a finite ( though small ) probability of occurence . this is because , if the traps were uniformly distributed , the diffusing particle would be trapped at a constant rate , which would result in an exponential decay of @xmath2 . it also means that the _ quenched _ disorder average must be carried out since an annealed averaging would smear out the trap positions .    in this paper we ask the question of whether and how this interesting behavior of the diffusing particle with traps is modified when we introduce long range correlations in their positions @xcite . in the equivalent scalar elasticity problem , the correlated traps map to the clamping of sites with a correlated distribution . in this sense , the present problem extends the so - called _ fractino _ problem @xcite where a fractal boundary of an otherwise nonfractal object is clamped to the case where the bulk of the substrate is itself a fractal . it is also an extension of the _ fracton _ problem @xcite as the traps ( or clamped boundaries ) are introduced into the scalar elasticity of fractals . calculational difficulties have prevented this problem from receiving its share of attention even though , often in real physical situations , the positions of traps are correlated at least within limited length scales . for the case of uncorrelated trap distribution , a poisson distribution is usually assumed , which simplifies the theoretical calculation @xcite ; however this assumption fails for correlated traps . as far as computational calculations , techniques based on exact enumeration are highly computer time and memory intensive because of the high sensitivity of the behavior of the diffusing particle to the actual positions of the traps in the sample , leading to large fluctuations in the measured quantity from sample to sample . this translates to requiring an ensemble average to be taken over a large number of disorder configurations for meaningful results . moreover the onset of the asymptotic regime is very slow in this type of problem , which adds to the computational difficulties @xcite .    in this paper we establish that @xmath2 has a stretched exponential form _ qualitatively _ similar to the case of the poissonian trap distribution but with a substantially different stretch exponent . we show this behavior from the analysis of the density of eigenvalues of the transition probability matrix @xmath3 which describes the diffusion process . we also study the finite size scaling of the largest eigenvalue of @xmath3 , which serves as an extremely powerful tool to extract the stretch exponent characterizing the long time behavior of @xmath2 . a preliminary account of this work was presented at the hayashibara forum ( 1995 ) and appears in its proceedings . we choose percolation cluster @xcite formed at critical probability of occupation @xmath4 as the substrate for diffusion . the sites on the external perimeter ( hull ) @xcite as well as on the internal perimeter of the percolation cluster are made absorbing and once the diffusing particle reaches the perimeter sites it is trapped permanently . the hull sites and the sites constituting the internal perimeter form fractals at @xmath4 @xcite . thus placing the traps along the perimeters induces spatial correlations among the traps because of the long range correlations among the sites of a fractal . while we do not pursue the distinction of external versus internal perimeters in this work , there _ are _ interesting effects when different boundary conditions are applied to them , at least in two dimensions @xcite .    the time evolution of the probability @xmath5 that the diffusing particle is at site @xmath6 at time @xmath1 is markovian and the process can be described in the continious time limit by the following master equation @xmath7 where @xmath8 is the hopping rate from site @xmath9 to @xmath6 . in this problem we consider only nearest neighbour hopping and so @xmath8 is non - zero only for the nearest neighbour pairs @xmath10 . we have chosen @xmath8 to be @xmath11 for all occupied nearest neighbours @xmath9 , where @xmath12 is the full coordination number of the underlying lattice ; otherwise @xmath8 is set equal to zero . what this amounts to is that once the particle hops to a trap site there is no further time evolution of that particular random walk and only those walks which have escaped getting trapped in the perimeter sites evolve further . thus the traps act as permanent particle absorbers . processes whose time evolution is governed by eq . ( [ eq : markov ] ) can be cast into an eigenvalue problem of the transition probability matrix @xmath3 . the matrix elements @xmath8 of @xmath3 control the dynamics of the random walk and the locations of the nonzero elements of @xmath3 have the information about the structure of the underlying fractal substrate responsible for the correlations among the trap sites on its perimeters . the density of normal modes of @xmath3 is related to the return to the starting point probability of the diffusing particle by the laplace transform @xcite . the probability distribution of the number of walks @xmath13 which return to their starting point after time @xmath1 , denoted as @xmath14 , was studied in @xcite and found to be a truncated log - normal distribution . thus @xmath15 in @xcite it was further found from the first moment of this distribution that @xmath16 in the asymptotic long time limit , where @xmath2 is obtained as @xmath17 , the bar above the quantity indicating the quenched disorder average . the exponent @xmath18 is the same one as that which appears in the long time behavior of the width of the log - normal distribution , @xmath19 , @xmath20    since @xmath2 has a stretched exponential behavior according to eq . ( [ eq : rtsp ] ) , we expect the density of normal modes @xmath21 ( which is the inverse laplace transform of @xmath2 ) to also have a stretched exponential form @xcite , @xmath22 in the limit of small @xmath23 , where @xmath24 denotes the eigenvalues of @xmath3 , and the exponents @xmath25 and @xmath26 are related to each other by @xmath27 . the results from @xcite gave support for the presence of such a behavior although their numerical estimates of @xmath26 need to be improved as they did not take into account the proper normalization of @xmath21 as well as the substantial non - asymptotic effects @xcite . we can also relate the behavior of @xmath28 to the finite size effects of the edge of the spectrum . first , note that , for any finite substrate however large , the largest eigenvalue of @xmath3 , which we denote by @xmath29 , must be less than one . this is because the eigenvalue of one would indicate the existence of a stationary state whereas there are always particles which are trapped , leading to a _ leakage _ of probability for the diffusing particle in contrast to the diffusion _ without _ traps . the value of @xmath30 then corresponds to the slowest time scale of the problem , and there is always a gap between @xmath29 and one , the latter being what the maximum eigenvalue would be for diffusion without traps . the interesting question is whether this gap is bounded from below by a nonzero constant or it approaches zero as the substrate size increases .    our argument for the behavior @xmath31 as the substrate size @xmath32 follows the observation of @xcite for the case of uncorrelated trap distribution . in that case , the return probability @xmath2 was shown to be a stretched exponential which is slower than a pure exponential ( although not as slow as a power law which would be the case in the absence of traps ) , and this behavior was attributed to the dominance of the trap - free regions although they occur relatively infrequently . in our case of critically correlated trap distribution , we also have a slow , stretched exponential for @xmath2 which we interpret as a similar dominance of the nearly trap - free regions . although there are more traps in the critically correlated case , there are also much greater fluctuations in their density . thus regions of relatively low concentrations of the traps are likely to be present , probably with a hierachical size distribution . since the leakage of probability is _ nearly _ zero in a _ nearly _ trap free region , we would expect @xmath29 to approach one as the size of the largest trap free region grows indefinitely . another argument is as follows : if there were indeed a nonzero lower bound for @xmath33 , then it would also give a bound for the slowest relaxation rate . thus @xmath2 would have to decay at least exponentially in time corresponding to this rate . since @xmath2 in fact decays as a stretched exponential with the stretch exponent less than one ( @xcite as well as this work ) , no such bound can exist for the rate ; rather @xmath34 as @xmath32 . we further propose a scaling relation between @xmath35 and @xmath36 via @xmath37 where a possible numerical prefactor has been neglected . this relation follows from the assumption that , say , the largest ( @xmath29 ) and second largest ( @xmath38 ) eigenvalues scale relative to the value one in the same way when @xmath32 . that is , if we assume @xmath39 where @xmath40 , then the difference must also scale in the same way , @xmath41 but this latter integral must scale as @xmath42 if the integral of the density of normal modes @xmath43 is normalized to unity . this means @xmath44 , thus eq . ( [ eq : eg ] ) follows . the assumption of the same asymptotic behavior for the two largest eigenvalues is plausible if we consider the difference between @xmath29 ( or @xmath38 ) and one to be the reflection of the finite size of the largest ( or the second largest ) trap free region ( respectively ) . as long as the large trap free regions are geometrically similar ( as would be the case in a hierarchical distribution of such regions ) , and as long as their sizes all go to infinity as @xmath32 , it seems reasonable that the gaps between @xmath45 ( @xmath46 ) and one behave in the same manner .    if we substitute a stretched exponential form of @xmath43 in eq . ( [ eq : eg ] ) we get an incomplete gamma function . on retaining only the leading term of the incomplete gamma function and taking natural log of both sides we get , @xmath47 where a general stretched exponential form of @xmath43 @xmath48 has been assumed and @xmath49 .    while the above gives a particular scaling prediction for the case of the critically correlated trap distribution , the argument leading to eq . ( [ eq : eg ] ) applies more generally . for example , a similar procedure should apply to the uncorrelated trap distribution of @xcite . also , for the trapless case ( or the _ ants _ @xcite ) , where the largest eigenvalue is actually one ( because there _ is _ a stationary state ) , a similar argument should work with , say , the second ( @xmath38 ) and third largest ( @xmath50 ) eigenvalues . of course , if there is no trap , a probability leakage , per se , does not occur . however , each mode with a large @xmath24 tends to be associated with a _ blob _ of high connectivity region , to which the same kind of argument can be applied . indeed , in the trapless case , eq .  ( [ eq : eg ] ) together with the power law density of states @xmath51 where @xmath52 is the spectral dimension of the substrate , leads to a power law relation between @xmath36 and @xmath53 , @xmath54 which is identical to the relation proposed on the basis of the finite size scaling of the largest nontrivial normal mode and numerically verified in @xcite . such finite size scaling relations provide a very powerful technique to obtain the quantitative characterization of the density of states as they reduce the computational effort drastically , requiring only the information on the highest ( or second highest ) mode . in this section we give numerical results for the exponent @xmath55 for the stretched exponential decay of the density of states as discussed above . we extract this exponent directly from the density of states ( cf . ( [ eq : spec ] ) ) and also independently from the finite size scaling of the largest eigenvalue of @xmath3 ( cf . ( [ eq : fs ] ) ) in two and three dimensions ( square and simple cubic lattice , respectively ) .    in order to reduce the computational time and memory requirements , we take advantage of the fact that we are interested only in the asymptotic long time relaxation of the system which is controlled by those normal modes with large eigenvalues @xmath24 . thus we use the arnoldi - saad algorithm @xcite to reduce the original @xmath3 into a smaller matrix which contains the approximate information about the highest normal modes . we analyze the density of states per site , @xmath43 . it is obtained from the eigenvalues of @xmath3 by binning them linearly in the @xmath56 space . the number of eigenvalues in each bin is divided by the bin width , the size of the cluster and also by the number of clusters over which the quenched disorder average is performed . the number of independent cluster realizations over which the disorder average was taken is of the order @xmath57 for cluster sizes @xmath58 @xmath59 , @xmath60 , and @xmath61 and of the order @xmath60 for @xmath58 @xmath62 . we retained in the final results only those bins which contained eigenvalues contributed from every substrate realization to avoid partial binning .    since @xmath43 has a stretched exponential decay we have plotted @xmath63 versus @xmath56 in fig . [ fig1 ] which is expected to have the form @xmath64 . there is an excellent data collapse for all the cluster sizes both in 2d and 3d . the solid curve drawn through the data is obtained by fitting the data to an expression of this form . the numerical estimates of @xmath55 and the numerical values of @xmath65 and @xmath66 corresponding to the central values of @xmath55 are tabulated in table  [ table1 ] . note that the value of @xmath55 can not be accurately obtained simply from the slope of the double log plot of @xmath67 with respect to @xmath56 as the value of the constant @xmath66 might be appreciable . the size of the symbols in fig . [ fig1 ] is larger than the cluster to cluster statistical fluctuations . the quoted error bars are obtained visually from changing the effective estimate of @xmath55 for the nonlinear fit until it no longer fits the data points . the small scaling regime which is typical for these kinds of problems ( i.e. , diffusion with traps ) @xcite makes the precise extraction of @xmath55 difficult . this accounts for the large error bars in our results as compared to the case of diffusion in a percolation cluster without absorbing sites @xcite . in the latter case the density of states has a power law form and the scaling regime increases appreciably with the size of the cluster , which resulted in much smaller error bars for the extracted exponents even with the same numerical technique .    in fig . [ fig2 ] we plot @xmath68 with respect to @xmath35 in 2d and 3d where @xmath36 is the size of the substrate . the solid curves are obtained by fitting the data to an expression @xmath69 , which is what we expect from the finite size scaling analysis of the largest eigenvalue . the size of the symbols are larger than the cluster to cluster fluctuations . the number of clusters over which the disorder average was performed is of the order 1000 for the larger clusters , same as for the density of states , and for smaller clusters of size @xmath70 , @xmath71 , @xmath57 and @xmath62 the average was taken over @xmath60 clusters . the estimates of @xmath55 , and the values of @xmath65 , and @xmath72 corresponding to the central values of @xmath55 which we obtain from the nonlinear fit are tabulated in table  [ table2 ] . the error bar for the value of @xmath55 is obtained in a similar way to that for the density of states . we note that the estimates of @xmath55 and @xmath65 obtained from the density of states and from the finite size scaling of the largest eigenvalue , given in table  [ table1 ] and [ table2 ] , respectively , are in good agreement with each other . the values of @xmath26 thus obtained , however , turn out to be about a factor of two larger than the corresponding stretch exponents for the uncorrelated distribution of traps @xcite . the estimates of @xmath55 also differ from those of @xmath73 as given by @xcite due , we believe , primarily to the failure of @xciteto take proper account of @xmath56 being not quite in the asymptotic region ( thus proper normalization and prefactors becoming important ) ( cf . @xcite ) . however , the exponent @xmath25 ( eq .  ( [ eq : chi ] ) ) is relatively insensitive to @xmath55 ( @xmath74 ) and the analysis of @xmath2 in @xcite is not affected by these problems . we also believe that their main conclusions remain valid . in summary , we have studied the problem of diffusion on a substrate with permanent traps which are distributed with critical correlation in their positions . the critical correlation has been modeled by placing the traps on the perimeters of critical percolation clusters , which are obtained by monte carlo simulation . the statistical behavior of the diffusing particle in the time domain is then mapped to the scalar elasticity problem with fixed , fractal boundaries , and a numerical analysis of the vibrational density of states of the latter problem is carried out using an approximate diagonalization algorithm of arnoldi and saad @xcite . this problem may be considered a generalization both of the _ fracton _ problem of alexander and orbach @xcite , where there are no traps but the substrate is fractal , and of the _ fractino _ problem of sapoval et al @xcite where a fractal boundary is clamped ( equivalent to traps ) but the substrate itself is not fractal . we have shown that introducing critical spatial correlations among the traps does not change the _ qualitative _ behavior of the density of states as compared to the case of uncorrelated traps , but that there are substantial effects of correlations _ quantitatively_. this strongly suggests that the diffusion process is dictated by the presence of nearly trapless regions similarly to the case of uncorrelated traps , but the long range correlations among the trap positions lead to a decrease in the number and size of these regions nearly free of traps , which induces a faster decay in the density of states and consequently a smaller probability of return to the starting point after time @xmath1 . this is reflected in a much larger exponent in the stretched exponential form of density of states than the case of uncorrelated traps .    on the other hand , as compared to the case of the same _ fractal _ substrate but with _ no _ traps , the density of states is _ qualitatively _ different , since the latter problem produces a power law density of states in the low energy limit known as _ fractons _ , accumulating to an infinite density toward the maximum eigenvalue of one @xcite . in contrast , with the perimeters acting as traps , the density of states becomes a stretched exponential , rapidly falling to zero toward the maximum eigenvalue . r. f. kayser and j. b. hubbard , phys . * 51 * , 79 ( 1983 ) . s. havlin and d. ben - avraham , advances in physics . * 36 * , 695 ( 1987 ) and references therein . e. courtens , r. vacher and e. stoll , in _ fractals in physics _ , a. aharony and j. feder , eds . ( north - holland , amsterdam , 1989 ) . m. donsker and s. r. s varadhan , commun . pure appl . math . * 28 * , 525 ( 1975 ) , and * 32 * , 721 ( 1979 ) . a. giacometti and a. maritan , phys . e. * 49 * , 227 ( 1994 ) . s. alexander and r. orbach , j. phys . ( paris ) * 43 * , 625 ( 1982 ) . b. sapoval , th . gobron , and a. margolina , phys . 67 * , 2974 ( 1991 ) . p. grassberger and i. procaccia , j. chem . phys . * 77 * , 6280 ( 1982 ) . s. havlin , m. dishon , j. e. kiefer , g. h. weiss , phys . 53 * , 407 ( 1984 ) . h. nakanishi and s. mukherjee , physica a ( 1995 ) , in press . d. stauffer and a. aharony , _ introduction to percolation theory _ , ( taylor and francis , london , 1992 ) . b. duplantier , j. de physique * xx * , xx ( 19xx ) . see , e.g. , j. feder , _ fractals _ ( plenum , new york , 1988 ) . s. mukherjee and h. nakanishi ( unpublished ) . s. mukherjee , h. nakanishi , and n. fuchs , phys . e * 49 * , 5032 ( 1994 ) , and references therein . a. giacometti and h. nakanishi , phys . e * 50 * , 1093 ( 1994 ) , see also a. giacometti , h. nakanishi , a. maritan and n. h. fuchs , j. phys . a * 25 * , l461 ( 1992 ) . a. giacometti and h. nakanishi ( unpublished ) . y. saad , linear algebra appl . * 34 * , 269 ( 1980 ) and w. e. arnoldi , quart . appl * 9 * , 17 ( 1951 ) .    .estimates of the exponent @xmath55 and the constants @xmath65 and @xmath66 from @xmath43 for the square lattice in two dimensions and simple cubic lattice in three dimensions . the data for @xmath63 have been fitted to the form @xmath75 . [ cols="^,^,^,^",options="header " , ]
Abstract:

Answer: we study diffusion on a substrate with permanent traps distributed with critical positional correlation , modeled by their placement on the perimeters of a critical percolation cluster . we perform a numerical analysis of the vibrational density of states and the largest eigenvalue of the equivalent scalar elasticity problem using the method of arnoldi and saad . we show that the critical trap correlation increases the exponent appearing in the stretched exponential behavior of the low frequency density of states by approximately a factor of two as compared to the case of no correlations . a finite size scaling hypothesis of the largest eigenvalue is proposed and its relation to the density of states is given . the numerical analysis of this scaling postulate leads to the estimation of the stretch exponent in good agreement with the density of states result .


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: real graphene has always defects , edges , point defects , chemical impurities  and grain boundaries . grain boundaries ( gbs ) are topological defects , trails of disorder that separate two pieces of pristine hexagonal graphene sheets . they reside practically in any graphitic material , graphite@xcite , soot@xcite , single- and multi - layer graphene@xcite , fullerenes@xcite , carbon nanotubes@xcite , with varying abundance . in graphene fabrication , for instance , chemically synthesized samples tend to have more gbs than mechanically cleaved samples@xcite . for electron mobility high purity is important and gbs are best avoided , while other properties may tolerate defects better . we refer here to strictly two - dimensional gbs ( albeit not necessarily planar ) which should not be confused with genuinely three - dimensional gbs in graphite . but not are defects , impurities and gbs always something bad ; they are interesting also in their own right  even useful@xcite . after all , being extended , gbs modify graphene more than point defects . they affect graphene s magnetic , electronic , structural , and mechanical properties@xcite .    considering the relevance and prevalence of graphene gbs , they ought to deserve more attention . most related studies are on gbs in graphite and moir patterns therein@xcite , or on intramolecular junctions in carbon nanotubes@xcite ; gbs in single- or few - layer graphene , whether experimental@xcite or theoretical@xcite , are still scarce . often graphitic grain boundaries are idealized by pentagon - heptagon pairs@xcite , where the pair distances determine lattice mismatch . while this ideal model works in certain occasions@xcite , it does not work for gbs that are rough , corrugated , have ridges@xcite , and  most importantly  are decorated by dangling bonds or reactive sites alike@xcite . in this paper we investigate graphene gbs beyond simple models . we explore computationally the trends in structural properties , chemical reactivity and vibrational properties for an ensemble of gbs that are _ not ideal _ , but have certain roughness . what this roughness really means will be clarified in sec.iii . the electronic structure method we use to simulate the gbs is density - functional tight - binding ( dftb)@xcite , and the ` hotbit ` code@xcite . dftb models well the covalent bonding in carbon@xcite , and suits fine for our simulations that concentrates on trends . we leave the details of the method to ref .   and only mention  this will be used later  that the total energy in dftb , including the band - structure energy , the coulomb energy , and the short - range repulsion , can be expressed as a sum of atomic binding contributions , @xmath0 here @xmath1 is the atom count and @xmath2 atom @xmath3 s contribution to cohesion energy ; in pristine graphene @xmath4  ev for all atoms ( dftb has some overbinding ) . the quantity @xmath2 allows us to calculate energies in a local fashion : for a group of carbon atoms in a given zone @xmath5 ( set of atoms ) , the sum @xmath6 measures how much the energy of zone @xmath5 is larger than the energy of equal - sized piece of graphene . in addition , since @xmath7  ev comes from @xmath8 equivalent bonds , the value @xmath9  ev infers a dangling bond for atom @xmath3 . hence , albeit less accurate than density - functional theory , dftb enables straightforward analysis of local quantities , in addition to fast calculations and extensive sample statistics . to investigate trends , we constructed an ensemble of @xmath10 graphene gbs using a procedure we now describe . ( i ) cut two ribbons out of perfect graphene in a given chirality direction , with random offset . chirality index @xmath11 means we cut graphene along the vector @xmath12 , where @xmath13 are the primitive lattice vectors . the cut direction can be expressed , equivalently , vy the chiral angle @xmath14 $ ] . ( ii ) passivate ribbons other edges with hydrogen , leaving carbon atoms on the other edges bare . ( iii ) thermalize the ribbons , still separated , to @xmath15  k using langevin molecular dynamics ( md ) . ( iv ) perform reflection operation on the other ribbon with respect to a plane along ribbon axis and perpendicular to ribbon plane . ( v ) perform random translations along the ribbon direction . ( vi ) merge the two ribbons bare carbon edges with md at @xmath15  k. ( vii ) cool the system to room temperature within @xmath16  ps . ( vii ) optimize the geometry using the fire method@xcite , and finally obtain structures like the one in fig.[fig1 ] . the total process takes about @xmath8  ps .    , made by merging two graphene edges with the same chirality , but different orientation and random translation along the boundary ; our samples are hence characterized by chirality . in this figure chiral indices are @xmath17 , chiral angle @xmath18 . dashed vertical lines stand for periodicity across the horizontal direction . ] our ensemble of @xmath10 gb samples are chosen to get representatives for all chiral angles from @xmath19 to @xmath20 . our limit for the width in cutting the ribbons is @xmath21   . ribbons have to be wide enough ( @xmath22   ) for faithful description of a semi - infinite graphene , and , by limiting the total number of atoms below @xmath23 , also the periodicity of gb becomes limited below @xmath24   . ( edges with @xmath25 and @xmath26 , that would have short periodicity , are scarce . ) our procedure to make gbs is neither the only nor the best one , but is suitable for trend - hunting . sure enough , it s a downside that the procedure excludes gbs from ribbons with different edge chirality , such as merged zigzag and armchair edges . different chirality , unfortunately , would mean different periodicity and practical problems . to have a single number , @xmath27 , to identify gbs is vital for finding the trends we concentrate on . for gbs identified by two numbers , @xmath28 and @xmath29 , trend - hunting would be harder ; those gbs can be investigated afterwards , using the insights we learn here . further , while the temperature @xmath15  k , motivated partly by experiments@xcite , allows rearrangements in merging , the heat of fusion renders initial temperatures irrelevant . we use canonical md in the merging to allow _ gb formation process to be steered by intrinsic driving forces and energetics_.    the time scale of the process is limited , as usual , in part by computational constraints . however , prolonging the process would cause no fundamental structural changes as the merging process itself is nearly instantaneous ; there is no diffusion after merging and saturation of the dangling bonds , there is only annealing of the worst local defects and stress - release that causes buckling ( sec.[sec : buckling ] ) . the usual time scales in graphene growth are hence not relevant to our gb formation process . only intrinsic roughness remains in the gb samples  roughness related to chiral angle and randomness in cutting offsets and translations . usually all atoms in ribbons were also bound to gbs , but in few structures ( near armchair chirality ) the merging process squeezed dimers out , and either inflected them out of plane or detached them completely . since ribbons found these own ways to optimize geometries  ways not anticipated beforehand  it suggests that the design of our construction process is not dominant . and armchair @xmath30 . free zigzag edge segments are metastable and prefer local reconstruction ( two adjacent hexagons reconstruct into a pentagon and a heptagon ) ; these so - called reczag edges are shown with gray connected symbols . the reczag segments are , however , important only for free edges and irrelevant for grain boundaries studied here . ]    before going into gbs themselves , let us look at the free bare edges of graphene . fig.[fig2 ] shows the edge energies as a function of the chiral angle , @xmath31 meaning zigzag and @xmath30 armchair edges . to ignore the effect of hydrogen passivated edge , carbons bound to hydrogen are neglected in edge energy calculation , @xmath32    the edge energy between zigzag and armchair varies linearly ; fluctuations in energy are due to random offset in the cut , occasionally producing pentagons . the edge energy in zigzag is high due to strong and unhappy dangling bonds ; in armchair the dangling bonds are weakened by the formation of triple bonds in the armrest parts@xcite . it was recently predicted theoretically@xcite and later confirmed experimentally@xcite that the zigzag edge is actually metastable , and prefers reconstruction into pentagons and heptagons at the edge , forming a so - called reczag edge . as it turned out , reczag is energetically even better than armchair . hence , if we reconstruct the zigzag segments in edges with small @xmath27 , we usually lower the edge energy , as seen in fig.[fig2 ] where reconstructions are added by hand . reczag edges are , however , irrelevant for our gbs where edges are _ not _ free , and are ignored because we want the dangling bonds to spontaneously find contact from the other merging edge . the edge energies were investigated and presented here for comparison with dft calculations@xcite . the accuracy in edge energies is better than @xmath21  % , and we expect same accuracy in gb energetics . the gb energy per unit length is @xmath33 measuring how much gb costs energy relative to the same number of carbon atoms in graphene . the energies for the whole ensemble of gbs are shown in fig.[fig3 ] , as a function of the chiral angle . the value @xmath34  ev /   appears only with @xmath19 and @xmath20meaning pristine graphene .    ) and armchair ( @xmath30 ) chiral angles . ] most gb energies are less than or equal to the edge energies of the free graphene edges , albeit with variation . this means that gbs regain , on average , other ribbon s edge energy during the merging , and the energy of fusion ranges from @xmath35  ev /  ; on average half of the free edges dangling bonds get passivated . the picture is not this simple , however , as md simulation creates different polygons that cause strain . the randomness of the polygon formation is manifested by the energy variations in fig.[fig3 ] . compared to ideal gbs with pentagons and heptagons only , our ensemble reveals the full complexity that rough gbs have . , is manifested by the appearance of mainly pentagons and heptagons , but also squares , octagons and nonagons . in b ) each dot represents one polygon in all the grain boundary polygons in all of our @xmath10 samples . abscissa shows the area of that polygon ( determined by triangulation ) and ordinate shows the polygon energy cost per atom , relative to an atom in pristine graphene . ] gb energies are the highest around @xmath36 , which can be understood from structural analysis , shown in fig.[fig4]a as polygon distribution within the gb zone . what we count into gb zone are all the polygons that were _ not _ part of the ribbons prior to merging . polygons larger than nonagons , which appear more like vacancies instead of polygons , are omitted here . it is around @xmath36 where the abundance of hexagons is at minimum and gbs are invaded by other polygons . the abundance of pentagons and heptagons is as high as the abundance of hexagons , but also squares , octagons and nonagons are found . close to @xmath25 and @xmath26 hexagons prevail . from this we may conclude that around @xmath36 the edge geometries have the largest mismatch , resulting in various polygons , consequent strains , and high energy .    to understand how much different polygons cost energy , we used the quantity @xmath37 measuring how much atoms , on average , cost more in @xmath38-gon relative to atoms in pristine graphene ; this quantity is for illustration only and considers polygons as separate items  the sum of @xmath39 s for given gb is not the total energy . fig.[fig4]b shows @xmath39 s for the polygons within gb zones as a function of the area of the polygon ( determined by triangulation ) . for simple geometrical reasons for small polygons the area distribution is narrow ; large polygons have more freedom to change their shape . for small polygons the distributions in @xmath39 , on the contrary , are wider ; this is partly due to smaller @xmath38 in eq.([eq : polye ] ) . clearly , the cost @xmath39 of any polygon depends on its environment , just as it also depends for hexagons , for which @xmath40 but note that fig.[fig4]b already contains the effect of the polygon environment and all potential cross - correlations ( such as pentagons often neighboring heptagons ) . it would be interesting to investigate polygon statistics also from transmission electron microscopy , now that aberration - corrected measurements can achieve atom accuracy  @xcite . in the generation process gbs are free to deform . the optimized gbs will typically end up having inflection angles , as the inset in fig.[fig5 ] illustrates . the data points in fig.[fig5 ] show the inflection angles for the gb samples as a function of chirality . flat gbs occur with @xmath19 and @xmath20 , that is , with pristine graphene alone . the notable trend in this scattered plot is the systematically small inflection angles , meaning flat gbs , around @xmath41 ; near @xmath42 inflection angles are more scattered meaning geometries that vary from flat to sharply kinked gbs . this trend can be understood by edge profiles : the stronger dangling bonds at zigzag edge , when brought into contact with another edge , can induce larger distortions than inert armchair edge ( dangling bonds at armchair are partly quenched by triple - bond formation ) . the steps in edge profiles near @xmath25 , moreover , are @xmath43   , whereas the steps in edge profiles near @xmath26 are only @xmath44   . this means that , to make bonds , the edge atoms near @xmath25 need more pulling than edge atoms near @xmath26 , causing the buckling . it is clear that the numbers in fig.[fig5 ] have no direct relevance as such , because our geometrically optimized gb samples are in vacuum , and measure only @xmath45    across the gb . we argue , however , that the inflection angle measures how much given gb would buckle in an experiment  large inflection angle meaning tendency to stick out ( sticking out would be bound by geometric constraints and by surface adhesion ) . for example , a gb with @xmath26 on a support will always remain smooth , whereas a gb with @xmath25 will either remain smooth or buckle , to be seen as a ridge or a mountain range in scanning probe experiments.@xcite indeed , this buckling tendency was seen in an experiment by ervenka _ at al . _ , measuring corrugation heights up to @xmath8    within gb regions@xcite . small corrugation heights and inflection angles are natural in samples on flat substrates , but also large inflection angles are realistic , for example in soot particles . @xcite     means planar structure ( occurring only with pristine graphene , @xmath19 and @xmath20 ) . the error bars remind of technical ambiguities in angle determination . the main trend here is that zigzag edges cause both large and small inflection , whereas armchair edges cause only small inflections . inset : inflection angle illustrated . ] now we pose the general question : how reactive are gbs ? we approach this question by examining hydrogen adsorption energies and dangling bonds ( db ) . to this end , we first develop a connection between hydrogen adsorption energy and the electronic structure given by dftb . chirality , including various polygons . b ) characterizing chemical reactivity from the dftb electronic structure directly : we calculated hydrogen adsorption energy for selected @xmath46 atoms as adsorption sites and plot for each atom adsorption energy and @xmath2 . given the correlation between the adsorption energy and @xmath2 , we access reactivity directly from the electronic structure , without additional calculations . the adsorption site numbers refer to panel a ) ; atom , having three @xmath47 angles has a small adsorption energy , atom has a dangling bond , and other atoms have strained environment . c ) using the correlation established in panels a ) and b ) , we calculate the averaged density of dangling bonds per unit length along grain boundary for all samples , as a function of @xmath27 . ] fig.[fig : reactivity]b shows the hydrogen adsorption energies for carbon atoms in @xmath48 representative gb samples , from total @xmath46 hydrogen adsorption calculations . it appears that adsorption to given carbon atom @xmath3 is strong if atom s cohesion @xmath49 decreases . carbon atoms part of regular hexagons have adsorption around @xmath50  ev ( in fig.[fig : reactivity]a ) but adsorption increases if atom is surrounded by other polygons and bonding angles deviate from @xmath51 ( , , and in fig.[fig : reactivity]a ) . the common denominator in these examples is that the change in hydrogen adsorption energy is caused by the strain in bond angles and in bond lengths . the strongest adsorptions around @xmath52  ev , in turn , are caused by dangling bonds that have @xmath53  ev ( in fig.[fig : reactivity ] ) ; such strong adsorption never occurs with three - coordinated atoms . ( we gave the argument about db energetics already in sec.[sec : dftb ] . ) we note that dftb hydrogen adsorption energy to zigzag ( @xmath54  ev ) , for example , agrees reasonably with dft energy ( @xmath55  ev)@xcite . therefore , we can characterize the reactivity _ directly _ by the dftb electronic structure using the quantities @xmath2 , without any adsorption calculations . surely , dangling bonds can be identified from geometry by defining criteria for coordination numbers , but this approach is prone to errors , especially in disordered regions that gbs are .    using @xmath53  ev as a criterion for a dangling bond , we then analyzed the reactivity for the whole gb ensemble . fig.[fig : reactivity]c shows the average number of dangling bonds per unit length in a gb with a given chiral angle . the highest density , one db per @xmath56   , occurs with @xmath57 . the lowest density , one db per @xmath58   , occurs with @xmath59 ; pristine graphene with no dbs becomes probable with @xmath25 and @xmath41 .    by analyzing the histogram in fig.[fig : reactivity]c another way , @xmath60  %  of the samples have no dbs , @xmath61  %  of the samples have one db every @xmath62   , and @xmath63  %  of the samples have one db every @xmath64   . these numbers agree with a recent scanning tunneling microscopy ( stm ) experiment . dangling bonds are highly localized states just below the fermi - level , and hence seen as a bump in constant - current stm with low bias . the stm images of ervenka _ et al . _ in ref . show periodic appearance of sharp peaks , with @xmath65    periodicity for @xmath66  %  of the samples and with @xmath67    periodicity for @xmath68  % of the samples . even if these periodicities should be caused by adsorbed impurities and not from bonds that dangle , it is still likely to be result from _ reactive _ sites within gbs  and answers the original question we posed in this section . the agreement with experiment gives confidence in the objectivity of the gb construction process . understanding the trends in defects and reactivity will hopefully help in the design of functionalized graphene compounds  @xcite . the structure of a gb , given its constituent polygons , affects directly on its vibrational spectrum , and gives experimentally complementary information . fig.[fig : vibr ] shows the projected vibrational density of states ( pvdos ) for the gbs as a function of chiral angle and wave number . pvdos was calculated by solving vibrations for the whole gb , by projecting the eigenmodes to atoms within gb area , and by renormalizing the modes  for all the @xmath10 gb samples . hence fig.[fig : vibr ] shows a lot of data in a compact form .    ) was calculated by solving vibrational spectrum from the dynamical matrix , and projecting the vibration eigenmodes on atoms within the grain boundary ; the contour plot is then gathered and smoothed from all the @xmath10 samples . intensity increases from blue to red . ] spectra show three main bands , two at low energies @xmath69  @xmath70 and @xmath71  @xmath70 , and one  the so - called g - band in raman spectroscopy  at high energy @xmath72  @xmath70 . the band at @xmath71  @xmath70 is steady across all @xmath27 s and can not be used for structural identification . the g - band , in turn , loses intensity and comes down some @xmath73  @xmath70 in energy around @xmath36 . this can be explained by structure . as discussed in sec.[sec : structure ] , around @xmath36 hexagons are at relative minimum , implying a non - uniform and less rigid structure , and causing floppiness in high - energy modes . the high - energy modes have bond stretching between nearest - neighbors and are hence sensitive to local structural changes . low - energy modes , again , are more collective and hence insensitive to local changes in structure , as long as they remain somehow `` graphitic '' . if fig.[fig : vibr ] and fig.[fig4 ] are compared carefully , one finds that the local intensity maxima of the g - band occur precisely for @xmath27 with local maxima in the abundance of hexagons . the intensity increase around @xmath74  @xmath70 and @xmath36 is mainly due to renormalization of pvdos . the observations above are qualitatively supported by earlier experiments . using raman spectroscopy , a g - band shift from @xmath75  @xmath70 to @xmath76  @xmath70 was reported by ferrari and robertson in ref .  ; the shift was identified to be a consequence of structural change from nanocrystalline graphite to amorphous phase . while still being mainly planar ( apart from inflection ) , gb zone around @xmath36 , such as fig.[fig : reactivity]a with @xmath77 , indeed appears amorphous . some defects , like singular stone - wales defect or reczag edge of graphene , have their own name because they are well defined . grain boundaries , on the contrary , are like snowflakes  there is no flake like another , but it s enough to know that they are roughly hexagonal , flat , small and cold . for the same reason it s enough to know how grain boundaries usually look and feel . knowing trends is valuable . we investigated gbs from different viewpoints , discovering trends with complementary information , measurable also experimentally . we investigated ( i ) geometry ( polygons and inflection angles ) measurable with transmission electron microscope or atomic probe microscope ; ( ii ) energy , that is manifested in geometry and thereby measurable ; ( iii ) reactivity , measurable with adsorption experiments ; ( iv ) vibrational properties , measurable with raman spectroscopy . we are confident that the trends are genuine , part because of agreements with earlier experiments , part because the trends all make intuitive sense : energy , inflection angles , reactivity and vibration trends make sense given the structure , and the structure trends make sense given the graphene edge mismatch . we hope the trends help to explain experiments  and also help simulations and experiments to design graphene structures for given functions . this work was supported by the academy of finland ( project @xmath78 the and the finnano mep consortium ) and finnish cultural foundation . computational resources were provided by the nanoscience center in the university of jyvskyl and by the finnish it center for science ( csc ) in espoo .
Abstract:

Answer: grain boundaries are topological defects that often have a disordered character . disorder implies that understanding general trends is more important than accurate investigations of individual grain boundaries . here we present trends in the grain boundaries of graphene . we use density - functional tight - binding method to calculate trends in energy , atomic structure ( polygon composition ) , chemical reactivity ( dangling bond density ) , corrugation heights ( inflection angles ) , and dynamical properties ( vibrations ) , as a function of lattice orientation mismatch . the observed trends and their mutual interrelations are plausibly explained by structure , and supported by past experiments .


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: our knowledge of the @xmath2-ray  sky has increased dramatically during this last decade , due principally to the @xmath2-ray  instruments on board the compton gamma ray observatory ( cgro)@xcite : egret , a spark chamber plus calorimeter instrument with @xmath2-ray  sensitivity in the energy interval 30 mev to 30 gev ; comptel , a compton telescope in the interval 0.1 to 30 mev ; batse , an omnidirectional x - ray and soft @xmath2-ray  `` burst '' detector consisting of large nai scintillators sensitive to 30 kev to 2 mev photons ( with smaller spectroscopic nai crystals for measurements up to 110 mev ) ; and osse , consisting of nai - csi phoswiches detecting photons of 0.1 to 10 mev . their observations have revolutionized our understanding of such extragalactic phenomena as blazars and gamma ray bursts ( grbs ) , as well those within our own galaxy , such as pulsars . a salient characteristic of these extragalactic sources is their extreme time variability . in the case of blazars , flare states are observed to occur in which the luminosities vary by over an order of magnitude within time - scales as small as hours . grbs are singular , extragalactic events due either to the the coalescence of massive compact objects or to the hypernovae of massive early stars , whose luminosities decay over periods of seconds to days depending upon the wavelength observed . for such objects , the importance of correlated , _ multiwavelength _ observations is paramount . by monitoring the relative time lags of blazar flare peaks at radio , optical , x - ray , and @xmath2-ray  wavelengths , for example , the properties of the source regions at the bases of agn ( active galactic nucleus ) jets are constrained . gamma ray observations of blazars by egret have been particularly crucial , as the luminosity of these objects peaks in the gev - to - tev part of the electromagnetic spectrum . the transient nature of these sources therefore argues for _ continuous _ coverage of the @xmath2-ray  sky .    until recently , observations with cgro and other space - based and ground - based telescopes have provided coverage of these and other sources up to the limiting sensitive energy of egret , approximately 30 gev . from there , a gap in our knowledge of the @xmath2-ray  sky spectrum has existed up to 200 - 300 gev , where ground - based @xmath2-ray  shower detectors presently have their energy thresholds . it is possible that within this gap there are novel features in the @xmath2-ray  sky , such as a gamma - ray line or continuum emission from postulated neutralino annihilation at the center of the galaxy . future instruments with sensivity in this unexplored region , such as agile @xcite , recently approved by the italian space agency ( asi ) , or the proposed glast @xcite payload , may uncover exciting new phenomena .    at the present time , however , our view of the @xmath2-ray  sky has diminished . with the effective turnoff of the egret @xmath2-ray  instrument on the cgro due to the nearly complete consumption of its spark chamber gas , there will soon be no operating instrument capable of observing the @xmath2-ray  sky in the energy interval @xmath3 to @xmath4 gev . ground - based @xmath2-ray  detectors , based on the atmospheric cerenkov technique ( act)@xcite , turn on at current energy thresholds of @xmath5200 - 500 gev . within the next several years energy thresholds for some act observatories are expected to go to as low as 20 to 50 gev , but lower thresholds than these are unlikely to be achieved due to the sizable effect of earth s magnetic field on the @xmath2-ray - induced air showers , and the lower cerenkov photon yield which must be detected against the night sky background . for @xmath2-ray  energies much lower than @xmath6 gev , then , spaced - based detectors are required . this observational gap for the energy window @xmath3 to @xmath6 gev will continue to exist for the next several years . eventually , this gap will be eliminated with the launch of a next - generation @xmath2-ray  satellite , such as glast @xcite , but such a mission is unlikely to occur before the year 2005 . also the agile satellite @xcite , to be launched in 2002 , will have a limited sensitivity above 50 gev .    in this paper we describe how the alpha magnetic spectrometer ( ams ) can largely fill this gap by acting as a @xmath2-ray  detector with sensitivity in the energy interval of 0.3 to 100 gev during its three - year mission on board the international space station alpha ( issa ) from 2003 to 2006 . ams , described in detail elsewhere , has as its primary mission the search for cosmic ray antinuclei , which are detected via the negative curvature of their trajectories through a magnetic spectrometer ( see sect . [ ams.sec ] ) . it is of interest , however , to analyze if a large particle detector like ams , operated on the space station for three years starting in 2003 , could be used to detect high energy @xmath2-rays  . there are at least two ways this could be done without significantly affecting the experiment sensitivity to antimatter : a ) adding a light ( e.g. @xmath7 ) converter at the entrance of the magnetic spectrometer , either passive ( e.g. , a high - z thin plate ) or active ( e.g. a multilayer tracking detector ) , or b ) implementing an high granularity imaging shower detector at the bottom of the experiment . option ( a ) consists in a minor modification to the baseline instrument , while option ( b ) , although would give a much higher @xmath2-ray  detection capabilities , requires the addition of a new detector . in this paper we present a study which has been perfomed on option ( a ) . we will show that in this option , ams can also detect @xmath2-rays  with performance characteristics similar to egret in the energy region of 0.3 to 20 gev , and with significantly enhanced capabilities between 20 and @xmath5 100 gev , a region which is not well unexplored . we refer to this modified instrument as `` ams/@xmath2 '' . we show that ams/@xmath2  can continue the valuable work of egret by providing continued monitoring of extragalactic and galactic @xmath2-ray  sources and by participating in multiwavelength observational campaigns . in addition , ams/@xmath2  will have unprecedented sensitivity to the @xmath2-ray  sky between the energies of 20 to @xmath8 gev ( albeit at a level somewhat lower than required for detection of known sources with power law spectra ) , so that ams/@xmath2  might provide us with unexpected discoveries in this region . the next section gives a brief description of the baseline ams instrument , and the modifications required for adding @xmath2-ray  detection capability in option ( a ) . section 3 describes the performance characteristics of ams/@xmath2  determined from monte carlo analyses , and section 4 uses these results to address the question of what @xmath2-ray  astrophysics can be done with ams . the results obtained analyzing option ( a ) give a hint of the physics potential which could be reached implementing option ( b ) for the flight of ams on the space station ; option ( b ) would gain a factor of about @xmath9 in @xmath2-raysstatistics over option ( a ) . detailed analysis of option ( b ) , however , wil be the subject of future work . the alpha magnetic spectrometer has been built by a large international collaboration of high energy physics institutions from the u.s . , italy , china , finland , france , germany , taiwan , russia , and switzerland . it recently had a successful test flight on the space shuttle mission sts-91 in june , 1998 @xcite , where it was carried in the cargo bay and observed for several days in both the zenith and nadir directions , the latter for measurement of albedo cosmic ray backgrounds . it is scheduled to be secured to an external payload attachment point on the international space station alpha in 2003 , where it will remain as a zenith - pointing instrument for a minimum of 3 years . the basic design of ams is shown in figure [ ams-nim.fig ] . the magnet spectrometer consists of a permanent ring dipole magnet made of very high grade nd - b - fe rare earth material whose magnetic energy product and residual induction are respectively @xmath10 g - oe and 14,500 g , yielding an analyzing power of @xmath11 t - m@xmath1 with less than 2 tonnes of magnet mass . four silicon strip detector tracking planes are located within the magnetic volume , with a fifth and sixth plane located just above and below the magnet . these are based on technology developed for the aleph and l3 vertex detectors at lep @xcite , and consist of @xmath12 m thick , double - sided , high purity silicon wafers . the total area of this tracker is about an order of magnitude larger than the size of previous microstrip silicon detectors . the momentum resolution of the silicon spectrometer is given in figure [ y94725.fig ] at low rigidities ( below 8 gv ) the resolution is dominated by multiple scattering ( @xmath13 ) , while the maximum detectable rigidity ( @xmath14 ) is about @xmath15 gv . the parameters of the silicon spectrometer are given in table 1 . in addition to measuring particle rigidity , the silicon planes will provide six independent measurements of @xmath16 for charge determination . .ams silicon tracker parameters ( in parenthesis the values used on the precursor flight ) . [ cols="^,^ " , ]     in models 2d and 4d , the density profile is fixed to follow the `` universal '' profile found by nfw , for which @xmath17 and @xmath18 . the last halo profile is an approximation to the true isothermal profile ( see binney & tremaine , 1987 ) : @xmath19 figure [ jpsi-nim.fig ] shows the resulting @xmath20 curves for the various halo profiles . the remaining two parameters that determine the @xmath2-ray  flux are , from eq . [ dfdo2.eq ] , the neutralino mass @xmath21 and the annihilation cross sections @xmath22 that result in @xmath2-ray  production . monochromatic photons from @xmath23 and @xmath24 are necessarily produced via box diagrams @xcite , and so have a smaller branching ratio that tree - level annihilation into @xmath25 and @xmath26 pairs . secondary @xmath2-rays  from the decay of tree - level @xmath25 s and @xmath26 s give a much larger continuum flux , which , depending on the susy parameters , may be well above the extragalactic @xmath2-ray  background flux ( section [ egrb.sec ] ) near the galactic center @xcite . monochromatic @xmath2-ray  detection is limited both by detector energy resolution ( fig . [ energyres.fig ] ) and aperture . should the expected event rate be significantly greater than unity , the line width must be small enough so that the signal is not lost in the galactic and extragalactic @xmath2-ray  backgrounds . the rapid deterioration of energy resolution above 100 gev limits ams detection of _ monochromatic _ @xmath2-rays  to neutralino masses @xmath27 gev . ( this limit does not apply to detection of continuum photons from @xmath25 and @xmath26 decays ; see section [ continuum.sec ] . ) however , recent analysis by the dama / nai collaboration has confirmed the presence of an annual modulation in their detector s event rate ( at 98.5% c.l . ) , suggesting the existence of a local wimp density with a @xmath28 mass range @xmath29 gev , and wimp - nucleon scalar elastic cross section @xmath30 nb , where @xmath31 is the fraction of the local dark matter density that is due to wimp neutralinos @xcite . within this mass range our rms energy resolution is flat at 24% . these mass and cross section ranges have also been shown to be fully compatible with a mssm neutralino being the major component of dark matter in the universe @xcite . figure [ mono-nim.fig ] illustrates the sensitivity of ams to detection of line photons from the annihilation channel @xmath32 as a function of halo model . plotted are integrated @xmath2-ray  counts @xmath33 as a function of angle @xmath34 from the galactic center for the four halo profile models , accumulated over a two - year period by ams , assuming @xmath35 gev and a velocity - cross section product @xmath36 @xmath37/s . with an energy acceptance width of 20 gev , approximately 80% of the line photons are captured . the dot - dashed line shows the galactic and extragalactic background count for this energy window ; the parameterization of bergstrm , ullio , and buckley @xcite was used for the galactic background contribution . the dashed and dotted lines correspond respectively to the `` berezinsky - like '' halo ( model 2f of dehnen and binney @xcite ) and the quasi - isothermal halo , while the upper and lower solid lines correspond to models 4d and 2d of dehnen and binney , which are the nfw - like models . bergstrm and ullio @xcite have recently performed a full one - loop calculation of the @xmath23 annihilation reaction , assuming the minimal supersymmetric extension of the standard model ( mssm ) . they performed extensive scans in susy parameter space to map out the range of values of the annihilation cross section versus neutralino mass . we have reproduced their figure 3a in figure [ twogamma-nim.fig ] , along with the ams sensitivity for the two most optimistic halo profiles , models 2f and 4d , expressed as a lower limit on @xmath38 versus @xmath21 for detection of neutralino annihilation . the detection criteria were arbitrarily chosen to be @xmath39 and @xmath40 , where @xmath41 is the total number of annihilation @xmath2-rays  detected over a 2-year operating period within the central few degrees of the galactic center , and @xmath42 is the total number of background photons from diffuse galactic emission detected within the energy window and the same solid angle about the galactic center . ullio and bergstrm @xcite also performed a mssm full one - loop calculation of neutralino annihilation via @xmath43 , assuming a pure higgsino state for the neutralino . for this channel the @xmath2-ray  has an energy @xmath44 . in figure [ onegamma-nim.fig ] we have reproduced figure 6 of ullio and bergstrm , which shows their @xmath45 versus higgsino mass ( solid line ) , along with the unitary lower bound on this cross section ( dashed line ) . superimposed on this plot is a line segment which shows the minimum @xmath45 required for detectability by ams over a 2-year period , assuming the `` berezinsky - like '' halo model 2f of dehnen and binney @xcite . none of the other halo models intersect the ullio and bergstrm curves . thus for this pure higgsino annihilation channel only the most singular halo profile can yield a detectable signal .    in general , then , only the more singular of the halo profile models , in which the dark matter density continues to increase as @xmath46 , @xmath47 , can give hope for line photon detection by ams , and then only in a restricted region of mssm parameter space . it should be kept in mind that a fraction of this restricted region may already have been eliminated by the absence of a detection by egret , but such an analysis has yet to be completed @xcite . continuum photons from the decay of @xmath25s and @xmath48s produced by tree - level annihilations @xmath49 are far more abundant than the line @xmath2-rays  considered above , and so possibly provide the greatest chance for detection , assuming the neutralino mass is above that of the @xmath25 and @xmath48 . berstrm , ullio , and buckley @xcite have used the pythia monte carlo code @xcite to determine the the continuum photon energy distribution from @xmath50 annihilation into @xmath51 , @xmath52 , and @xmath53 pairs . they find that the hardest distributions come from @xmath51 and @xmath52 pairs , and that the annihilation rate is dominated by @xmath51 and @xmath52 production at higher ( higgsino - like ) neutralino masses . from a scan of mssm parameter space , they find that the models with the highest rate of @xmath51 and @xmath52 production yield an observed differential photon flux of @xmath54 where @xmath55 for @xmath56 gev . the @xmath20 factor is given by eq . [ jpsi.eq ] .    using the @xmath20 functions of the four halo models described above ( fig . [ jpsi-nim.fig ] ) , the differential flux coming from a solid angle @xmath57 about the galactic center can be obtained and compared to the flux from the galactic and extragalactic backgrounds . using the collection area*time product @xmath58 ( fig . [ areatime.fig ] ) for ams , the differential event count is obtained . figures [ cont.300.1-nim.fig ] to [ cont.300.32-nim.fig ] show these event counts after integration over @xmath34 out to three different solid angles about the galactic center , @xmath59 , @xmath60 , and 1.0 steradians . the dot - dashed line is the diffuse galactic and extragalactic background photon count ; the upper and lower solid lines correspond to halos for models 4d and 2d , respectively ; the dashed line corresponds to the `` berezinsky - like '' halo model 2f ; and the dotted line comes from the quasi - isothermal halo model . the faux - data points shown in the figures take into account the effect of detector energy resolution ; these are the counts obtained in each energy bin after convolution with the monte carlo energy resolution functions ( fig . [ energy.fig ] ) . there are a number of salient features , the primary one being the strength of the signal compared to the diffuse background for essentially all of the halo models considered . for this particular ( optimized ) mssm model it is likely that egret would have already provided a detection given any of the halo models considered here . although there is a large variation in signal between the various halo models , it appears nevertheless that for _ any _ of these halo models a substantial region of mssm parameter space can be constrained by the absence of a continuum @xmath2-ray  signal . these figures also show that the differing angular dependences of the background and signal result in there being an optimal solid angle @xmath57 that maximizes signal to noise , with more singular profiles having a stronger signal to noise at smaller solid angles . these differing angular dependences can also be used to discriminate the signal from the diffuse backgrounds . finally , we point out that continuum photons result from @xmath61 as well as from @xmath51 and @xmath52 pair production , producing scaled spectra similar to those from @xmath25 and @xmath48 decays . this implies that ams/@xmath2  sensitivity to continuum @xmath2-rays   may extend down to the current lep limits on the neutralino mass @xcite . we have shown that with minor modifications the alpha magnetic spectrometer can become a powerful @xmath2-ray  detector as well , with overall performance characteristics being comparable , if not superior , to those of egret . with @xmath2-ray  energy resolution extending past 100 gev , and with an aperture that is nearly flat above @xmath62 gev , ams/@xmath2  can address a number of outstanding issues in @xmath2-ray  astrophysics that relate to the relatively unexplored region of @xmath63 gev . for one , ams/@xmath2  will likely confirm or refute the hypothesis that unresolved blazars are responsible for the bulk of the extragalactic @xmath2-ray  background ; ams/@xmath2  will also extend the spectrum of the diffuse galactic background to above 100 gev , helping to resolve current difficulties in interpreting the egret diffuse galactic background measurement @xcite . ams/@xmath2  should roughly double the total number of blazars detected in @xmath2-rays , and will be enable multiwavelength observational campaigns to include the gev region of blazar spectra during the flight years of 2003 - 2006 . there is an additional possibility that an indirect detection of the cosmic uv and optical photon background can be made through the detection of extinctions in high - redshift blazars above @xmath64 gev . ams/@xmath2  will also likely observe gev @xmath2-ray  emission from one or more gamma ray bursts during its operational lifetime . ams/@xmath2  will also search for both line and continuum emission of @xmath2-rays  from the region of the galactic center created by the annihilation of dark matter neutralinos . although the sensitivity to line emission appears marginal , there is nevertheless a finite , though small , region of halo / mssm phase space which allows a detection by ams . however , a much larger region of dark - matter halo / mssm parameter space can be constrained in a search for continuum @xmath2-rays  by ams/@xmath2 . we thank b. dingus , s. hunter , and p. ullio for helpful conversations , and p. sreekumar for the same and for supplying us with additional , unpublished information on the egret measurement of the extragalactic @xmath2-ray  background . we also thank the many members of the ams collaboration who have indirectly contributed to this paper , and especially s.c.c . ting for his critical comments . kurfess , d.l . bertsch , g.j . fishman , and v. schfelder ( 1997 ) aip conf . 410 , 509 . d. engovatov et al . ( 1997 ) ieee nucl . ( albuquerque ) . agile , phase a report , asi small program scientific mission , october 1998 , ar.das.98.501 . a.m. hillas and j.r . patterson ( 1990 ) j. phys . 16 , 1271 . m. acciarri et al . ( 1995 ) nucl . instr . meth . a360 , 103 ; r. battiston ( 1995 ) nucl . s. ahlen et al . ( 1994 ) nucl . instr . meth . a350 , 351 . g. ambrosi et al . ( 1998 ) , submitted to nucl . wind ( 1974 ) nucl . 115 , 431 . r. mukherjee , _ et al . _  ( 1997 ) ap.j . 490 , 116 . p. sreekumar , _ et al . _  ( 1998 ) ap.j . 494 , 523 . et al . _  ( 1997 ) ap.j . 481 , 205 . thompson , _ et al . _  ( 1993 ) ap.j.suppl . jim scheib ( nasa / jsc ) ( 1995 ) , personal communication . t.p . li and y.q . ma ( 1983 ) 272 , 317 . swanenberg , _ et al . _  ( 1978 ) nature 275 , 298 . g. ghisellini , _ et al . _  ( 1998 ) astro - ph/980731 . g. henri and g. pelletier ( 1991 ) ap.j . r.d . blandford and a. knigl ( 1979 ) ap.j . 232 , 34 . c. von montigny , _ et al . _  ( 1995 ) ap.j . 440 , 525 . r. sambruna , l. maraschi , and c.m . urry ( 1996 ) ap.j . 463 , 444 . stecker , o.c . de jager , and m.h . salamon ( 1996 ) ap . j. 473 , l75 . stecker , o.c . de jager , and m.h . salamon ( 1992 ) ap.j . 390 , l49 . mattox , _ et al . _  ( 1997 ) ap.j . hartman , w. collmar , c. von montigny , and c.d . dermer ( 1997 ) aip conf . 410 , 307 . thompson , _ et al . _  ( 1995 ) ap.j.suppl . 101 , 259 . stecker and m.h . salamon ( 1996 ) ap.j . 464 , 600 . p. padovani , g. ghisellini , a.c . fabian , and a. celotti ( 1993 ) mnras 260 , l21 . j. chiang and r. mukherjee ( 1998 ) ap.j . 496 , 752 . dunlop and j.a . peacock ( 1990 ) mnras 247 , 19 . dermer ( 1995 ) ap.j . 446 , l86 . kerrick , _ et al . _  ( 1995 ) ap.j . 452 , 588 . m.h . salamon and f.w . stecker ( 1998 ) ap.j . 493 , 547 . m. sikora , m.c . begelman , and m.j . rees ( 1994 ) ap.j . 421 , 153 . gould and g.p . schreder ( 1966 ) phys . m.h . salamon and f.w . stecker ( 1997 ) ap.j . 493 , 547 . salamon , f.w . stecker , and o.c . de jager ( 1994 ) ap.j . k. hurley ( 1998 ) astro - ph/9812052 . t. piran ( 1998 ) astro - ph/9801001 . k. hurley et al . ( 1994 ) nature 372 , 652 . b. dingus et al . ( 1997 ) in _ proc . 4th huntsville symposium on grbs_. e. waxman ( 1995 ) ap.j . m.g . baring and a.k . harding ( 1997 ) ap.j . 491 , 663 . r. plaga ( 1995 ) nature 374 , 430 . thompson and c.e . fichtel ( 1982 ) a&a 109 , 352 . a. dar and n shariv ( 1995 ) phys . 75 , 3052 . f.w . stecker and m.h . salamon ( 1996 ) phys . 76 , 3878 . bignami , c.e . fichtel , r.c . hartman , and d.j . thompson ( 1979 ) ap.j . 232 , 649 . d. kazanas and j.p . protheroe ( 1983 ) nature 302 , 228 . d. mller and k.k . tang ( 1987 ) ap.j . 312 , 183 . et al . _  ( 1994 ) ap.j . 436 , 769 . barwick , _ et al . _  ( 1997 ) ap.j . . w. menn , _ et al . _  ( 1997 ) 25th icrc 3 , 409 ( dublin ) . g. alverson , _ et al . _  ( 1993 ) phys . m. mori ( 1997 ) ap.j . 478 , 225 . c. alcock _ et al . _  ( 1997 ) ap.j . 486 , 697 . w. sutherland ( 1998 ) astro - ph/9811185 . m. fukugita , c.j . hogan , and p.j.e . peebles ( 1998 ) ap.j . 503 , 518 . g. jungman , m. kamionkowski , and k. griest ( 1996 ) phys . 267 , 195 . m.s . turner ( 1986 ) phys . d34 , 1921 . reid ( 1993 ) ann . m. fich and s. tremaine ( 1991 ) ann . rev . w. dehnen and j. binney ( 1996 ) astro - ph/9612059 . j. binney and s. tremaine ( 1987 ) _ galactic dynamics _ ( princeton ) . navarro , c.s . frenk , and s.d.m . white ( 1996 ) ap.j . 462 , 563 . navarro , c.s . frenk , and s.d.m . white ( 1997 ) ap.j . 490 , 493 . s. cole and c. lacey ( 1995 ) astro - ph/9510147 . berezinsky , a.v . gurevich , and k.p . zybin ( 1992 ) phys . b294 , 221 . flores and j.r . primack ( 1994 ) ap.j . a.v.kratsov , a.a . klypin , j.s . bullock , and j.r . primack ( 1998 ) ap.j . p. ullio and l. bergstrm ( 1997 ) hep - ph/9707333 . l. bergstrm and p. ullio ( 1997 ) hep - ph/9706232 . r. bernabei et al . ( 1998 ) roma ii university preprints roma2f/98/27 and roma2f/98/34 . a. bottino , f. donato , n. fornengo , and s. scopel ( 1998 ) hep - ph/9808456 . l. bergstrm , p. ullio , and j.h . buckley ( 1997 ) astro - ph/9712318 . b.l . dingus and m.h . salamon ( 1998 ) in preparation . t. sjstrand ( 1994 ) comp . m. biasini ( 1999 ) in proc . xiii rencontres de physique ( la thuile ) . m. pohl and j.a . esposito ( 1998 ) ap . j. 507 , 327 . r. battiston ( 1999 ) in proc . xiii rencontres de physique ( la thuile ) ; also as astro - ph/9907152 . begeman , a.h . broeils , and r.h . sanders ( 1991 ) mnras 249 , 523 . e.d . bloom ( 1996 ) sp . sci .
Abstract:

Answer: the modeled performance of the alpha magnetic spectrometer ( ams ) as a high - energy ( 0.3 to 100 gev ) gamma - ray detector is described , and its gamma - ray astrophysics objectives are discussed .    * the gamma ray detection capabilities of the alpha magnetic spectrometer *    * r. battiston@xmath0 , m. biasini@xmath0 , e. fiandrini@xmath0 , + j. petrakis@xmath1 , and m.h . salamon@xmath1 + *    @xmath0sezione infn and dipartimento di fisica , universit di perugia + @xmath1department of physics , university of utah +    * pacs , keywords : * 95.55.k , 98.54.c , 95.35.+d , 98.35.g ; high - energy gamma rays , blazars , neutralinos , gamma - ray background


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: astrophysical discs are observed to transport angular momentum . it has been hypothesised that such discs are transporting angular momentum through a turbulent process @xcite . despite decades of theoretical research we still lack a sufficient understanding of turbulence to make quantitative predictions . several candidate processes exist to sustain turbulence in these discs : the magneto - rotational instability ( mri ) ( e.g. * ? ? ? * ) , gravitational instability ( e.g. lodato & rice 2004 , 2005 ) the vertical shear instability @xcite and the baroclinc instability ( e.g. * ? ? ? the mri is still the leading candidate for angular momentum transport in accretion discs ( see @xcite for a recent review ) . however , at the low temperatures and ionization fractions expected in protoplanetary discs , non - ideal mhd effects become important , qualitatively changing the nature of the turbulence and its associated transport properties , or rendering it ineffective resulting in ` dead - zones ' ( e.g. * ? ? ? a corollary to the angular momentum transport problem in protoplanetary discs is the transport of dust particles . there is large amounts of astrophysical ( e.g * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) and cosmochemical ( e.g * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) evidence to suggest large scale radial transport of dust particles occurs . in particular , the level of diffusion that results from the turbulence is an unknown parameter . the turbulent diffusion coefficient ( @xmath2 ) is often parametrised in terms of the turbulent kinematic viscosity ( @xmath3 ) as @xmath4 where @xmath5 is the schmidt number . assumptions of isotropic , kolmogorov - like turbulence lead to the inference that @xmath6 @xcite ; however , a large range of reported experimental / theoretical / numerical values exist that span approximately two orders of magnitude in the range @xmath7 @xcite . direct observational measurements of the strength of turbulent diffusion in the continuum will be complicated by dust - drag , grain growth / fragmentation and optical depth effects . however , it is not only dust particles that will experience turbulent diffusion ; any rare , gas - phase species ( tracer species ) where changes in temperature results in a concentration gradient will also experience turbulent diffusion . if , for example , a tracer species is predominately produced at a given radius , the tracer will then diffuse away from this radius , with a distribution that is strongly dependent on the value of the schmidt number @xcite . @xcite considered the case of deuterated water and argued the d / h water distribution in chondrites implies the schmidt number is smaller than unity . snow - lines , where a gas phase tracer ( for example h@xmath0o , co@xmath0 , co ) condenses out of the gas to form ice below some temperature represents a scenario where a sharp concentration gradient can occur . species such as h@xmath0o , co@xmath0 & co are relatively abundant , such that the surface layers can be optically thick . additionally , model degeneracies make it difficult to constrain snow - lines directly ( even with optically thin isotopologues ) . however , if the production / destruction of yet rarer tracers are regulated by the presence or absence of gas phase h@xmath0o , co@xmath0 or co then these rare tracers could be used as proxies to detect snow - lines . for the co snow - line it is expected that n@xmath0h@xmath1 and h@xmath0co will only be abundant when co freezes out @xcite . such an expectation has been born out in observations of star - forming cores @xcite , and similar results were obtained in the discs _ sma _ survey of several nearby protoplanetary discs @xcite . recently , @xcite imaged a hole in n@xmath0h@xmath1 using _ alma _ at a radius of @xmath8 au , co - incident with the expected location of the co snow - line ( based on a free - out temperature of @xmath9  k ) . the sharpness of the inner edge will depend strongly on the strength of the turbulent diffusion , with weaker diffusion resulting in a sharper hole .    in this letter we demonstrate how the distribution of a tracer species that is only abundant outside a snow - line ( being destroyed inside ) is strongly dependant on the schmidt number , and that _ alma _ observations could be able to constrain the schmidt number in protoplanetary discs . we consider a 1d axis - symmetric disc , where the evolution of the gas surface density ( @xmath10 ) and the surface density of any ( gas - phase ) tracer species ( @xmath11 ) is given by ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) : @xmath12\label{eqn : gas_disc}\\ \frac{\partial \sigma_i}{\partial t}&+&\frac{1}{r}\frac{\partial}{\partial r}\left [ r\sigma_i u_g - d^g_r r\sigma_g\frac{\partial x_i}{\partial r}\right]=\nonumber\\ & & \sum_js_i(\sigma_g,\sigma_j , r , t)\label{eqn : tracer_disc}\end{aligned}\ ] ] where @xmath13 is the concentration of the tracer species , @xmath14 is the net radial gas velocity , @xmath15 is the radial gas turbulent diffusion co - efficient and @xmath16 is a source / sink term that represents the production and destruction of the tracer species .      the following general model can be applied to any snow - line ( h@xmath0o , co@xmath0 , co , etc . ) where the chemical and transport time - scales decouple . here we specifically consider the case of n@xmath0h@xmath1 destruction at the co snow - line as observed in tw hya @xcite . inside the co snow - line n@xmath0h@xmath1 is destroyed by gas phase co , outside the co snow line this destruction channel is no - longer dominant and it instead is destroyed at a much slower rate by dissociative recombination @xcite . simulations without transport suggest the n@xmath0h@xmath1 abundance drops by several orders of magnitude inside the co snow - line ( e.g. * ? ? ? * ) with an abundance that depends on the square of the gas - phase co abundance @xcite .      in order for chemical tracers at the snow - line to be useful in terms of probing the strength of the turbulent diffusion , we must de - couple the chemical and dynamical time - scales . namely , the desorption time - scale must be faster than the transport time - scales in - order to create a sharp snow - line ; furthermore , the destruction time - scale of the tracer species inside the snow - line must also be fast . the transport time - scales of interest are the time with which to move a radial distance @xmath17 ( where @xmath17 is the disc s scale height , not @xmath18 . ] which is of order the radial scale length ) . therefore , the advection time - scale ( @xmath19 ) is : @xmath20 and the diffusive time - scale ( @xmath21 ) is : @xmath22 it is well known ( e.g. * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) for a steady disc model , but it is a more general result - see clarke & pringle ( 1988 ) . ] the logarithmic concentration gradient rapidly approaches @xmath23 , thus the diffusive time - scale is identical to the advection time - scale ( this is somewhat unsurprising as they are both driven by the same process ) . thus , the relevant time - scale for the movement of a individual tracer molecule over a radial scale @xmath17 is ( @xmath24 ) @xmath25 years . we want to compare this transport time - scale to the time - scale for the desorption of the snow - line species , along with the destruction of the tracer species inside the snow - line . considering our example of n@xmath0h@xmath1 and the co snow - line then the desorption time - scale is obtained by balancing desorption with absorption ( with rate constant @xmath26 ) , so the desorption time - scale becomes @xcite : @xmath27 where @xmath28 is the dust - grain size , @xmath29 is the dust - to - gas mass ratio and @xmath30 is the co abundance . additionally , the destruction time - scale for n@xmath0h@xmath1 is ( from the @xcite simplified network ) : @xmath31 where @xmath32 is the destruction rate constant . therefore , the co / n@xmath0h@xmath1 system clearly satisfies @xmath33 , for conditions experienced in protoplanetary discs . as such we may ignore the details of the chemical rate equations and simply model the destruction of any remaining n@xmath0h@xmath1 to occur instantaneously at the snow - line radius ; although , we emphasise that the following analysis can be applied to any tracer species with similar destruction time - scales . therefore , in this situation the source function @xmath34 is drastically simplified to : @xmath35 where @xmath36 is the mass - accretion rate , @xmath37 is the concentration at large radius , @xmath38 is the dirac delta - function and @xmath39 is the radius of the snow - line . this source function represents the instantaneous destruction of any remaining tracer species at @xmath40 . we note it ignores any possible vertical structure of the snow - line , which will necessarily spread out the destruction region of the tracer species @xcite and we discuss the implications of this caveat in section  4 . we will now restrict ourselves to a steady disc problem . in that case equation [ eqn : gas_disc ] becomes @xmath41 : where @xmath42 is the accretion rate , and we have neglected the very small contribution due to an inner boundary at finite radius . furthermore , the equation for the gas tracer becomes : @xmath43=r s_i(r)\label{eqn : tracer_steady}\ ] ] where @xmath44 is the radial schmidt number . therefore , radially integrating equation  [ eqn : tracer_steady ] and using our expression for @xmath36 we find : @xmath45 using equation  [ eqn : source ] we can integrate equation  [ eqn : conc_steady_intm ] to find the radial concentration distribution : @xmath46 where @xmath47 is the heaviside step function . setting @xmath48 at @xmath40 we find the solution : @xmath49 therefore , we see that the radial profile of the concentration is strongly sensitive to the value of the schmidt number . in figure [ fig : conc ] we show how the concentration varies with radius and schmidt number . it is important to emphasise that concentration distribution is independent of assumptions of the ( unknown ) viscosity , and mass - accretion rate and can in principle provide a ` clean ' measurement of the schmidt number .     in the range 0.1 - 10 . ] unfortunately it is not possible to directly observe the concentration gradient . what is directly observed is the surface - brightness distribution of the relevant species . the surface - brightness distribution is sensitive to the surface density of the tracer species rather than its concentration . thus , we must multiply the concentration of our tracer species by the gas surface density . adopting a power - law gas distribution of the form @xmath50 with a cut - off radius @xmath51 , the surface density of the tracer species is : @xmath52 \label{eqn : tracer_surf}\ ] ] in the range @xmath53 and @xmath54 elsewhere . the actual brightness distribution from rotational emission lines ( such as n@xmath0h@xmath1 j=4 - 3 ) , as well as being sensitive to the surface density distribution is also sensitive to the background gas temperature and excitation temperature of the molecule ( which in turn is a function of density and temperature ) along with the optical depth . however , in the case that the tracer species is optically thin , and the gas density is far above the critical density ( @xmath55 ) for the rotational transition , then the molecular line is thermalised . thus , if the temperature gradient is weak compared to the scales of interest and @xmath56 then we can approximate the surface brightness as being directly proportional to the surface density of the tracer species , provided the line remains optically thin .    for the gas at 30  au the gas density is typically : @xmath57 where @xmath58 is the mean molecular weight of the gas . comparing this to the critical density of the @xmath59 transition of n@xmath0h@xmath1 which has a critical density of @xmath60 @xmath61 ( e.g. * ? ? ? * ) , we see that @xmath56 . we note , since temperature and density are expected to be power - laws with radius ( e.g. * ? ? ? * ; * ? ? ? * ) then the small additional corrections due to temperature and density effects of converting @xmath11 to the surface brightness distribution ( @xmath62 ) will manifest themselves as changes in the power - law index @xmath63 in equation  [ eqn : tracer_brightness ] .    assuming the disc to be axis - symmetric and observed face - on , we may write the brightness distribution on the sky as : @xmath64 \label{eqn : tracer_brightness}\ ] ] in the range @xmath65 and @xmath66 elsewhere , where @xmath67 is a constant , @xmath68 is the angular size on the sky and @xmath69 is the angular size of the snow - line given by : @xmath70    thus , we see that the brightness distribution still retains the strong sensitivity to the schmidt number . since the angular resolution required to probe the value of the schmidt number is only available through mm - interferometry , we can use our brightness distribution to calculate synthetic visibilities . the visibilities can obtained by a hankel transform of the axis - symmetric brightness - distribution such that : @xmath71 where @xmath72 is a radial baseline co - ordinate defined as @xmath73 where @xmath74 & @xmath75 are the usual baseline co - ordinates . following the observed n@xmath0h@xmath1 @xmath59 emission by @xcite we calculate our synthetic observations at 372ghz with our standard model adopting the best fit parameters from @xcite of @xmath76 au , @xmath77 au and @xmath78 . we calculate our synthetic observations assuming the source is observed face - on at a distance of 150pc . in figure [ fig : brightness_vis ] we show our simulated visibility curves in the left - hand panels and the surface density profiles in the right hand panel . in the top panels we vary the schmidt number between 0.1 - 10 , in the middle panels we vary @xmath63 from 1 - 3 and in the bottom panels we vary @xmath39 from 15 - 60 au . the visibility curves in figure  [ fig : brightness_vis ] clearly show that variations in the schmidt number give rise to significant differences . furthermore , comparisons between varying the index of the gas surface density @xmath63 , snow - line radius and the schmidt number are not too degenerate . the simplified model presented here contains four free parameters @xmath79 , which would all need to be constrained by fitting the visibilities . inspection of figure  [ fig : brightness_vis ] suggests that sensitivities @xmath80 % at base - lines @xmath81 km would allow the theoretically suggested range of schmidt numbers ( 0.1 - 10 ) to be observationally constrained . we have shown that observations of snow - lines in protoplanetary discs using a tracer species ( for example n@xmath0h@xmath1 in the case of the co snow line @xcite ) can be used to probe the schmidt number , a unknown parameter in studies of turbulent transport in accretion discs , where current estimates span a range of two - orders of magnitude ( @xmath82 0.1 - 10 ) . the _ alma _ telescope is a sub - mm / mm interferometer ; once completed it will have @xmath83 individual antennas , offering @xmath84 baselines with separations upto 12  km . at 372 ghz ( @xmath85 n@xmath0h@xmath1 j=4 - 3 line ) this provides a maximum spatial resolution of @xmath86 arcsec . figure  [ fig : brightness_vis ] clearly shows that _ alma _ posses the required number of base - lines ( @xmath87 ) with separations in the range 0.1 - 1 km to constrain schmidt values within the current range of uncertainty , provided the observations are sensitive enough . taking the tw hya n@xmath0h@xmath1 observations as reference @xcite - a source brightness @xmath88 mjy beam@xmath89 km s@xmath89 with an rms noise of 8.1 mjy beam@xmath89 km s@xmath89 and beam size @xmath90 arcsec made using 23 - 26 antennas then using the current observation to constrain the schmidt number seems unlikely ; however , using the velocity channels separately can be used to increase the effective resolution ( e.g. * ? ? ? ] - a similar level of sensitivity , but with a beam size of @xmath91 arcsec , could be obtained assuming a fully operational _ alma _ of @xmath8550 antennas . thus , such observations are feasible and sufficiently high resolution to allow constraints to be placed on the schmidt number at the distance to tw hya . longer integration times and larger base - lines would be required to reach similar levels of sensitivity at distances of 150pc . we have argued that the sharpness of the hole in tracer species at the snow - line probes the value of the schmidt number independent of the assumed properties of the turbulence ( e.g. assumed value of the viscous ` @xmath92 ' parameter ) . since several snow - lines are expected to occur at different radii , then measurements at different snow - lines would allow the radial dependence of the schmidt number to be probed . in particular , recent simulation suggest that different non - ideal mhd effects ( which dominate at different radii @xcite ) lead to different schmidt numbers @xcite . thus , comparing the simulation predictions of the schmidt number for various turbulent driving mechanisms with the observed value would allow inferences about the nature of the turbulence to be made . furthermore , independently measuring the rms turbulent velocity ( @xmath93 ) ( e.g. * ? ? ? * ) then combining it with a measurement of the schmidt number would allow the viscosity ( including estimates of @xmath94 ) to be calculated , since @xmath95 . we have constructed an idealised model to investigate whether snow - lines could begin to probe the strength of turbulent diffusion . as such there are several model improvements that must be made before fitting to real data . therefore , our model presented in this letter is a ` proof - of - concept ' rather than a road map for observational modelling . for example , a real protoplanetary disc is not one - dimensional . as such the vertical temperature structure is not constant and passively heated discs cool as one approaches the mid - plane @xcite . therefore , the snow - line is unlikely to occur exactly at a fixed radius , but is more - likely to be an extended structure with a scale variation of @xmath96 , with is time - varying position @xcite . additionally the conversion of the gas phase to ice particles at the snow - line will result in turbulent diffusion of the gas and ice particles away from the snow - line ( in identical manner to that discussed for the snow - line tracer discussed here ) . as such , there is unlikely to be a very sharp change in the gas abundance at the snow - line but rather a smoother change . furthermore , the chemical time - scales may not fully decouple from the transport time - scales . in the n@xmath0h@xmath1 case considered here we have argued that the time - scales are likely to be decoupled ; this may not be the case for all snow - line tracer species , thus dynamical modelling which includes turbulent diffusion is needed to determine the importance of this effect . the model presented here should provide a stringent upper limit of the schmidt number , and good measurement if it is small ( @xmath97 ) ; however , if it is large , and the sharpness of the passive tracer has width of @xmath96 , then the 1d model would only provide an order of magnitude estimate and a better model is need to constrain the schmidt number . finally , if the snow - line resides in a dead - zone , where there is limited or no turbulence then it is unlikely this method can be used to cleanly probe the schmidt number ; but , dead - zones are not expected at the large radius of the co / n@xmath0h@xmath1 system discussed here . in this letter we have shown that the recent observations of snow - lines through tracer species ( e.g. n@xmath0h@xmath1 or h@xmath0co in the case of the co snow - line @xcite ) could allow direct observational measurements of the schmidt number in astrophysical accretion discs . in the case that the chemical time - scale is suitably de - coupled from the transport time - scales then the concentration gradient of the tracer outside the snow - line directly depends on the schmidt number in a power - law fashion ( @xmath98 ) independent of the choice of the turbulent @xmath92 parameter . we argue that the effect of turbulent diffusion on surface brightness distribution of such a snow line tracer is detectable with _ alma _ observations of discs in nearby star - forming regions which can possess high enough angular resolution to constrain the current theoretically / numerically estimated values of the schmidt number ( @xmath99 ) . observations of different snow - lines ( e.g. h@xmath0o , co@xmath0 & co ) at different radii in the discs would allow the radial dependence of the schmidt number to be probed . coupling these snow - line observations with observational estimates of the gas - surface density and turbulent line - widths would allow direct estimates of the strength and nature of the turbulence in astrophysical accretion discs .
Abstract:

Answer: sharp chemical discontinuities can occur in protoplanetary discs , particularly at ` snow - lines ' where a gas - phase species freezes out to form ice grains . such sharp discontinuities will diffuse out due to the turbulence suspected to drive angular momentum transport in accretion discs . we demonstrate that the concentration gradient - in the vicinity of the snow - line - of a species present outside a snow - line but destroyed inside is strongly sensitive to the level of turbulent diffusion ( provided the chemical and transport time - scales are decoupled ) and provides a direct measurement of the radial ` schmidt number ' ( the ratio of the angular momentum transport to radial turbulent diffusion ) . taking as an example the tracer species n@xmath0h@xmath1 , which is expected to be destroyed inside the co snow - line ( as recently observed in tw hya ) we show that _ alma _ observations possess significant angular resolution to constrain the schmidt number . since different turbulent driving mechanisms predict different schmidt numbers , a direct measurement of the schmidt number in accretion discs would allow inferences about the nature of the turbulence to be made .