Instruction: Write an abstract for the following acadamic paper.

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

Paper: the two - dimensional ( 2d ) ferromagnetic xy model has had a long history of studies concerning the dynamics of spinwaves at low temperatures . it was first analyzed by villain@xcite and by moussa and villain@xcite , for low temperatures , in the harmonic approximation . they found that the dynamic structure function @xmath1 has a spinwave peak of the form @xmath3 where @xmath4 is the critical exponent describing the decay of the static spin - spin correlation function , and @xmath5 is the spinwave frequency at wavevector @xmath6 . nelson and fisher@xcite treated the model in a fixed - length hydrodynamic description , obtaining the following expression around the spinwave peak , @xmath7 where for small @xmath8 , @xmath9 , with @xmath10 being the spinwave velocity . these theoretical approaches led to spinwave peaks diverging unphysically at the spinwave frequency . here we describe an alternative theoretical description in which the divergence does not occur .    at higher temperatures , the spinwave or harmonic approximation begins to fail , and xy magnets are known to undergo a topological phase transition@xcite associated with the unbinding of vortices . we generally consider temperatures below this unbinding temperature @xmath11 , where a spinwave approximation will be valid . at very low temperatures any nonlinear effects should be very weak , do not lead to vortex production , and can be included as the damping or scattering of the spinwave modes , that results in a spinwave lifetime . thus our approach will be to start from the harmonic approximation to the spin hamiltonian , @xmath12 , and then include the damping that results from the terms fourth order in spin operators ( @xmath13 ) in an expansion of the original full spin hamiltonian , @xmath14 . the resulting damping strength @xmath0 will imply a finite lifetime @xmath15 for the spinwave modes . we are interested to investigate whether the presence of such lifetime effects can eliminate the divergence of the dynamic structure function @xmath1 at the spinwave frequency . we calculate the spinwave damping by a feynmann diagram technique for finite temperatures , applied to the fourth - order hamiltonian @xmath13 , using the leading diagrams to second order in @xmath13 . the only diagrams we include are those that give a contribution to the _ imaginary part _ of the spinwave energy , i.e. , to the damping or inverse lifetime . these are the terms that involve a _ dynamical _ effect on the spinwave energies . the associated integral expressions for these diagrams require the evaluation of four - dimensional integrals involving a delta - function . methods for evaluating this type of integral in three dimensions have been developed by gilat and co - workers ; here we show how to extend their work and apply the same technique in four dimensions . in principle , we could also use the feynmann technique to determine the leading changes in the _ real part _ of the spinwave energy , @xmath16 , due to the perturbation @xmath13 . this real part would result as the usual effect due to finite temperature : a shift in the spinwave frequencies to lower values as the temperature increases . it is essentially a _ static _ effect , primarily due to the renormalization of the nearest neighbor exchange interaction strength @xmath17 with temperature . calculation of these spinwave frequency shifts by the feynmann technique would require including diagrams first order in @xmath13 ; the so - called single - loop diagrams . in this paper , however , we are primarily concerned with accounting for dynamical effects . we will include the renormalization of the spinwave energies using the memory function formalism , combined with some input information from classical monte carlo ( mc ) data . in the memory function approach , the position of the spinwave peak in @xmath1 does not occur exactly at the zero - temperature dispersion relation , @xmath5 , but rather , close to the value @xmath18 , where @xmath19 is the temperature - dependent second moment . @xmath19 can be estimated from static quantities , either using low - temperature analytic expressions or results from monte carlo , the latter being necessary for temperatures approaching @xmath11 . the memory function approach already includes to leading order the exchange strength renormalization , as can be seen in the shift in the position of the spinwave peaks . thus we use it rather than making any renormalization of the exchange @xmath17 . a few comments can be made about how the damping @xmath0 is included into the memory function calculation . note that the damping calculation sketched above is completely independent of the memory function calculation ; the latter is a formalism for calculating dynamic correlations for a given hamiltonian and operators . due to the xy symmetry of the model , the natural coordinates of the spins are the in - plane angle @xmath20 and out - of - plane spin component @xmath21 . fourier transforming and using a standard canonical transformation , the hamiltonian is written in terms of creation and annihilation operators @xmath22 and @xmath23 , with time dependencies varying as @xmath24 and @xmath25 , respectively . once damping is included , the time dependencies become @xmath26 and @xmath27 , respectively . henceforth the calculation proceeds as usual in the memory function approach , except for this small modification of the operator time dependencies via the shifted operator frequencies , @xmath28 . in the memory function calculation , spinwave sum and difference processes enter . these are characterized by combinations of frequencies , @xmath29 , where @xmath6 is the wavevector at which the dynamic correlations are measured , and @xmath30 is summed over to include all possible such processes . based on the above modifications of the operator frequencies , there also results the modification for these processes , @xmath31 . the damping rates for different wavevectors always add for both types of processes . the paper is organized as follows : first , in sec . [ model ] , we present the model and develop the expansion of the hamiltonian to quartic order in creation / annihilation operators . in sec . [ diagrams ] , the expressions for the damping via the feynmann technique are presented . then we show how these are evaluated using the gilat procedure . in sec . [ memory ] , the key aspects of the memory function calculation of @xmath1 are reviewed , together with a focus on the modification due to the damping . some details of the monte carlo and spin dynamics simulations are summarized in sec . [ mc+sd ] . this is followed in sec . [ results ] by presentation of the results both for very large and relatively smaller systems , in order to discuss the dynamical finite size effects . the memory function results also are compared with standard monte carlo  spin - dynamics calculations of the dynamic correlation functions in order to evaluate the validity and accuracy of the method . our conclusions are presented in sec . [ conclusions ] . the ferromagnetic spin hamiltonian under consideration is @xmath32 where the first sum is over all sites @xmath33 of a 2d square lattice , the second is over all nearest neighbor displacements @xmath34 , @xmath17 is the nearest neighbor coupling strength , and each @xmath35 is a three - component classical spin - vector . because each bond is counted twice , a factor of @xmath36 is included on @xmath17 . we primarily concentrate on the xy model , with anisotropy parameter @xmath37 , however , for completeness , formulas covering the easy - plane anisotropy range , @xmath38 , are presented .      in order to develop the damping as a perturbation of the harmonic approximation , we need to expand the full spin hamiltonian , including the harmonic ( quadratic ) and quartic deviations from the ground state . the ground state corresponds to spins aligned in any direction within the xy plane . writing the spins in hamiltonian ( [ hamil ] ) using coordinates @xmath39 and @xmath40 , @xmath41 the hamiltonian becomes @xmath42 ~ , \ ] ] where we use a simplified notation for the nearest neighbors , @xmath43 . at low temperature , it is reasonable to expand for small @xmath40 and @xmath44 , keeping up to quartic order in these , obtaining , @xmath45 where the harmonic and quartic hamiltonians are @xmath46,\end{aligned}\ ] ] @xmath47 ~ .\end{aligned}\ ] ] it is convenient to carry out the calculations in wavevector space by defining fourier transforms @xmath48 together with a similar definition for @xmath49 . here @xmath6 denotes a two - dimensional wavevector , the sums are over the first brillouin zone ( bz ) , and @xmath50 is the total number of lattice sites . then we rewrite @xmath12 and @xmath13 as @xmath51 ~ , \ ] ] and @xmath52 for the square lattice , with coordination number @xmath53 , the @xmath54 function @xmath55 is an even function of @xmath6 , i.e. , @xmath56 . the sums in ( [ h1 ] ) are constrained to @xmath57 . using the following canonical transformation to creation / annihilation operators , @xmath58 @xmath59^{1/4 } , \quad\quad \beta_{\bq}=\left[\frac{1-\gamma_{\bq } } { 4s^2(1-\lambda\gamma_{\bq})}\right]^{1/4 } ~ , \ ] ] the quadratic hamiltonian is diagonalized , @xmath60 ~ , \ ] ] @xmath61 the diagonal parts of @xmath13 can then be written as @xmath62 where @xmath63 we have dropped terms in @xmath13 whose expectation values in the ground state are zero . note also that if a particular @xmath64 is zero , then @xmath65 diverges , however , when this occurs the corresponding factors multiplying @xmath65 go to zero faster , canceling the divergence . the time evolution of the creation / annihilation operators as given from the harmonic hamiltonian is @xmath66 as discussed in the introduction , the effect of the @xmath13 perturbation is to shift the spinwave frequencies , in the most general case , by some complex self - energy @xmath67 . for the moment , we ignore the real part of @xmath67 , whose effect will be included by the memory function approach , and concentrate instead on the contributions to the imaginary part of @xmath67 , which determine the spinwave damping . to find @xmath67 , the contributions of @xmath13 to the single particle green s function will be evaluated , keeping only the diagrams that have a nonzero imaginary part . the lowest order diagrams that satisfy this requirement are found to be second order in @xmath13 ; there is no contribution from first order diagrams .    following the standard procedure,@xcite the time - ordered zero - temperature green s function , in second order , is @xmath68 here @xmath69 is the time ordering operator , and @xmath70 is the ground state of @xmath12 . physically , this green function represents the amplitude for a quasiparticle ( i.e. , spinwave ) created at time @xmath71 in state @xmath72 to be destroyed at time @xmath73 in state @xmath74 , after undergoing two scattering events due to @xmath13 . the free - particle green function , @xmath75 involves no scattering . due to momentum conservation , these green functions are proportional to the kronecker delta , @xmath76 , so we can drop the double momentum indices . also , because of time translational invariance , they depend only on the time difference @xmath77 . introducing the time fourier transforms , @xmath78 the free propagator becomes the well - known expression , @xmath79 where @xmath80 . we apply the rules for feynmann diagrams to evaluate the full propagator , @xmath81 , in terms of @xmath82 and @xmath83 . the free propagator @xmath84 is represented as a straight line with an associated momentum @xmath30 and frequency @xmath85 . the interaction term in ( [ h_1 ] ) is represented as a wavy line as in fig . [ aq1q2q3q4 ] , with an outgoing straight line ( free propagator ) for each creation operator and an incoming straight line for each destruction operator , both with assigned momentum and frequency . the two leading diagrams that contribute to @xmath86 are shown in fig . [ diagrams ] . following the standard feynmann rules for these diagrams , with summations over the intermediate frequencies and wavevectors , and using momentum and frequency conservation at each junction , @xmath86 can be written as @xmath87 where the self - energy function due to the scattering parts of the diagrams is @xmath88 ~ .\end{aligned}\ ] ] @xmath89 is a sum over the `` direct '' and the `` exchanged '' diagrams ; the factor of 4 is the multiplicity of each of these under rearranging the exterior lines . @xmath89 is seen to be @xmath86 with its incoming and outgoing lines removed . the full interacting propagator @xmath81 can be approximated as a sum over repeated applications of @xmath90 , i.e. , @xmath91 which , after summing , leads to the well - known dyson equation , @xmath92 the location of the pole gives the perturbed quasiparticle energy , @xmath93 , and the damping @xmath94 is given by @xmath95    to apply the result to the system at finite temperature @xmath69 , and obtain the self - energy associated with the temperature green s function , the integrals over frequency are transformed to sums according to @xmath96 @xmath97 thus the frequency integrals in ( [ sigma ] ) become @xmath98 [ n(\omega_{\bq}+\omega_{\bk+\bp-\bq})-n(\omega_{\bp } ) ] } { \omega-(\omega_{\bk+\bp-\bq}+\omega_{\bq}-\omega_{\bp } ) } ~ , \end{aligned}\ ] ] where the bose factors are @xmath99 the latter form applies to the classical limit . the imaginary part of @xmath100 is derived from the analytic relation , @xmath101 where @xmath102 takes the principal part and @xmath103 is the usual delta - function . because @xmath6 and @xmath104 each have two components , the sum in ( [ sigma ] ) is equivalent to a four - dimensional integration , once we take the continuum limit . each dimension converts as , for example , @xmath105 , where @xmath106 is the system size , in terms of the lattice constant @xmath107 , leading to @xmath108 the factors of @xmath50 here are canceled by those in the expression for @xmath109 , leading to a size - independent result . thus there finally results for the damping function , @xmath110 where the kernel @xmath111 is given by @xmath112 \times \nonumber \\ & & [ n(\omega_{\bq})-n(-\omega_{\bk+\bp-\bq } ) ]   [ n(\omega_{\bq}+\omega_{\bk+\bp-\bq})-n(\omega_{\bp } ) ] ~ .\end{aligned}\ ] ] we use this result to calculate the damping at the spinwave frequency : @xmath113 , in the thermodynamic limit @xmath114 . note that in the classical limit , the bose factors are simply proportional to temperature . therefore , we only need to calculate @xmath115 at one temperature ( @xmath116 ) and then scale with @xmath117 to get its value at other temperatures . the summation / integration in ( [ gamma ] ) involving a delta function presents an interesting mathematical problem in its evaluation . general integrals of this form ( i.e. , with one delta - function ) in @xmath118 dimensions can be considered , @xmath119 where @xmath120 is a @xmath118-dimensional vector . in two dimensions , for example , the delta - function will pick off a particular contour that contributes to the integral ; in three dimensions that contour will actually be a surface , and in four dimensions it will be a three - dimensional volume . depending on the form of the function @xmath121 , the contributing region need not be simply connected , but may consist of separated distinct parts . for the higher dimensional cases , it is difficult to visualize the region or regions that contribute to the integral . thus it is necessary to have a way to evaluate the integral without appealing directly to geometrical properties . in a series of papers , gilat and co - workers@xcite have considered how to evaluate these integrals in three dimensions , in work related the calculation of spectral density functions . we have used the `` linear analytic '' ( la ) method@xcite of gilat and raubenheimer , extended to four dimensions , which gives good precision for few sample points , at little cost in complexity . suppose we calculate @xmath122 over a small range of @xmath85 and then take the limit @xmath123 . to do so , a 4d cubic grid of small cells of size @xmath124 is set in the @xmath120-space . then we consider the contribution of each cell to @xmath122 , thereafter summing over the cells . suppose the cells are small enough so that the kernel @xmath125 has only a small variation within a cell , and , therefore , we can use its value at the cell center , @xmath126 . then we can write @xmath127 @xmath128 the idea of the la method is to expand @xmath121 around the cell center , @xmath129 shifting coordinates to @xmath130 , and writing the gradient at the cell center as @xmath131 , where the unit vector @xmath132 points along the gradient evaluated at the cell center , we have @xmath133 @xmath134 the variable @xmath135 is a rescaling of @xmath85 . if we use the integral representation of the delta function , in the form , @xmath136 then the contribution from one cell is found to be @xmath137 where the @xmath138 are the values of parameter @xmath135 at the 16 corners of the cell , i.e. , @xmath139 with all 16 possible choices of the signs , and @xmath140 is the number of minus signs used in each @xmath138 . one can note that similar simple formulas result for this type of integral in other dimensions as well.@xcite    in the actual application of eq . ( [ dcell ] ) , inside the cell , the parameter @xmath135 will range between minimum and maximum values : @xmath141 these values hold at two opposite corners of the cell . by eq . ( [ w - omega ] ) , there will be a corresponding range of @xmath85 over which this cell contributes to @xmath122 . then for the general problem where @xmath142 is desired over a large range of @xmath85 , what is done is to set a uniform grid of @xmath85-values , finding the contribution of each cell to @xmath122 for each value on the grid , say , @xmath143 . typically it works well to have the grid fine enough so that each cell contributes to @xmath142 at many different @xmath144 .      for the damping problem , eq . ( [ gamma ] ) , we have the replacements , @xmath145 with the 4d variable @xmath146 . we want to find @xmath147 right at @xmath148 for some particular wavevector , @xmath30 . to accomplish this , we set a range @xmath149 , and define the grid of values @xmath150 , using @xmath151 . for the integrations over @xmath6 and @xmath104 , the cell sizes actually can be quite large : we used the cell width @xmath152 , where @xmath153 is a `` system size , '' which gave results consistent with those using much smaller cells ( larger @xmath154 ) . that is , each axis , @xmath155 , @xmath156 , @xmath157 , @xmath158 , was partitioned into 41 cells . first , to get good precision , a narrow range symmetrical around the spinwave frequency was used : @xmath159 , @xmath160 . second , for comparison , the behavior of the damping function was calculated over the whole range of possible spinwave frequencies , i.e. , @xmath161 , @xmath162 . in both of these approaches , interpolation on the @xmath85-grid gives @xmath163 . the scheme was found to be very robust , with both approaches agreeing to better than one percent precision . we should also note one other important calculational detail . because eq . ( [ dcell ] ) depends on the reciprocal gradient of @xmath164 , it is necessary to avoid any points in the @xmath165-space where this gradient vanishes . we accomplished this by taking wavevectors @xmath166 , ( at @xmath153 ) , and then shifting the @xmath165-cell center positions slightly from the usual values : @xmath167 , @xmath168 . here @xmath169 , etc . , are integers ranging from 0 to @xmath170 , and we used @xmath171 to shift the @xmath165-cells away from high symmetry points . the choice of an odd number for @xmath154 also helps to avoid certain high symmetry points . some slight improvement in precision can be achieved with a larger value of @xmath154 , however , at the expense of a tremendous increase in cpu time , which varies as @xmath172 . at @xmath153 , the typical relative errors in @xmath147 were found to be much less than one percent , when compared to calculations using @xmath173 . some results for the damping function , which is proportional to @xmath117 , are shown in fig . [ gam_omega - fig ] for @xmath37 . in general , it is found that @xmath147 either has a cusp or a peak at the spinwave frequency ( arrows in fig .  [ gam_omega - fig ] ) . the damping strength increases with @xmath174 , as seen in fig . [ gam_k - fig ] . along the ( 11 ) direction , however , @xmath115 tends to saturate halfway out in the zone . keeping in mind that @xmath175 , the damping strength we find here is quite weak , especially when we consider its effects at @xmath176 where the spinwave approximation is valid . at small nonzero values of the anisotropy parameter @xmath177 , there will be only a very slight change in @xmath115 , because @xmath177 does not appear explicitly in @xmath13 , but only enters the calculation via the dispersion relation and via the @xmath178 and @xmath179 functions of ( [ canonical ] ) . to demonstrate this effect , fig . [ gam_lambda - fig ] shows @xmath115 at various values of @xmath177 for @xmath180 along the ( 10 ) direction . in general , the damping decreases with increasing @xmath177 , i.e. as the system becomes more isotropic . the calculation of the memory function requires the knowledge of @xmath115 over the entire bz . we can reduce the region of calculation to @xmath181 , because , by symmetry , @xmath182 , and @xmath183 . in this region @xmath115 was calculated on a uniform square grid with spacing @xmath184 . linear interpolation was used to get values between these grid points . for the calculation of the dynamic correlation function @xmath185 , we use the memory function formalism due to reiter and sjlander,@xcite who applied the procedure to the isotropic heisenberg model . the technique is based on a low - temperature expansion of the memory function.@xcite the formalism was applied to the 1d isotropic antiferromagnet in an external field by de raedt et al@xcite and to a 1d easy - plane antiferromagnet by gouva and pires.@xcite more recently , it was also applied to 2d easy - plane ferromagnets by menezes et al,@xcite although some numerical difficulties inhibited numerical evaluation of the analytic formulas . some of these problems were surpassed in the work of pereira et al.@xcite in all of these works , no mode damping was included . the presence of damping actually relieves the numerical difficulty of the problem , as we see below . the dynamic correlation function is defined in terms of the fourier transformed spin - fields as @xmath186 where @xmath178 is any component of @xmath187 , and by symmetry , @xmath188 . in this and following formulas for correlation functions , we note the distinction between continuum ( functions of @xmath6 ) and discrete ( subscripts @xmath6 ) fourier transforms , which are related by factors of @xmath2 for each dimension.@xcite this is necessary so that we can make a careful comparison between the theory and the mc - sd simulations , _ including the absolute magnitudes_. specifically , the continuum correlation function in ( [ sqw ] ) , which gives a thermodynamic weight over a _ wavevector range _ @xmath189 , is related to the discrete correlation at a _ specific _ wavevector @xmath6 , using discrete fourier transforms as defined in ( [ fourier ] ) , by@xcite @xmath190 in the analytic calculations , many intermediate quantities are discrete functions , however , for the comparison with the mc - sd data , both results will be converted and displayed as continuum quantities .    due to the symmetry in the @xmath191-plane , the calculation is carried out by making averages of the rotationally symmetric quantity , @xmath192 , with @xmath193 . then in the memory function calculation , the continuum dynamical structure function can be written in terms of a second - order memory function , @xmath194 , to be described below , according to @xmath195 ^ 2 + \left [ \omega ~{\rm i m } ~ \sigma_{\bq}^{\perp}(\omega ) \right]^2 } ~ , \end{aligned}\ ] ] where @xmath196 is the liouville operator , and the second moment @xmath197 is determined from a ratio , @xmath198 the location of the spinwave peak that results from eq . ( [ sqw - memory ] ) is very sensitive to the value of @xmath197 , therefore we are especially concerned to get an accurate evaluation of it . the denominator of ( [ moment2 ] ) is a static correlation function , whose value will be taken from the monte carlo simulations , although for low temperatures it can be evaluated approximately as    [ sqperp ] @xmath199 @xmath200    the numerator of ( [ moment2 ] ) has been shown@xcite to have a leading low - temperature dependence , also related to static nearest neighbor correlations , @xmath201 the correlation @xmath202 , therefore its contribution can be ignored for small @xmath177 . there being two bonds per lattice site , the in - plane near - neighbor correlation is related to the average energy per bond by @xmath203 this correlation can also be obtained from the mc simulations , although a good low - t approximation is @xmath204 a comparison of the analytic and mc static parameters will be given below . to the leading order in temperature , the time - dependent memory function has been shown to be given from    @xmath205    @xmath206    where @xmath207 is projection onto the non - secular variables , see ref . , whose action is equivalent to @xmath208 using expressions ( [ mapping ] ) expanded to second order in @xmath39 , @xmath40 , and ( [ canonical ] ) , it is straightforward to write the spins @xmath209 in terms of the creation / annihilation operators . in the leading order in temperature , the resulting expressions involve expectation values of products of four creation / annihilation operators . the damping is included in the time evolution according to , @xmath210 then , there results , after some manipulations , expressions to give the memory function , @xmath211 the two terms correspond to spinwave sum ( @xmath212 ) and difference ( @xmath213 ) processes , with weights @xmath214 where    @xmath215 \bigr\ } ~ , \end{aligned}\ ] ]    @xmath216 \bigr\ } ~ . \end{aligned}\ ] ]    the frequencies and damping rates of these processes are    @xmath217    @xmath218    finally , the frequency - dependent memory function is given by the one - sided fourier - laplace transform , @xmath219 from this expression , eq . ( [ m - memory ] ) , and equations following , we obtain a sum over all the sum / difference processes , @xmath220   \nonumber \\ & &   + w_{- } \bigl [    \frac{\omega+\omega_{-}-i\gamma_{-}}{(\omega+\omega_{-})^{2}+\gamma_{-}^{2 } } + \frac{\omega-\omega_{-}-i\gamma_{-}}{(\omega-\omega_{-})^{2}+\gamma_{-}^{2 } } \bigr ] \bigr\ } ~ .\end{aligned}\ ] ] where the @xmath221 dependencies of @xmath222 , and @xmath223 have been suppressed . this is the expression used to evaluate @xmath224 , and thus the dynamic structure function , via eq . ( [ sqw - memory ] ) . to evaluate ( [ sigma_sum ] ) on a discrete @xmath225 lattice , the wavevector @xmath6 must be one of the allowed values , @xmath226 where @xmath227 and @xmath228 are integers . it is also necessary to make the summations over @xmath104 such that the vectors @xmath229 are both given by integers times @xmath2 . thus , it is convenient to change from the summation over @xmath104 to a summation over a new variable @xmath30 defined by @xmath230 where @xmath231 and @xmath232 are integers . then , we also have @xmath233 equal to an allowed wavevector . this insures that the sums will be over spinwaves that are allowed on the given finite lattice . this approach is essential for calculation of the finite size dynamical effects . on the other hand , for the infinite sized system , no variable change is needed , and there is the usual replacement , @xmath234 . note that in the limit of zero damping , eq .  ( [ sigma_sum ] ) is inapplicable , except for the continuum limit ( i.e. , @xmath114 ) . the resulting continuum integrals for the infinite system are singular , and the analytic relation ( [ delta - relation ] ) must be applied , leading to @xmath235   \nonumber \\ & &   + p \int d^{2}p ~ w_{- }   \bigl [ \frac{1}{\omega+\omega_{- } } + \frac{1}{\omega-\omega_{- } } \bigr ]   \nonumber \\ & &   -i \pi   \int d^{2}p ~ w_{+ }   \bigl [ \delta(\omega+\omega_{+ } ) + \delta(\omega-\omega_{+ } ) \bigr ] \nonumber \\ & &   -i \pi   \int d^{2}p ~ w_{- }   \bigl [ \delta(\omega+\omega_{- } ) + \delta(\omega-\omega_{- } ) \bigr ] \bigr\ } ~ .\end{aligned}\ ] ] indeed , in the very low temperature range ( @xmath236 ) where the damping is extremely weak , it may be more efficient to use this expression rather than eq . ( [ sigma_sum ] ) , which has poorer convergence with @xmath154 . ( [ sigma_int ] ) was applied by pereira et al.@xcite ; the imaginary part can be evaluated by techniques described here ( sec .  [ integration ] ) , and the principal - valued real part can be calculated by an extension@xcite of techniques used in 1d.@xcite in ref .   the last expression was found to give finite amplitude spinwave peaks , without damping present .    before presenting the results obtained via the memory function approach , the key aspects of the monte carlo and spin - dynamics simulations are summarized . to test the validity of the memory function technique , we compare it with numerical simulations on @xmath237 square lattices ( @xmath238 , @xmath37 ) using monte - carlo and spin - dynamics ( sd ) simulations , which include effects due to all thermodynamically allowed excitations . as mentioned above , we also use some data for static correlations from the monte carlo calculations as input to the memory function dynamical calculations . the techniques used here have been described in ref .   and are based on the simulation methods developed in ref .  . in the mc calculation , we applied a combination of metropolis single - spin moves and over - relaxation moves that modify all three spin components , and in addition , wolff single - cluster operations@xcite that modify _ only _ the @xmath191 spin components . the over - relaxation and cluster moves are important at low temperatures , where the @xmath191 spin components tend to freeze and single spin moves become inefficient . in the combined mc - sd simulations , a set of mc generated states from thermodynamic equilibrium at temperature @xmath69 are used as initial conditions for numerical integrations in time via a fourth order runge - kutta scheme . specifically , the first 4000 mc steps were used for equilibrating the system . then , every 400 mc steps a state was taken to initiate the sd integration . this was repeated 500 times ; the displayed @xmath1 data are averages over the 500 time integrations . the total mc sequence was then @xmath239 steps ; every 25 steps a state was used to calculate averages of the static correlation functions . in the time integration , the basic runge - kutta time step was @xmath240 , and every 11 time steps , data samples ( for @xmath241 , etc . ) were saved for @xmath6 s in the ( 10 ) , ( 01 ) , and ( 11 ) directions only ( due to machine memory restrictions ) . a total of @xmath242 data samples were saved , followed by applying a fast - fourier - transform to obtain @xmath243 . the total time of integration of @xmath244 is large enough so that no smoothing window is necessary for the data generation . then the frequency resolution of these results is good , @xmath245 , so that any fine details in the low temperature results should be visible . for the lowest temperatures ( @xmath246 ) , where we wanted to see the complex fine details due to the finite size effects , some simulations were made saving data samples every 28 times steps , leading to a total integration time @xmath247 , and frequency resolution @xmath248 . first , a comparison of the low-@xmath69 formulas and static mc results needed as input to the memory function formalism is presented . then dynamical results from the memory function and from md - sd will be presented for large systems , and finally , for smaller systems , where the finite size features will be analyzed . the continuum static structure function , @xmath249 was measured directly in the mc simulations . a comparison with the results based on the low - temperature expression ( [ sqperp ] ) is given in fig . [ sq - fig ] . there we also show a fit of the mc data to a form similar to the analytic expression , but with a fitting coefficient ( @xmath250 ) and power ( @xmath251 ) , @xmath252 from the fits , the parameter @xmath250 is somewhat larger than the low - temperature limit , @xmath253 . this can be expected due to the renormalization of the effective @xmath17 to lower values with temperature . the power @xmath251 falls with increasing @xmath69 , which must also be due to anharmonic effects . there is little difference between the low - t expression and the mc data at the lowest temperature . the deviations between the mc and low - t theory increase substantially by @xmath254 , especially at higher wavevectors . note that the fits are only shown for demonstration ; we used the mc data itself as inputs to the memory function calculations . the other static input to the memory function calculation , for determining the second moment , is the quantity @xmath255 , which is related to the nearest neighbor correlation @xmath256 by eq . ( [ lsls ] ) . estimates of that latter from the mc evaluation of the average energy per spin [ eq . ( [ s0sa - en ] ) ] are compared to the theoretical expression ( [ s0sa - theory ] ) in fig . [ s0sa - fig ] . the analytical expression is seen to give results consistent with the mc data up to @xmath257 . when combined to give the second moment , the mc data produces results as shown in fig . [ moment2-fig ] . it is essentially this quantity that has the strongest influence on the locations of the spinwave peaks in @xmath1 . these are the data used as input to the memory function calculation . we evaluate @xmath194 by using the sum over discrete wavevectors @xmath104 as displayed in eq . ( [ sigma_sum ] ) , and then apply eq . ( [ sqw - memory ] ) to determine @xmath193 . if @xmath154 is taken large enough , then the sum should be a good approximation to the infinite size limit . in general , we can see whether larger @xmath154 is necessary by the presence of finite size features in the data for @xmath224 and @xmath243 . typically we used @xmath258 to give a reasonable approximation to the continuum limit for @xmath259 and larger @xmath6-values ; at lower temperatures or wavevectors even larger @xmath154 is necessary to smooth out all the finite size features , due to the very weak damping . in fig . [ sigt3qpi4pi4-fig ] , some typical results for the real and imaginary parts of @xmath194 are shown for @xmath37 , @xmath260 , and @xmath261 , where @xmath262 was used in the sum . the resulting dynamical structure function is compared to the corresponding mc - sd data for @xmath238 in fig . [ sxxt3qpi4pi4-fig ] . at the ( zero-@xmath69 ) spinwave frequency @xmath263 , i m @xmath194 achieves its minimum value , which is not zero due to the presence of the damping , while re @xmath194 has a small nonzero value there as well . the spinwave peak in the memory data occurs at @xmath264 , which is just slightly greater than the associated value of the second moment , @xmath265 . the spinwave peak position will occur where the denominator of eq . ( [ sqw - memory ] ) achieves its minimum value , which gives , to a good approximation , @xmath266 for this case , the memory function data give @xmath267 , which is consistent with this relation . the peak position is very sensitive to the value of the second moment , which is why we have tried to estimate it as accurately as possible from the static mc data . still , the mc - sd data show the spinwave peak closer to @xmath268 , slightly lower than the memory function prediction . using eq . ( [ sqw - memory ] ) , the spinwave peak height is proportional to @xmath269^{-1}$ ] , while the peak width is approximately equal to @xmath270 . when compared to the mc - sd data , the memory function calculation usually overestimates the peak height and underestimates the peak width , although both are in reasonable agreement with the mc - sd data . in fact , these errors tend to be larger at lower temperature .    in figs . [ sxxt1-fig][sxxt5-fig ] , we display some further comparisons of the memory function calculations and mc - sd results at other values of @xmath6 for temperatures ranging from @xmath271 to @xmath272 . for the lower temperatures there is very good agreement between the memory theory and mc - sd peak positions , while the peak heights and widths do not agree as well . on the other hand , at higher temperatures the opposite is true : the memory data show greater errors in the peak positions , especially at large @xmath8 , but the peak heights and widths are predicted very well . probably these discrepancies are caused by a fairly large _ underestimate _ of the damping strength at low temperature , as well as errors in the estimate of the second moment at higher temperatures . however , in general , the tails of the spinwave peaks far from @xmath5 are described very well by the memory function calculations . when we consider the results using relatively smaller values of wavevector , @xmath243 shows not only a spinwave peak , but additional features ( small peaks ) which depend on the system size . these features are obvious in the mc - sd @xmath1 data shown above . similar size - dependent dynamical features have been noticed in earlier mc - sd simulations of xy - symmetry@xcite and easy - axis symmetry@xcite magnets , but without a clear theoretical explanation . we now consider a closer look at the details of these features , making a comparison with the memory function calculations at smaller values of @xmath154 . some examples showing the resulting memory function data for @xmath243 for @xmath238 systems are given in fig . [ sxxt3q16]a [ @xmath273 , @xmath274 , fig . [ sxxt3q9]a [ @xmath275 , @xmath274 , and fig . [ sxxt1q5]a [ @xmath276 , @xmath277 . in fact , there appear various sequences of small peaks . the number of separate sequences appears to depend on the wavevector in integer units of @xmath2 . note that the logarithmic scale exaggerates the size of these subpeaks , they are typically a few orders of magnitude smaller than the spinwave peak . for sufficiently low temperatures , these effects also appear clearly in the mc - sd results : see figs . [ sxxt3q16]b , [ sxxt3q9]b and [ sxxt1q5]b , although the corresponding features there are smeared out compared to the memory function calculation . the effects are also present at larger wavevectors , however , in that case there are many more of these finite size subpeaks , which tend to merge together into a smoother background . these extra peaks can each be associated with particular spinwave sum or difference process . to see that this is the case , consider first the difference processes , whose frequencies @xmath278 are always less than the spinwave frequency @xmath5 . a plot of the weights @xmath279 [ eq .  ( [ weights ] ) ] that determine the contributions to the memory function versus the corresponding frequencies @xmath280 , can show how strongly each process [ i.e. , each @xmath30 as defined in ( [ k - discrete ] ) ] contributes to i m @xmath224 and hence to @xmath1 . in fig . [ omminust3q16 ] , we made such a plot for an even integer @xmath6 , namely , @xmath281 , for @xmath238 and @xmath260 . each dot there corresponds to one allowed wavevector @xmath30 ; some of the values of @xmath30 in units of @xmath2 are indicated there . the sequence of @xmath30 s , ( 8,0 ) , ( 8,1 ) , ( 8,2 ) , etc . , leads to multiple contributions at zero frequency , that is , a weak central peak , which is a general result for any even integer @xmath6 . the sequence , ( 9,0 ) , ( 9,1 ) , ( 9,2 ) , etc . , produces a sequence of small peaks in @xmath1 out to @xmath282 , and , the next sequence , ( 10,0 ) , ( 10,1 ) , ( 10,2 ) , etc . , gives a sequence of peaks out to @xmath283 , and so on . each of these @xmath30-sequences can be identified with corresponding peak - sequences in the memory function calculation of @xmath1 , as seen in fig . . indeed , for the wavevector @xmath281 , there are @xmath284 of these sequences below the spinwave peak , including the ones that contribute only at @xmath285 . the higher sequences are more spread out , and merge with the other ones , making a complicated fine structure in @xmath1 . similar features are seen in the mc - sd data , fig . [ sxxt3q16]b , although it is clear that the individual peaks of each sequence are smeared out , resulting in one wider feature from the entire sequence . another example of these finite size features , for an odd integer wavevector , @xmath286 is shown in fig . [ ompmt3q9 ] , again for @xmath238 and @xmath260 . for this smaller integer wavevector there are now only 5 sequences of fine - structure difference process peaks . it is significant that there is no such peak at zero frequency . this is clearly caused by the absence of any allowed spinwave difference process exactly at @xmath287 , although there are many such processes with very weak weights _ close _ to @xmath287 . in general , for any odd integer wavevector , @xmath1 will have a relative minimum at @xmath285 , because of the absence of allowed difference processes there .    in fig . [ ompmt3q9 ] we also show some of the weights of the sum - processes , @xmath288 , indicated by plus signs . for either sum or difference processes , numbers in parenthesis refer to the values of @xmath30 in units of @xmath2 . the sum frequencies must fall above the spinwave frequency @xmath5 . comparing figures [ ompmt3q9 ] and [ sxxt3q9 ] , we see that the finite size features above @xmath5 in @xmath1 can each be identified with a particular spinwave sum process . the identification , however , is somewhat more complex than that for the difference processes , because the different sequences of processes overlap in frequency . the match between the sum peak frequencies in the memory data and the mc - sd data is not as good as for the differences processes below @xmath5 . we can note that at these higher frequencies , any relative errors appear larger on an absolute scale , and , the spinwave sum / difference frequencies do not include any exchange softening effects . finally we show another example , for @xmath289 at @xmath290 , @xmath238 in fig . [ ompmt1q5 ] , where both the difference and sum processes are shown . in this case there are only 3 dominant sequences of difference processes leading to finite size features . although the memory function formalism applied above , strictly speaking , applies to low temperatures , it is interesting to consider what kind of predictions it makes if applied at higher temperatures . it is clear from the above comparisons between mc - sd data and the memory function data that at low temperatures , where the spinwave approximation should be valid , the damping of the modes has been underestimated . this is seen especially in the @xmath271 memory data where all the finite size features , as well as the spinwave peak , are too narrow . thus it is interesting to see whether the spinwave width and shape are better described by this theory at higher temperature , because the theoretical damping rates grow as @xmath117 . the other reason to look at higher temperature is to consider how the background intensity in @xmath1 for frequencies @xmath291 , depends on the temperature . for one thing , the increase in temperature should spread all of the finite - size effect peaks and smooth the spinwave background , but can be expected to raise its intensity as well . for comparison , as @xmath69 crosses the vortex unbinding temperature @xmath11 , previous mc - sd simulations@xcite show that this background grows with @xmath69 , until it dominates over the spinwave peak . the effect has been described by a theory employing a gas of freely moving vortices,@xcite although simulations have given some indication that vortex lifetimes@xcite are rather too short to take this theory literally , and , vortices tend to be pinned by the lattice and also have complex motions not resembling free translation . thus it is interesting to compare the memory approach to the mc - sd data at higher temperatures , where the latter will include effects due to all possible types of excitations , beyond linear spinwaves .    for @xmath37 , the vortex unbinding temperature has been estimated to be @xmath292 . in fig . [ hi - t - sxx ] , we show some comparisons of the continuum ( @xmath293 ) memory and mc - sd ( @xmath238 ) data for @xmath1 , at temperatures @xmath294 for two particular values of @xmath6 . as already seen , the spinwave peak locations in the mc - sd data here also fall lower than the memory function theory prediction . also , there is a significant growth of the low - frequency background ( @xmath291 ) in the mc - sd data and in the memory function calculations . for @xmath295 the low - frequency background already dominates over the spinwave peak , if the latter can be seen at all . the effect is stronger in the mc - sd data : there is a clear excess of intensity around zero - frequency in @xmath1 in the md - sd data compared to the memory function data . this effect is absent at @xmath296 , however , where the mc - sd and memory data coincide for low frequencies ; in this low temperature region the spinwave approximations used in the memory approach are valid . for temperatures near and above @xmath11 , this extra central peak intensity , not described within the spinwave picture of the memory function approach , may possibly be attributed to nonlinear excitations , including specifically , vortices . this is consistent with the phenomenological vortex gas description of mertens et al:@xcite the excess central peak component is stronger for lower wavevectors , because the vortex pairs being generated above @xmath11 disorder the spins on long length scales . the presence of spinwave damping included into the memory function approach has been suggested here as a way in which finite amplitude spinwave peaks result in the 2d xy model , canceling the divergences in previous theories . a rough calculation of the damping strength has been made using the feynmann diagram technique applied to the quartic corrections of the harmonic hamiltonian . the diagrams used involve only two scattering events ; clearly this gives an underestimate of the damping , and it is likely that including the third order diagrams could make a substantial modification to the results . generally the damping @xmath115 obtained increases away from the center of the brilluoin zone , being close to zero near the center and being greatest along the ( 11 ) direction . in this classical calculation , the damping varies as the square of the temperature ; this results in extremely weak damping strength ( @xmath297 ) for temperatures below @xmath298 . as the temperature approaches and passes @xmath11 , the damping becomes comparable to the spinwave frequency . the comparison of the memory function data and mc - sd data for @xmath1 is good in several aspects : the theory gives finite spinwave peaks , they are positioned correctly for low temperatures , and the multiple finite - size - effect peaks are reproduced as in the mc - sd data . the most significant problem in the memory function approach is the extremely narrow peak widths  clearly caused by an underestimate of the true damping strength , especially at lower temperatures and wavevectors . for very low temperatures in the present theory , the longest wavelength spinwaves will not scatter even once before crossing the entire system . of course , this is a quantum picture , applied to the classical system , so perhaps a better approach in this low - t limit would be to consider some kind of adiabatic scattering in the nonuniform potential due to the spinwave population . this would essentially be the equivalent of the quantum calculation at very high order diagrams rather than at the lowest order . when the memory calculation is carried out for large systems , that is , in the continuum limit , the finite - size peaks merge into a smooth background . indeed , in this limit , the spinwave peak itself has contributions on the high frequency side from spinwave sum processes , and on the low frequency side from spinwave difference processes . the combination of all of these contributions , essentially an infinite number of sum and difference processes , leads to an overall smooth and nondivergent spinwave peak . even in the limit of _ zero damping _ , the spinwave peak are nondivergent , according to continuum calculations [ using eqs . [ sigma_int ] ] by pereira et al.@xcite however , when damping is included , the heights of the spinwave peaks are smaller , and we expect the results to be more realistic .    in the memory function theory , the spinwave peak height was seen to be proportional to @xmath269^{-1}$ ] . the function @xmath299 , to a leading approximation , is proportional to the temperature @xmath69 . from eq . ( [ sqw - memory ] ) , in regions away from the spinwave peak @xmath1 should be proportional to @xmath299 . thus we see that the memory function theory demonstrates how the background grows with temperature while the spinwave peak relatively tends to fall . for high enough temperature , the distribution of spinwave difference processes near zero frequency dominates the spectrum of @xmath1 , leading to a central peak stronger than the spinwave peak . clearly this means applying the theory well above the limiting temperature where we expect spinwave approximations to be valid . however , this gives results very similar to those in the mc - sd simulations , although the latter can be expected to include nonlinear effects such as vortices . in fact , we have seen a strong indication for nonlinear effects ( beyond spinwave difference processes ) above @xmath11 as extra intensity in @xmath1 near zero frequency in the mc - sd data .    in conclusion , application of the memory function formalism , including a leading approximation to the mode damping , has been shown to accomplish three things for the xy model : i ) produce nondivergent spinwave peaks ; ii ) explain the finite size features in @xmath1 both above and below @xmath5 ; iii ) explain the background in @xmath1 below @xmath5 for the infinite system . note that strictly speaking , damping is necessary only for item ii ) , because @xmath1 can not be calculated in the memory function formalism for a finite system without damping . the continuum memory function approach gives finite size spinwave peaks and a low - frequency background even in the absence of damping . the theory also suggests that a sizable contribution to the growth of the central peak for @xmath69 above @xmath11 may be due to spinwave difference processes . the excess central peak intensity in the mc - sd data , when compared to the memory function data , must be due to nonlinear excitations , and gives evidence for dynamical effects due to vortices . gmw thanks universidade federal de minas gerais for their hospitality and support on this project . the authors thank the support by nsf - int-9502781 , nsf cda-9724289 , fapemig , and cnpq . permanent address : department of physics , kansas state university , manhattan , ks 66506 - 2601 . [ email : _wysin@phys.ksu.edu_ ] . j.  villain , j.  phys . ( paris ) * 35 * , 27 ( 1974 ) . f.  moussa and j.  villain , j.  phys . c * 9 * , 4433 ( 1976 ) . d.  r.  nelson and d.  s.  fisher , phys . b * 16 * , 4945 ( 1977 ) . v. l. berezinskii , sov . jetp * 32 * , 493 ( 1970 ) ; * 34 * , 610 ( 1972 ) . j. m. kosterlitz and d. j. thouless , j. phys . c*6 * , 1181 ( 1973 ) r.  d.  mattuck , _ a guide to feynmann diagrams in the many - body problem _ , mcgraw - hill international book company ( 1976 ) . a.  l.  fetter and j.  d.  walecka , _ quantum theory of many - particle systems _ , mcgraw - hill , new york ( 1971 ) . a.  a.  abrikosov , l.  p.  gorkov and i.  e.  dzyaloshinski , _ quantum field theoretical methods in statistical physics _ , pergamon press , oxford ( 1965 ) . g. gilat and j. raubenheimer , phys .  rev . * 144 * , 390 ( 1966 ) . g. gilat and z.  kam , phys . rev . lett . * 22 * , 715 ( 1969 ) . g. gilat , j.  comp . * 10 * , 432 ( 1972 ) . in two - dimensions , ( [ dcell ] ) becomes @xmath300 and in three dimensions , @xmath301 for arbitrary dimension @xmath118 , the generalization is @xmath302 g.  reiter and a sjlander , phys . lett .  * 39 * , 1047 ( 1977 ) ; j.  phys . c : solid state phys . * 13 * 3027 ( 1980 ) . h.  mori , prog . * 33 * , 423 ( 1965 ) . b.  de raedt , h.  de raedt and j.  fivez , phys . b * 23 * , 4597 ( 1981 ) . m.  e.  gouva and a.  s.  t.  pires , j.  phys . c : solid state phys .  * 20 * , 2431 ( 1987 ) . s.  menezes , a.  s.  t.  pires and m.  e.  gouva , phys . b * 47 * , 12280 ( 1993 ) . a.  f.  a.  pereira , m.  s.   thesis , universidade federal de minas gerais , brazil ( 1999 ) . a.  f.  a.  pereira , m.  e.  gouva and g.  m.  wysin , unpublished . the usual conversion from discrete to continuum fourier weights is to divide by the volume @xmath303 per allowed point in @xmath6-space : @xmath304 the @xmath305 normalization of the fourier transform definition ( [ fourier ] ) simplifies many formulas in this paper , including the @xmath50-independent discrete correlation ( [ sqperp ] ) . distributing this weight over the @xmath6-space volume @xmath306 , leads to an @xmath50-_dependent _ continuum correlation @xmath307 . to avoid this @xmath50-dependence in the continuum quantity , it is only necessary to use the alternative fourier transform definitions , @xmath308 this would complicate many intermediate formulas , but in the end the continuum @xmath1 will be independent of @xmath50 and also have a prefactor of @xmath309 . g.  m.  wysin , unpublished . m.  e.  gouva , unpublished . m.  e.  gouva and g.  m.  wysin , phys . b * 56 * , 14192 ( 1997 ) . c. kawabata , m. takeuchi and a.r . bishop , j. magn . mater . * 54 - 57 * , 871 ( 1986 ) ; j. stat . phys . * 43 * , 869 ( 1986 ) ; g.m . wysin and a.r . bishop , phys . b * 42 * , 810 ( 1990 ) ; u. wolff , nucl . phys . * b300 * , 501 ( 1988 ) ; phys . lett . * 62 * , 361 ( 1989 ) . evertz and d.p . landau , phys . b * 54 * , 12302 ( 1996 ) . d.  p.  landau , j.  mag . mater .  * 200 * , 231 ( 1999 ) ; m.  krech and d.  p.  landau , phys . b * 60 * , 3375 ( 1999 ) . m.  e.   gouva , g.  m.   wysin , s.   a.   leonel , a.  s.  t.   pires , t.   kamppeter and f.  g.   mertens , phys . b * 59 * , 6229 ( 1999 ) . mertens , a.r . bishop , g.m . wysin and c. kawabata , phys . b * 39 * , 591 ( 1989 ) . d.a . dimitrov and g.m . wysin , phys . b * 53 * , 8539 ( 1996 ) ; j. phys . c * 10 * , 7453 - 65 ( 1998 ) .
Abstract:

Answer: the effect of damping of spinwaves in a two - dimensional classical ferromagnetic xy model is considered . the damping rate @xmath0 is calculated using the leading diagrams due to the quartic - order deviations from the harmonic spin hamiltonian . the resulting four - dimensional integrals are evaluated by extending the techniques developed by gilat and others for spectral density types of integrals . @xmath0 is included into the memory function formalism due to reiter and solander , and menezes , to determine the dynamic structure function @xmath1 . for the infinite sized system , the memory function approach is found to give non - divergent spinwave peaks , and a smooth nonzero background intensity ( `` plateau '' or distributed intensity ) for the whole range of frequencies below the spinwave peak . the background amplitude relative to the spinwave peak rises with temperature , and eventually becomes higher than the spinwave peak , where it appears as a central peak . for finite - sized systems , there are multiple sequences of weak peaks on both sides of the spinwave peaks whose number and positions depend on the system size and wavevector in integer units of @xmath2 . these dynamical finite size effects are explained in the memory function analysis as due to either spinwave difference processes below the spinwave peak or sum processes above the spinwave peak . these features are also found in classical monte carlo  spin - dynamics simulations . = 1 = 1 = 250 = 350