Instruction: Write an abstract for the following acadamic paper.

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

Paper: the origin of crs remains an enduring question in astrophysics . the arrival directions of these particles are highly isotropic . this is expected from the propagation of charged particles in the interstellar medium where the directions of the particle momenta are randomized over time by the effective scattering in the encountered magnetic fields . due to the scrambling action of these fields , small anisotropy contrats are expected to be imprinted upon the distribution of arrival directions of crs as observed on earth . only for the population of ultra - high energy crs , magnetic deflections could be small enough to allow for mirroring to some extent the distribution of sources and the observed patterns , but the small flux combined to the potential absence of particles with low electric charge at these energies still prevents such a charged - particle astronomy with current data .    during the past decade , multiple observatories located in both hemispheres have reported significant observations of large - scale and small - scale anisotropies in the tev - pev energy band . these results have challenged the long - standing description of cr propagation in terms of a typical spatial diffusion process from stationary sources located preferentially in the disk of the galaxy , leading to a dipole moment only in the direction of the cr gradient and with an amplitude steadily increasing with the energy . the current picture , much more complex and elaborated , is discussed in section  [ sec : tev - pev ] .    in the pev - eev energy range , the possible increase of anisotropy amplitudes does not allow , yet , for compensating the decrease of the flux and of the collected statistics . there are essentially upper limits on the amplitude and phase measurements for the first harmonic in right ascension , which are presented in section  [ sec : pev - eev ] .    at eev energies and above , the study of ultra - high energy crs has experienced a jump in statistics as well as improved instrumentation in the past decade . this has allowed a better sensitivity in searching for anisotropies . currently , the stronger constraints come from measuring the first harmonic modulation in right ascension . at the highest energies , searches for clustering of arrival directions as well as searches for correlations between cr arrival directions and the positions of nearby extragalactic objects are the best suited tools in attempting to decipher the origin of these particles . these searches are summarized in section  [ sec : eev - zev ] .    prior to reviewing these results , and since most of the analysis techniques aim at extracting the moments of the observed angular distributions to characterize the anisotropies , general reminders on the formalism of the moment reconstruction are first provided in section  [ sec : analysis ] . at energies @xmath0tev , measurements of crs are indirect due to too low primary fluxes and are usually performed through extensive air shower arrays . these arrays operate almost uniformly with respect to sidereal time thanks to the rotation of the earth : the zenith - angle - dependent shower detection is then a function of declination but not a function of right ascension . thus , the most commonly used technique is the analysis in right ascension only , through harmonic analysis of the counting rate within the declination band defined by the detector field of view . considered as a function of the right ascension only , the flux of crs can be decomposed in terms of a harmonic expansion : @xmath1 the customary recipe to extract each harmonic coefficient makes use of the orthogonality of the trigonometric functions  @xcite . modelling any observed arrival direction distribution , @xmath2 , as a sum of @xmath3 dirac functions over the circle , @xmath4 , the coefficients can be estimated from the discrete sums : @xmath5 here , the re - calibrated harmonic coefficients @xmath6 and @xmath7 are directly considered , as it is traditionally the case in measuring _ relative _ anisotropies . in addition , over - lined symbols are used to indicate the _ estimator _ of any quantity . the statistical properties of the estimators @xmath8 can be derived from the poissonian nature of the sampling of @xmath3 points over the circle distributed according to the underlying angular distribution @xmath9 . in case of small anisotropies ( _ i.e. _ @xmath10 and @xmath11 ) , the harmonic coefficients are recovered with an uncertainty such that @xmath12 . for an isotropic realization , @xmath13 and @xmath14 are random variables whose joint p.d.f . , @xmath15 , can be factorized in the limit of large number of events in terms of two gaussian distributions whose variances are thus @xmath16 . for any @xmath17 , the joint p.d.f . of the estimated amplitude , @xmath18 , and phase , @xmath19 , is then obtained through the jacobian transformation  : @xmath20 from this expression , it is straightforward to recover the rayleigh distribution for the p.d.f . of the amplitude , @xmath21 , and the uniform distribution between 0 and @xmath22 for the p.d.f . of the phase , @xmath23 overall , this formalism provides the amplitude of the different harmonics , the corresponding phase ( right ascension of the maximum intensity ) , and the probability of detecting a signal due to fluctuations of an isotropic distribution with an amplitude equal or larger than the observed one as @xmath24 . note that the formalism aforementioned can be applied off the shelf only in the case of a purely uniform directional exposure , condition which is generally not fulfilled . at the _ sidereal _ time scale , the directional exposure of most observatories operating with high duty cycle ( _ e.g. _ surface detector arrays ) is however only moderately non - uniform . different approaches are then adopted in the literature , and the non - uniform directional exposures presented in next sections are carefully accounted for .      in general , and in contrast to the simplified approach presented in the last subsection , the flux of crs @xmath25 can depend on both the right ascension and the declination and thus be decomposed in terms of a multipolar expansion onto the spherical harmonics @xmath26 : @xmath27 non - zero amplitudes in the @xmath28 modes arise from variations of the flux on an angular scale @xmath29 radians . the partial - sky coverage of ground - based observatories prevents the multipolar moments @xmath30 to be recovered in a direct way through the customary recipe making use of the orthogonality of the spherical harmonics basis  @xcite . indirect procedures have to be used , one of them consisting in considering first the pseudo-multipolar moments @xmath31 and then the system of linear relationships relating these pseudo moments to the real ones : @xmath32 assuming a bound @xmath33 beyond which @xmath34 , these relations can be inverted allowing the recovering of the moments @xmath30 . however , the obtained resolution on each moment does not behave as @xmath35 ( with @xmath36 a constant factor depending on @xmath37 ) as expected from naive statistical arguments , but increases exponentially with @xmath33  @xcite .    it is to be noted that , in many practical cases where the declination dependence of the directional exposure @xmath38 is too difficult to estimate down to the required accuracy , analyses are performed by reproducing the observed declination dependence of the event counting rate . in such cases , anisotropies are captured in right ascension only , be it in several declination bands . hence , in general , the small values of the energy - dependent @xmath30 coefficients combined to the available statistics in the different energy ranges does not allow for an estimation of the individual coefficients with a relevant resolution as soon as @xmath39 . however , based on analysis techniques previously developed in the cmb community , it is possible , under some restrictions detailed and discussed in  @xcite , to reconstruct the angular power spectrum coefficients @xmath40 within a statistical resolution independent of the bound @xmath33 . the starting point is to consider any observed distribution of arrival directions as a particular realization of an underlying gaussian process . the simplest non - trivial situation to describe the underlying process is then to consider that the anisotropies cancel in ensemble average and produce a second order moment that does not depend on the position on the sphere but only on the angular separation between @xmath41 and @xmath42 . in this case , the underlying @xmath30 coefficients vanish in average and are not correlated to each other ( _ i.e. _ diagonal covariance : @xmath43 ) so that the @xmath44 coefficients can be viewed as a measure of the variance of the @xmath30 coefficients . in this situation , it can then be shown that the pseudo - power spectrum @xmath45 ( which is directly measurable ) is related to the real power spectrum through @xmath46 where the operator @xmath47 describing the cross - talk induced by the non - uniform exposure between genuine modes is entirely determined by the knowledge of the exposure function . various techniques are then available to estimate the @xmath44 coefficients from the measured @xmath48 ones . caution should be kept in mind when interpreting angular power spectra measured with partial - sky coverages : invoking a random nature for the arrival direction distributions of crs allows for a stochastic modelling of non - well known or unknown source positions and propagation mediums ; in opposition to cmb studies where the power spectrum results from primordial and fundamental fluctuations . as a result , the angular distributions in the uncovered regions of the sky could show patterns that would lead to different power spectra . in the tev - pev energy range , complex patterns have been revealed in the arrival directions of crs thanks to the large statistics collected in the last decade by several experiments , and anisotropy contrats at the @xmath50 level are now established at large scales . consistent measurements from experiments located in both hemispheres were reported : tibet as@xmath51  @xcite , super - kamiokande  @xcite , milagro  @xcite , eas - top  @xcite , minos  @xcite , argo - ybj  @xcite , and hawc  @xcite in the northern hemisphere , and icecube  @xcite and icetop  @xcite in the southern hemisphere . a collection of amplitude and phase measurements from  @xcite is shown in figure  [ fig : ampphase - tev ] including a series of former experiment results : the amplitude is observed to increase with energy up to @xmath52tev before flattening , and the phase is observed to be smoothly evolving before undergoing a sudden flip at @xmath53pev . note that small corrections to the amplitudes related to the different latitudes and zenithal - dependent efficiencies of the different experiments are not applied here . the first harmonic parameters as derived in figure  [ fig : ampphase - tev ] are generally considered with special interest due to their relationships with the dipole moment . in the context of spatial diffusion , the dipole moment is naively expected to provide a way to probe the particle density gradient shaped by the diffusion in interstellar magnetic fields on scales of the scattering diffusion length . in this picture , for stationary sources smoothly distributed in the galaxy , the dipole vector should align roughly with the direction of the galactic center with an amplitude increasing with energy in the same way as the diffusion coefficient , typically @xmath54 . however , this simple picture is not confirmed by the measurements , showing that the dipole amplitude is not described by a single power law and that the dipole phase does not align with the galactic center and undergoes a rapid flip at an energy of 0.1 - 0.3 pev . recent studies have put this picture into question by considering spatial and temporal stochastic distributions of sources : even for isotropic diffusion governed by a single diffusion coefficient growing as @xmath54 with energy , fluctuations of the dipole parameters are then naturally induced by local and young sources even if their contribution to the total flux is only subdominant  @xcite . on top of the background diffuse emission , the contribution of such sources within some energy intervals can lead , depending on the source population and the diffusion region , to phase rotations and amplitude modulations as a result of the summation of competing dipole vectors associated to each source .        within the framework of isotropic diffusion , a phase flip with vanishing amplitude requires , nevertheless , fine tuning of the contribution of local sources on top of the diffuse emission . additional effects may enter into play to shape the dipole parameters actually observed , such as the anisotropic diffusion induced by the local ordered magnetic field whose strength is such that circular motions predominate over random scattering . in this case , the observed dipole could result from the projection of the density gradient of crs onto the direction of the ordered magnetic field  @xcite . assuming that the measured first harmonic amplitude @xmath55 and phase in right ascension @xmath56 are predominantly determined by the dipole vector @xmath57 , the equatorial dipole components can be derived as @xmath58 , with @xmath59 a parameter depending on each experiment and @xmath60 and @xmath61 the vector components in the equatorial plane associated to the unit vectors pointing towards local sidereal times 0  h and 6  h. assuming on the other hand that the local standard of rest ( lsr ) corresponds to the local plasma frame in which crs are isotropic , the contribution to the observed dipole components @xmath62 arising from the relative motion of the solar system through the lsr , known as the compton - getting effect  @xcite , can be subtracted . doing so , the summary plot of the reconstructed tev - pev dipole components in the equatorial plane obtained recently in  @xcite is shown in figure  [ fig : ahlers ] . the numbers attached to the data indicate the median energy of the bins as @xmath63 . the dashed line and gray - shaded area indicate the magnetic field direction and its uncertainty ( projected onto the equatorial plane ) inferred from ibex observations  @xcite . a close alignment of the inferred dipole components with the local magnetic field direction is observed . overall , within this framework of anisotropic diffusion , although the projection of the dipole onto the magnetic field axis does not allow one to reconstruct the cr gradient , the magnetic hemisphere of the cr gradient can be determined by the dipole phase and aligns with galactic longitudes between 120@xmath64 and 300@xmath64 below 0.1 - 0.3 pev . in this regard , and among other requirements , the vela snr is shown to be a plausible source candidate in  @xcite .          beyond the large - scale anisotropies captured by the first harmonic in right ascension , smaller but significant anisotropy contrasts at the @xmath65 level at intermediate and smaller scales have also been reported in the last decade  @xcite . the angular power spectrum as derived from icecube data with a median energy of 20  tev , shown as the black points in figure  [ fig : powerspectrum ] , provides a summary plot to visualize the relatively large amount of power at low multipole moments @xmath66 as well as the small but still significant power at moments up to @xmath67 . the expected power spectrum fluctuations in 90% of cases of an isotropic angular distribution is indicated by the light gray band . these results have challenged the long - standing picture of cr propagation in terms of a spatial diffusion that would lead predominantly to a dipole moment only . several interpretations have been put forward . the most conservative one is based on the fact that the observed complex patterns , whose positions appear randomly distributed , are mainly the consequence of the particular structure of the turbulent magnetic field within the last sphere of diffusion encountered by crs  @xcite . with a density gradient of crs at the entrance of this last sphere of diffusion , the turbulent magnetic field is frozen on the time scale of cr propagation and thus plays the same role as a structured one . the field is then expected to connect regions of higher density outside from the last sphere of diffusion to regions of lower density as seen from earth , and _ vice versa_. this mechanism was shown to produce intermediate- and small - scale anisotropies through monte - carlo studies  @xcite , and more recently through an analytical approach  @xcite . in the former approach , numerical integration of test particles in low @xmath68-compressible magnetohydrodynamic turbulence with the external mean magnetic field as the only free parameter is shown to reproduce reasonably well the power spectrum , as seen in figure  [ fig : powerspectrum ] through the different solid color lines standing for different primaries  @xcite . in the latter approach , an initially dipolar distribution is shown to be distorded into a much complex pattern described by higher - order multipoles when introducing second order corrections to the diffusion approximation . these corrections aim at describing the propagation regime between the diffusive particle transport and the deterministic flow of particles . the asymptotic strength of the high - order multipoles at times longer than the diffusion relaxation time is observed to reproduce the data in a reasonable way , as shown by the dotted line in figure  [ fig : powerspectrum ]  @xcite . overall , thanks to the contemporary measurements of cr arrival directions , the information encompassed in the angular distributions appear today as a tool allowing a possible probe of the local magnetic field environment .    besides , it is also to be noted that the power spectrum at small scales measured at tev energies may also provide relevant signatures to probe and study the electric field induced by the motion of the heliosphere relative to the plasma rest frame where the electromagnetic field can be considered as purely magnetic due to the high conductivity of the medium  @xcite . the anisotropy would then be the consequence of the changes of energy associated to the acceleration / deceleration of the particles , which , given the size and speed of the heliosphere , are expected to produce changes of intensity within small angular scales at the required level for tev crs . at higher energy , in the pev - eev energy range , the expected increase of anisotropy contrasts does not compensate , yet , the decrease in the collected statistics with increasing energy . the most constraining data , collected in the top panel of figure  [ fig : peveev ] in terms of the equatorial component of the dipole , are provided by the icetop  @xcite , kascade - grande  @xcite and auger  @xcite experiments . except for the icetop amplitude , only upper limits are currently provided in this energy range . meanwhile , it is interesting to note that an apparent constancy of phase , even though the significances of the amplitudes are relatively small , has been pointed out previously in surveys of measurements with ion chambers and counter telescope made in the range @xmath69  @xcite . a clear tendency for maxima to occur around 20 hours l.s.t . was observed . as already pointed out by linsley long time ago , this is potentially indicative of a real underlying anisotropy , because a consistency of the phase measurements in ordered energy intervals is indeed expected to be revealed with a smaller number of events than needed to detect the amplitude with high statistical significance  @xcite . interestingly , the phases reported by contemporary experiments collected in the bottom panel of figure  [ fig : peveev ] show a consistent tendency to align in the general right ascension of the galactic center .    above @xmath70eev , a change of phase is observed towards , roughly , the opposite of the one at energies below @xmath70eev . the percent limits to the amplitude of the anisotropy exclude the presence of a large fraction of galactic protons at eev energies  @xcite . accounting for the inference from @xmath71 data from both the pierre auger observatory and the telescope array that protons are in fact abundant at those energies , this might indicate that this component is extragalactic , gradually taking over a galactic one . the low level of anisotropy would then be the sum of two vectors with opposite directions , naturally reducing the amplitudes . this scenario is to be explored with additional data . increased statistics is thus necessary to probe the anisotropy contrast levels that may exist in this energy range and contain valuable information about the old - age question on the way the transition between galactic and extragalactic crs occurs . also , a current limitation of the measurements is that neither spectra nor anisotropies can yet be studied as a function of the mass of the particles with adequate statistical precision , measurements that would allow a distinction between galactic and extragalactic angular distributions . above @xmath52eev , the flux of ultra - high energy crs ( uhecrs ) is expected to be of extragalactic origin . although the actual sources of uhecrs are still to be identified , their distribution in the sky is expected to follow , to some extent , the large - scale structure of the matter in the universe . it is thus interesting to highlight that in addition to the upper limit on the equatorial dipole component above 8  eev , the measured amplitude is also shown in figure  [ fig : peveev ] as the corresponding @xmath72value is as low as @xmath73  @xcite . assuming that the only significant contribution to the anisotropy is from a dipolar pattern , the amplitude of this signal converts into a @xmath74 dipole amplitude  @xcite . this hint may constitute in the near future the first detectable signature of extragalactic crs observed on earth .    to characterize further the angular distribution above 10  eev , the dipole moment on the sphere is of special interest . an unambiguous measurement of this moment as well as of the full set of spherical harmonic coefficients requires full - sky coverage . currently , this can be achieved by combining data from observatories located in both the northern and southern hemispheres . to this end , a joint analysis using data recorded at the pierre auger observatory and the telescope array above @xmath75  eev has been performed in  @xcite . thanks to the full - sky coverage , the measurement of the dipole moment reported in these studies does not rely on any assumption on the underlying flux of crs . note that , in contrast to the results discussed in section  [ sec : tev - pev ] , the directional exposure derived in these studies is a truly function of both right ascension and declination , so that the reconstructed anisotropy patterns are not absorbed along the declination . eev smoothed out at a 60@xmath64 angular scale , in km@xmath76  yr@xmath77  sr@xmath77 units  @xcite.,width=283 ]    the resulting entire mapping of the celestial sphere has revealed a dipole moment with an amplitude @xmath78 , captured with a chance probability of @xmath79 . no other deviation from isotropy can be observed at smaller angular scales . the recovered moment can be visualized in figure  [ fig : augerta ] , where the average flux smoothed out at an angular scale @xmath80 per solid angle unit is displayed using the mollweide projection , in km@xmath76yr@xmath77sr@xmath77 units . this map is drawn in equatorial coordinates . the direction of the reconstructed dipole is shown as the white star . large - scale anisotropies of crs with energies in excess of @xmath81eev are closely connected to the sources and the propagation mode of extragalactic uhecrs , see e.g.  @xcite . due to scattering in the extragalactic magnetic fields , large deflections are expected even at such high energies for field amplitudes ranging in few nanogauss and extended over coherence lengths of the order of one megaparsec , or even for lower amplitudes if the electric charge of uhecrs is large . for sources distributed in a similar way to the matter in the universe , the angular distribution of uhecrs is then expected to be influenced by the contribution of nearby sources , so that the milky way should be embedded into a density gradient of crs that should lead to at least a dipole moment . the contribution of nearby sources is even expected to become dominant as the energy of crs increases due to the reduction of the horizon of uhecrs induced by energy losses more important at higher energies . once folded through the galactic magnetic field , the dipole pattern expected from this mechanism is transformed into a more complex structure presumably described by a lower dipole amplitude and higher - order multipoles . however , in these scenarios , the dipole moment could remain the only one at reach within the sensitivity of the current generation of experiments . on the other hand , the detection of significant multipole moments beyond the dipole one could be suggestive of non - diffusive propagation of uhecrs from sources distributed in a non - isotropic way .    while it has a small influence      -radius windows for the auger events with @xmath82  eev . also indicated are the super - galactic plane ( dashed line ) and centaurus  a ( white star )  @xcite.,width=283 ]    the energy losses of uhecrs limit the horizon of the highest - energy particles . for small - enough magnetic deflections , the distribution of the arrival directions of uhecrs above @xmath83eev could mirror the inhomogeneous distribution of the nearby extragalactic matter . the search for anisotropy at small and intermediate angular scale at the highest energies is thus potentially the most powerful way to infer the sources of uhecrs .    at the pierre auger observatory , data have been subjected to comprehensive anisotropy searches for different energy thresholds between @xmath84 and @xmath85eev , and within different angular windows , between @xmath86 and @xmath87  @xcite . searches for significant excesses anywhere in the sky have been performed , as well as searches for correlations with known astrophysical structures and with objects that are considered plausible candidates for uhecr sources . out of all the searches performed , none of the analyses provides any statistically significant evidence of anisotropy . the two largest departures from isotropy , both with post - trial probability @xmath88 , are found for energies in excess of @xmath89eev when looking , on the one hand , within @xmath90 of the direction of centaurus a - the closest radio - loud active galactic nucleus ( agn ) , and when looking , on the other hand , within @xmath91 of the objects collected in the flux - limited catalog of agns observed in x - rays ( swift - bat-70  @xcite ) closer than 130  mpc and brighter than @xmath92erg / s . note that the most significant excess ( 4.3@xmath93 pre - trial ) observed in the blind search corresponds to a region close to the super - galactic plane and to the direction of cen a , at a similar energy threshold , 54  eev , and similar angular scale , @xmath94 . a sky map in galactic coordinates of the li - ma significances of overdensities in 12@xmath64-radius windows for the events with @xmath82  eev is shown in figure  [ fig : cena ] . -radius windows for the telescope array events with @xmath95  eev  @xcite . ]    at the telescope array , similar searches for anisotropies have also been performed above 57  eev . using data collected over a 5-year period , a cluster of events has been found by oversampling the sky map using 20@xmath64-radius circles  @xcite . the cluster of events , centered at about 19@xmath64 off of the super - galactic plane but on top of no known specific object , has a li - ma statistical significance of 5.1@xmath93 ( pre - trial ) . once penalized for the angular scan performed each 5@xmath64 , the significance of the excess is 3.4@xmath93 . a confirmation of this excess requires more statistics . overall , the absence of strong anisotropies at ultra - high energies is amenable to different interpretations . if the bulk of particles is made of light nuclei , that might indicate a large number of sources . if in turn the bulk of particles is made by large-@xmath96 nuclei , then the lack of anisotropy might be caused by large deflections . information on the mass of the primaries at these energies is thus also of relevance in the study of the distribution of the arrival directions .      in the context of testing the distribution of the arrival directions of the highest - energy crs , it is worth mentioning the result of a full - sky study , conducted in a collaboration among auger , telescope array and icecube  @xcite . it consists in the search for correlations between the directions of uhecrs observed at the auger observatory and at the telescope array , and those of very high - energy neutrino candidates detected by icecube . although no indications of correlations at discovery level are found , it is to be noted that the smallest post - trial @xmath72values ( corresponding to @xmath97 ) are obtained when considering the correlations between the directions of cascade events observed by icecube and those of the crs . the excess of correlations , found at angular scales of @xmath98 , arises mostly from pairs of events in the region of the sky where the telescope array has reported an excess of events and in regions close to the super - galactic plane in correspondence with the largest excess observed in auger data . further insight will potentially arise from increased statistics and eventually with the inclusion of cr composition information that may become available so as to better model possible effects of magnetic field deflections . this may help to understand if there is a contribution in the astrophysical neutrino signal observed by icecube correlated to the sources of the observed uhecrs . during the past decade , important observational results have been reported on the angular distributions of tev - pev crs . while only dipolar excesses were expected , the myriad of reported anisotropies has led to important progresses on the understanding of the propagation regime of low - energy galactic crs .    in contrast , the quest for finding uhecr sources is more difficult than expected a decade ago . future work will profit from the increased statistics and ability to perform anisotropy searches with distinction based on the mass of the primaries as anticipated with the upgraded instrumentation at the pierre auger observatory  @xcite . however , another jump in statistics appears necessary , keeping similar observable resolutions . m. amenomori et al . , apj 626 ( 2005 ) l29 g. guillian et al . , phys . rev . d75 ( 2007 ) 062003 a.a . abdo et al . , phys . ( 2008 ) 221101 m. aglietta et al . , apj 692 ( 2009 ) l130 j. de jong , proc . 32nd icrc , beijing , china , 2011 b. bartoli et al . d88 ( 2013 ) 082001 s. benzvi et al . 33rd icrc , rio de janeiro , brazil , 2013 r. u. abbasi et al . , apj 718 ( 2010 ) l194 m. g. aartsen et al . , apj 765 ( 2012 ) 55 g. di sciascio and r. iuppa ,  homage to the discovery of cosmic rays  , chapter 9 , nova science publishers , new york , 2013 a. d. erlykin and a. wolfendale , astropart.phys . 25 ( 2006 ) 183 p. blasi and e. amato , jcap 1201 ( 2012 ) 010 l. g. sveshnikova , o. n. strelnikova and v. s. ptuskin , astropart . phys . 50 - 52 ( 2013 ) 33 p. mertsch and s. funk , phys . ( 2015 ) 021101 m. ahlers , phys . lett . 117 ( 2016 ) 151103 a. h. compton and i. a. getting , phys . 47 ( 1935 ) 817 h. o. funsten et al . , apj 776 ( 2013 ) 30 r. u. abbasi et al . , apj 740 ( 2011 ) 16 b. bartoli et al . d88 ( 2013 ) 082001 a. u. abeysekara et al . , apj 796 ( 2014 ) 108 m. g. aartsen et al . , apj 826 ( 2016 ) 200 g. giacinti and g. sigl , phys . lett . 109 ( 2012 ) 071101 m. ahlers , phys . 112 ( 2014 ) 021101 m. ahlers and p. mertsch , apj 815 ( 2015 ) l2 v. lpez - barquero et al . , apj 830 ( 2016 ) 19 l. drury , proc . 33rd icrc , rio de janeiro , brazil , 2013 , arxiv:1305.6752    a. chiavassa for the kascade - grande collaboration , proc . 34th icrc , the hague , the nederlands , 2015 i. al samarai for the pierre auger collaboration , proc . 34th icrc , the hague , the nederlands , 2015 j. delvaille , f. kendziorski , and k. greisen , j. phys . japan 17 ( suppl . a - iii ) ( 1962 ) 76 d. edge et al . , j. phys . g 4 ( 1978 ) 133 the pierre auger collaboration , apjs 203 ( 2012 ) 34 r. u. abbasi et al . , astropart . 86 ( 2017 ) 21 the pierre auger collaboration , apj 802 ( 2015 ) 111    the pierre auger and telescope array collaborations , apj 794 ( 2014 ) 172 o. deligny for the pierre auger and telescope array collaborations , proc . 34th icrc , the hague , the nederlands , 2015 d. harari , s. mollerach and e. roulet , phys . d89 ( 2014 ) 123001 p. tinyakov and f. urban , j. exp . 120 3 ( 2015 ) 533 the pierre auger collaboration , apj 804 ( 2015 ) 15 w. h. baumgartner et al . , apj 207 ( 2013 ) 19 r. u. abbasi et al . , apj 790 ( 2014 ) l21 the pierre auger , telescope array and icecube collaborations , jcap 01 ( 2016 ) 037 r. engel for the pierre auger collaboration , proc . 34th icrc , the hague , the nederlands , 2015
Abstract:

