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: according to the instanton liquid model ( ilm ) instantons are largely responsible for the spontaneous breaking of chiral symmetry and the low energy properties of light hadrons @xcite . this has been demonstrated using phenomenological instanton liquid models . the two main tasks involved in the construction are    * producing an equilibrium instanton ensemble assuming some interaction between the pseudoparticles . this is typically done with a variational calculation @xcite . * computing hadronic observables using the above instanton ensemble as the gauge background . this is a formidable task and the computation necessarily involves several approximations . it would thus be very desirable to check the physical picture emerging from the ilm , starting from first principles .    in the last few years there has been a considerable effort in this direction , using lattice simulations of qcd as a tool . the properties of the instanton ensemble have been extracted from monte carlo generated lattice configurations using several different `` smoothing '' techniques . the most important parameters , the topological susceptibility , the instanton density , and the average instanton size were found to be in qualitative agreement with the ilm . on the other hand  although there is a considerable amount of indirect evidence @xcite  so far the lattice has not given any direct results concerning the second part of the above programme .    in ref . @xcite cooled lattice configurations were used to show that low eigenmodes of the dirac operator saturate the point - to - point hadronic correlators . more recently a strong correlation was found between the chiral density of low modes of the staggered dirac operator and the topological charge density @xcite . these results suggest that instantons probably really play an important role in qcd . all this work however has involved monte carlo generated or cooled lattices that contained other gauge field fluctuations in addition to the instantons . therefore , based solely on these studies it is impossible to separate the part of the effects due to the instantons and the part ( if any ) that is due to the rest of the gauge field fluctuations .    to my knowledge , the only lattice study so far that makes such a separation possible has appeared in ref . @xcite . there , the instanton content of smoothed gauge configurations was reconstructed . hadron spectroscopy done on these artificial instanton backgrounds was compared to the spectroscopy on the corresponding smoothed configurations . while the smoothed ensemble yielded physically reasonable results for the pion and the rho mass , the instantons failed to reproduce this . in fact wilson spectroscopy on the instanton ensemble was hardly different from free field theory ( trivial gauge field background ) . in view of the positive results cited above , this was highly unexpected , especially that at first sight the only difference between the compared two ensembles appeared to be the presence versus the lack of confinement on the smoothed and the artificial instanton ensembles respectively . according to the ilm , confinement is not essential for the low energy properties of light hadrons @xcite . the findings of ref . @xcite indeed suggest that the physics in the chiral limit is rather insensitive to the string tension . with less than half of the original string tension , the low fermion modes on the cooled configurations still produced physically sensible hadron correlators .    in the present paper i address the question of how to reconcile the insensitivity of light hadrons to the string tension ( also expected from the ilm ) with the failure of the instantons alone to reproduce the correct physics ( found in @xcite ) . i will exclusively consider the su(2 ) gauge group . along the way i will also obtain some useful insight into the structure of the low end of the spectrum of lattice dirac operators and in particular i will show the importance of using a ( close to ) chiral lattice dirac operator . this is essential for a faithful reproduction of continuum - like quark ( quasi)zero modes . the plan of the paper is as follows . in section [ se : zmz ] i briefly discuss the instanton liquid picture of physics in the chiral limit , and in particular the role of the zero modes and the so called `` zero mode zone '' . in section [ se : cs ] i point out the role of the explicit chiral symmetry breaking of the lattice dirac operator and the importance of minimising it . i show that the optimised clover dirac operator introduced in @xcite performs much better in this respect than the wilson operator . i explicitly study the zero mode zone in the simplest nontrivial case , that of an instanton antiinstanton pair and find that the mixing of zero modes closely follows the pattern expected in the continuum .    in section [ se : spectrum ] i find that the optimised clover operator produces spectroscopy results markedly different from the wilson operator but it still falls short of reproducing the correct physics on the instanton backgrounds . instead , it yields a peculiar mixture of qcd and free field theory . i demonstrate that this is due to the presence of a large number of low lying free field modes in the spectrum of the dirac operator with the instanton backgrounds that produce artificially light hadrons . the free modes are completely absent from the smoothed configurations . they can be also eliminated from the instanton configurations by superimposing an ensemble of zero momentum gauge field fluctuations ( torons ) over the instantons . these gauge field fluctuations survive even deep in the perturbative regime . i argue that they are probably needed to compensate for the anomalously large small - volume effects present in the instanton ensemble due to the lack of confinement . by comparing the density of eigenvalues around zero on the physical ( smoothed ) , the toron , and toron plus instanton ensembles , i obtain direct evidence that the instantons are responsible for chiral symmetry breaking . after the elimination of the free quark modes , the instanton backgrounds are shown to reproduce the correct physical pion and rho correlators . section [ se : con ] contains my conclusions along with some speculations on still unanswered questions . finally , in the appendix , for reference i collected some general properties of wilson type lattice dirac operators that i used in the main text . some of these can be also found in the literature . in this section i describe how instantons can produce low lying modes of the dirac operator and why this is important for the physics in the chiral limit . for simplicity i start with the continuum theory and in the second part of this section i discuss the complications that appear on the lattice . let @xmath0 be the continuum massless dirac operator which implicitly depends on a gauge field background . the basic building block of the physics in the chiral limit is the @xmath1 limit of the massive quark propagator , @xmath2 . from the spectral decomposition in terms of eigenmodes of @xmath0 @xmath3 it is easily seen that in the chiral limit the most important modes are the ones with eigenvalues @xmath4 close to 0 . the crucial assumption of the ilm is that it is the instantons that are responsible for generating the bulk of these lowest eigenmodes . this can be qualitatively understood as follows . according to the atiyah - singer index theorem , in the background of an instanton ( antiinstanton ) @xmath0 has at least one negative ( positive ) chirality zero mode . in the presence of an infinitely separated instanton antiinstanton pair we still expect to find opposite chirality zero modes localised on the two objects . if the members of the pair are brought closer to one another , the two degenerate zero eigenvalues will in general split into two complex ones still close to the origin . a very simplistic but qualitatively correct description of this can be given as follows . let @xmath5 and @xmath6 be the zero modes of the dirac operators @xmath7 and @xmath8 respectively . since we consider @xmath0 in different gauge backgrounds , we explicitly indicate this ; @xmath9 and @xmath10 refer to an instanton and an antiinstanton . we now want to describe the spectrum of @xmath11 where @xmath12 is some superposition of the gauge fields @xmath9 and @xmath10 . since @xmath5 and @xmath6 are of opposite chirality , @xmath13 . let us arbitrarily complete the set of these two vectors into an orthonormal basis and construct the matrix of @xmath11 in this basis . the elements of this matrix in the subspace spanned by @xmath5 and @xmath6 are @xmath14 where @xmath15 the diagonal matrix elements vanish because @xmath0 maps left handed vectors into right handed ones and vice versa . by a suitable choice of the phases , the off - diagonal elements can always be made real and due to the anti - hermiticity of @xmath0 they are of equal magnitude and have opposite signs . if we now assume that @xmath16 and @xmath17 lie approximately in this two - dimensional subspace , i.e. the matrix of @xmath11 contains the above @xmath18 block - diagonal part , then @xmath11 can be easily seen to have two complex conjugate eigenvalues @xmath19 with the corresponding ( unnormalised ) eigenvectors being @xmath20 . this mechanism can be generalised to gauge fields which are superpositions of several instantons and antiinstantons . the basis in this case consists of the zero modes of the individual ( anti)instantons . the only additional approximation involved is the assumption that the zero modes corresponding to different ( anti)instantons are orthogonal . if the pseudoparticles are well separated , this is approximately true . one can thus consider the matrix of @xmath0 restricted to the subspace of zero modes  called the `` zero mode zone ''  in a general instanton background . if the off diagonal elements of this matrix are small then it will have complex eigenvalues close to the origin , in addition to a number of exact zero eigenvalues corresponding to the total topological charge . moreover , the complex eigenmodes will be approximately linear combinations of several instanton and antiinstanton zero modes and thus become highly delocalised , making it possible to propagate quarks to large distances .    the most important quantities in this construction are the off diagonal instanton antiinstanton matrix elements , the @xmath21 s . in the simplest case of only one pair , @xmath21 is known to depend on the distance and the relative orientation of the members of the pair as @xmath22 where @xmath23 describes the relative orientation in group space , @xmath24 and @xmath25 is the relative position of the instanton and the antiinstanton @xcite . @xmath26 is a function of the instanton and antiinstanton scale parameters and the distance between them . at small distances it depends on the ansatz used to combine the field of the instanton and the antiinstanton , at large separation it is expected to fall off as @xmath27 . we have seen how instanton zero modes can give rise to the low - lying eigenmodes that are the most important ones in the chiral limit . close to the chiral limit the number of modes giving a substantial contribution to the quark propagator becomes smaller and smaller . it is thus crucial to reproduce these low eigenmodes and eigenvalues as faithfully as possible on the lattice if we want to test the ilm . the most important obstruction to this is that due to the nielsen - ninomiya theorem it is impossible to construct a local lattice dirac operator describing one fermion species with exact chiral symmetry @xcite . in the present paper i shall limit the discussion to wilson type lattice dirac operators that describe one fermion species but explicitly break chiral symmetry . in the absence of chiral symmetry the fermion zero modes corresponding to the topology of the gauge field are not protected ; any small fluctuation of the gauge field can shift the zero eigenvalues .    in order to understand how this happens , a useful starting point is the following property shared by the lattice and the continuum dirac operator . let @xmath28 be an eigenvector of @xmath0 corresponding to the eigenvalue @xmath29 . if @xmath0 has no degeneracies then @xmath30 if and only if @xmath29 is real . this is a simple consequence of the so - called `` @xmath31 hermiticity '' , @xmath32 , a common property of the continuum and the lattice dirac operator . a proof of the above statement is given in the appendix .    in the continuum , @xmath0 is anti - hermitian and has eigenvalues on the imaginary axis . thus the only real eigenvalue it can have is 0 . chirality makes a sharp distinction between zero and non - zero modes in the continuum . zero modes have a chirality of @xmath33 while all the non - zero eigenmodes have zero chirality . since the eigenmodes and their chiralities depend continuously on the gauge field , zero modes are protected , they can not be shifted away from zero by continuous deformations of the gauge field . on the other hand , the lattice dirac operator is not anti - hermitian , and its eigenvalues are not constrained to be on the imaginary axis . while lattice zero modes are still protected from being shifted off the real axis by smooth deformations of the gauge field ( this would cause their chirality to jump to zero ) , there is nothing preventing them from moving continuously along the real axis . indeed , a lattice dirac operator will in general not have any exact zero modes unless the gauge field is fine tuned . it is now clear that the lattice analogues of the continuum zero modes are real modes of the dirac operator . these have non - vanishing chiralities but in general their magnitude is smaller than 1 . the crucial role chiral symmetry plays in this discussion can be understood by noting that chiral symmetry of @xmath0 ( i.e. @xmath34 ) together with @xmath31 hermiticity would imply that @xmath0 is anti - hermitian and has protected zero modes . the smaller the explicit chiral symmetry breaking of @xmath0 is , the more its low lying spectrum will resemble that of an anti - hermitian operator . recently it has been shown that the closest a lattice dirac operator can come to being chirally symmetric is by obeying the ginsparg - wilson relation and thus having an ultralocal chirality breaking in the quark propagator @xcite . in this case the spectrum of @xmath0 lies on a circle in the complex plane passing through the origin . indeed the low end of this spectrum is almost on the imaginary axis and there are protected zero modes . i shall now discuss how the continuum description of the zero mode zone presented in the previous subsection has to be modified on the lattice . assume that there is an instanton on the lattice . in the limit when the instanton size goes to infinity ( in units of the lattice spacing ) , the instanton `` does not see '' the lattice at all , and the continuum description applies . the lattice dirac operator , @xmath7 has an exact zero eigenvalue with a corresponding @xmath35 chirality eigenmode @xmath5 . if the instanton is not infinitely large , the real eigenvalue gets shifted away from zero , its magnitude depends on the instanton size and the chirality breaking of @xmath7 @xcite . the corresponding eigenmode also becomes only an approximate eigenvector of @xmath31 and thus @xmath36 . let us now consider the dirac operator in the presence of an instanton antiinstanton pair . following the continuum discussion in the previous subsection , its matrix elements between the ( lattice approximate ) `` zero modes '' @xmath5 and @xmath6 can be parametrised as @xmath37 where @xmath38 is real , @xmath39 and @xmath40 are non - negative real numbers . in the su(2 ) case the diagonal elements are automatically real and the off - diagonal ones can be also made real by a suitable choice of the phases of the eigenvectors . this will be proved in the appendix . due to the @xmath31 hermiticity of @xmath0 , since @xmath5 and @xmath6 are approximate eigenvectors of @xmath31 , the off - diagonal terms have opposite signs and about the same magnitude . the eigenvalues of this matrix are @xmath41 with the corresponding ( unnormalised ) eigenvectors being @xmath42 the most important feature of these eigenvectors is that  as in the continuum  they are mixtures of the instanton and antiinstanton zero modes with roughly the same magnitude ( the ratio being @xmath43 ) and with a relative phase @xmath44 , up to a correction proportional to @xmath38 . this mechanism produces highly delocalised eigenmodes in the zero mode zone and thus facilitates the propagation of quarks by jumping from instanton to antiinstanton . maintaining this feature on the lattice is thus crucial for the description of the zero mode zone . the requirement for the proper mixing of zero modes is the inequality @xmath45 for producing qualitatively continuum - like eigenmodes , the difference of the diagonal elements has to be much smaller than the magnitude of the off - diagonal mixing matrix elements of @xmath11 . otherwise the two eigenvalues will be real and the corresponding eigenmodes rapidly become localised on the instanton and the antiinstanton respectively . recall that in the continuum , the mixing matrix elements @xmath46 are proportional to @xmath47 , where @xmath48 is the invariant angle of the instanton antiinstanton relative orientation ( see eq . ( [ eq : relor ] ) ) . @xmath38 ( and @xmath49 ) is expected to be of the order of the real eigenvalues of @xmath7 and @xmath8 , both small but nonzero numbers . this means that inequality ( [ eq : ineq ] ) will be violated for @xmath48 sufficiently close to @xmath50 . we can minimise the range in relative orientation , where this happens by making @xmath38 small , i.e. having the would - be zero modes close to zero . in this section i study how the general features of the zero mode zone described in the previous section are realised by two commonly used lattice dirac operators , the wilson and the clover operator . the gauge background i use for this test contains an instanton of size @xmath51 and an antiinstanton of size @xmath52 , separated by a distance @xmath53 in the time direction . the lattice size is @xmath54 and i always use antiperiodic boundary conditions in the time direction , periodic in all other directions . these parameters are typical in the spectroscopy calculations that i shall present in the next section . let me start with the wilson operator . in the presence of only the ( anti)instanton with the above parameters , @xmath55 has a real eigenvalue at ( 0.17 ) 0.51 with chirality ( 0.39 ) -0.53 . it is clear that for these small instantons the wilson operator does not even come close to satisfying the conditions of the previous section : the `` zero modes '' are far away from zero with chiralities of magnitude much smaller than 1 . the simplified description of the instanton - antiinstanton configuration in terms of the @xmath18 matrix is then inadequate . if the instanton and antiinstanton are close to being parallel @xmath56 then the real part of the ( anti)instanton related eigenvalues is around 0.3 . moreover already at @xmath57 the two complex eigenvalues become real and the corresponding modes are rather localised at the instanton and the antiinstanton respectively . a much better description of the zero mode zone is provided by the optimised clover action proposed in ref . ( see also @xcite for more discussion and a recent scaling test of this action . ) the clover coefficient was chosen to minimise the range in which the real eigenvalues occur in the physical branch of the spectrum on locally smooth gauge backgrounds . as explained in the previous section , this is exactly what is needed for a good description of the mixing of zero modes . the optimal value for @xmath58 was found to be around 1.2 on a set of ape smeared @xmath59 quenched su(3 ) configurations . physically these are quite similar to the su(2 ) configurations appearing in the present study , therefore in what follows i shall always choose @xmath60 . the real modes occurring in the single ( anti)instanton background already reveal that the clover action is much closer to the continuum than the wilson action . the eigenvalues are located at -0.014 and -0.002 and the chiralities of the corresponding eigenmodes are 0.993 and -0.995 respectively . in the instanton plus antiinstanton background the lowest two modes are always complex for almost any relative orientation of the pair , even at @xmath61 . we can hope that the simple description in terms of the mixing matrix works well in this case . indeed , the magnitude of the diagonal elements is small , neither @xmath38 nor @xmath29 goes above 0.02 in magnitude for any relative orientation . in fig . [ fig : t12 ] i plotted the two off - diagonal matrix elements , @xmath62 versus @xmath48 along with the one - parameter fits of the form const.@xmath63 . as expected , the mixing matrix elements are very well described by the continuum ansatz of eq.([eq : relor ] ) . a comparison of the eigenvalues of the @xmath18 mixing matrix and the explicitly computed eigenvalues of the lattice dirac operator in the corresponding gauge backgrounds shows agreement to within @xmath64 in the whole range of @xmath65 . we can conclude that the @xmath60 clover operator should give a satisfactory treatment of the zero mode zone as long as the instanton size does not drop below around @xmath51 . this is of course true only if the gauge configuration is locally smooth . otherwise the fermion mass acquires an additive renormalisation which completely destroys the above features important for the description of the zero mode zone . having a fermion action that is expected to describe the zero modes and their mixing reasonably well , we can discuss the main topic of the paper , namely , what part of qcd is reproduced by instantons alone . in this section i present quenched hadron spectroscopy results obtained with the @xmath60 clover fermion action . the gauge field backgrounds on which the spectroscopy is done are the ones used in ref . the starting point is a set of 28 @xmath54 su(2 ) configurations generated with a fixed point action . the lattice spacing is @xmath66fm , as fixed by the sommer parameter of the heavy quark potential . these monte carlo generated configurations are not yet suitable for our purposes since they are not locally `` smooth '' . both the instantons and the details of the zero mode zone are obscured by short distance fluctuations . cycling , a smoothing technique based on the renormalisation group , has been shown to preserve the main physical features of the gauge configurations . after 9 cycling steps the string tension @xcite and @xmath67 versus @xmath68 are essentially unchanged . moreover , these lattices are smooth enough that their instanton content can be unambiguously identified . in fact , about @xmath69 of the action is accounted for by instantons ( assuming there is no interaction between them ) . we can create artificial lattice configurations that have the same instanton content as the smoothed ones . the instanton sizes and locations are reproduced but the relative orientation in group space is distributed according to the su(2 ) haar measure . recently a study of the relative orientation appeared in @xcite , however reproducing the relative orientations would make our procedure much more cumbersome . i compare the spectroscopy done on the smoothed ( cycled ) lattices and the corresponding artificial instanton configurations with the same instanton content . these ensembles are both locally smooth enough that the optimised clover action should give a satisfactory description of their zero mode zone . in fig . [ fig : pirho ] the vector meson mass is plotted as a function of the pseudoscalar mass for both ensembles . the correlators were always fitted with the assumption that there is only one lightest particle dominating both the pseudoscalar and the vector channels . i also show the @xmath70 free field line . in fact , the wilson fermion action produces exactly degenerate free field like pions and rhos on the instanton backgrounds . compared to the wilson action , the clover action shows a marked improvement , but it still fails to reproduce the @xmath67 vs.  @xmath68 observed on the smoothed lattices . [ fig : pi ] shows the pion mass sqared vs.  the bare quark mass for both ensembles . on the smoothed ensemble @xmath71 , as expected from pcac . this is clearly not the case on the instanton ensemble , where a best fit to the form @xmath72 gives @xmath73 which is between @xmath74 and the free field value , @xmath75 . we also note that since both ensembles are locally very smooth , the additive mass renormalisations are very close to zero . therefore , it is also meaningful to compare hadron masses obtained at the same bare quark mass on the two gauge ensembles . this comparison reveals that the instanton ensemble typically produces much lighter hadrons than the smoothed ensemble .      in spite of its substantially improved chiral properties , even the clover dirac operator fails to reproduce the correct physics of qcd on the instanton backgrounds . instead , it yields a peculiar mixture of qcd and free field theory . is there a problem with the instanton liquid model or the physics of these instanton lattices is still not close enough to the continuum ? to answer this question we look at the low lying eigenmodes of the dirac operator in more detail .    in fig . [ fig : dspec ] i plot the lowest 30 eigenvalues in the complex plane for the instanton and the smoothed ensemble . the eigenvalues of all 28 configurations are superimposed . the boundary condition is antiperiodic in the time direction , periodic in all other directions . in both cases most of the eigenvalues lie very close to a circle passing through the origin . this is a sign of the approximate chiral symmetry of the action . the most striking difference between the two ensembles is the appearance of a gap at around im@xmath29=0.5 on the instanton ensemble . it is also instructive to compare the density of eigenvalues . since the low eigenvalues are almost imaginary i project them on the imaginary axis and plot the densities as a function of the imaginary part of the eigenvalues . this corresponds to the density of eigenvalues around zero in the continuum .    in fig . [ fig : ev_hist1 ] we can compare the densities corresponding to the two ensembles . the real modes have been removed , these would show up as spikes at the origin . the densities at zero seem to agree quite well . according to the banks - casher relation the chiral condensate is proportional to the density of modes around zero @xcite ( excluding exact zero modes , which do not contribute in the infinite volume limit ) . this shows that on both the instanton and the smoothed ensemble chiral symmetry is broken and the value of the chiral condensate is approximately the same . apart from the vicinity of the origin , however , the two distributions are quite different . on the instanton ensemble there is a substantial `` piling up '' of modes above im(@xmath76 and then a `` thinning '' of modes at 0.5 . what is the reason of this huge difference between the two distributions ? the answer can be easily given by looking at the eigenvectors of the dirac operator . it turns out that the quark density of the eigenvectors on the instanton ensemble corresponding to the peak of the distribution are much more delocalised than all the other eigenvectors occurring in either ensemble . in fact , these modes are essentially free quark modes . this can be seen as follows . on a trivial gauge field configuration ( all links are 1 ) with periodic bondary conditions in all directions , there are @xmath77 trivial ( constant ) zero modes , this is the number of c - number degrees of freedom corresponding to one fermion flavour . in the presence of antiperiodic boundary condition in the time direction , the lowest eigenmodes are shifted away from zero , to @xmath78 since in our case @xmath79 . both eigenvalues are @xmath77-fold degenerate . the mixing of the free field modes into a particular eigenmode @xmath28 can be characterised by @xmath80 , where @xmath81 is the projection of the normalised eigenmode @xmath28 onto the eigenspace corresponding to the eigenvalue @xmath82 . this quantity is not gauge invariant therefore i work in landau gauge . any generic fermion mode on a non - trivial but locally smooth configuration will have a nonzero projection on the @xmath82 eigenspaces . the question is whether this is just an accidental mixing or the given ensemble really contains close to free field modes . to decide this , a good quantity to look at is @xmath83 . if the mixing is accidental , we expect @xmath84 to fluctuate around 1 , independently of the corresponding eigenvalue @xmath29 . on the other hand , if there are free field like modes on a given configuration then @xmath84 will increase substantially when @xmath29 approaches the free field eigenvalue @xmath85 .    in fig . [ fig : proj ] i plotted @xmath84 versus the distance of the given eigenvalue from the free field mode @xmath85 . the figure shows a random selection of eigenvalues with imaginary parts between 0.0 and 0.2 . in the smoothed ensemble @xmath84 fluctuates around 1 everywhere , there is no trace of the free field modes . on the other hand , in the instanton ensemble , the modes close to @xmath85 have a substantially larger projection on the @xmath85 eigenspace than on the @xmath86 subspace . in fact , these modes , close to @xmath85 have @xmath87 , so they are essentially free field modes . a detailed study of the chiral density @xmath88 reveals that all the modes with @xmath89 look like mixtures of instanton zero modes with the density concentrated in several lumps . this is to be contrasted with the chiral density of the modes with @xmath90 that spreads almost homogenously over the whole lattice , as expected of the lowest free field eigenmodes . now we can understand why the instanton ensemble gave substantially smaller hadron masses than the smoothed ensemble . the reason is that the former had a large number of free quark modes . even if the zero mode zones of the two ensembles are similar , the free field modes provide a very efficient way to propagate quarks to large distances . this is why the instanton ensemble shows a peculiar mixture of qcd and free field characteristics . in the smoothed ensemble confinement completely eliminates the freely propagating quarks .      according to the instanton liquid model , instantons alone , without confinement , can reproduce most of the properties of the light hadrons . our result shows that in the presence of instantons only , without confinement , free quark modes provide a very effective way to propagate quarks and they contaminate the hadron correlators producing unusually light hadrons .    one possible way to suppress the free quark modes is to add to the instantons some locally very smooth background that contains only long wavelength fluctuations . there has to be enough long distance structure to suppress the free propagation of quarks but on the other hand they have to be locally smooth enough not to distort the instantons . cycling is a very efficient way to create such backgrounds . to this end i started with an ensemble of wilson @xmath91 ( @xmath92fm ) @xmath93 su(2 ) gauge configurations , performed 4 cycling steps and finally inverse blocked them to size @xmath54 . this produced a set of lattices extremely smooth on the few lattice spacing scale but containing the longest wavelength fluctuations characteristic to a @xmath94@xmath95 lattice . their spatial size was a bit smaller , their temporal size a bit larger than the confinement scale . moreover these configurations had a small enough physical size that they almost never contained instantons . after superimposing these lattices on the instanton configurations by simply multiplying the corresponding links ( in landau gauge ) , i checked that the original topological charge of the instanton configurations changed very little .    to show how the addition of the smooth background affects the hadron correlators , in fig . [ fig : pirho_b2.4 ] i plotted typical pseudoscalar and vector correlators ( at @xmath96 ) . the change is dramatic . while on the instanton configurations the pion and the rho are very light ( crosses ) , the addition of the smooth background brings the correlators ( bursts ) very close to the physical ones obtained on the smoothed lattices ( boxes ) . the good agreement of the physical correlators and the instanton plus smooth background correlators persists in the whole range of masses where i tested it ( @xmath97 ) .    it would be tempting to interpret this as a consequence of confinement ; the superimposed smooth backgrounds still contained the longest scale fluctuations of the confining lattices and these removed the free field modes that contaminated the spectroscopy before . but is confinement really needed for that ? this can be easily tested by replacing the superimposed backgrounds with ones of much smaller physical size . the lattices that i used for this purpose were generated in exactly the same way as before , except that the starting @xmath98 lattices were produced at wilson @xmath99 . these lattices have a tiny physical size , they are expected to be completely perturbative . in fig . [ fig : pirho_b3.0 ] i plotted again the pion and rho correlators at @xmath96 . quite surprisingly , the correlators have absolutely no dependence on the @xmath100 at which the background was created , they are the same on the instantons plus confining and instantons plus perturbative configurations .    apparently , there is still something nontrivial on these physically very tiny perturbative lattices that has the same effect on the correlators as a confining background . what could this be ? the answer can be given by noting that the @xmath101 limit of yang - mills theory on the torus is not completely trivial because there is not only one gauge field configuration with zero action . the set of flat ( minimal action ) gauge field configurations can be characterised by the su(2 ) conjugacy classes represented by four mutually commuting gauge group elements corresponding to the ( constant ) polyakov loops in the four perpendicular directions . these are essentially zero momentum gauge field configurations , with a constant gauge potential ( link ) , sometimes referred to as torons . their contribution to physical quantities is a finite size effect , and has been estimated in ref . @xcite .    to explicitly check that it is really the torons that affect the correlators so strongly , i extracted the `` toron content '' of the perturbative background fields . this was done in the following way . after landau gauge fixing i averaged the polyakov loops in all four directions . since the polyakov loops along any given direction were almost constant ( up to small perturbative fluctuations ) , the averages were close to being su(2 ) elements . i projected the averages back onto su(2 ) and then distributed the average polyakov loop evenly among the links in the given direction . this resulted in configurations with constant gauge fields ( links in any direction ) that carried the average polyakov loops of the original perturbative configurations . i shall refer to these lattices as `` toron '' configurations . we can now compare the correlators in the instantons plus toron background and the instantons plus the full perturbative background . the result  again in the typical case of @xmath96  is shown in fig . [ fig : toron ] . the instantons and the torons together fully reproduce the correlators , demonstrating that the most important fluctuations on the perturbative lattices are indeed the torons . it is also quite instructive to compare the distribution of the lowest eigenvalues of the dirac operator on the different types of gauge configurations . [ fig : bc_hist ] shows the densities on the smoothed lattices ( solid lines ) , the perturbative ones ( dotted line ) and the lattices containing both the instantons and the perturbative backgrounds ( dashed line ) . the three distributions agree quite well away from zero . however around zero , the perturbative ensemble produces only a very low density ( compatible with zero ) , while the two other ensembles yield qualitatively similar densities , both nonzero . this demonstrates very clearly that it is the instantons that are responsible for creating a nonzero eigenvalue density around zero which through the banks - casher relation @xcite implies the spontaneous breaking of chiral symmetry . we also note that the peak corresponding to the free field modes ( see fig .  [ fig : ev_hist1 ] ) has been completely eliminated by the perturbative backgrounds . this again can be attributed to the torons , as can be checked by looking at the eigenvalue distribution on the toron lattices ( not shown ) which is very similar to the dotted line in fig . [ fig : bc_hist ] . in the present paper i addressed the question of how much of qcd in the chiral limit is reproduced by restricting the fluctuations of the gauge field to instantons only . i reconstructed the instanton content of smoothed monte carlo lattice configurations and compared hadron spectroscopy on this instanton ensemble to the spectroscopy on the original `` physical '' smoothed configurations . the fermion action used for the spectroscopy was a `` chirally optimised '' clover action . i explicitly studied the fermion zero modes in the presence of an instanton and also the mixing of the zero modes corresponding to an instanton and an antiinstanton . i found that in these simple cases the optimised action described the continuum features of the zero mode zone rather well . in spite of this , the optimised action still failed to reproduce the physical hadron spectrum . it yielded anomalously light hadrons with a mixture of qcd and free - field like features .    a closer look at the spectrum of the dirac operator revealed that on the instanton ensemble it had a large number of free quark modes . these were absent on the physical smoothed configurations . they provided a very efficient way of propagating quarks to large distances thereby substantially reduced the hadron masses and contaminated the spectroscopy . superimposing an ensemble of very smooth , essentially perturbative gauge field fluctuations on the instantons , was enough to eliminate the free quark modes and to restore the physical hadron correlators . by comparing the density of low eigenvalues of the dirac operator i also obtained direct evidence that it is the instantons that are responsible for creating a nonzero density of modes around zero and by the banks - casher relation also for chiral symmetry breaking . the important fluctuations on the superimposed perturbative lattices , needed to recover the physical hadron correlators , turned out to be zero momentum ( constant link or gauge potential ) configurations , torons . i would like to emphasise here that this toron ensemble was a perturbative one , produced at large @xmath100 i.e. small physical lattice size . the meson correlators appeared to be largely independent of the @xmath100 at which the superimposed background was created , as long as it was above the finite temperature phase transition on the given lattice size . it seems rather surprising that such `` mild '' gauge field fluctuations as the torons , can have so profound an impact on physical quantities . after all , torons exist only on the torus , they are absent in an infinite space - time and thus constitute only finite size effects . on the other hand we know that the instantons by themselves do not confine @xcite . this means that the instanton ensemble does not produce a mass gap and it is expected to exhibit much stronger finite size effects than qcd , from which the instantons have been extracted . in fact , the free field modes that are eliminated by the torons are also absent in the infinite volume limit , since they are not normalisable . in free field theory the density of modes around zero is proportional to the linear size of the box @xmath102 , whereas the density of instanton related modes grows proportionally to the volume . it is thus not inconceivable that in increasingly bigger volumes , although the torons have less and less influence , at the same time they become less and less needed since the relative importance of the free field modes also dies out . the torons might be needed only for compensating the unnaturally large finite size effects due to the absence of confinement in the instanton ensemble .    what is the picture emerging from this lattice study concerning the instanton liquid model ? the qcd vacuum contains an ensemble of instantons that is generated by the non - perturbative gauge field dynamics . the instantons by themselves break chiral symmetry but hadron correlators in the instanton backgrounds are strongly contaminated by freely propagating quarks . this yields anomalously light mesons and a small splitting between the pseudoscalar and the vector channel . the free quarks can be eliminated and the physical hadron correlators can be restored by superimposing a perturbative ensemble of zero momentum gauge configurations ( torons ) on the instantons . this suggests that the contamination by free quarks is most likely a finite volume effect which , in qcd , is completely suppressed by confinement .    more evidence in favour of this scenario could be gathered by studying the volume dependence of these effects on larger volumes . another important question , not answered by the present study is why it is exactly the perturbative toron ensemble that is needed to restore the physical meson correlators . it would be also very desirable to extend this investigation to a larger class of observables , in particular baryon correlators , for which the su(3 ) case need be considered . finally , for more precise tests a better chiral action , such as the overlap @xcite could be used . in the appendix , for reference , i collected some properties of wilson type lattice dirac operators that were used in the main body of the paper . some of these results can be found in the literature , in particular in refs . @xcite . i assume that the dirac operator satisfies @xmath31 hermiticity , i.e. @xmath103    let @xmath0 be a @xmath31-hermitian dirac operator with no degeneracy in its spectrum , i.e.  having as many different eigenvalues as the dimension of the space it acts on . let @xmath29 , and @xmath49 be two eigenvalues of @xmath0 with the corresponding eigenvectors being @xmath28 and @xmath104 . then @xmath105 if and only if @xmath106 . _ proof : _ a simple consequence of @xmath31 hermiticity is that @xmath107 which implies @xmath108 it immediately follows that if @xmath109 then @xmath106 . to prove the converse we first note that since @xmath0 has as many different eigenvalues as the dimension of the space it acts on , its eigenvectors form a basis . we have already seen that the only eigenvector on which @xmath110 can have a nonzero projection , is @xmath111 . if this projection were also zero , then @xmath110 would be orthogonal to all the vectors in a complete set and consequently it would be zero . this is impossible since @xmath31 has a trivial kernel .    for proving that @xmath112 it was essential that the eigenvalues of @xmath0 be non - degenerate and its eigenvectors form a complete set . generally this is expected to be the case unless the gauge field is fine tuned . there is however an important exception . when topology  as seen by the fermions  changes through a smooth deformation of the gauge fields , two opposite chirality real eigenvalues collide and leave the real axis as a complex conjugate pair . it can be shown that when the two eigenvalues coincide , the corresponding eigenspace is one - dimensional , the eigenvector has zero chirality , @xmath0 does not possess a complete set of eigenvectors and therefore has no spectral decomposition . a simple consequence of statement 1 is that if @xmath0 has a complete non - degenerate spectrum , then the chirality of an eigenvector , @xmath88 is nonzero if and only if the corresponding eigenvalue , @xmath29 , is real .    in the remainder of the appendix i shall prove the following property of wilson type su(2 ) lattice dirac operators . let @xmath113 and @xmath114 be real eigenvectors of two ( not necessarily the same ) su(2 ) lattice dirac operators such that the corresponding eigenvalues and their complex conjugates are non - degenerate . it is then possible to choose the phases of @xmath113 and @xmath114 such @xmath115 , @xmath116 , @xmath117 , and @xmath118 are all real for any ( third ) lattice dirac operator @xmath0 . _ proof : _ besides @xmath31 hermiticity , this property depends on an additional symmetry of the dirac operator relating @xmath0 to its complex conjugate @xcite . this is specific to the su(2 ) case . let @xmath119 be the charge conjugation operator and @xmath120 , where @xmath121 is a pauli matrix acting in colour space . using that @xmath122 ( where the last equality follows because we work in a representation with all @xmath123 s hermitian ) , and that for any @xmath124 @xmath125 , it is not hard to prove that @xmath126 now if @xmath28 is an eigenvector of @xmath0 with eigenvalue @xmath29 then @xmath127 which means that @xmath128 is an eigenvector of @xmath0 with eigenvalue @xmath129 . due to the non - degeneracy of the eigenvalues and that @xmath130 , it follows that @xmath131 in the special case when @xmath29 is real , the phase choice @xmath132 eliminates the extra phase factor on the r.h.s . ( [ eq : phase ] ) .    now writing @xmath133 , and making use of eqs . ( [ eq : dstar ] ) and ( [ eq : phase ] ) with the above phase choice , the non - diagonal matrix elements can be easily seen to be real . the same argument also shows that the diagonal matrix elements are automatically real , regardless of the phase choice . this completes the proof . we would like to note that while any wilson type dirac operator has a spectrum symmetric with respect to the real axis ( this follows from @xmath31 hermiticity ) , in general there is no simple relation between the eigenvectors corresponding to complex conjugate eigenvalues . the only exception is the su(2 ) case , when they are related by eq . ( [ eq : phase ] ) . i thank pierre van baal for discussions , tom degrand and anna hasenfratz for correspondence . the computations were performed using a computer code based on the milc @xcite code and the arpack @xcite package . i thank the milc collaboration and the authors of the arpack package for making their code publicly available . this work has been supported by fom . de  forcrand , m. garca prez , and i .- o .  stamatescu , nucl . phys .  * b499 * ( 1997 ) 409 ; c.  michael , and p.s . spencer , phys . * d52 * ( 1995 ) 4691 ; d.a . smith , and m.j . teper , phys . * d58 * ( 1998 ) 014505 ; b.  alles , m.  campostrini , a.  di  giacomo , y.  gunduc , e.  vicari , phys . * d48 * ( 1993 ) 2284 ; a.  hasenfratz , and c.  nieter , phys . * b439 * ( 1998 ) 366 . edwards , u.m . heller , j.  kiskis , and r.  narayanan , phys . * 82 * ( 1999 ) 4188 ; r.g . edwards , u.m . heller , j.  kiskis , and r.  narayanan , chiral condensate in the deconfined phase of quenched gauge theories , hep - lat/9919941 .
Abstract:

Answer: i address the question of how much of qcd in the chiral limit is reproduced by instantons . after reconstructing the instanton content of smoothed monte carlo lattice configurations , i compare hadron spectroscopy on this instanton ensemble to the spectroscopy on the original `` physical '' smoothed configurations using a chirally optimised clover fermion action . by studying the zero mode zone in simple instances i find that the optimised action gives a satisfactory description of it . through the banks - casher formula , instantons by themselves are shown to break chiral symmetry but hadron correlators on the instanton backgrounds are strongly influenced by free quark propagation . this results in unnaturally light hadrons and a small splitting between the vector and the pseudoscalar meson channels . superimposing a perturbative ensemble of zero momentum gauge field fluctuations ( torons ) on the instantons is found to be enough to eliminate the free quarks and restore the physical hadron correlators . i argue that the torons that are present only in finite volumes , are probably needed to compensate the unnaturally large finite size effects due to the lack of confinement in the instanton ensemble . = 1ex plus0.5ex minus0.2ex    * instantons and chiral symmetry on the lattice * + * tams g. kovcs * + _ instituut - lorentz for theoretical physics , p.o.box 9506 + 2300 ra , leiden , the netherlands _ + e - mail : kovacs@lorentz.leidenuniv.nl +


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: stochastic geometry has recently been used for the analysis and performance evaluation of wireless  ( ad hoc as well as cellular ) networks ; in this approach , one models node locations as a spatial point process , e.g. , homogeneous poisson point processes , and one computes various network statistics , e.g. , interference , successful transmission probability , coverage  ( or , outage ) probability etc . as spatial averages . this often leads to tractable performance metrics that are amenable to parametric optimization with respect to network parameters ( node density , protocol parameters , etc . ) . more precisely , this approach yields spatial averages of the performance metrics for given network parameters ; then the parameters can be chosen to optimize performance . this approach takes a macroscopic view of the network with the underlying assumption that all nodes in the network have identical statistical characteristics .    in practice , due to randomness and heterogeneity in networks , nodes need to adapt to local spatial and temporal conditions  ( e.g. , channel conditions and topology ) to reach optimum network wide performance . for example , nodes in wireless lans adjust their window sizes based on acknowledgment feedback ; in cellular networks nodes are scheduled based on channel conditions and adapt their transmit powers based on the measured sinrs , which in turn depend on the transmit powers set by other nodes . in all such scenarios , distributed adaptive algorithms are used to reach a desired network wide operating point e.g. that maximizing some utility . while the behavior of such distributed optimization protocols is often well understood on a given topology , there are usually no analytical characterizations of the statistical properties of the optimal state in large random and heterogeneous networks . the main aim of this work is to use stochastic geometry to study spatial adaptations of medium access control in aloha that aim at optimizing certain utilities . while we identify a utility for which stochastic geometry can be used to compute the spatial distribution of map and the expected utility , we are far from being able to do so for all types of utilities within the @xmath0-fair class and we discuss the difficulties to be faced . let us start with a review of the state of the art on aloha . wireless spectrum is well known to be a precious and scarce shared resource . medium access control  ( mac ) algorithms are employed to coordinate access to the shared wireless medium . an efficient mac protocol should ensure high system throughput , and should also distribute the available bandwidth fairly among the competing nodes . the simplest of the mac protocols , aloha and slotted aloha , with a `` random access '' spirit , were introduced by abramson  @xcite and roberts  @xcite respectively . in these protocols , only one node could successfully transmit at a time . reference  @xcite modeled node locations as spatial point processes , and also modeled channel fadings , interferences and sinr based reception . this allowed for spatial reuse and multiple simultaneous successful transmissions depending on sinr levels at the corresponding receivers . all the above protocols prescribe identical attempt probabilities for all the nodes . reference  @xcite further proposed opportunistic aloha in which nodes transmission attempts are modulated by their channel conditions .    among the initial attempts of map adaptation in aloha , reference @xcite analyzed protocol model and proposed stochastic approximation based strategies that were based on receiver feedback and were aimed at stabilizing the network . references  @xcite also optimized nodes attempt probabilities  ( or thresholds ) in order to maximize the spatial density of successful transmissions . reference  @xcite analyzed both plain and opportunistic aloha in a network where all the nodes communicate to one access point . they assumed statistically identical rayleigh faded channels with no dependence on geometry  ( i.e. , no path loss components ) . they demonstrated a paradoxical behavior where plain aloha yields better aggregate throughput than the opportunistic one . reference  @xcite also studied optimal random access with sinr based reception . however , they considered constant channel gains . they developed a centralized algorithm that maximizes the network throughput , and also an algorithm that leads to max - min fair operation . reference  @xcite modeled network as an undirected graph and studied aloha under the protocol model . they designed distributed algorithms that are either proportional fair or max - min fair . reference  @xcite built upon the model of  @xcite , and formulated the channel access problem as a non - cooperative game among users . they considered throughput and delay as performance metrics and proposed pricing schemes that induce socially optimum behavior at equilibrium . however , they set time average quantities  ( e.g. , throughput , delay ) as utilities  ( or costs ) , and concentrated on symmetric nash equilibria . consequently , in their analysis , dependence on local conditions vanishes .    in none of the above aloha models , nodes account for both wireless channel randomness and local topology for making their random access decisions , as we do in the present paper . there is a vast literature on the modeling of csma by stochastic geometry which will not be reviewed in detail here . the very nature of this mac protocol is adaptive as each node senses the network and acts in order to ensure that certain exclusion rules are satisfied , namely that neighboring nodes do not access the channel simultaneously . however , csma as such is designed to guarantee a reasonable scheduling , not to optimize any utility of the throughput . the closest reference to our work is probably  @xcite where the authors study an adaptation of the exclusion range and of the transmit power of a csma node to the location of the closest interferer . this adaptation aims at maximizing the mean number of nodes transmitting per unit time and space  ( while respecting the above exclusion rules ) . this mean number is however only a surrogate of the rate . in addition , the adaptation is only w.r.t . to the location of the nearest interferer . we study spatial adaptation of aloha in ad hoc networks . the network setting is described in section  [ sec : network - model ] . we consider quasi - static networks in which mobiles learn the topology , and incorporate this information in their medium access probability  ( map ) selection . we consider the cases where nodes are benevolent and cooperate in a distributed way to maximize the global network throughput or to reach either a proportional fair or max - min fair sharing of the network resources . we analyze the case where nodes account only for their closest interferers , for all nodes in a given ball around them , or even all nodes in the network . section  [ sec : aggr - interf ] is focused on the distributed algorithms that maximize the aggregate throughput or lead to max - min fairness in such networks . in the proportional fair case , we show that nodes can compute the optimal maps as solutions to certain fixed point equations . in the maximum throughput case , the optimal maps are obtained through a gibbs sampling based algorithm . in the max min case , the optimal maps are obtained as the solution of a convex optimization problem . section  [ ss : stocg ] contains the stochastic geometry results . the model features nodes forming a realization of a homogeneous poisson point process in the euclidean plane . we compute the map distribution in such a network in the proportional fair case using shot noise field theory . to the best of our knowledge , this distribution is the first example of successful combination of stochastic geometry and adaptive protocol design aimed at optimizing ceratin utility function within this aloha setting . we also show that the mean value of the logarithm of the throughput obtained by a typical node can be derived from this distribution . finally , we discuss the difficulties to be faced in order to extend the result to other types of utilities . the numerical results are gathered in section  [ ss : numerics ] . the aim of this section is two - fold : 1 ) check the analytical results against simulation and 2 ) quantify the gains brought by adaption within this setting . we model the ad - hoc wireless network as a set of transmitters and their corresponding receivers , all located in the euclidean plane . this is often referred to as `` bipole model ''  ( * ? ? ? * chapter  16 ) . there are @xmath1 transmitter - receiver pairs communicating over a shared channel . the transmitters follow the slotted version of the aloha medium access control  ( mac ) protocol  ( see figure  [ fig : bipole - model ] ) . a transmitter , in each transmission attempt , sends one packet which occupies one slot . each transmitter uses unit transmission power . we assume that each node has an infinite backlog of packets to transmit to its receiver . the euclidean distance between transmitter @xmath2 and receiver @xmath3 is @xmath4 , and the path - loss exponent is @xmath0  ( @xmath5 ) . we also assume rayleigh faded channels with @xmath6 being the random fading between transmitter @xmath2 and receiver @xmath3 . moreover , we assume that the random variables @xmath7 are independent and identically distributed with mean @xmath8 . thus all @xmath6s have cumulative distribution function  ( cdf ) @xmath9 with @xmath10 . all the receivers are also subjected to white gaussian thermal noise with variance @xmath11 , which is also constant across slots . we assume that a receiver successfully receives the packet of the corresponding transmitter if the received sinr exceeds a threshold @xmath12 . let @xmath13 be the indicator variable indicating whether transmitter @xmath3 transmits in a slot , and @xmath14 be @xmath3 s medium access probability . thus @xmath15 when node @xmath3 transmits , the received sinr at the corresponding receiver is @xmath16 then the probability of successful reception @xmath17 can be calculated as follows . @xmath18 thus @xmath19 where @xmath20 further simplifying , @xmath21 then , the rate or throughput of transmitter @xmath3 is given by @xmath22 . the thermal noise appears merely as a constant multiplicative factor in the expression for the successful transmission probability  ( see  ) . moreover , in interference limited networks , the impact of thermal noise is negligible as compared to interference . we focus on such networks , and thus we ignore the thermal noise factor throughout . in this section , we analyze adaptations of spatial aloha that maximize aggregate throughput or achieve proportional fairness or max - min fairness .      the throughput maximizing medium access probabilities solve the following optimization problem . we first argue that the optimum in the above optimization problem is attained at one of the vertices of the hypercube formed by the constraint set . to see this , suppose @xmath24^{\mathcal{n}}$ ] is an optimal solution , and @xmath25 for some @xmath26 . clearly , @xmath27 since the partial derivative is independent of @xmath14 , @xmath14 can be set to either @xmath28 or @xmath29 without reducing the value of the objective function . this proves our claim . in the following we focus only on such extreme solutions . then the above problem is equivalent to finding an @xmath30 such that @xmath31 if and only if @xmath32 is an optimal solution . thus we are interested in @xmath33    [ [ an - iterative - solution ] ] an iterative solution + + + + + + + + + + + + + + + + + + + + +    we can pose this problem as a strategic form game with the users as players  @xcite . for each player its action @xmath34 lies in @xmath35 , and the utility function @xmath36 is given by @xmath37 this is a potential game with the above objective function as the potential function  @xcite . thus the best response dynamics converges to a nash equilibrium . this algorithm can be implemented in a distributed fashion if each node @xmath3 knows @xmath38 for all @xmath2 , and also @xmath39 and @xmath40 for all @xmath41 after each iteration . however , a nash equilibrium can be a suboptimal solution to the above optimization problem . to alleviate this problem , we propose a gibbs sampler based distributed algorithm , wherein each node @xmath3 chooses action @xmath29 with probability @xmath42 the parameter @xmath43 is called the temperature . the gibbs sampler dynamics converges to a steady state which is the gibbs distribution associated with the aggregate throughput and the temperature @xmath43  @xcite . in other words , we are led to the following distribution on the action profiles : @xmath44 where @xmath45 is a normalizing constant . when @xmath43 goes to @xmath28 in an appropriate way  ( i.e. , as @xmath46 , where @xmath47 is the time ) , the distribution @xmath48 converges to a dirac mass at the action profile @xmath49 with maximum aggregate utility if it is unique . notice that the aggregate utility @xmath50 equals the aggregate throughput . thus the action profile @xmath49 is a solution to the original throughput optimization problem . the first two terms in the utility function @xmath51 can be seen as `` selfish '' part of user @xmath3 , whereas the last summation term is `` altruistic '' part . the user makes a decision based on whether the `` selfish '' part dominates or viceversa .    in a quasi - static network where topology continuously changes , although at a slower time scale , different sets of nodes are likely to be scheduled to transmit under different topologies . thus , in terms of long term performance , maximum throughput medium access is not grossly unfair . the proportional fair medium access problem can be formulated as follows . @xmath52 the objective function can be rewritten as @xmath53 we thus have a convex separable optimization problem . the partial derivative of the objective function with respect to @xmath14 is @xmath54 which is continuous and decreasing in @xmath14 over @xmath55 $ ] . we conclude that at optimality , for each @xmath56 , @xmath57 observe that user @xmath3 s optimal attempt probability is independent of the attempt probabilities of other users . in particular , if @xmath58 , user @xmath3 can perform iterations @xmath59 autonomously . furthermore , @xmath60 clearly , @xmath61 , i.e. , @xmath62 is a contraction . thus the fixed point iterations converge to the optimal @xmath14 starting from any @xmath63 $ ] . the characterization of the optimal attempt probabilities reflects the altruistic behavior of users . more precisely , user @xmath3 s attempt probability is a function of @xmath64 which are measures of @xmath3 s interference to all other users . in particular , if @xmath65 , i.e. , if @xmath3 s transmission does not cause significant interference to the other users , then @xmath3 transmits in all the slots . unlike the throughput maximization problem , there is no `` selfish '' component in the decision making rule .    as the target sinr @xmath66 , @xmath67 for all @xmath68 , and the proportional fair attempt probabilities satisfy @xmath69 for all @xmath26 . this yields @xmath70 for all @xmath26 . this is expected , because in the limiting case a transmission can succeed if and only if there is no other concurrent transmission . this is hence aloha without spatial reuse , and it is well known that in this case , the optimal access probability is @xmath71 asymptotically  @xcite . our analysis in this section follows  @xcite . the max - min fair medium access problem can be formulated as @xmath72 where constraint functions are defined for all @xmath73^{\mathcal{n}}$ ] . the following is an equivalent convex optimization problem  ( see  @xcite for details ) : @xmath74    the lagrange function of this problem is given by  @xcite @xmath75 with @xmath76 being the lagrange multipliers . minimization of the lagrange function  ( which is concave in @xmath77 and @xmath78 ) gives @xmath79 wang and kar  @xcite suggest that the lagrange multipliers be updated using the gradient projection method . more precisely , for all @xmath26 , @xmath80^+ , \label{eqn : aggr - interf - lambda - update}\end{aligned}\ ] ] where @xmath81 is the step size at the @xmath82th iteration . further more ,  ( * ? ? ? * theorem  2 ) implies that a solution arbitrary close to an optimal solution can be reached via appropriate choice of step sizes . however , all the users need to exchange variables in order to perform updates .    finally , the _ directed link graph _ corresponding to our network is a directed graph in which each vertex stands for a user  ( i.e. , a transmitter - receiver pair ) in the network . there is an edge from vertex @xmath3 to vertex @xmath2 in the directed link graph if transmission of user @xmath3 affects the success of transmission of user @xmath2 . two vertices @xmath3 and @xmath2 are said to be connected if either of the following two conditions hold :    1 . there is an edge from @xmath3 to @xmath2 or viceversa , 2 . there are vertices @xmath83 such that @xmath84 and @xmath85 are connected for @xmath86 . clearly , the directed link graph for our network model is a complete graph ; for any pair of vertices @xmath3 and @xmath2 there is an edge from @xmath3 to @xmath2 and also from @xmath2 to @xmath3 . in particular , the directed link graph is a single _ strongly connected component _  @xcite . thus  ( * ? ? ? * corollary  1 ) implies that the above optimization also obtains the lexicographic max - min fair medium access probabilities that yield identical rates for all the users . note that a user needs to know the entire topology , and in a few cases , also needs to communicate with all the nodes to implement the adaptation rules developed in sections  [ sec : aggr - interf - max - thput]-[sec : aggr - interf - max - min - fair ] . in this section , we carry out analysis under the simplifying assumption that the aggregate interference at a receiver is dominated by the transmission from the closest interferer . this is a reasonable approximation in a moderately dense network , specifically when the path loss attenuations are high . throughout this section , we use the notation @xmath87 for all @xmath56 . in words , @xmath88 is the strongest interferer of node @xmath3 , and @xmath89 is the set of nodes to which node @xmath3 is the strongest interferer . we assume that there is always a unique @xmath88 for each @xmath3 . then , accounting only for the nearest interferer , the approximate probability of successful transmission for node @xmath3 is @xmath90 the analysis of sections  [ sec : aggr - interf ] can be adapted to this simplified scenario . the throughput maximization problem can now be posed as follows . @xmath91    as in section  [ sec : aggr - interf - max - thput ] , we can argue that some @xmath92 attains the optimal throughput . again , an equivalent optimization problem is @xmath93 or alternatively , @xmath94    we now formulate a strategic form game among users , with action sets @xmath35 and utility functions given by @xmath95 again a gibbs sampler based algorithm yields the optimal set of the transmitting users . also , user @xmath3 only needs to know the distances of user @xmath88 and all the receivers in @xmath89 and their actions to make its decision . notice that user @xmath3 must choose @xmath96 if @xmath97 such users can set their actions to @xmath29 , and need not undergo gibbs sampler updates . [ [ discussion ] ] discussion + + + + + + + + + +    consider a scenario where a node s closest interferer does not transmit , i.e. , has zero attempt probability . nonetheless , this node always has an active closest interferer  ( unless there are no other nodes in the network ) . a better approximation of the success probabilities , and hence of the throughput , is obtained by always accounting for the closest active interferer . towards this , let us define @xmath98 for all @xmath32 . we are now faced with the following optimization problem . @xmath99 we can now define users utility functions as follows . @xmath100 where @xmath101 . the analogous distributed algorithm  ( gibbs sampler based ) can again be shown to lead to the optimal solution . we now aim to solve the following optimization problem . @xmath102 following the discussion in section  [ sec : aggr - interf - prop - fair ] , we obtain @xmath103 again , if @xmath104 , iterations @xmath59 converge to the optimal @xmath14 starting from any @xmath63 $ ] . if @xmath105 , i.e. , if @xmath3 s transmission does not cause significant interference to the users for whom @xmath3 is closest interferer , then @xmath3 transmits in all the slots . the same user may not transmit  ( in any slot ) under the maximum throughput objective if @xmath106 and @xmath88 and users in @xmath89 transmit . we can have explicit formulae for the attempt probabilities in a few special cases .    1 . suppose @xmath89 is singleton for each @xmath3 . if @xmath107 , then @xmath108 2 . _ linear network topology :  _ we now consider a scenario where nodes are placed along a line , say @xmath109 , and are indexed sequentially . we also assume that for any node @xmath3 the potential interferers are the two immediate neighbors @xmath110 and @xmath111 . this also amounts to assuming @xmath112 for all @xmath3 . then , assuming @xmath113 , the proportional fair attempt probability of node @xmath3 satisfies @xmath114 this is quadratic equation in @xmath14 , which on solving gives @xmath115      similar to section  [ sec : aggr - interf - max - min - fair ] , the max - min fair medium access problem is @xmath116 again , the constraint functions are defined for all @xmath73^{\mathcal{n}}$ ] . let us recall the definition of the directed link graph associated with the network . in this section , we only account for the interference due the closest interferer . thus there is an edge from vertex @xmath3 to vertex @xmath2 if and only if @xmath117 . we assume that the directed link graph is connected  ( i.e. , all the vertices in the graph are connected with each other ) . if it is not connected , the max - min fair medium access problem on the entire graph decomposes into separate max - min fair medium access problems on each of the connected subgraphs , which can be solved independently . we now pursue the following convex optimization problem which is equivalent to the above max - min fair medium access optimization problem . @xmath118 the last set of constraints along with the connected assumption  ( of the directed link graph ) forces @xmath119 to be equal for all @xmath26 . this confirms equivalence to the initial optimization problem . now the lagrange function is @xmath120 where @xmath121 , are the lagrange multipliers . an approach similar to section  [ sec : aggr - interf - max - min - fair ] prescribes the following update rules . for all @xmath122 , @xmath123^+ , \label{eqn : close - interf - lambda - update } \\ \mu_{ij}(n+1 )   = & [ \mu_{ij}(n ) + \beta(n)(\theta_i - \theta_j)]^+ , \label{eqn : close - interf - mu - update}\end{aligned}\ ] ] where @xmath81 is the step size at the @xmath82th iteration as before . observe that any user @xmath3 can perform updates  - via local information exchange . more precisely , it needs to communicate only with user @xmath88 and the users in @xmath89 . the aim of the present paper is primarily of theoretical nature and it is beyond our scope to discuss implementation issues . let us however stress that the discussed adaptations of the map are implementable . we will focus on the proportional fair case in view the main focus of the paper . assume each receiver has a distinctive pilot signal with fixed power . since we assumed a quasi - static network in which the nodes move at a slower time scale , each node can then learn the distance that separates it from a given receiver by listening to its pilot signal and by performing a time average over  ( so as to smooth out fading ) . once this data is available for all receivers , a given transmitter can then solve the key fixed point equation that characterizes its optimal map . distinctive pilot signals can be obtained by a collection of orthogonal codes chosen at random by the receivers . in practice , it is enough for a transmitter to detect the `` dominant '' receivers  ( i.e. , those within a certain distance to it ) , so that the scheme will work in an infinite network with a finite  ( properly tuned ) number of such codes . we now assume that the transmitting nodes are scattered on the euclidian plane according to a homogeneous poisson point process of intensity @xmath124 . for each transmitter , its corresponding receiver is at distance @xmath125 in a random direction . the traffic and channel models are the same as in section  [ sec : network - model ] . as before the transmitters use slotted aloha to access the channel , and a receiver successfully receives the packet from its transmitter if the received sinr exceeds a threshold @xmath12 . finally , transmitters adapt their attempt probabilities as described in section  [ sec : aggr - interf ] . each transmitter is associated with a multi dimensional mark that carries information about the adaptive transmission probability and the transmission status . let @xmath126 denote a marked poisson point , where    * @xmath127 denotes the poisson point process of intensity @xmath124 , representing the location of transmitters in the euclidean plane . * @xmath128 denote the marks of the poisson point process @xmath129 , which consist of three components : * * @xmath130 denote the angles from transmitters to receivers . these angles are i.i.d . and uniform on @xmath131 $ ] and independent of @xmath132 . we will call them the primary marks . * * @xmath133 denote the maps of the nodes ; @xmath134 is a secondary mark ( i.e. a functionals of @xmath132 and its primary marks , see below ) . * * @xmath135 are indicator functions that take value one if a given node decides to transmit in a given time slot , and zero otherwise . clearly , @xmath136 . in particular , given @xmath134 , @xmath137 is independent of everything else including @xmath138 . the locations of the receivers will be denoted by @xmath139 with @xmath140 . it follows from the displacement theorem @xcite that @xmath141 is also a homogeneous poisson point process of intensity @xmath124 . the above assumptions will be referred to as the poisson model . we will also consider below a more general case where the above marked point process is just stationary . let us consider response functions @xmath142 defined for each @xmath143 as follows @xmath144 for all @xmath145 , the shot noise field @xmath146 associated with the above response function and the marked point process @xmath129 is @xmath147 notice that this shot noise is not that representing the interference at @xmath148 . it rather measures the effect of the presence of a transmitter at @xmath148 on the whole set of receivers . consider a typical node at the origin , @xmath149 , with marks @xmath150 . let @xmath151 denote the palm distribution of the stationary marked point process @xmath129  ( * ? ? ? * chapter  1 ) . the fixed point equation determining the map of node @xmath152 reads  ( see  ) @xmath153 we have a similar equation for each node and the sequence @xmath133 is readily seen to be a sequence of marks of @xmath132 and its primary marks . it follows directly from monotonicity arguments that @xmath154 notice that we have not used the specific assumptions on the point process so far . hence we have the following general connection between the optimal map distribution and the shot noise @xmath155 :    for all stationary marked point processes @xmath156  ( not necessarily poisson ) , for all @xmath157 , @xmath158 with @xmath151 the palm distribution of @xmath129 . we now use the fact that @xmath129 is an independently marked poisson point process  ( * ? ? ? * definition  2.1 ) . from slivnyak s theorem  ( * ? ? ? * theorem  1.13 ) , @xmath159 for all @xmath160 . consequently , @xmath161    it follows from  ( * ? ? ? * proposition  2.6 ) and from the fact that @xmath162 is a homogeneous poisson point process that one can write the laplace transform @xmath163 of the shot noise @xmath164 as @xmath165 where @xmath166 . [ thm : map - distribution ] under the above poisson assumptions , the attempt probability of the typical node has the distribution @xmath167 with @xmath168 given by ( [ eq : lapl ] ) . let @xmath169 denote the density of the shot noise field @xmath170 . then @xmath171 where @xmath172 if @xmath173 and 0 otherwise . now using parseval s theorem @xmath174 with @xmath175 the fourier transform of the real valued random variable @xmath176 and @xmath177 the complex conjugate of @xmath178 . the claim follows after substituting @xmath179 and @xmath180 . [ lem : laplace - transform ] for @xmath181 , the laplace transform @xmath182 can be simplified as @xmath183      this subsection is devoted to the analysis of the mean value of the logarithm of the throughput of the typical node : @xmath184 with @xmath185 defined in  . since we know the cdf @xmath186 of @xmath187 , the first term poses no problem . the second term can be rewritten as  ( see  ) @xmath188.\ ] ] under the law @xmath151 , the points @xmath189 of @xmath132 form a homogeneous poisson point process of intensity @xmath124 . however , the marks @xmath190 do not have the law identified in the last section . in fact , the mark @xmath134 of a point @xmath191 satisfies the following modified fixed point equation : @xmath192 with the convention that @xmath193 if there is no solution in @xmath55 $ ] . we can use the same argument as above to conclude that @xmath194 iff @xmath195    conditioned on there being two nodes at @xmath28 and @xmath148 , the other points form a homogeneous poisson point process of intensity @xmath124 . this allows one to prove the following . [ thm : map - typical - receiver ] under the above poisson assumptions , given that there is a node at 0 and a node at @xmath196 , the attempt probability of the node at @xmath148 has the distribution @xmath197 with @xmath198 and @xmath140 . due to the circular symmetry , the first integral in the expression of @xmath199 depends on @xmath148 only through @xmath200 . thus the density of @xmath201 also depends on @xmath200 only , and it will be denoted by @xmath202 when @xmath203 . the density of @xmath187 identified in the last subsection will be denoted by @xmath186 . the main result of this section is :    [ thm : mean - utility ] under the above poisson assumptions , in the proportional fair case , the mean utility of a typical node is @xmath204    see appendix  [ proof : thm - mean - utility ] . a preliminary concern when trying to use euclidean stochastic geometry of the infinite plane in the maximum throughput and max - min fairness cases is that it is not clear whether the associated infinite dimensional optimization problems make sense in the first place . in the proportional fair case , each node computes its optimizing map in one step as the solution of a fixed point equation that is almost surely well defined  ( in terms of a shot noise ) even in the infinite poisson population case . unfortunately , this does not extend to the other two cases . this does not mean that there is no hope at all . consider for instance the maximum throughput case accounting only for the closest interferer , and further simplify it by measuring interference at the transmitter rather than at the receiver . in other words , consider the same optimization problem as in section  [ sec : close - interf - max - thput ] but with @xmath88 being the closest transmitter to transmitter @xmath3  ( rather than to receiver @xmath3 ) . then , in the poisson case , the infinite dimensional optimization problem can be shown to reduce to a _ countable collection of finite optimization problems_. this follows from the fact that there are no infinite `` descending chains '' in a homogeneous poisson point process  ( * ? ? ? * chapter  2 )  ( a _ descending chain _ is a sequence of nodes @xmath205 such that @xmath206 for all @xmath82 ) . as a result , the poisson point process can be decomposed into a countable collection of finite `` descending trees '' , where each path from the leaves to the root is a descending chain . the associated optimization is hence well defined and the problem can be reduced to evaluating the solution of the optimization problem of section  [ sec : close - interf - max - thput ] on the typical descending tree . hence , there is hope to progress on this and on related cases . this will however not be pursued here and is left for future research . in this section , we study the proposed adaptive spatial aloha schemes quantitatively . we compute various metrics formulated through the stochastic geometry based analysis , and we also perform simulation . the simulation not only validates the analytical model , but also illustrates the performance of schemes for which we do not have an analytical characterization . we used maple and matlab to evaluate the integrals of section [ ss : stocg ] . the infinite integral that shows up in the expression of the laplace transform ( [ eq : lapl ] ) is handled without truncation by maple and matlab . the singularity at @xmath207 in the contour integrals ( [ eq : parsev ] ) leveraging parseval s theorem is a false singularity and it is also handled without further work by either maple or matlab . the matlab code is particularly efficient and is used throughout the analytical evaluations described below . we consider a two dimensional square plane with side length @xmath208 , and @xmath1 nodes placed independently over the plane according to the uniform distribution ; this corresponds to @xmath209 in the stochastic geometry model . , the poisson random variable with mean @xmath210 is highly concentrated around its mean . thus we can use @xmath210 nodes for all the realizations in our simulation . ] each node has its receiver randomly located on the unit circle around it , again as per the uniform distribution . thus @xmath211 for all @xmath3 . we set @xmath181 and @xmath212 . to nullify the edge effect , we take into account only the nodes falling in the @xmath213 square around the center while computing various metrics . while all other parameters remain , we vary @xmath208 and @xmath1 for different simulations . for each parameter set we calculate the average of the performance metric of interest over @xmath214 independent network realizations . we validate the analytical expression against the simulation for the case of proportional fair medium access . for illustration , we plot the cumulative distribution function  ( c.d.f . ) of the map in figure  [ fig : map - cdf - analysis ] . here we set @xmath215 and consider two values of @xmath1 , @xmath216 and @xmath217 , which correspond to @xmath218 and @xmath219 respectively . the plots show that the stochastic geometry based formula  ( see theorem  [ thm : map - distribution ] ) quite accurately predicts the nodes behavior in simulation . we also study the distributions @xmath220  ( also referred to as @xmath202 for @xmath203 ) defined in theorem  [ thm : map - typical - receiver ] . figure  [ fig : map - cdf - receiver ] shows their plots for @xmath219 and two values of @xmath200 , @xmath221 and @xmath222 . again the plots based on the analytical expression and those based on simulation closely match . clearly , under @xmath151 , nodes closer to the origin are more likely to be inactive .    , the map distribution of a node at distance @xmath223 from origin , under @xmath151.,width=288,height=240 ]      in this section we illustrate the performance of various adaptive schemes and their benefit over plain aloha . we compute the performance metrics via simulation and also through analytical expressions whenever the latter are available . in such cases the analytical results and the simulation validate each other . first we set @xmath224 and @xmath225 . we consider the maximum throughput medium access , however , only accounting for the closest interferers . in figure [ fig : gibbs - closei ] , we show steady state behavior of our gibbs sampling based algorithms ; we have set the temperature @xmath226 . as expected , the improved maximum throughput medium access  ( see the discussion at the end of section  [ sec : close - interf - max - thput ] ) insures a better exclusion behavior . under this scheme , a lesser number of nodes transmit , and neighboring nodes are unlikely to transmit simultaneously . so , this is expected to deliver better aggregate throughput . now we keep @xmath208 fixed at @xmath227 , but vary @xmath1 from @xmath228 to @xmath229 ; this corresponds to varying @xmath124 from @xmath230 to @xmath231 in the analytical expressions . we evaluate the aggregate throughputs of various aloha schemes including plain aloha . the average throughputs are plotted in figure  [ fig : throughput ] . _ although in some of the schemes we derive the attempt probabilities only considering the closest interferers , we always take into account the aggregate interference while calculating the throughput_. when the number of nodes is small , both the throughput maximizing medium access and plain aloha have identical performance ; both prescribe attempt probabilities close to one for all the nodes . when the number of nodes increases beyond @xmath232 , the throughput maximizing medium access significantly underestimates the interference , and thus its performance deteriorates . on the other hand , the aggregate interference based proportional fair scheme significantly outperforms plain aloha in terms of aggregate throughput also . this benefit is sustained even as the number of nodes increases . we also notice that the improved version of maximum throughput medium access  ( see the discussion at the end of section  [ sec : close - interf - max - thput ] ) yields best performance among all the schemes , and its performance does not deteriorate until a much higher number of nodes .        in figure [ fig : log - throughput ] , we plot the logarithms of the aggregate throughputs corresponding to the two proportional fair medium access schemes . the figure illustrates that as the number of nodes increases , the performance of the closest interferer based medium access worsens in comparison to the performance of the aggregate interference based scheme - this is not visible merely looking at the corresponding aggregate throughputs  ( see figure  [ fig : throughput ] ) ) .        in figure [ fig : map - cdf ] , we plot the c.d.f . of the map for the proportional fair case . our objective is to compare the case when nodes account for the aggregate interference with when they account for the closest interferer only . we set @xmath215 and plot map distributions corresponding to two values of @xmath1 , @xmath233 and @xmath234  ( corresponding to @xmath219 and @xmath218 respectively ) . as expected , the nodes attempt more aggressively when they account for the closest interferer only . while there is a favored probability interval in the aggregate interference case , there are more than one such intervals in the closest interferer case . also , in the latter case , about @xmath235 of nodes attempt in almost all slots , irrespective of @xmath1 , the total number of nodes . this can be understood by noticing that the probability that a node is not the closest interferer to any other node is not sensitive to @xmath1 . we have shown the feasibility of the performance analysis of distributed adaptive protocols that aim at maximizing some global utility in a large random network using stochastic geometry . more precisely , the most natural distributed adaptation of the medium access probability of aloha that aims at proportional fairness optimization was shown to have a tractable optimal map distribution . this distribution is obtained from the law of a certain shot noise field that describes the interference created by a typical node to all receivers but his . in the poisson case , the distribution of the optimal map is obtained as a non - singular contour integral which is amenable to an efficient evaluation using classical numerical tools . the network performance at optimum can in turn be deduced from the latter using campbell s formula . this approach is shown to provide an analytic way of quantifying the gains brought by this proportionally fair adaptive version of aloha compared to plain aloha . this line of thoughts opens several research directions . the first one is the extension to other types of fairness , still in the framework of aloha . we indicate that this is possible at least under certain simplifications of the interference model . the second and broader question is whether this approach can be extended to mac protocols other than aloha . an example would be an adaptation of the exclusion radius of csma / ca to the full environment of a node aiming at maximizing some utility of the throughput . a third general question concerns the evaluation of the `` price of decentralization '' . when the discussed protocols are suboptimal because of their greedy / distributed nature , is it possible to use stochastic geometry to evaluate the typical discrepancy between the performance of the distributed scheme and the optimal centralized one ? this work was carried out at laboratory of information , networking and communication sciences  ( lincs ) paris , and was supported by inria - alcatel lucent bell labs joint research center .    10    n.  abramson . : another alternative for computer communications . in _ proceedings of the fall joint computer conference  ( afips ) _ , montvale , new jersey , november 1970 . f.  baccelli and b.  blaszczyszyn . , volume 3  issue  3 - 4 of _ foundations and trends in networking_. now publishers , 2009 . f.  baccelli and b.  blaszczyszyn . , volume 4  issue  1 - 2 of _ foundations and trends in networking_. now publishers , 2009 . f.  baccelli , b.  blaszczyszyn , and p.  muhlethaler . an aloha protocol for multihop mobile wireless networks . , 52(2):421436 , february 2006 . f.  baccelli , b.  blaszczyszyn , and p.  muhlethaler . stochastic analysis of spatial and opportunistic aloha . , 27(7):11051119 , september 2009 . f.  baccelli , j.  li , t.  richardson , s.  shakkottai , s.  subramanian , and x.  wu . on optimizing csma for wide area ad hoc networks . , 72(1 - 2):3168 , 2012 . d.  p. bertsekas and r.  g. gallager . . prentice hall , 1992 . s.  boyd and l.  vandenberghe . . cambridge university press , 2004 . p.  bremaud . . number  31 in texts in applied mathematics . springer - verlag , 1999 . m.  franceschetti and r.  meester . . cambridge university press , 2008 . b.  hajek . stochastic approximation methods for decentralized control of multiaccess communications . , 31(2):176184 , march 1985 . m.  k. hanawal , e.  altman , and f.  baccelli . stochastic geometry based medium access games . in _ ieee infocom _ , orlando , florida , usa , march 2012 . hsu and h .- j . channel - aware aloha with sinr capture : when is the knowledge of channel not helpful ? in _ proceedings of ieee pimrc _ , toronto , canada , september 2011 . mohsenian - rad , v.  w. wong , and r.  schober . optimal sinr - based random access . in _ proceedings of ieee infocom _ , san diego , ca , usa , march 2010 . d.  monderer and l.  s. shapley . potential games . , 14(1):124143 , may 1996 . m.  osborne and a.  rubinstein . . the mit press , 1999 . l.  g. roberts . packet system with and without slots and capture . , 5(2):2842 , april 1975 .    x.  wang and k.  kar . distributed algorithms for max - min fair rate allocation in aloha networks . in _ 42nd annual allerton conference on communication , control , and computing _ , urbana , il , usa , september - october 2004 .    x.  wang , k.  kar , and j .- s . lexicographic max - min fair rate allocation in random access wireless networks . in _ 45th ieee conference on decision and control _ , san diego , ca , usa , december 2006 . let @xmath236 be a point process with marks and @xmath237 be the point process @xmath238 shifted by @xmath239 . given that @xmath240 , @xmath241 where the mapping @xmath242 associates with @xmath243 and the real number @xmath45 the solution of @xmath244 with the usual convention if there is no solution in @xmath55 $ ] . it follows from slivnyak s theorem that @xmath245=\\ & & \hspace{-.7 cm } \mathbb{e}^0 \left[\sum_{n\ne 0 } \log \left(1-\frac{h(||x_n - y_0||,\phi^r\setminus \{y_0,y_n\}\circ\theta_{x_n})}{\frac{\vert x_n - y_0\vert^{\alpha}}{t r_0^{\alpha } } + 1}\right ) \right]=\\ & & \hspace{-.7cm}\frac 1 { 2\pi } \int\limits_0^{2\pi } \mathbb{e } \sum_{n } \log \left(1-\frac{h(||x_n - ( r_0,\phi)||,\phi^r \setminus \{y_n\}\circ\theta_{x_n } ) } { \frac{\vert x_n-(r_0,\phi)\vert^{\alpha}}{t r_0^{\alpha } } + 1}\right ) { \rm d}\phi.\end{aligned}\ ] ] it now follows from campbell s formula that @xmath246\\ & & \hspace{-.6 cm } = \int_{\mathbb{r}^2 } \int_v   \log\left(1-\frac{v\overline r_0}{||x-(r_0,\phi)||^\alpha+\overline r_0}\right ) \lambda { \rm d}x\\ & & \hspace{3 cm } \mathbb{p}^0 ( h(||x-(r_0,\phi)||,\phi^r\setminus\{y_0\})= { \rm d}v)\\ & & \hspace{-.6 cm } = \int_{\mathbb{r}^2 } \int_v   \log\left(1-\frac{v\overline r_0}{||x-(r_0,\phi)||^\alpha+\overline r_0}\right ) ) \lambda { \rm d}x \\ & & \hspace{3 cm } \mathbb{p } ( h(||x-(r_0,\phi)||,\phi^r)= { \rm d}v),\end{aligned}\ ] ] where the last relation follows from slivnyak s theorem .
Abstract:

Answer: this work aims at combining adaptive protocol design , utility maximization and stochastic geometry . we focus on a spatial adaptation of aloha within the framework of ad hoc networks . we consider quasi - static networks in which mobiles learn the local topology and incorporate this information to adapt their medium access probability  ( map ) selection to their local environment . we consider the cases where nodes cooperate in a distributed way to maximize the global throughput or to achieve either proportional fair or max - min fair medium access . in the proportional fair case , we show that nodes can compute their optimal maps as solutions to certain fixed point equations . in the maximum throughput case , the optimal maps are obtained through a gibbs sampling based algorithm . in the max min case , these are obtained as the solution of a convex optimization problem . the main performance analysis result of the paper is that this type of distributed adaptation can be analyzed using stochastic geometry in the proportional fair case . in this case , we show that , when the nodes form a homogeneous poisson point process in the euclidean plane , the distribution of the optimal map can be obtained from that of a certain shot noise process w.r.t . the node poisson point process and that the mean utility can also be derived from this distribution . we discuss the difficulties to be faced for analyzing the performance of the other cases ( maximal throughput and max - min fairness ) . numerical results illustrate our findings and quantify the gains brought by spatial adaptation in such networks .