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 process of oxidation of metal surfaces is of considerable fundamental scientific interest as well of paramount technological importance.@xcite corrosion and passivation are just two examples of either destructive or useful processes linked directly to this phenomenon , which are well - known from everyday life . oxidation of metal surfaces begins with dissociative chemisorption of oxygen on a clean surface and is followed by the formation of a film of metal oxide . in - between these events there are many elementary processes . the study of al and mg oxidation is particularly important because al and mg belong to the group of so - called simple metals and are considered model systems for studies of oxidation of transition metals . both metals exhibit high reactivity with oxygen and oxidize rapidly . aluminum oxides ( alumina ) form several different phases where the structure of some of them was only recently identified by combined _ ab initio _ density functional theory ( dft ) and experimental studies.@xcite magnesium oxide is known to eventually form crystals of the rock - salt structure,@xcite but can also experience complex reconstructions at partial coverages.@xcite both mg and al oxides play an important role in catalyst support and find many other useful applications . the number of experimental@xcite studies and theoretical@xcite treatments devoted to the oxidation of magnesium is relatively limited compared to the more extensively studied oxidation of aluminum . in this paper we present a systematic study of the initial oxygen incorporation and island formation by first - principles theory calculations and we use the results to interpret high - quality x - ray photo - electron diffraction ( xpd ) measurements for low o@xmath1 doses .    the existing picture of the mg(0001 ) surface oxidation processes dates back to the early 1980 s . based on the extensive low - energy electron diffraction ( leed ) , auger electron spectroscopy , electron energy loss spectroscopy , and work function measurements namba et al . proposed  @xcite a four - stage model consisting of : ( 1 ) dissociative oxygen chemisorption at random sites followed by oxygen incorporation , ( 2 ) assembly of incorporated oxygen atoms into islands and lateral growth , ( 3 ) oxide formation from surface mg atoms and subsurface o atoms which starts at o@xmath1 exposures around 23 l ( langmuir ) or 0.51 monolayer of total coverage , and ( 4 ) oxide thickening . several other experimental studies provided arguments for a three - stage model of oxidation that was proposed at approximately the same time . in this model @xcite atomic oxygen is directly incorporated into magnesium ( below the top mg layer ) right after on - surface dissociation . in the next step there is simultaneous formation of an oxide layer and a decrease in the oxygen on the surface.@xcite finally , there is oxide thickening and transformation into a rocksalt structure . this model of immediate incorporation of oxygen gained support from the x - ray photo - electron spectroscopy study of ghijsen et al.,@xcite and particularly from the measured sharp decrease in the work function of the mg(0001 ) surface upon initial o exposure.@xcite the formation of an o-(@xmath2 ) subsurface layer up to monolayer coverage and its subsequent transformation into epitaxial oxide was also reported.@xcite the positions of the subsurface oxygen atoms were deduced to be octahedrally coordinated where the oxygen has an environment similar to that in mgo.@xcite this is in agreement with a more direct experimental evidence for oxygen incorporation in the octahedral interstitial sites of the first two interlayer spacings of mg during the initial stages of magnesium oxidation,@xcite and with later inelastic ion scattering experiments on a polycrystalline mg surface.@xcite a very recent scanning tunneling microscopy ( stm ) study of the oxidation of mg(0001 ) suggests that at low oxygen exposures ( up to 2 l ) the incorporated oxygen atoms form a single layer underneath the top layer of mg.@xcite the idea of immediately populated subsurface sites seems also to be supported by recent measurements by mitrovic et al.,@xcite however , they show that most of the oxygen ( 90% ) remains over the surface and only a small fraction goes below the surface . thus , the question regarding the most favored adsorption sites is still debated and requires further analysis .    on the theoretical side , the local density - functional theory ( dft ) calculations of bungaro et al.@xcite show that oxygen is incorporated below the surface forming a @xmath3 subsurface lattice . a recent dft - based lattice gas model study@xcite of island formation in the early mg(0001 ) oxidation stages has shown that at very low filling the oxygen atoms adsorb in the top - most subsurface ( tetrahedral ) b - sites in an abab ... stacking of magnesium . at higher oxygen load also the deeper sites are filled , and the adsorbed oxygen atoms form dense clusters just below the surface and in further subsurface locations.@xcite for larger clusters the oxygen - oxygen interaction effects drive the oxygen atoms further into the magnesium subsurface layers and bulk . the rest of this paper is organized as follows . in section ii systematic dft calculations in the range of coverages between 1/16 and 3 monolayers ( ml ) are presented and used for the study of the initial oxygen incorporation and island formation . the previous experimental and theory results are the starting base for our first - principles calculations of the very beginning of o adsorption and oxide formation at the mg(0001 ) surface . in section iii we describe our high - quality xpd measurements for low o@xmath1 doses corresponding approximately to global coverages between 0.1 and 2 ml of atomic oxygen , and in section iv our dft results are used to interpret the experimental xpd data . the first - principles calculations are carried out using the plane - wave density - functional dacapo code@xcite with the generalized gradient approximation ( gga ) for the exchange - correlation energy functional@xcite and with ultrasoft pseudopotentials@xcite to represent the ionic cores . the clean mg(0001 ) surface is modeled by periodic slabs consisting of six to seven magnesium layers separated by 21  of vacuum . the plane - wave basis set with 25 ry energy cutoff the o atoms are adsorbed on one side of the slab only and the electric field arising due to the asymmetry of the system is compensated for by a consistent dipole correction.@xcite a @xmath4 mesh of monkhorst - pack special @xmath5-points for the @xmath6 atom surface unit cell and a fermi - surface smearing of 0.2 ev are applied to the brillouin - zone integrations . for smaller cells the number of k - points are increased accordingly up to a @xmath7 mesh for the @xmath2 surface unit cell . the positions of atoms in the three to four topmost magnesium layers , and of all the oxygen atoms , are fully optimized until the sum of the hellmann - feynman forces on all unconstrained atoms converges to less than 0.05 ev /  . the forces acting on the ions in the unit cell are derived from the converged charge densities , and the atom dynamics is determined using a preconditioned quasi - newton method based on the broyden - fletcher - goldfarb - shanno algorithm.@xcite    the study involves very low coverages of oxygen , from one o atom in a large surface cell up to coverages of three monolayers . we define the coverage @xmath8 as the ratio of the number of adsorbed atoms to the number of atoms in an ideal substrate layer . the binding energy per adsorbate is calculated relative to the energy of an isolated , spin - polarized o@xmath1 molecule and the clean , relaxed mg - surface .      before studying the effect of oxygen chemisorption on the mg(0001 ) surface , we determined the bulk and bare - metal - surface structures . the calculated lattice constants for hcp mg is @xmath9   and @xmath10 , which agree very well with the experimental values @xmath11   and 1.623 ( ref .  ) and other gga calculations.@xcite the bulk modulus , ignoring lattice vibrations , is found to be 34.1 gpa,@xcite in good agreement with the measured value 35.5 gpa ( obtained from the elastic constants of ref .  ) and with the gga calculations of ref .  . the relaxations of the surface interlayer spacing with respect to the bulk spacing show 1.5% expansion , 0.4% contraction , and 1.3% contraction of the first , second and third interlayer distance , respectively . the expansion of the first interlayer distance compares well with the experimental@xcite value of 1.7% , and with other gga calculations.@xcite the small contraction of the second interlayer spacing disagrees with results of experiment and other calculations which predict a small expansion . since the magnitude of this relaxation is very small this discrepancy should not have any significant effect on the results for o - adsorption . the calculated work function is 3.72 ev for the relaxed surface and agrees well with the experimental value 3.84 ev . @xcite      having obtained reliable results for the clean , relaxed mg(0001 ) surface we calculate the binding of oxygen in on - surface and subsurface sites . the locations of the possible adsorption sites of mg(0001 ) are sketched with crosses in fig . [ figsketch ] . in our calculations the results for single or nearest - neighbor pairs of o atoms adsorbed at the @xmath6 , @xmath13 , and @xmath14 surface cells were used as 1/16 , 1/8 , 1/4 , and 1/2 ml coverage data .    in fig . [ f2 ] we present the calculated low - coverage binding energies with respect to the energy of an isolated o@xmath1 molecule . the adsorption of oxygen atoms on the mg(0001 ) surface induces further changes to the surface structure and the electronic properties . however , the positional changes of the mg atoms are in most cases not very large . for the @xmath13 supercell with a single o - atom slightly below the surface ( site b2 in fig .  [ figsketch ] ) or in first - sublayer adsorption sites ( sites a4 or a5 ) the mg atoms move less than 3% vertically and 0.6% laterally from their clean - surface positions . to estimate the sensitivity of the binding energy to positional relaxations of the atoms we also calculated an even distribution of four oxygen atoms in one of the tetrahedral subsurface sites ( sketched as site b2 in fig . [ figsketch ] ) in the @xmath6 system , and compared the binding energy to the binding energy of one site - b2 atom in the @xmath13 system . the two configurations relax to slightly different atomic positions ( but both converged within the force requirement ) and we found the difference in binding energy per oxygen atom to be 5 mev . we therefore estimate the relaxation - imposed inaccuracy of the binding energies to be approximately 5 mev .    _ on - surface adsorption_. it has previously been suggested that for low coverages the o atoms are preferentially adsorbed into _ on - surface _ sites,@xcite as is the case for , e.g. , oxygen adsorption on al(111),@xcite or that on - surface adsorption sites exist although having lower binding energies than the lower - lying adsorption sites.@xcite for the al(111 ) surface an additional o is incorporated only when 1 ml of on - surface o is complete.@xcite a similar situation is found for a transition metal , the ru(0001 ) surface.@xcite for this reason on - surface adsorption was closely studied . in summary , however , we find the initial oxygen adsorption to be a subsurface process .    of the on - surface adsorption sites the natural candidates are the c1 and b2top hollows ( fig . [ figsketch ] ) . however , our calculated low - coverage binding energies , plotted in figure [ f2 ] , show that for single atom adsorption at coverages @xmath15 ( corresponding to a single o atom in the @xmath6 , @xmath13 , and @xmath14 surface unit cells ) the only stable on - surface adsorption site is the one denoted as c1 , and the binding energy of this site is significantly less than that of any of the subsurface sites . we find _ no _ single - adsorbant energy minimum near the position where we would expect the on - surface b2top - site . instead , when relaxing dilute systems with o in the b2top sites the surrounding mg atoms move laterally by a small amount , enough to let the o - atoms sink into the subsurface b2 ( tetrahedral ) site . only if a pair of o atoms are adsorbed as nearest neighbors , or at increased o loading ( @xmath16 ) , can one of the o atoms be stabilized on the surface in the b2top site , albeit at smaller binding energy than in the slightly lower - lying b2-site . thus , for low coverages our dft results point to the c1 site as the only possible , but energetically unfavorable , on - surface adsorption site . these results are in agreement with the theory results of bungaro et al.@xcite as to subsurface adsorption sites being more favorable than on - surface ones , but contradict their prediction of a stable on - surface b2top - site for low coverages . it should be noted , however , that for coverages @xmath17 bungaro et al.@xcite relaxed only the o atoms but not the mg - lattice , and that their @xmath5-point sampling was limited to the @xmath18 point only . thus our results underline the importance of relaxing the substrate . _ the subsurface adsorption_. by extensive searches for subsurface adsorption sites we find that within the two top mg - layers the tetrahedral sites ( sites b2 , a4 , and a5 ) have higher binding energies than the on - surface c1 and the subsurface octahedral c3-site for all coverages @xmath19 ( compare fig .  [ f2 ] and table  [ t1 ] ) . thus , the first o atom adsorbed in the oxidation process binds in a subsurface site of type b2 , its binding energy being only 0.09 ev and 0.15 ev higher ( fig .  [ f2 ] ) than that of o in subsurface site a4 and a5 , respectively . the preference for the b2 site agrees with the results of ref .  , but we find a binding energy about 0.6 ev larger . in addition , they find that the on - surface b2top site is much more favored than the subsurface a4 site , in clear disagreement with our results . these difference are likely due to the lack of mg - atoms relaxation and the inferior @xmath5-points sampling used in that study . .[t1 ] the binding energy @xmath20 per o - atom for a double - site o occupancy of different chemisorption sites ( see fig .  [ figsketch ] ) and for different surface unit cells . two o atoms in the @xmath6 and @xmath13 surface unit cells correspond to the coverages @xmath21 and 1/2 , respectively . the exact vertical position of the sites b2 are distinguished by  top " for a position ( slightly ) above the top mg - layer , and  mix " for a position in which some of the top mg - atoms are above and some are below the oxygen atom , whereas b2 denotes adsorption below all atoms of the top - most mg - layer . all pairs of sites are nearest neighbors . [ cols="<,^,<,^ " , ]     in order to check the possibility of formation of ionic complexes ( which would represent a seed for further oxide formation ) consisting of the mg atoms of the topmost layer and subsurface and on - surface oxygen , we calculated the most preferred sites for an additional o atom in the presence of an already adsorbed one . in studying possible combinations of the oxygen adsorption sites we concentrated on the ones including the b2 sites . as is evident from table  [ t1 ] , even at low coverages with one o - atom isolated from other o - atoms the addition of another o atom leads to an increased o - population in the b2 sites . the two ( neighboring ) atoms in the b2 sites are by @xmath220.030.04 ev more favored than one atom in the b2 site and another one in either the a4 or a5 site .      for 1 ml coverage with occupation in only one kind of site the energetic ordering of the most stable sites is different than for lower coverages ( fig . [ f2 ] and table  [ thigh ] ) . the b2-type sites become less favorable than the a4 site and the actual b2 site is shifted to the b2top position . however , allowing for simultaneous occupancy of the b2 and b2top sites at 1 ml in the @xmath14 unit cell ( table  [ thigh ] ) we find that this combined b - site structure is by 0.06 ev more favored than the single a4 occupancy shown in fig . [ f2 ] . lcrccc sites & @xmath20 & @xmath24 & @xmath25 & @xmath25 + & [ ev ] & [ ev ] & ( 1.4 l)&(9.7 l ) +   + b2 & a5 & 5.25 & -0.64 & 0.41 & 0.30 + b2mix & a4 & 5.21 & -0.62 & 0.34 & 0.37 + b2top & c3 & 5.03 & -0.49 & 0.51 & 0.43 +   + b2top & b2 & 5.39 & -0.46 & 0.46 & 0.46 + a4 & 5.33 & 0.00 & 0.42 & 0.48 + b2top & 5.27 & -0.99 & 0.44 & 0.46 + a5 & 5.23 & 0.11 & 0.53 & 0.42 + c1 & 5.15 & -0.12 & 0.44 & 0.46 + c3 & 5.04 & 0.62 & 0.56 & 0.51 +   + a4 + b7 & 5.61 & -0.07 & 0.41 & 0.50 + a5 + b8 & 5.55 & 0.27 & 0.56 & 0.45 + b2 + a5 & 5.54 & -0.03 & 0.44 & 0.33 + rocksalt 2ml & 5.45 & 0.72 & 0.54 & 0.45 + flat 2ml & 5.42 & -1.14 & 0.29 & 0.36 +   + a4 + b7 + a10 & 5.70 & 0.04 & 0.45 & 0.50 + a5 + b8 + a11 & 5.67 & 0.09 & 0.53 & 0.43 + b2 + a5 + b8 & 5.66 & 0.04 & 0.46 & 0.35 + rocksalt 3ml & 5.61 & 0.78 & 0.55 & 0.50 + flat 3ml & 5.47 & -1.14 & 0.30 & 0.37 +    as more oxygen is adsorbed into the mg surface the possible adsorbate configurations grow in number , partly due to simple combinatorics , and partly because new adsorption sites appear within an increasingly distorted mg background lattice . this is already apparent for nearest - neighbor adsorption in the low - coverage region ( table  [ t1 ] ) where a slight vertical distortion of the top mg - layer in some cases leads to adsorption in the b2top site and in a ` mixed ' b2-b2top position . in table [ thigh ] we name the 1 ml configurations by the adsorbate position , but for higher coverages the notation is merely an approximate description of the structure , and we refer the reader to fig .  [ figmediumstructures ] for a sketch of the multitude of similar but not identical configurations at coverages 23 mls . the stable structures found in the 2 and 3 ml coverage region fall in two general classes : layered structures with oxygen and magnesium in buckled or flat layers having a hexagonal structure [ fig . [ figmediumstructures ] ( a)(c ) , and ( e ) ] , and rock - salt structures [ fig .  [ figmediumstructures ] ( d ) ] on top of clean magnesium , strained with respect to the clean - surface lattice constant . the buckled structures are labeled according to the approximate positions of the oxygen atoms ( comparing to the low - coverage adsorption sites of fig . [ figsketch ] ) . the buckled structures involve subsurface oxygen only . the vertical separation of o and mg within a buckled layer is approximately 0.6  , and the layers can exist both on top of the clean mg(0001 ) [ fig . [ figmediumstructures](c ) ] and as subsurface layers [ fig .  [ figmediumstructures](a ) and ( b ) ] , and at 2 or 3 ml coverage . thus one can discriminate the structures buckled up and down , depending on the oxygen atoms staggering over or below the mg layer , respectively . the flat structures [ fig . [ figmediumstructures](e ) ] have vertical o - mg separation 0.010.06  , and exist for both 2 and 3 ml coverage , but do not exist as subsurface layers . in the flat structures the top - most oxygen atom sits slightly above the top mg layer , corresponding approximately to the position of site b2top in fig . [ figsketch ] . the binding energies of the 13 ml structures are given in table  [ thigh ] . note that for several other configurations with 2 ml or 3 ml coverage of the tetrahedral adsorption sites no energy minimum is found . this includes any subsurface flat phase , phases of mixed flat and buckled layers , and mixed buckle - up and buckle - down layers . further , we tested structures that include an octahedral ( c3 , c6 , c9 ) site in the 2 ml and 3 ml coverage and found binding energies per atom at least 0.30.4 ev smaller than the energetically best structures at the same coverage , similar to the results of the low coverage results ( fig . [ f2 ] ) . the preference for subsurface adsorption at the mg(0001 ) surface was originally suggested by a sharp work function decrease upon oxidation.@xcite previous local - density - approximation dft calculations@xcite have shown that only oxygen incorporation into b2 sites lowers the work function and for one monolayer load this lowering was found to be 0.3 ev , which is much smaller than measured experimentally . we do not find 1 ml adsorption in the b2 sites , but for the neighboring b2top sites our results presented in table [ thigh ] confirm the lowering of the work function . however , this lowering is three times larger than that reported in ref . and we note that our b2top sites are slightly above and nearly coplanar with the topmost mg layer . this larger decrease of the work function upon oxygen adsorption agrees with recent measurements.@xcite    it is generally expected that an adlayer of negative o ions will increase the metal - surface dipole layer and the work function . for example , this is clearly the case for oxygen on the al(111 ) surface . @xcite our results suggest that 1 ml of on - surface ( b2top ) oxygen on the mg(0001 ) surface leads to a decrease in the work function . this decrease for on - surface o is consistent with experimental observations.@xcite on the other hand , for 1 ml coverage the b2top sites are not the energetically most favorable . we conjecture that the following scenario takes place . at low coverage the only on - surface sites are the fcc hollows ( c1 ) , but also the lower - lying subsurface b2 sites are exposed to o atoms arriving at the surface . thus an o - atom that arrives to the surface will find the b2 site having a higher adsorption energy than the c1 site and will immediately diffuse to the former . the b2-sites are most favorable up to 0.5 ml coverage ( fig . [ f2 ] ) . with growing o - coverage the energetic situation changes and the deeper tetrahedral sites become more favorable . however , there is an energy barrier to overcome in order to get into these lower - lying a4 and a5 sites so for 1 ml coverage , instead of moving further into the surface , the oxygen atoms stay in the neighborhood of the b2 sites , now shifted into the b2top position . the b2top site is energetically slightly less favored than the tetrahedral subsurface sites ( by about 0.06 ev ) but requires no crossing of significant energy barriers . with an increased oxygen load ( around 2 ml ) the energy barrier is lowered and some of the oxygen atoms can move into the a5 sites . again , the occupation of the b2top and a5 sites ( the  flat 2ml " structure ) is not the most favorable energetically ( table  [ thigh ] ) but presumably this metastable state is most easily available for the o atoms . the work function change @xmath26 for this configuration is similar to @xmath26 for the energetically favorable a4+b7 subsurface structure . with a further increase of o load ( 3 ml ) the energy barriers that separate the lower lying sites are lowered and oxygen populates the b8 sites . it seems that this scenario is supported by our xpd experiments presented in sec . iii .    _ island formation_. in a previous theory study based on a lattice - gas - model@xcite it was shown that oxygen atoms that are adsorbed in the mg(0001 ) surface form dense clusters in tetragonal subsurface sites . for example , at a dosing corresponding to a global coverage of 1 ml , the oxygen atoms were found to be distributed , on average , with 55% in site b7 , 30% in site a4 , and 15% in site a10 and with local coverages of 23 mls within the clusters .    in the present , more detailed investigation we find that past the adsorption of the first , isolated o - atoms in the oxidation process the o - atoms show a tendency for island formation by pairing up and showing higher binding energies in combinations with subsurface sites . an inspection of the binding energies presented in fig . [ f2 ] and in table  [ t1 ] shows that the higher the coverage , the larger the binding energy of the lower - lying sites ( sites a4 and a5 ) . our results show that in general the binding energy is higher in more close - packed systems , thus showing that clusters or islands of o atoms are preferred to isolated or pairs of o atoms . in fact , for 1 ml coverage , having all o atoms in site a4 is preferred by a small amount to all other subsurface sites and to the on - surface b2top site . this agrees well with the previous lattice - gas - model study.@xcite with coverage increasing above 0.5 ml , the o - atoms from the b2 sites are pushed up into the surface b2top sites . thus they end up , at 1 ml load , in the b2top - sites located almost coplanarly with the mg atoms of the topmost layer , with the binding energy comparable to that of o in sites a4 and a5 . the preference for the formation of subsurface islands in tetrahedral a4 sites or lower agrees both with the lattice - gas model simulations and with the results of bungaro et al . @xcite but disagrees with the results of thiry et al . @xcite who found adsorption in octahedral sites . a new feature that was not observed by bungaro et al . , and which is probably connected with a more careful treatment of the lattice relaxation effect in our work , is the appearance of the on - surface o - islands for 1 ml coverage , with binding energy slightly lower than in subsurface islands of a4-type . also worth noting are the relatively small binding energy differences for oxygen in different sites for @xmath27 . the tendency to form islands continues with further increase of coverage , showing some quenching in the adsorption energy at 3 ml coverage . a similar behavior was found in the lattice - gas model simulations . @xcite    _ layered structures_. as the coverage is increased beyond a single monolayer the layered structures of fig . [ figmediumstructures](a)(c ) and ( e ) appear , along with a rocksalt surface structure , fig .  [ figmediumstructures](d ) . the layered structures can be described as a mg(0001 ) surface with slightly changed atomic positions to accommodate o atoms within the top layers of the mg surface , thereby creating a layered structure . each of the layers has a honeycomb structure with o at every second vertex and mg at the remaining vertices . the layers are stacked in an aa@xmath28 stacking sequence , with o on top of mg atoms and _ vice versa _ at a layer separation of approximately 2.32.7  . again , we find ( fig . [ figmediumstructures ] ) that subsurface layered structures are preferred . comparing the adsorption binding energies listed in fig . [ figmediumstructures ] it is evident that with increased o coverage , the difference in binding energies for the subsurface layered and rocksalt structures is diminishing , indicating that eventually , for still higher coverage , the rocksalt structure may become the most favorable one . _ directional bonding_. figure  [ figdens ] displays the electron charge redistribution for three of the 2 ml structures . whereas the rocksalt structure ( fig .  [ figdens](b ) ) has an almost isotropic charge redistribution around each atom , as is typical for ionic binding , the two layered structures show more complicated , anisotropic charge redistributions . the flat layered structure ( fig .  [ figdens](c ) ) shows most clearly the anisotropy , with more charge redistribution along the in - plane nearest - neighbor mg - o lines than in other directions . to supplement the dft calculations and to provide a reference ground we have also carried out xpd measurements of the clean mg(0001 ) surface and after low dosing of oxygen . xpd has been chosen because of its chemical sensitivity and its sensitivity to local real space order . it is a powerful technique for surface structural investigations,@xcite and it has been shown that full - hemispherical xpd patterns provide very direct information about the near - surface structure . at photoelectron kinetic energies above about 500 ev , the strongly anisotropic scattering by the ion cores leads to a forward focusing of the electron flux along the emitter - scatterer axis . the photoelectron angular distribution , therefore , is to a first approximation a forward - projected image of the atomic structure around the photoemitters . analysis of the symmetry and positions of forward - focusing maxima thus permits a very straightforward structural interpretation of xpd data . furthermore , detailed structural parameters can be determined by comparing the experimental xpd patterns to calculated ones . the relatively simple and efficient single - scattering cluster ( ssc ) formalism@xcite has proven adequate in most cases . the agreement between ssc calculations and experimental xpd pattern can be quantified using a reliability factor such as the _ r_-factor @xmath25 defined previously.@xcite    the experiments have been performed in the university of fribourg s vg escalab mk ii spectrometer modified for motorized sequential angle - scanning data acquisition . clean mg(0001 ) surfaces have been prepared by cycles of sputtering ( 500 ev ar@xmath29 ) and annealing ( 130  @xmath30c ) . before o@xmath1 exposure , the surface displayed a sharp @xmath3 leed pattern with little background . after o@xmath1 exposure with the sample held at room temperature , the coverage was determined from the relative intensities of the o 1__s _ _ and mg 2__p _ _ photoelectron peaks . experimental o @xmath31 xpd patterns were obtained after exposure of the mg(0001 ) surface to 0.15 , 0.7 , 1.4 and 9.7  l of o@xmath1 , corresponding to _ global _ o coverages of 0.1 , 0.4 , 1 and 1.7 ml , respectively . the mg(0001 ) surface with various coverages of oxygen atoms was studied with leed and xpd . experimental xpd patterns for o@xmath1 doses of 0.15 and 9.7  l are shown in fig . [ figexpssc ] . the most prominent feature of the 0.15 l pattern is a strong intensity maximum at normal emission . this observation is a clear and direct evidence for o adsorption below the topmost mg layer : only an atom located directly above the oxygen photoelectron emitter can give rise to a forward - focusing maximum at normal emission , i.e. , in the center of the plot . a further conclusion can immediately be drawn from the evolution of the xpd pattern with o@xmath1 exposure . apart from differences due to counting statistics , the experimental xpd patterns are strikingly similar , not only regarding the intensity maximum at normal emission but regarding all the prominent diffraction features . it must therefore be concluded that the local atomic geometry is the same over the entire exposure range . the same conclusion is obtained from the evolution of the leed pattern with increasing o@xmath1 exposure : the leed pattern stays ( @xmath32 ) , and only the background rises slightly , which indicates an increasing amount of disorder with increasing coverage . since no superstructures are observed , the presence of any other islands than islands of ( @xmath32 ) periodicity that are larger than about 50100   can therefore be excluded . in other words , both the leed and xpd measurements indicate that locally the geometry is the same irrespective of the global coverage ( at least below 1.7 ml ) . already at 0.1 ml global coverage this local structure thus contains at least 1 ml oxygen , and with increasing coverage it will grow in domain size without significant changes in geometrical structure . therefore , determining the local geometrical structure at 1ml local coverage will determine the geometries even at the lowest global coverages accessible to the experiments .    in the further analysis we compare simulated ssc diffraction pattern to experimental xpd images of the mg(0001 ) surface both with dosing 1.4 l ( approx . 1 ml global coverage ) and 9.7  l ( approx .  1.7 ml global coverage ) . that the o / mg(0001 ) system does not simply form a mgo(111 ) rocksalt structure is clearly seen by comparing a simulated diffraction pattern for mgo(111 ) [ fig . [ figexpssc](c ) ] to the experimental 9.7  l xpd pattern shown in fig .  [ figexpssc](b ) : the simulation does not reproduce any of the diffraction features satisfactorily . therefore , the most relevant atomic structures , as determined by the dft calculations , were used to simulate xpd patterns by means of ssc calculations . the resulting ssc calculations were compared to the experimental 1.4  l and 9.7  l xpd patterns , the latter of which is shown in fig . [ figexpssc](b ) . none of the ssc calculations for an oxygen atom in an isolated site gave satisfactory agreement with experiment . the calculations , however , confirm the conclusion drawn from the experimental patterns that one oxygen atom is adsorbed directly below a mg atom : only the simulations for oxygen in the a4 and a5 sites exhibit the characteristic intensity maximum at normal emission that is seen in experiment [ fig . [ figexpssc](a),(b ) ] . the fact that the experimental xpd patterns do not agree with simulated diffraction patterns of disperse , low local coverages even for low overall coverage in the experiment is consistent with the leed measurements discussed above , indicating again that the oxygen will form island with a local coverage of 1 ml or more . accordingly , further ssc calculations concentrated on the atomic positions obtained from dft calculations using a @xmath2 surface cell .    among the simulated diffraction patterns for 1 ml local coverage only the site a4 and the site a5 ssc calculations exhibit a central maximum . in contrast , the ssc calculations for each of the c3 , b2top , and c1 sites ( leaving all other sites empty ) do not show such maximum . this means that the comparison between ssc calculations and experimental xpd patterns thus disqualifies the site b2top - only occupation at 1 ml local coverage , in agreement with the theoretical binding energy being smaller for b2top than for a4 at 1 ml coverage . along with the binding energy obtained from the dft calculations the comparison to the xpd patterns also rules out octahedral subsurface ( site c3 ) occupation , in disagreement with the results of thiry et al . @xcite    consequently , ssc calculations were also performed for the atomic positions obtained from dft calculations considering 2 and 3 ml coverage . in table [ thigh ] we summarize the 1 ml , 2 ml and 3 ml local coverage ssc simulations [ one to three o atoms per @xmath3 cell ] by listing their reliability factors with respect to the 1.4  l and 9.7  l exposure xpd data . best agreements between experiment and ssc calculations ( lowest reliability factors ) are found for the layered structures `` b2mix & a4 '' and `` b2 & a5 '' at 1 ml local coverage , `` flat 2ml '' [ fig . [ figexpssc](d ) ] and `` b2 + a5 '' at 2 ml local coverage , and `` flat 3ml '' and `` b2 + a5 + b8 '' for 3 ml local coverage . these results are consistent with the work function changes discussed above . the general trend is that the rocksalt structures ( 2 ml and 3 ml ) do not agree very well with the experimental xpd patterns ( high values of @xmath25 ) , whereas some of the layered structures agree better . in particular we notice that among the 2 ml structures the flat structure ( `` flat 2ml '' ) has the best agreement with the 1.4  l experiment , whereas one of the buckled layer structures ( `` b2 + a5 '' ) fits best to the 9.7  l experiment . this indicates that as the global coverage is increased by higher o@xmath1 dosing , more of the surface becomes covered with buckled mg and o layers and less with the flat surface . for the 1.4  l dosing experiment the flat surface is favored , in disagreement with the dft results . 1.4  l dosing corresponds to @xmath221 ml of global coverage , but since the distribution of oxygen is nonuniform due to the tendency of the atoms to form islands the experimentally observed surface is patched . for small surface oxide patches of 2 mls ( or more ) local oxygen coverage the strain can get released over most of the island and thus , experimentally , mostly the flat structure is seen at low dosing . this gives a low @xmath25-value for the flat structure , but still keeps the order of the buckled and rocksalt structure found in dft calculations : buckled is also at 1.4 l more favored experimentally than rocksalt . thus , the relevant dosing ( among the measured ones ) to use for comparing dft and the xpd experimental structures is 9.7  l. this structure probably does not fully cover the surface ( if clusters are at least 2 ml dense ) , but certainly does so to a larger extent than for dosing 1.4  l. at dosing 9.7  l the dft calculations and the xpd experiment agree on the buckled surface .    based on both the calculated and experimental structures discussed above we conclude that at relatively low dosage ( corresponding to 23 ml coverage ) o / mg(0001 ) forms the layered oxide structure . the rocksalt structure typical of mgo starts to grow only at higher o dosing . we have performed extended first principles calculations of oxygen adsorption and of the initial stages of mg(0001 ) oxidation . a variety of configurations and a wide range of coverages were considered in order to determine the most stable structures . these were compared with x - ray photoelectron diffraction experiment and simulations . at low coverages ( @xmath12 ml ) both our dft calculations and experiment show that oxygen adsorbs in subsurface sites . our dft calculations show that the first o atom chemisorbed in the oxidation process binds in a subsurface tetrahedral site of the b2 type . the importance of substrate lattice relaxation in accurate determination of the most stable sites is demonstrated . at higher coverage o adsorbed in subsurface tetrahedral sites shows a tendency to form subsurface islands which results in an increased binding energy . for the 2 ml ( 3 ml ) coverages we find some rather unanticipated surface oxide structures , consisting of two ( three ) mixed oxygen - magnesium layers on top of an almost undistorted mg(0001 ) surface . these layered oxide structures have hexagonal symmetry and can be flat or buckled . for 23 ml coverage the rocksalt structure is found to be unfavorable compared to the buckled layer structures . however , our dft results show that the rocksalt structure may become energetically competitive at an increased coverage . this work was supported in part by the swedish foundation for strategic research ( ssf ) , the swedish research council ( vr ) , the swedish foundation for international cooperation in research and higher education ( stint ) , the carl tryggers foundation , and the polish state committee for scientific research ( kbn ) , project 5 p03b 066 21 . r.f .  would like to thank p.  aebi for continuous support and help with the xpd experiments . the allocation of computer time at the unicc facility at chalmers and gteborg university is gratefully acknowledged . note that the multipole _ r_-factor @xmath25 is not a normalized quantity , therefore its absolute value depends on many parameters and only relative differences are significant for structural analysis . so far , no rigorous theory for error estimation has been developed for @xmath25 . it is our experience , however , that _ r_-factor differences of 10 % are significant and allow to safely rule out structures with higher @xmath25 values .
Abstract:

Answer: first principles density functional calculations are used to study the early oxidation stages of the mg(0001 ) surface for oxygen coverages @xmath0 monolayers . it is found that at very low coverages o is incorporated below the topmost mg layer in tetrahedral sites . at higher oxygen - load the binding in on - surface sites is increased but at one monolayer coverage the on - surface binding is still about 60 mev weaker than for subsurface sites . the subsurface octahedral sites are found to be unfavorable compared to subsurface tetrahedral sites and to on - surface sites . at higher coverages oxygen adsorbs both under the surface and up . our calculations predict island formation and clustering of incorporated and adsorbed oxygen in agreement with previous calculations . the calculated configurations are compared with the angle - scanned x - ray photoelectron diffraction experiment to determine the geometrical structure of the oxidized mg(0001 ) surface .


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: smoothing of two - dimensional event distributions is a procedure routinely used in many fields of data analysis . in practice , smoothing means the convolution @xmath0 of the measured data @xmath1 with a kernel function @xmath2 ( often also called ` filter ' or ` window function ' ) . although the raw data may be an image in the term s common meaning [ i.e. the data set can be represented as a function @xmath3 where @xmath4 is some intensity , and @xmath5 and @xmath6 are spatial coordinates ] , the two coordinates @xmath5 and @xmath6 can , in principle , describe any two - dimensional parameter space . the coordinates @xmath5 and @xmath6 are assumed to take only discrete values , i.e.  the events are binned into @xmath7 intervals . the only requirement on @xmath4 that we shall assume in all of the following is that @xmath4 is the result of a counting process in some detector , such that @xmath8 . an image , as defined above , is a two - dimensional histogram and is thus often a coarse representation of the underlying probability density distribution ( e.g.  merrit & tremblay 1994 , vio et al . . however , for certain experiments , an unbinned event distribution may not even exist  for instance if the @xmath5 and @xmath6 values correspond to discrete pha ( spectral energy ) channels . also , _ some _ binning can be desirable , for instance , in cases where the dynamic range of the data under consideration is large . if the bin size is sufficiently small , the unavoidable loss of spatial resolution introduced by binning the raw event distribution may be a small price to be paid for a data array of manageable size .    smoothing of high - resolution data is of interest whenever the signal ( defined as the number of counts per pixel above the expected background ) in the region of interest in @xmath9 space is low , i.e.  is less than or of the order of 10 , after the raw event distribution has been sorted into intervals whose size matches approximately or exceeds the instrumental resolution . it is crucial in this context that the observed count statistics are not taken at face value but are corrected for background , which may be internal , i.e.  originating from the detector ( more general : the instrumental setup ) , or external . if this correction is not applied , the observed intensity ( counts ) distribution @xmath3 may be high across the region of interest , suggesting good count statistics , even if the signal above the background that we are interested in is actually low and poorly sampled . the statistics of the observed counts alone can thus be a poor indicator of the need for image smoothing . rebinning the data set into larger , and thus fewer , intervals improves the count statistics per pixel and reduces the need for smoothing . this is also the basic idea behind smoothing with a kernel of the form @xmath10 ( circular ` top - hat ' or ` box - car ' filter of radial size @xmath11 ) , the only difference being that smoothing occurs semi - continuously ( the step size being given by the bin size of the original data ) whereas rebinning requires an additional phase information [ the offset of the boundaries of the first bin with respect to some point of reference such as the origin of the @xmath7 coordinate system ] . however , when starting from an image binned at about the instrumental resolution , both rebinning and conventional smoothing share a well known drawback , namely that any improvement in the count statistics occurs at the expense of spatial resolution . although conventional smoothing algorithms usually employ more sophisticated functional forms for the kernel than the above ` top - hat ' filter ( the most popular probably being a gaussian ) , the problem remains that a kernel of fixed size is ill - suited for images that feature real structure on various scales , some of which may be much smaller or much larger than the kernel size . in such a situation , small - scale features tend to get over - smoothed while large - scale structure remains under - smoothed . _ adaptive - kernel smoothing _ ( aks ) is the generic term for an approach developed to overcome this intrinsic limitation by allowing the kernel to vary over the image and adopt a position - dependent ` natural ' size . aks is closely related to the problem of finding the optimal adaptive kernel estimator of the probability density distribution underlying a measured , unbinned event distribution . the advantages of adaptive kernel estimators for the analysis of discrete , and in particular one - dimensional , astronomical data have been discussed by various authors ( e.g.  thompson 1990 ; pisani 1993 , 1996 ; merritt and tremblay 1994 ; vio et al . an overview of adaptive filtering techniques in two dimensions is given by lorenz et al . ( 1993 ) . a common feature of all non - parametric adaptive kernel algorithms is that the ` natural ' smoothing scale for any given position is determined from the number of counts accumulated in its immediate environment . following the aforementioned principle , smoothing occurs over a large scale where few counts have been recorded , and over a small scale where count statistics are good . aks algorithms differ , however , in the prescription that defines how the amplitude of the local signal is to be translated into a smoothing scale . a criterion widely used for discrete data is that of silverman ( 1986 ) . it determines the size , @xmath11 , of the local kernels relative to that of some global ( i.e.  non - adaptive , fixed ) kernel ( @xmath12 ) by introducing a scaling factor which is the inverse square root of the ratio of the globally smoothed data to their logarithmic mean . for images , and using the same notation as before , this means @xmath13 where @xmath14 , and @xmath15 represents the convolution of the measured data with a kernel of fixed size @xmath12 . however , whether or not this approach yields satisfactory results depends strongly on the choice of the global smoothing scale @xmath12 ( vio et al . 1994 ) . in the context of discrete data sets , pisani ( 1993 ) suggested a least - squares cross - validation procedure to determine an optimal global kernel size in an iterative loop . however , for binned data covering a large dynamical range ( see section  [ example ] for an example ) , the dependence of the result on the size of the global kernel becomes very sensitive indeed , and the iteration becomes very time - consuming . also the dependence on the somewhat arbitrary scaling law ( eq .  [ sman ] ) remains . other adaptive filtering techniques discussed recently in the literature include the hfilter algorithm for square images ( richter et al . 1991 , see also lorenz et al . 1993 ) and the akis algorithm of huang & sarazin ( 1996 ) . closely related are image decomposition techniques including wavelet - based algorithms ( starck & pierre 1998 , and references therein ) and adaptive binning ( e.g. , sanders & fabian 2001 , cappellari & copin 2003 , diehl & statler 2005 ) .    in the following , we present asmooth , an aks algorithm for images , i.e.  binned , two - dimensional datasets of any size , which determines the local smoothing scale from the requirement that the _ above the background - corrected _ signal - ro - noise ratio ( s.n.r . ) of any signal enclosed by the kernel must exceed a certain , preset value . the algorithm is similar to akis ( huang & sarazin 1996 ) in that it employs a signal - to - noise ratio ( s.n.r . ) criterion to determine the smoothing scale . however , other than akis , asmooth does not require any initial fixed - kernel smoothing but determines the size of the adaptive kernel directly and unambiguously from the unsmoothed input data . asmooth also goes beyond existing aks algorithms in that its s.n.r . criterion takes the background ( instrumental or other ) of the raw image into account . this leads to significantly improved noise suppression in the case of large - scale features embedded in high background . our approach yields smoothed images which feature a near - constant ( or , alternatively , minimal ) signal - to - noise ratio above the background in all pixels containing a sufficient number of counts . in contrast to most other algorithms which require threshold values to be set ( e.g. , for the h coefficients in the case of the hfilter technique ) , asmooth is intrinsically non - parametric . the only external parameters that need to be specified are the minimal and , optionally , maximal signal - to - noise ratios ( above the background ) required under the kernel . the simplicity of the determination of the local smoothing scale from the counts under the kernel and an estimate of the local background greatly facilitates the translation of the smoothing prescription into a simple and robust computer algorithm , and also allows a straightforward interpretation of the resulting smoothed image . asmooth adjusts the smoothing scale ( i.e.  the size of the smoothing kernel ) such that , at every position in the image , the resulting smoothed data values share the same signal - to - noise ratio with respect to the background ; one may call this the ` uniform significance ' approach . the only external parameter required by asmooth is the desired minimal s.n.r . , @xmath16 .    in order to ensure that statistically significant structure is not over - smoothed to an s.n.r . level much higher than @xmath16 , an s.n.r . range can be specified as a pair of @xmath16 , @xmath17 values . note , however , that the maximal s.n.r . criterion is a soft one and , also , is applied only at scales larger than the instrumental resolution ( which is assumed to be similar to , or larger than , the pixel scale ) ; under no circumstances will asmooth blur significant pointlike features ( pixels whose s.n.r .  in the unsmoothed image exceeds @xmath16 ) in order to bring their s.n.r .  value down below the @xmath17 threshold . the background - corrected s.n.r .  of any features is computed according to one of two definitions . the weaker requirement is given by the definition of the `` significance of detection '' above the background , @xmath18 where @xmath19 and @xmath20 are the number of counts under the smoothing kernel and the background kernel , respectively , and @xmath21 is the @xmath22 uncertainty of the background counts . alternatively , and by default , the more stringent definition of the `` significance of the source strength measurement '' can be used , @xmath23 with @xmath24 being the @xmath22 uncertainty of the source counts accumulated under the smoothing kernel . either of these definitions assume gaussian statistics by implying that the @xmath25 error of a measurement is equal to @xmath26 times the @xmath22 error . in addition to the desired s.n.r .  limits @xmath27 , estimates of the background @xmath28 and the associated background error @xmath29 are optional additional parameter . to allow background variations across the image to be taken into account , @xmath28 and @xmath29 can be supplied as images of the same dimensions as the raw image ; in the case of a flat background @xmath28 and @xmath30 reduce to global estimates of the background and background error per pixel , i.e. , single numbers . note that , more often than not , @xmath31 as the background estimate will originate from model predictions rather than being the result of another counting experiment . if no background information is supplied , asmooth determines a local background from an annular region around the adaptive smoothing kernel , extending from 3 - 4@xmath11 for a gaussian kernel , and from 1 - 4/3@xmath11 for a top - hat kernel . internally , the threshold s.n.r . values @xmath16 , @xmath17 are translated into a minimal and a maximal integral number of counts , @xmath32 , @xmath33 , to be covered by the kernel . more precisely , the criterion is that @xmath34 where @xmath35 is the characteristic , position - dependent scale of the respective kernel . @xmath36 in eq . [ mc_crit ] are determined from the definition of the minimal and maximal s.n.r . value @xmath37 , @xmath38 where , in analogy to the definition of @xmath36 ( cf . eqs .  [ convol],[mc_crit ] ) , @xmath39 and @xmath40 are the integral number of background counts under the respective kernel and the associated error . from eq . [ taudef ] follows @xmath41 for an adaptive circular top - hat kernel of size @xmath35 ( cf . eq .  [ tophat ] ) , eq . [ mc_crit ] translates into @xmath42 , and the interpretation is straightforward : at least @xmath32 , but no more than @xmath43 , counts are required to lie within the area @xmath44 that the smoothing occurs over . in the case of a uniform background , the value of @xmath45 in eq .  [ n_min ] is simply given by @xmath46 where @xmath47 is the global background level per pixel in the input image .    for any given pair of @xmath48 values , a gaussian kernel @xmath49 will yield considerably larger effective smoothing scales than a top - hat , as , in two dimensions , more than 60 per cent of the integral weight fall outside a @xmath50 radius , whereas , in the case of a circular top - hat kernel , all of @xmath32 needs to be accumulated within a @xmath51 radius . ( note that , according to eq . [ gaussian ] , it is the weights _ per unit area _ that follow a gaussian distribution . the weights per radial annulus _ do not _ , which is why , for the kernel defined in eq . [ gaussian ] , the fraction of the integral weight that falls outside the @xmath51 radius is much larger than the 32 per cent found for a one - dimensional gaussian . ) which kernel to use is up to the user : asmooth offers a choice of gaussian ( default ) and circular top - hat . the algorithm is coded such that the adaptively smoothed image is accumulated in discrete steps as the smoothing scale increases gradually , i.e. @xmath52 where @xmath53 starts from an initial value @xmath54 which is matched to the intrinsic resolution of the raw image ( i.e. , the pixel size ) , and @xmath55 is given by @xmath56 the adaptively smoothed image is thus accumulated in a `` top - down '' fashion with respect to the observed intensities as asmooth starts at small kernel sizes to smooth the vicinity of the brightest pixels , and then increases the kernel size until , eventually , only background pixels contribute . note that condition [ maskcrit ] ensures that pixels found to contain sufficient signal at a scale @xmath53 will not contribute to the image layers @xmath57 subsequently produced with smoothing scales @xmath58 . consequently , each feature is smoothed at the smallest scale at which it reaches the required background - corrected s.n.r . ( see eq . [ mc_crit ] ) , and low - s.n.r . regions are smoothed at an appropriately large scale even in the immediate vicinity of image areas with very high s.n.r .    in order to take full advantage of the resolution of the unbinned image , the size @xmath54 of the smallest kernel is chosen such that the area enclosed by @xmath59 is about one pixel . for the circular top - hat filter of eq . [ tophat ] this means @xmath60 ; for the gaussian kernel of eq . [ gaussian ] we have @xmath61 . subsequent values of @xmath62 are determined from the requirement that eq . [ mc_crit ] be true . if a near - constant s.n.r . value is aimed at with high accuracy , i.e. , if a @xmath17 value very close to @xmath16 is chosen , the smoothing scale @xmath53 will grow in very small increments , and the smoothing will proceed only slowly . in all our applications we found values of @xmath63 to yield a good compromise between cpu time considerations and a sufficiently constant signal - to - noise ratio of the smoothed image . if no value for the optional input parameter @xmath16 is supplied by the user , the code therefore assumes a default value of @xmath64 which meets the above requirement for all reasonable values of @xmath16 .    while the intrinsic resolution of the raw image ( i.e. the pixel size ) determines the smallest kernel size @xmath54 , the size of the image as a whole represents an upper limit to the size of the kernel . although the convolution can be carried out until the numerical array representing the kernel is as large as the image itself , this process becomes very cpu time intensive as @xmath53 increases . once the smoothing scale has exceeded that of the largest structure in the image , the criterion of eq . [ mc_crit ] can never be met as only background pixels contribute . since the only features left unsmoothed at this stage are insignificant background fluctuations , the algorithm then smoothes the remaining pixels with the largest possible kernel . unavoidably , the signal - to - noise ratio of these last background pixels to be smoothed does not meet the condition of eq . [ mc_crit ] . note that in the case of the background being determined locally from the data themselves ( the default ) , it is the size of the background kernel ( an annulus around the smoothing kernel ) that reaches the limit first . hence , for a gaussian , the largest smoothing scale ( 1@xmath11 size of the kernel ) is 1/8 of the image size , for a top - hat it is 3/8 of the image size . since the algorithm uses fast - fourier transformation ( fft ) to compute the many convolutions required , the overall performance of the code is significantly improved if the image size in pixels is an integer power of two . the smoothed image obtained from the above procedure strictly conserves total counts ( within the limitations set by the computational accuracy ) and provides a fair representation of the original data at all positions . we first demonstrate the performance of our algorithm by applying it to an image of x - ray emission from a massive cluster of galaxies taken with the acis - s detector on board the chandra x - ray observatory . then we use simulated data to test how faithfully asmooth reproduces the true count distribution of the input image . because of the large range of scales at which features are detected in the selected chandra observation , this x - ray image is ideally suited for a demonstration of the advantages of aks . if photon noise is to be suppressed efficiently , the very extended emission from the gaseous intracluster medium needs to be smoothed at a rather large scale . at the same time , a small smoothing scale , or no smoothing at all , is required in high - s.n.r . regions in order to retain the spatial resolution in the vicinity of bright point - like sources ( stars , qsos , agn ) superimposed on the diffuse cluster emission . fixed - kernel smoothing can meet , at most , one of these requirements at a time . our choice of data set has the additional advantage that the selected image was taken with an x - ray detector that is also very sensitive to charged particles which makes the image particularly well suited to emphasize the importance of a proper treatment of the background . +   +   +    fig . [ ms1054 ] ( top left panel ) shows the raw counts detected with acis - s in the 0.5 to 7 kev energy band in an 89 ks observation of the galaxy cluster ms1054.4@xmath650321 . the diffuse emission originates from an electron - ion plasma trapped in the gravitational potential well of the cluster and heated to temperatures of typically @xmath66 k ( corresponding to @xmath67 kev ) . a detailed analysis of this observation is given by jeltema et al . the image shown here ( @xmath68 pixel@xmath69 ) covers a subsection of the detector spanning a @xmath70 arcmin@xmath69 region ; the pixel size of 0.492 arcsec corresponds roughly to one on - axis resolution element of the telescope - detector configuration . note the high background throughout the image , caused by high - energy cosmic rays . [ ms1054 ] summarizes asmooth results obtained with a gaussian kernel , and for s.n.r . target values of @xmath71 and @xmath72 , in a three - by - three array of plots below the image of the raw data . the left - most column of fig .  [ ms1054 ] shows the adaptively smoothed images for the three different @xmath16 values . asmooth fully preserves the high information content of the raw data in the high - s.n.r . regions corresponding to bright , small - scale features , while at the same time heavily smoothing the low - s.n.r . regions of the image where the signal approaches the background value . note , however , that for small values of @xmath16 ( top row ) the goodness of the local estimation of the signal above the background is relatively low and noise is not removed efficiently on all scales . figure  [ plotlog ] illustrates the gradual assembly of the adaptively smoothed image by showing the change in various asmooth parameters as a function of smoothing step .    at the highest pixel intensities above the s.n.r . threshold , asmooth occasionally returns smoothed count values that _ exceed _ the counts in the unsmoothed . since the final image is accumulated from the partial images resulting from the successive convolution with the whole set of differently sized kernels ( cf . .  [ aks_acc ] ) , the total smoothed intensity can become larger than the actually observed counts at pixel positions where @xmath73 . this artifact is caused by the limited resolution of the images and can be noticeable when a large bin size is chosen for the original data . the integral number of counts in the image is always conserved . = 0.5      as illustrated in fig . [ plotlog ] asmooth applies a wide range of smoothing scales to the input data as the algorithm attempts to fulfill the requirement given by eq . [ mc_crit ] throughout the image . in addition to the adaptively smoothed image , asmooth also returns , in an idl data structure , maps of the background - corrected s.n.r .  of the pixel values in the smoothed image and the kernel size used in the smoothing process , respectively ,    the second and third columns of the three - py - three panel in fig . [ ms1054 ] show both of these maps for the asmooth images shown in the first , left - most column of plots . the maps of asmooth kernel sizes as well as fig  [ plotlog ] demonstrate how very small kernel sizes of less than or about one pixel ( @xmath22 radius ) are assigned to very few bright pixels ; accordingly , these pixels remain essentially unsmoothed . for @xmath74 ( middle row ) , for instance , some 27 per cent of the image pixels are found to satisfy the criterion of eq . [ mc_crit ] at smoothing scales between one and 62 pixels ( radius ) , while the majority of the remaining pixels ( more than 72 per cent ) do not contain enough signal to reach the required s.n.r .  even at the largest permissible smoothing scale . the majority of the image is smoothed at the largest possible scale of 63.8 pixels at which the dimensions of the background kernel array ( a 1@xmath11 wide annulus surrounding the two - dimensional gaussian smoothing kernel computed out to @xmath75 radius ) equals that of the image itself . the s.n.r . maps ( fig . [ ms1054 ] , right - most column ) , finally , give evidence of the near - constant s.n.r .  of all regions smoothed with kernels other than the largest one . = 0.53    the qualitative demonstration of the performance of asmooth shown in fig . [ ms1054 ] can be made more quantitative by comparing it with the results obtained with fixed - kernel smoothing . to this end we examine the properties of the residual images obtained by subtracting the respective smoothed image from the observed raw image . following pia & puetter ( 1992 ) who introduced this criterion in the context of bayesian image reconstruction , we state that , for an ideal smoothing algorithm , this residual image should contain only uncorrelated poissonian noise around a zero mean . a _ global _ mean of zero is guaranteed  within the numerical accuracy of the convolution code  by the requirement that any smoothing algorithm conserve counts . the requirements that the mean also be zero _ locally _ , that the residual signal have poissonian variance , and that the residual signal be uncorrelated across the image are much harder to meet . in the following we shall examine how well adaptive and fixed kernel smoothing fulfill these requirements . _ is the residual signal consistent with poissonian noise shifted to zero mean as expected if shot noise dominates ? _ if so , the square of the noise ( given by the residual signal ) , @xmath76 , should equal the mean ( given by the smoothed signal @xmath77 ) . since the smoothed image contains a wide range of mean values , we can test for this condition only within bins of similar mean . figure  [ asm_csm2 ] shows the ratio @xmath78 as a function of @xmath77 for both asmooth and a number of fixed kernels of various sizes . note how the residual signal obtained with asmooth exhibits a near - poissonian variance over a larger range of mean values than does the residual produced by smoothing with a fixed kernel of essentially any size . only the smallest fixed kernel size tested here , @xmath79 pixels , yields comparable results  a kernel of such small size , however , also provides essentially no smoothing . = 0.53    _ is the signal in the residual image spatially uncorrelated ? _ figure  [ asm_csm1 ] shows the autocorrelation function of the residual images obtained with asmooth and fixed kernels of various sizes . here we use the standard definition of the autocorrelation @xmath80 as a function of radial lag @xmath81 : @xmath82 where the angular brackets signify averaging over the position @xmath83 within the image ( e.g.  peebles 1980 ) . in the normalized representation of @xmath84 shown in fig . [ asm_csm1 ] an ideal smoothing algorithm would produce zero signal at all lags except for a singular value of unity at @xmath85 . asmooth comes close to this ideal : the signal in the residual images generated with asmooth is essentially uncorrelated at all scales except for a weak ( @xmath86 ) positive signal at scales smaller than about the maximum adaptive kernel size . note the very different result when fixed gaussian kernels of size @xmath12 are used : strong spatial correlations are observed in the residual images at all scales smaller than about the kernel size . only for very small kernel sizes does non - adaptive smoothing come close to meeting the requirement that any residual signal after smoothing be spatially uncorrelated . however , such small kernel sizes perform very poorly in the presence of significant structure on a large range of scales ( cf . [ ms1054 ] , center column ) . the example shown in the previous section has the advantage of using real data but , for this very reason , does not allow the user to assess quantitatively how the asmoothed image compares to the true counts distribution underlying the noisy input image . we therefore present in this section results obtained for a simulated data set that contains multiple extended and point - like features , as well as a constant background component . figure  [ simplot ] summarizes the characteristics of the model used for our simulation ( top row ) . the same figure shows , in the bottom row , results obtained with asmooth for @xmath87 , and with the adaptive _ binning _ algorithm wvt ( diehl & statler 2005 ) for a target s.n.r . value of 5 . note how the requirement of a minimal s.n.r . _ above the local background _ .  4 and 5 ) causes asmooth to reduce noise much more aggressively than wvt in spite of a nominally lower s.n.r .  target value . as demonstrated by the error distribution depicted in the final panel of fig . [ simplot ] the asmooth residual image exhibits a near gaussian error distribution around zero mean , in addition to featuring poissonian variance and being spatially uncorrelated ( section  [ ms1054 ] ) . asmoothed images can thus be considered to be fair representations of the ` true ' input image . we describe asmooth , an efficient algorithm for adaptive kernel smoothing of two - dimensional image data . asmooth determines the local size of the kernel from the requirement that , at each position within the image , the signal - to - noise ratio of the counts under the kernel , and above the background , must reach ( but not exceed greatly ) a certain preset minimum . qualitatively , this could be called the ` uniform significance approach ' . as a consequence of this criterion , noise is heavily suppressed and real structure is enhanced without loss of spatial resolution . due to the choice of boundary conditions , the algorithm preserves the total number of counts in the raw data . the asmoothed image is a fair representation of the input data in the sense that the residual image is consistent with pure noise , i.e. , it the residual possesses poissonian variance and is spatially uncorrelated . as demonstrated by the results obtained for simulated input , the adaptively smoothed images created by asmooth are fair representations of the true counts distribution underlying the noisy input data .    since the background is accounted for in the assessment of any structure s s.n.r . , a feature that distinguishes asmooth from most other adaptive smoothing or adaptive binning algorithms , the interpretation of the adaptively smoothed image is straightforward : all features in the asmoothed image are equally significant at the scale of the kernel used in the smoothing . a map of these smoothing scales is returned by asmooth together with the smoothed image . note that these are the _ smallest _ scales at which real features reach the selected s.n.r . consequently , background regions never meet the s.n.r . criterion and are smoothed at the largest possible scale for which the kernel size equals the size of the image . such regions of insufficient signal are easily flagged using an s.n.r . map which is also provided by our algorithm . we emphasize that our definition of `` signal '' as meaning `` signal _ above the ( local ) background _ '' was chosen because many , if not most , real - life appplications do not provide the user with a priori knowledge of the intensity of the background or any spatial background variations . also , superpositions of features ( such as the point sources on top of diffuse emission in our chandra acis - s example ) require a _ local _ background estimate if the intrinsic signal of features on very different scales is to be assessed reliably . the approach taken by asmooth is thus intrinsically very different from the one implemented in , for instance , the adaptive binning algorithm wvt ( diehl & statler 2005 ) and the difference in the resulting output is accordingly dramatic ( see section  [ simsec ] ) . asmooth is being used extensively in the analysis of astronomical x - ray imaging data gathered in a wide range of missions from rosat to xmm - newton , and has become the analysis tool of choice for the display of high - resolution images obtained with the chandra x - ray observatory . recent examples illustrating asmooth s performance can be found in ebeling , mendes de oliveira & white ( 1995 ) , brandt , halpern & iwasawa ( 1996 ) , hamana et al . ( 1997 ) , ebeling et al . ( 2000 ) , fabbiano , zezas & murray ( 2001 ) , krishnamurthi et al . ( 2001 ) , karovska et al . ( 2002 ) , bauer et al . ( 2002 ) , gil - merino & schindler ( 2003 ) , heike et al . ( 2003 ) , rasmussen , stevens & ponman ( 2004 ) , clarke , blanton & sarazin ( 2004 ) , ebeling , barrett & donovan ( 2004 ) , or pratt & arnaud ( 2005 ) , as well as in many chandra press releases ( ) .    under development is an improved version of the algorithm which accounts for poisson statistics using the analytic approximations of ebeling ( 2003 ) to allow proper treatment of significant negative features , such as absorbed regions , detector chip gaps , or instrument elements ( e.g. , the window support structure of the rosat pspc ) obscuring part of the image .    asmooth is written in the idl programming language . the source code is available on request from ebeling@ifa.hawaii.edu . a @xmath88 version of an early version of the code , called csmooth , is part of ciao , the official suite of data analysis tools for the chandra x - ray observatory , and is not identical to the algorithm described here . we thank all asmooth users in the community who took the time to provide comments and criticism that have resulted in a dramatic improvement of the algorithm and its implementation since its original inception in 1995 . discussions with barry labonte about the statistical properties of residual shot noise . thanks also to our referee , herv bourdin , for constructive criticism and many useful suggestions that helped to improve this paper significantly . he gratefully acknowledges partial financial support from many sources over the years during which asmooth was developed and improved , among them a european union eara fellowship , sao contract sv4 - 64008 , as well as nasa grants nag 5 - 8253 and nag 5 - 9238 . daw acknowledges support from pparc . 9 bauer f.e . 2002 , aj , 123 , 1163 brandt w.n . , halpern j.p .  & iwasawa k. 1996 , mnras , 281 , 687 ebeling h. , mendes de oliveira c.  & white d.a . 1995 , mnras , 277 , 1006 cappellari m.  & copin y. 2003 , mnras , 342 , 345 clarke t.e . , blanton e.l .  & sarazin c.l . 2004 , 616 , 178 diehl s.  & statler t.s . 2005 , mnras , submitted ebeling h.  et al . 2000 , apj , 534 , 133 ebeling h.  2003 , mnras , 340 , 1269 ebeling h. , barrett e.  & donovan , d. 2004 , apj , 609 , l49 fabbiano g. , zezas a.  & murray s.s . 2001 , apj , 554 , 1035 gil - merino r.  & schindler s. 2003 , a&a , 408 , 51 hamana t. , hattori m. , ebeling h. , henry j.p . , futamase t. , shioya y. 1997 , apj , 484 , 574 heike k. , awaki h. , misao y. , hayashida k. , weaver k.a . 2003 apj , 591 , 99l huang z.  & sarazin c. 1996 , apj , 461 , 622 jeltema t.e . , canizares c.r . , bautz m.w . , malm m.r . , donahue m. , garmire g.p . 2001 , apj , 562 , 124 karovska m. , fabbiano g. , nicastro f. , elvis m. , kraft r.p . , murray s.s . 2002 , apj , 2002 , 577 , 114 krishnamurthi a. , reynolds c.s . , linsky j.l . , martin e. , gagn m. 2001 , aj , 121 , 337 lorenz h. , richter g.m . , capaccioli m. , longo g. 1993 , a&a , 277 , 321 merritt d.  & tremblay b. 1994 , aj , 108 , 514 peebles p.j.e . 1980 , _ the large - scale structure of the universe _ , princeton , princeton university press pia r.k .  & puetter r.c . 1992 , pasp , 104 , 1096 pisani a. 1993 , mnras , 265 , 706 pisani a. 1996 , mnras , 278 , 697 pratt g.w . & arnaud m. 2005 , a&a , 429 , 791 rasmussen j. , stevens i.r .  & ponman t.j . 2004 , mnras , 354 , 259 richter g.m . , bohm p. , lorenz h. , priebe a. 1991 , astron . nachr . , 312 , 346 sanders j.s .  & fabian a.c . 2001 , mnras , 325 , 178 silverman b.w . 1986 , density estimation for statistics and data analysis , chapman & hall , london starck j .- l .  & pierre m. 1998 , a&as , 128 , 397 thompson a.m. 1990 , a&a , 240 , 209 vio r. , fasano , g. , lazzarin m. , lessi , o. 1994 , a&a , 289 , 640
Abstract:

Answer: an efficient algorithm for adaptive kernel smoothing ( aks ) of two - dimensional imaging data has been developed and implemented using the interactive data language ( idl ) . the functional form of the kernel can be varied ( top - hat , gaussian etc . ) to allow different weighting of the event counts registered within the smoothing region . for each individual pixel the algorithm increases the smoothing scale until the signal - to - noise ratio ( s.n.r . ) within the kernel reaches a preset value . thus , noise is suppressed very efficiently , while at the same time real structure , i.e.  signal that is locally significant at the selected s.n.r .  level , is preserved on all scales . in particular , extended features in noise - dominated regions are visually enhanced . the asmooth algorithm differs from other aks routines in that it allows a quantitative assessment of the goodness of the local signal estimation by producing adaptively smoothed images in which all pixel values share the same signal - to - noise ratio _ above the background_.    we apply asmooth to both real observational data ( an x - ray image of clusters of galaxies obtained with the chandra x - ray observatory ) and to a simulated data set . we find the asmoothed images to be fair representations of the input data in the sense that the residuals are consistent with pure noise , i.e.  they possess poissonian variance and a near - gaussian distribution around a mean of zero , and are spatially uncorrelated . techniques : image processing  methods : data analysis  methods : statistical