Answer: important observational results have been recently reported on the angular distributions of cosmic rays ( crs ) at all energies , calling into question the perception of crs a decade ago . these results together with their in - progress interpretations are summarized in this short overview paper , following the contribution presented at the european cosmic ray symposium .


Instruction: Write an abstract for the following acadamic paper.

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

Paper: the interpretation of astrophysical spectra rests heavily upon the availability of relevant atomic and molecular data . of critical importance are ( 1 ) the energies of the various accessible quantum states , ( 2 ) the rates of spontaneous radiative decay between states , and ( 3 ) the state - to - state rate coefficients for excitation and deexcitation in inelastic collisions . the last of these have proven quite elusive ; in many circumstances , our ability to interpret astrophysical spectra is severely limited by the availability of collisional rate coefficients . as a result of new spectroscopic facilities such as the _ spitzer space telescope _ , the _ herschel space observatory _ , and the planned atacama large millimetric array ( alma ) , the need for fundamental molecular data has become all the more acute . in modeling the emission from interstellar molecular clouds , in particular , there is a critical need for rate coefficients for the rotational excitation of molecules in collisions with h@xmath3 . the theoretical calculation of such data is extremely time consuming , involving the computation of a potential energy surface for the system in question , followed by extensive scattering calculations . even for those molecules for which such calculations have been carried out , the available data are often limited to collisionally - excited collisions involving a relatively small set of low - lying states . submillimeter and infrared spectra of molecular clouds often reveal the presence of higher - lying transitions for which collisional rate coefficients have not been computed ; here , some kind of extrapolation method is needed . for the case of linear molecules , the infinite order sudden ( ios ) and related approximations provide an obvious method of extrapolation ; for asymmetric molecules , there is no simple analogous method .    given a limited set of transitions of an asymmetric molecule for which collisional rate coefficients have been computed , previous efforts at extrapolation ( e.g.  neufeld & melnick 1987 , faure & josselin 2008 ) have typically adopted an expression with one or more adjustable parameters for how the collisional rate coefficient depends upon the relevant parameters for each transition ( energy difference , change in rotational quantum number , etc . ) the adjustable parameters are then chosen so as to optimize the fit to those transitions for which collisional rate coefficients are known , and then the expression is used to estimate the values for higher - lying transitions where calculated results are unavailable . for example , the recent study by faure & josselin ( 2008 ) considered the rate coefficients for the excitation of h@xmath3o by h@xmath3 , which will be of crucial importance to the interpretation of many _ spitzer _ and _ herschel _ spectra . here , faure & josselin extrapolated a set of rate coefficients obtained from quasi - classical trajectory calculations for transitions among the lowest 45 rotational states of water . each such state is characterized by an energy , @xmath4 , and rotational quantum numbers @xmath5 , @xmath6 and @xmath7 that describe the total angular momentum and its projection onto the principal axes of the molecule . for transitions with @xmath8 , @xmath9 and @xmath10  termed  high - propensity " transitions  faure & josselin assumed that rate coefficients depend solely upon @xmath11 , @xmath12 and @xmath13 . for all other transitions (  low propensity transitions " ) , they assumed the rate coefficients to be _ independent _ of @xmath11 , @xmath12 and @xmath13 and proportional to @xmath14 , where @xmath15 is adjusted to optimize the fit to the available data . as an alternative to the methods adopted previously , i have explored the use of an artificial neural network ( ann ) to accomplish the extrapolation . in contrast to previous techniques , the use of an ann does not make any _ a priori _ assumption about how the transitions might be categorized or how the rate coefficients might depend upon the parameters for each transition . the calculational method is discussed in  2 ( with the more technical details given in the appendix ) . the performance of the ann in reproducing a known set of collisional rate coefficients in presented in  3 . a discussion follows in  4 . since the mid-1980 s , artificial neural networks have been widely used to solve a wide variety of problems , particularly in pattern recognition . previous astrophysical applications have centered primarily on the automatic classification of astronomical objects based on morphology , photometry or spectroscopy ( e.g.  ball et al . 2004 , collister & lahav 2004 ) , although ann have also been used previously as approximation tools ( e.g. asensio ramos & socas - navarro 2005 ) . anns consist of a set of interconnected artificial neurons . such neurons are of three types : input neurons , output neurons , and hidden neurons . given a set of input values that are fed into the input neurons , a unique set of output values emerges from the output neurons . in the context of the application considered here , the input neurons are fed the relevant parameters for each transition ( e.g. the energies and relevant quantum numbers for the upper and lower states ) , and the collisional rate coefficient emerges from a single output neuron . the hidden neurons mediate the dependence of the ann s output upon its input . each individual neuron accepts a particular input value , and outputs a value that is  in general  a non - linear function of that input . this function is called the activation function , and is typically sigmoidal in shape . ( the inverse hyperbolic tangent is an often - used activation function in ann . ) with the exception of the input neurons , for which the inputs are a set of input parameters , the input for any neuron is the weighted sum of the outputs of a set of other neurons to which it is connected . during the  learning process `` , the ann is shown a ' ' training set " of input parameters , together with the appropriate outputs ; the relevant weights of the connections are varied until the outputs ( from the output neurons ) match the appropriate values to within the desired accuracy . the arrangement of the connections between the neurons is called the  architecture " of the network . in  feed - forward " networks such as the one used in this study , the structure is non - recursive : this means that , when considered in the appropriate order , the input value for any neuron reflects only output values for neurons that have already been considered . in some ann , the architecture is fixed , and the learning process simply involves varying the weights of the connections . in other schemes , such as the cascade - correlation network introduced by fahlman and lebiere ( 1990 ) and used in the present study , the architecture evolves during the learning process , with hidden neurons trained and added one at a time to improve the performance of the ann . figure 1 shows a simple example of a cascade - correlation network , with 3 input neurons ( blue circles numbered 1  3 ) , 2 hidden neurons ( numbered 4 and 5 ) , and one output neuron ( number 6 ) . the @xmath16th neuron has an input value @xmath17 ( shown in red ) and an output value @xmath18 = @xmath19 ( green ) , where @xmath20 is an activation function ( which may , in principle , be different for each neuron  see appendix ) . information propagates from left to right , starting with the input values @xmath21 , @xmath22 and @xmath23 , and ending with the output value @xmath24 . the @xmath17 and @xmath18 are related by a matrix of weights , * w * , with @xmath25 . in the example shown here , * w * has 12 non - zero elements . at the end of training , the weights have assumed values that optimize the performance on the training set : in other words , @xmath24 is close to the desired output value for each given set of inputs , @xmath26 $ ] .    in this study , i used the fast artificial neural network ( fann ) library , an open source software package developed by steffen nissen and evan nemerson . full details of my use of this software are presented in the appendix ; my intent there is to provide enough specificity to allow the reader to reproduce the calculational method exactly . a copy of the relevant computer code , written in c , will be provided upon request to the author .    in evaluating the performance of an ann , it is essential to train the network on a  training set " of data and then test its performance on a  test set " _ that it has never seen . given sufficient complexity in the network ( i.e. enough hidden neurons ) , it is always possible to fit the training set to any desired degree of accuracy . the key question is how well the ann performs on new , unseen data . in this study , i considered the collisional rate coefficients computed recently by dubernet et al . ( 2009 ) for the excitation of ortho - h@xmath3o in collisions with para - h@xmath3 . the close - coupling calculations that led to these results represent the state - of - the - art in scattering calculations and are extremely expensive computationally . this dataset provides the rate coefficients for all transitions among the 45 lowest rotational states of ortho - h@xmath3o , computed for several temperatures and several values of the initial and final rotational states of the colliding para - h@xmath3 molecule . for simplicity , i considered first the values obtained at a single temperature , 800  k , and with the initial and final h@xmath3 states being @xmath27 . this dataset then contains @xmath28 independent values , where @xmath29 is the total number of water rotational states . the factor of @xmath30 arises because the upwards and downwards rates for any transition are related by considerations of detailed balance and are therefore not independent of each other . as the training set , i adopted a subset of these 990 data points ; because the ultimate goal of this method is to estimate collisional rate coefficients for transitions that are _ higher _ in energy than those calculated previously , i chose a training set consisting of those collisional - induced transitions between the _ lowest _ @xmath31 rotational states , where @xmath32 . the training set then consists of @xmath33 data points , with the test set comprising the remaining @xmath34 data points . one key methodological choice concerns the input parameters that will be provided to the ann ; this determines the number of input neurons . each rotational state of ortho - h@xmath3o is characterized by an energy @xmath4 , and three rotational quantum numbers : @xmath5 , which describes the total rotational angular momentum , @xmath6 , which describes its projection on the principal axis for which the moment of inertia is smallest , and @xmath7 , which describes its projection on the principal axis for which the moment of inertia is largest . however , once @xmath5 and @xmath6 are specified , @xmath7 is uniquely determined . thus it is convenient to define a parameter @xmath35 , which allows each rotational state of ortho - h@xmath3o to be specified unambiguously by just two quantum numbers , @xmath5 and @xmath36 . the integer @xmath36 takes all odd values between @xmath37 and @xmath38 . in the present study , i experimented with three choices for the input parameter set . the three cases considered were : ( 1 ) the ann is fed just the energy difference between the initial and final state , @xmath39 , and therefore has one input neuron ; ( 2 ) the ann is fed @xmath39 , @xmath13 , and @xmath40 , and therefore has 3 input neurons ; ( 3 ) the ann is fed @xmath4 , @xmath5 , and @xmath36 for both the initial and final states , requiring a total of 6 input neurons . in  3 below , the performance of the ann is described for each of these three cases . in all cases , the ann has single output neuron which provides the rate coefficient for the downward collisionally - induced transition . because the rate coefficients cover several orders of magnitude , i trained the network to output the logarithm of the rate coefficient . two other choices , more technical in nature , determine the performance of the ann . the final number of hidden neurons , @xmath41 , is discussed in the appendix . based upon the experimentation described there , i adopted @xmath42 as the optimum value . with too few hidden neurons , the ann fails to fit even the training data set well ; with too many hidden neurons , the network fits the training set perfectly but fails to generalize well to the test set . because the initial weights on the connections are randomly assigned , the results vary slightly every time the ann is trained . accordingly , i took the average of the outputs of several ann in obtaining results for the test set . as discussed in appendix a , the performance of the ann fails to improve once the number of realizations , @xmath43 , exceeds @xmath44 ; this behavior indicates that the remaining errors in the ann predictions are not a consequence of the particular choice of initial weights , but are instead a fundamental limitation of the method . the results presented in  3 below all apply to an ann with @xmath42 and @xmath45 the top right , bottom left , and bottom right panels in figure 2 present the performance of an ann with @xmath46 for the cases with 1 , 3 , and 6 input neurons discussed in  2 above . in each panel , the horizontal axis represents the logarithm of the rate coefficient ( in units of @xmath47 ) obtained by dubernet et al . ( 2009 ) , and the vertical axis represents the corresponding output of the ann . red points apply to the 435 data points in the training set , while blue points apply to the 555 data points in the test set ( i.e.  to those transitions that were not seen by the ann during training ) . as expected from the discussion of high and low propensity transitions presented by faure & josselin ( 2008 ) , the ann performs much better for the cases with 3 and 6 input neurons ; only in these cases , is it informed about the rotational quantum numbers along with the transition energy . indeed , the performance in the 1-input - neuron is strikingly similar to the results obtained if the data are all fit with an expression of the form @xmath48 . adopting best - fit values of @xmath49 and @xmath15 determined from a least squares fit to the data yields an expression which is compared to the actual data in the top left panel of figure 2 . more detailed statistical information is provided in figures 3 and 4 . in figure 3 , the frequency distribution of the errors , @xmath50 , is shown for each case . black histograms apply to the entire test set , while blue , green and red histograms apply to transitions with low , medium and high actual rate coefficients ; for this purpose ,  low " applies to value smaller than @xmath51 ,  medium " to value in the range @xmath52 , and  high "  to a value greater than @xmath53 . clearly , the cases with 3 and 6 neurons yield histograms that are more strongly peaked than that obtained with 1 input neuron , indicating that the typical error in the ann output is considerably smaller . while the probability distributions for the 3- and 6-input - neuron cases are very similar when all the test data are considered together ( black histogram ) , the 6-input - neuron ann systematically underpredicts the results for transitions with high rate coefficients ( red histogram ; this behavior is also apparent from a careful inspection of figure 2 . ) paradoxically , the ann performs worse when it is provided with additional information ( i.e. the individual @xmath4 , @xmath5 , @xmath36 values for the upper and lower states rather than just the differences between those values , @xmath39 , @xmath13 , @xmath40 ) . in the 6-input - neuron case , the ann is evidently fitting a trend with energy , @xmath4 , that is present for the lower energy transitions ( i.e.  for the training set ) , but extrapolates poorly to the higher - energy transitions encountered in the test set . based upon figures 2 and 3 , i adopted the 3-input - neuron ann as the optimum architecture , along with the values @xmath42 and @xmath54 discussed in the appendix . for this architecture , the distribution of the output errors is presented in an alternative form in figure 4 . here i show the fraction of the data for which @xmath55 exceeds a given value , @xmath56 . the fraction of rate coefficients is shown on the vertical axis ( logarithmic scale ) , with the same color coding adopted in figure 3 , as a function of @xmath56 . for example , the fact that the red curve passes through the point ( 0.19 , 0.30 ) implies that one - half ( i.e.  10@xmath57 ) of those rate coefficients with high actual values ( i.e. greater than @xmath53 ) are in error by more than a factor of 1.5 ( i.e.  @xmath58 ) once again , results apply only to the test set . clearly , the ann output is less reliable for transitions with rate coefficients smaller than @xmath51 ( blue curve ) . even in this case , however , more than half of the rate coefficients are correct to within a factor 2.1    the results shown in figures 2  4 were all obtained with @xmath59 , for which the training set comprises the lowest 435 transitions and the test set contains the next 555 transitions . i have also investigated the performance of the ann as a function of the size of the training set , i.e.  as a function of @xmath31 . in figure 5 , i show the rms value of @xmath60 , as a function of n@xmath61 , for the entire set of transitions ( black points ) and for low , medium and high values of the actual rate ( blue , dark green and red as before ) . results are shown for @xmath62 , 15 , 20 , 25 , 30 , 35 , and 40 . these correspond to training sets containing 45 , 105 , 190 , 300 , 435 , 595 , and 780 data points , respectively , and test sets containing 945 , 785 , 800 , 690 , 555 , 395 , and 210 data points . the overall performance varies little once @xmath63 , i.e.  once the training set encompasses @xmath64 of the total available data .    one key application of the present study , of course , will be to train the ann on all 990 transitions for which dubernet et al . have computed collisional rate coefficients , and to then use the ann to estimate the rates for higher - lying transitions for which no calculations have been performed . the results of figure 5 provide a gauge of the likely success of that procedure , and of the number of additional rate coefficients that can be estimated reliably . if the accuracy of the extrapolation depends upon the _ relative _ size of the training set ( i.e.  the _ ratio _ of the number of transitions in the training set to the total number for which ann estimates are to be obtained ) , we may expect that an additional @xmath65 rate coefficients can be predicted with the use of an ann . on the other hand , it is possible that the _ absolute _ size of the training set is the relevant parameter ; in that case , there is no obvious limit to how far the extrapolation can be carried . a conservative approach would be to use an ann to extrapolate no more than an additional @xmath65 rate coefficients ; the resultant data set would include all transitions among the lowest @xmath66 rotational states of o - h@xmath3o . for higher - lying transitions , the method of faure & josselin ( 2008 ) promises to be more robust : in that method , the results for high - propensity transitions are obviously well - behaved , and those for low - propensity transitions tend gracefully to zero as @xmath39 tends to infinity . with the ann , by contrast , the behavior could be unstable if one extrapolates to a @xmath39 ( or @xmath13 or @xmath40 ) that is much larger than anything it saw in the training set . while the collisional excitation of water has received a great deal of theoretical attention recently , many other molecules will be observed by alma and herschel . figure 6 shows results obtained for three other molecules for which there is a large set of collisional rate coefficients in the basecol database : so@xmath3 ( from green 1995 ) , ortho - c - c@xmath67h@xmath3 ( chandra & kegel 2000 ) , and ortho - h@xmath3co ( green 1991 ) . like water , these are asymmetric top molecules with rotational states that can be characterized by the quantum numbers @xmath5 and @xmath36 , although in these cases , the available data apply to excitation in collisions with he . i have also obtained additional results for the collisional excitation of ortho - water with para - h@xmath3 , again using rate coefficients from dubernet et al . ( 2009 ) but now for cases in which the initial and/or final state of h@xmath3 is @xmath68 . in each case , the training set consisted of the lower - lying half of the entire available data set . the rms value of @xmath60 for the test set is printed on each panel ; as in the case of ortho - h@xmath3o colliding with h@xmath3 ( @xmath69 ) , typical rms values are roughly 0.35 \1 . artificial neural networks ( ann ) can be successfully trained to determine rate coefficients for the collisional excitation of astrophysical molecules . the type of ann favored here is a cascade correlation network which creates 16 hidden neurons during the course of training . there are three input neurons that inform the network of @xmath39 , @xmath13 and @xmath40 for a given transition , and one output neuron that provides the logarithm of the rate coefficient for collisional de - excitation . the performance of such an ann can be evaluated by training it on a subset of the set of transitions for which collisional rate coefficients are known , and then testing it on the remainder of that set . since the motivation for this study is to develop a method for estimating the rate coefficients for high - lying transitions for which values have not been computed , the training set involved the lowest - lying transitions among the available data . \5 . in the primary test case considered here , the excitation of ortho - h@xmath3 by h@xmath3 ( @xmath27 ) , the available data  consisting of 990 rate coefficients  can be modeled equally well provided that the network is trained on at least @xmath70 of the available data . the performance of the ann is summarized in table 1 , which lists the typical output errors in every case considered . here the rms and median values of @xmath71 value ) are listed for the entire test set and separately for low , medium and high values of the actual rate . overall , for ann with 16 hidden neurons and 3 input neurons that are trained on at least 20@xmath2 of the available data , the typical rms value of @xmath60 is @xmath72 ; for transitions with rate coefficients @xmath73 , the typical rms value is @xmath74 . the corresponding median values of @xmath75 are @xmath76 and @xmath77 . these median values imply that one - half of the ann outputs are typically accurate to within factors of 1.8 ( or @xmath0 in the case of rate coefficients @xmath73 . ) however , the results shown in figure 4 indicate that some of the ann outputs are highly discrepant ( with 4@xmath2 in error by more than a factor 10 ) . it should be noted that the errors given here apply to individual state - to - state rate coefficients . when these are used to model molecular line emissions , the errors in the resultant line strengths can be expected to be smaller than in the state - to - state rate coefficients because a given radiative transitions may be pumped by a combination of several different collisionally - induced transitions . the calculation made use of the fast artificial neural network ( fann ) library , version 2.0.0 , an open source software package developed by steffen nissen and evan nemerson and available at http://fann.sourceforge.net . the cascade training option was adopted for this problem , and was found to result in rapid training . in this algorithm ( fahlman & lebiere 1990 ) , the training starts with an empty network devoid of hidden neurons . hidden neurons are trained and added to the network one at a time ; each new hidden neuron is chosen from a selection of candidate neurons with different initial weights and one of several possible activation functions .    in using the fann library , the input neuron values were rescaled to small values ( relative to unity ) to ensure that the input neurons did not saturate : energies ( given as wavenumber in units of @xmath78 ) were therefore divided by @xmath79 and rotational quantum numbers were divided by 100 . no other preprocessing was performed . with these rescalings , the input values lay in the ranges [ 0 , 0.14 ] , [ 0.03 , 0.10 ] , and [ 0.16 , 0.18 ] respectively for @xmath39 , @xmath13 , and @xmath40 . the desired output neuron values were the logarithm to base 10 of the rate coefficient in units of @xmath80 , and covered the range [ 15.3 , 9.9 ]    the cascade algorithm was used with all the default settings except for one : a linear function was selected as the activation function for the output layer . this choice was necessitated by the fact that the output neuron values lay outside the range of the other activation functions . the desired error was set to a very small value to ensure that the ann made use of the maximum number of hidden neurons specified for each test . the default settings for the fann software make six activation functions available to the algorithm , all of which are continuous functions with a domain of [ @xmath81 and a range of either [ 0 , 1 ] or [ 1 , 1 ] . the network is optimized using the i - rprop algorithm of igel & hsken ( 2000 ) , a variant of the rprop algorithm introduced by riedmiller & braun ( 1993 ) . even if the settings and training set are unchanged , the results are slightly different every time the ann is run . this is because the initial weights are randomly assigned at the start of training . i therefore averaged the output values over a set of @xmath43 runs . figure 7 shows how the performance of the network depends upon @xmath43 , for a network with 3 input neurons and 16 hidden neurons and with @xmath59 . as in figures 3  5 , black curves apply to the entire test set , while blue , green and red histograms apply to transitions with low , medium and high actual rate coefficients ; for this purpose ,  low " applies to value smaller than @xmath51 ,  medium " to value in the range @xmath82 , and  high "  to a value greater than @xmath53 . the magenta curve refers to the training set . clearly , the performance on the test set is improved if multiple runs are averaged , but for @xmath43 greater than a few , no further improvements are realized . based upon the behavior exhibited in figure 7 , i selected @xmath45 as the standard parameter for the calculations that yielded the results discussed in  3 . figures 8 and 9 show how the performance depends upon the final number of hidden neurons . the color coding in figure 9 is the same as that in figure 7 . clearly , if enough hidden neurons are provided , the fit to the training set becomes extremely accurate ( lower right panel in figure 8) , with the rms error for 64 neurons being less than 3@xmath2 . however , the performance on the _ test _ set worsens somewhat once the number of hidden neurons exceeds 16 . this behavior indicates that the network is  overfitting " the data ( e.g.  smith 1993 ) . the following analogy is presented by the use of polynomial fitting functions . given a function of one variable that has been evaluated at @xmath83 values , a perfect polynomial fit to those values can always be obtained with a polynomial of order @xmath84 . however , such a polynomial may show strong wiggles between ( and beyond ) the points where the function was evaluated ( particularly if the function evaluations have errors ) . a lower order polynomial may fit the data points imperfectly but may be preferable for the purposes of interpolation and particularly extrapolation . based upon the behavior exhibited in figure 8 and 9 , i selected @xmath85 as the standard parameter for the calculations that yielded the results discussed in  3 . i thank george fekete for helpful discussions about ann , and alexandre faure for several valuable comments on an earlier version of the manuscript . i am grateful to the anonymous referee for several useful suggestions . i gratefully acknowledge the support of funding from research support agreement 1349625 issued by nasa s jet propulsion laboratory . dubernet , m. , grosjean , a. , daniel , f. , et al . 2006 , in ro - vibrational collisional excitation database : basecol  http://basecol.obspm.fr ( japan : journal of plasma and fusion research series , series 7 )                        & & & & & & & all & low@xmath86 & med . @xmath86 & high@xmath86 & all & low@xmath86 & med . @xmath86 & high@xmath86 + 2 & 3 & 435 & 555 & o - h@xmath3o & h@xmath3 00@xmath87 & 800 & 0.483 & 0.652 & 0.394 & 0.425 & 0.284 & 0.374 & 0.259 & 0.292 + 4 & 3 & 435 & 555 & o - h@xmath3o & h@xmath3 00 & 800 & 0.461 & 0.634 & 0.373 & 0.359 & 0.266 & 0.386 & 0.237 & 0.208 + 8 & 3 & 435 & 555 & o - h@xmath3o & h@xmath3 00 & 800 & 0.417 & 0.527 & 0.366 & 0.352 & 0.239 & 0.279 & 0.227 & 0.178 + 16 & 3 & 435 & 555 & o - h@xmath3o & h@xmath3 00 & 800 & 0.417 & 0.517 & 0.371 & 0.350 & 0.247 & 0.327 & 0.227 & 0.186 + 32 & 3 & 435 & 555 & o - h@xmath3o & h@xmath3 00 & 800 & 0.446 & 0.559 & 0.391 & 0.411 & 0.259 & 0.345 & 0.235 & 0.207 + 64 & 3 & 435 & 555 & o - h@xmath3o & h@xmath3 00 & 800 & 0.471 & 0.569 & 0.422 & 0.473 & 0.288 & 0.364 & 0.235 & 0.176 +   + 16 & 3 & 45 & 945 & o - h@xmath3o & h@xmath3 00 & 800 & 0.915 & 1.322 & 0.809 & 0.538 & 0.458 & 0.784 & 0.433 & 0.286 + 16 & 3 & 105 & 885 & o - h@xmath3o & h@xmath3 00 & 800 & 0.675 & 1.007 & 0.561 & 0.491 & 0.358 & 0.473 & 0.314 & 0.423 + 16 & 3 & 190 & 800 & o - h@xmath3o & h@xmath3 00 & 800 & 0.474 & 0.653 & 0.416 & 0.299 & 0.291 & 0.402 & 0.279 & 0.143 + 16 & 3 & 300 & 690 & o - h@xmath3o & h@xmath3 00 & 800 & 0.409 & 0.531 & 0.364 & 0.305 & 0.237 & 0.306 & 0.225 & 0.152 + 16 & 3 & 435 & 555 & o - h@xmath3o & h@xmath3 00 & 800 & 0.417 & 0.517 & 0.371 & 0.350 & 0.247 & 0.327 & 0.227 & 0.186 + 16 & 3 & 595 & 395 & o - h@xmath3o & h@xmath3 00 & 800 & 0.391 & 0.430 & 0.377 & 0.307 & 0.227 & 0.264 & 0.222 & 0.139 + 16 & 3 & 780 & 210 & o - h@xmath3o & h@xmath3 00 & 800 & 0.381 & 0.447 & 0.341 & 0.328 & 0.214 & 0.258 & 0.199 & 0.123 +   + 16 & 3 & 435 & 555 & o - h@xmath3o & h@xmath3 00 & 800 & 0.417 & 0.517 & 0.371 & 0.350 & 0.247 & 0.327 & 0.227 & 0.186 + 16 & 6 & 435 & 555 & o - h@xmath3o & h@xmath3 00 & 800 & 0.406 & 0.441 & 0.353 & 0.691 & 0.293 & 0.352 & 0.249 & 0.615 + 16 & 1 & 435 & 555 & o - h@xmath3o & h@xmath3 00 & 800 & 0.725 & 0.906 & 0.613 & 0.889 & 0.420 & 0.466 & 0.393 & 0.829 +   +   +   +   + & & & & & & & all & low@xmath86 & med . @xmath86 & high@xmath86 & all & low@xmath86 & med . @xmath86 & high@xmath86 + 16 & 3 & 495 & 495 & o - h@xmath3o & h@xmath3 00@xmath87 & 800 & 0.433 & 0.577 & 0.353 & 0.372 & 0.234 & 0.314 & 0.209 & 0.239 + 16 & 3 & 495 & 495 & o - h@xmath3o & h@xmath3 02 & 800 & 0.338 & 0.675 & 0.276 & 0.361 & 0.195 & 0.450 & 0.174 & 0.224 + 16 & 3 & 495 & 495 & o - h@xmath3o & h@xmath3 20 & 800 & 0.325 & 0.368 & 0.274 & n / a & 0.188 & 0.229 & 0.163 & n / a + 16 & 3 & 495 & 495 & o - h@xmath3o & h@xmath3 22 & 800 & 0.342 & 0.457 & 0.329 & 0.368 & 0.210 & 0.319 & 0.202 & 0.221 +   + 16 & 3 & 540 & 541 & o - c - c@xmath67h@xmath3 & he & 30 & 0.325 & n / a & 0.269 & 0.414 & 0.208 & n / a & 0.171 & 0.316 + 16 & 3 & 602 & 603 & so@xmath3 & he & 25 & 0.341 & 0.361 & 0.327 & 0.369 & 0.205 & 0.220 & 0.187 & 0.368 + 16 & 3 & 346 & 346 & o - h@xmath3co & he & 30 & 0.288 & 0.636 & 0.270 & 0.221 & 0.168 & 0.646 & 0.155 & 0.163 +        fig .  1  simple example of a cascade - correlation network , with 3 input neurons ( blue circles numbered 1  3 ) , 2 hidden neurons ( numbered 4 and 5 ) , and one output neuron ( number 6 ) . the @xmath16th neuron has an input value @xmath17 ( red ) and an output value @xmath18 = @xmath19 ( green ) , where @xmath20 is an activation function ( which may , in principle , be different for each neuron  see appendix ) . information propagates from left to right , starting with the input values @xmath21 , @xmath22 and @xmath23 , and ending with the output value @xmath24 . the @xmath17 and @xmath18 are related by a matrix of weights , * w * , with @xmath25 . the algebraic expressions shown on the figure show that the quantities @xmath88 , @xmath89 , @xmath90 , @xmath91 , @xmath92 , @xmath93 , @xmath94 , @xmath95 , and @xmath24 can be computed non - recursively ( i.e.  explicitly , in the above order ) , given input values of @xmath21 , @xmath22 , and @xmath23 .        fig . 2  the top right , bottom left , and bottom right panels present the performance of an ann with @xmath46 for the cases with 1 , 3 , and 6 input neurons discussed in  2 . in each panel , the horizontal axis represents the logarithm of the rate coefficient ( in units of @xmath47 ) obtained by dubernet et al . ( 2009 ) , and the vertical axis represents the corresponding output of the ann . red points apply to the 435 data points in the training set , while blue points apply to the 555 data points in the test set ( i.e.  to those transitions that were not seen by the ann during training ) . upper right panel : analogous results if the data are all fit with an expression of the form @xmath48 . 3  the frequency distribution of the errors , @xmath50 , for each of the 4 cases shown in fig .  1 . black histograms apply to the entire test set , while blue , green and red histograms apply to transitions with low , medium and high actual rate coefficients ; for this purpose ,  low " applies to value smaller than @xmath51 ,  medium " to value in the range @xmath82 , and  high "  to a value greater than @xmath53 . fig .  4  fraction of the test data for which @xmath55 exceeds a given value , @xmath56 . the fraction of rate coefficients is shown on the vertical axis ( logarithmic scale ) , with the same color coding adopted in figure 2 , as a function of @xmath56 . the left panel is simply a zoomed - in version of the right panel . fig .  5  rms value of @xmath60 , as a function of n@xmath61 , for the entire set of transitions ( black points ) and for low , medium and high values of the actual rate ( blue , dark green and red as in figs .  2 and 3 ) . 6  results obtained for excitation of so@xmath3 ( from green 1995 ) , ortho - c - c@xmath67h@xmath3 ( chandra & kegel 2000 ) , and ortho - h@xmath3co ( green 1991 ) in collisions with he ; and for excitation of ortho - water in collisions with para - h@xmath3 ( dubernet et al . ( 2009 ) for cases in which the initial and/or final state of h@xmath3 is @xmath68 . fig .  7  dependence of the performance of the ann upon @xmath43 , for a network with 3 input neurons and 16 hidden neurons and with @xmath59 . black curves apply to the entire test set , while blue , green and red histograms apply to transitions with low , medium and high actual rate coefficients ; for this purpose ,  low " applies to value smaller than @xmath51 ,  medium " to value in the range @xmath82 , and  high "  to a value greater than @xmath53 . the magenta curve refers to the training set .  performance of an ann with 3 input neurons and @xmath59 . four panels shows the results for 8 , 16 , 32 , and 64 hidden neurons . in each panel , the horizontal axis represents the logarithm of the rate coefficient ( in units of @xmath47 ) obtained by dubernet et al . ( 2009 ) , and the vertical axis represents the corresponding output of the ann . red points apply to the 435 data points in the training set , while blue points apply to the 555 data points in the test set ( i.e.  to those transitions that were not seen by the ann during training ) . fig .  9  dependence of the performance of the ann upon the number of hidden neurons , for a network with 3 input neurons and @xmath59 . black curves apply to the entire test set , while blue , green and red histograms apply to transitions with low , medium and high actual rate coefficients ; for this purpose ,  low " applies to value smaller than @xmath51 ,  medium " to value in the range @xmath82 , and  high "  to a value greater than @xmath53 . the magenta curve refers to the training set .
Abstract:

Answer: an artificial neural network ( ann ) is investigated as a tool for estimating rate coefficients for the collisional excitation of molecules . the performance of such a tool can be evaluated by testing it on a dataset of collisionally - induced transitions for which rate coefficients are already known : the network is trained on a subset of that dataset and tested on the remainder . results obtained by this method are typically accurate to within a factor @xmath0 ( median value ) for transitions with low excitation rates and @xmath1 for those with medium or high excitation rates , although 4@xmath2 of the ann outputs are discrepant by a factor of 10 more . the results suggest that anns will be valuable in extrapolating a dataset of collisional rate coefficients to include high - lying transitions that have not yet been calculated . for the asymmetric top molecules considered in this paper , the favored architecture is a cascade - correlation network that creates 16 hidden neurons during the course of training , with 3 input neurons to characterize the nature of the transition and one output neuron to provide the logarithm of the rate coefficient .