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: it is well known that the equilibrium momentum distribution of free noninteracting fermions is described by the fermi - dirac distribution function . when interactions are turned on , the momentum distribution is changed due to these interactions . similarly , if the noninteracting gas is confined in a harmonic trap , the distribution can also be calculated explicitly @xcite , however in the limit of a large number of fermions it is convenient to use the semiclassical thomas  fermi approximation @xcite or other semiclassical approaches @xcite . the situation becomes much more complicated if the interaction between fermions can not be neglected . then , numerical methods usually need to be applied , especially when the system is in a trap . on the other hand , the momentum distribution function is a quantity that can be directly accessed in experiments with ultracold atomic gases and therefore it is of a particular importance in trying to make contact between theory and experiment . we examine the problem of a mixture of different mass fermionic atoms on an optical lattice at low temperature . since most stable heavy isotopes of alkalies are bosonic , we envision mixtures of li@xmath0 with k@xmath1 , or light li@xmath0 or k@xmath1 mixed with heavy fermionic isotopes of sr or yb ( it turns out that if the heavy particle is a boson with strong enough intraspecies repulsion , and one is at a low enough temperature , then the statistics of the heavy particle does not enter into our model , as it appears effectively like a hardcore object ) . we have already performed numerical calculations on this system and shown that one can see interesting phenomena similar to viscous fingering as one tunes the trap curvature for the light species and moves from a phase separated state with the heavies on the outside to one with the heavies on the inside  @xcite . interesting ordered density - wave patterns appear at low temperature in these mixtures . the question is , how do we detect the presence of such order with current experimentally available techniques ? in situations where _ in situ _ imaging with single - site precision is available  @xcite , one would simply look for the ordered phases directly , just as they appear in the monte carlo snapshots of a particular configuration of the atoms . but can one see effects of the density - wave ordering in a time - of - flight image , or via direct bragg scattering of light off of the density - wave pattern ? we answer these questions here .    in section ii , we introduce the model and the techniques used to solve for the equilibrium properties of the mixtures of fermions in a harmonic trap . in section iii , we present our numerical results . conclusions follow in section iv . a mixture of light and heavy fermionic atoms in a harmonic trap ( each in one and only one hyperfine atomic state and hence acting like a spinless fermion ) , under the assumption that the quantum - mechanical effects of the hopping of the heavy atoms can be neglected , is described by the falicov  kimball hamiltonian @xcite @xmath2 where @xmath3 is an operator that creates a light fermionic atom at site @xmath4 and @xmath5 is the on  site interspecies interaction potential . the symbol @xmath6 is the light particle trap potential , and @xmath7 is its chemical potential . @xmath8 is the trap for the heavy atoms and @xmath9 is its respective chemical potential . the symbol @xmath10 or 1 is the number operator of the heavy particles , which can be treated as a classical variable since the heavy particles do not hop . the same hamiltonian can also describe fermi - bose mixtures , provided the bosonic atoms are the heavy ( localized ) ones and a strong repulsion between them prevents them from occupying one site by two bosons ( hard core bosons ) . this model can easily be extended to describe a system with soft - core bosons . in such a case @xmath11 and an additional term describing the boson - boson on - site interaction @xmath12 has to be added to the hamiltonian in eq . ( [ ham ] )  @xcite . the trap potentials for the light ( @xmath6 ) and heavy ( @xmath13 ) atoms are given by @xmath14 , v_i=\frac{j}{r^2}\left(x_i^2+y_i^2\right),\ \ \ v^h_i=\frac{j}{\left(r^h\right)^2}\left(x_i^2+y_i^2\right ) \label{trap}\ ] ] where @xmath15 is the position of site @xmath4 ( @xmath16 is the lattice constant ) . the steepness of the potential confining the light atoms varies from @xmath17 to @xmath18 , whereas it is fixed for the heavy atoms with @xmath19 . the chemical potentials @xmath7 and @xmath9 have been introduced in order to control the number of the light and heavy atoms , respectively . the model is solved by means of a variation of the monte carlo ( mc ) method . the method is based on the classical metropolis algorithm modified in such a way that systems with both quantum and classical degrees of freedom can be simulated @xcite . in each mc step , a new configuration of the heavy atoms is generated . then , the hamiltonian in eq . ( [ ham ] ) is numerically diagonalized to yield all the eigenenergies and eigenstates of the light atoms for the trial configuration of the heavy atoms . the new configuration of the heavy atoms is accepted according to the same rules as in the metropolis algorithm , but with the free energy of the light atoms used for comparing energies instead of the internal energy . the results are then averaged over all the configurations generated during the entire mc run . real - space configurations of the atoms for different model parameters ( interaction , shape of the trapping potential ) at different temperatures have been studied in ref . . in the present paper we use a similar method to investigate the momentum distribution and bragg scattering spectra . since diagonalization of the hamiltonian in eq . ( [ ham ] ) for a given configuration of the heavy atoms gives all the eigenstates of the light atoms , it allows one to also calculate the momentum distribution of the light atoms . under the assumption that the optical lattice potential is deep enough that we can restrict to a single - band model , then the field operator of these atoms @xmath20 ( expanded in terms of the bloch wavefunctions for the lowest band ) is given by @xmath21 where @xmath22 and @xmath23 is the bloch wavefunction . expanding the bloch wavefunction in terms of the wannier wavefunctions of the lowest band @xmath24 which are localized about lattice site @xmath25 yields @xmath26 . the summation runs over the first brilliouin zone . then , the free space momentum distribution can be calculated as @xcite @xmath27 where @xmath28 is the fourier transform of the wannier function . hence the momentum distribution can be approximated by @xmath29 in the single - band limit , with @xmath30 the number of lattice sites .    here , the quantum mechanical expectation value @xmath31 is calculated for a given configuration of the heavy atoms and therefore @xmath32 has to be averaged over the configurations . we ignore dynamic interactions between heavy and light atoms due to scattering in the course of the expansion in the time - of - flight experiments . this is a reasonable approximation since the density of the clouds is very small once the trap and the optical lattice potentials are dropped . since those potentials are dropped over a finite period of time , there can be some effects as the potentials are lowered , but these tend to be smaller in these systems because the repulsive interspecies interaction keeps the light and heavy particles away from each other in their initial distributions on the lattice prior to the time - of - flight experiment . since in our approximation the heavy atoms are localized , we are interested in the momentum distribution of only the light ones . the heavy atoms , however , may display some kind of density - wave ordering , which has been demonstrated in ref . . in particular , for some regimes of parameters they form a checkerboard pattern where the heavy atoms occupy only , _ e.g. _ , black squares . in another regime , superpositions of vertical and horizontal stripes or phase separation , with the heavy and light atoms occupying different regions of space have been observed . of course , for any nonvanishing interaction between both species of atoms the distribution of the light atoms is up to some degree correlated with the distribution of the heavy ones . as a result , also the light atoms are ordered , but unless the interaction is strong enough , the magnitude of the density - wave order is much smaller . this follows from the fact that the light atoms , in contrast to the heavy ones , gain kinetic energy while delocalized . the ordering of the heavy atoms can be analyzed by use of scattering of light ( bragg scattering , see _ e.g. _ , ref . ) . this kind of experiment is an equivalent to neutron or x - ray diffraction for the solid state . however , due to the difference of lattice constants between solid state crystals and optical lattices , the required wavelength corresponds to that of visible light . the observation of well - defined bragg peaks has been used to confirm a crystalline structure formed by atoms in an optical lattice @xcite . since this method allows one to determine the distance between atoms forming the crystalline structure , it can be applied in order to detect the checkerboard pattern as well . in the case of strong repulsive interaction between the light and heavy atoms , checkerboard squares of different color are occupied by different kinds of atoms and the problem is similar to the detection of antiferromagnetic ( af ) order . it has recently been proposed to use bragg diffraction of light to detect such an order @xcite in a hubbard system . in the case of af order , the spin  dependent scattering is achieved by using the probe light frequency near atomic resonance , where the interaction between light and atoms is neither purely diffractive nor purely absorptive . in a similar way the probe light can be tuned to be scattered in a different way by the light and heavy atoms . as a result , we can expect a @xmath33 peak in the case of a two - dimensional checkerboard pattern with the light scattered by the heavy atoms . moreover , the intensity of this peak can give information about the fraction of the system occupied by ordered atoms . this method can also be used to detect other types of correlations . if the `` labyrinthine '' patterns obtained in monte carlo simulations @xcite are superpositions or mixtures of horizontal and vertical stripes , the bragg scattering in this case should reveal @xmath34 and @xmath35 peaks . the bragg spectra integrated over the frequency gives the static structure factor @xmath36 , a key quantity that can be expressed as @xmath37 in a practical realization we get a new configuration @xmath38 in each mc step and in each step @xmath36 is calculated . then , @xmath36 is averaged over the entire mc run . simulations have been carried out on a @xmath39 square lattice with 625 heavy and 625 light atoms . in contrast to solid state simulations , there is no translational symmetry in the real system and therefore we use hard - wall boundary conditions at the edges . for steep harmonic trap potentials ( and at low @xmath40 ) most of the atoms are in the center of the system , however , if the potentials are shallow the results may be affected by the finite size of the cluster . [ fig1 ] illustrates the real - space distributions of the light and heavy atoms for different shapes of the trap for the light atoms [ see eq . ( [ trap ] ) for the trap potential parametrization ] . the trapping potential for the heavy atoms is kept at @xmath19 for all simulations .    from fig . [ fig1 ] one can see that if the light atoms trap potential is sufficiently steep and the repulsion between the two species of atoms is relatively strong ( @xmath41 ) , the heavy atoms are pushed out from the center of the trap . this is an intuitive result . however , when the trap potential for the light atoms becomes shallower they start to spread and occupy the periphery of the cluster even before the trap curvatures are set to be equal . for identical potentials , the heavy atoms are concentrated in the central part of the trap , surrounded by the light atoms . in such a configuration the light atoms can gain kinetic energy at the expense of the potential energy in the harmonic trap . [ fig2 ] shows the momentum distribution of the light atoms for the same parameters as in fig . since the strong repulsion leads to phase separation , at least in the limiting cases , the light atoms are able to move freely within the region not occupied by the heavy atoms . nevertheless , one can notice that there is no sign of the fermi surface in the momentum distribution , which results from the inhomogeneity of the system . the upper row of panels presents a comparison of the crossections of the actual distribution with that of a noninteracting fermionic gas in a harmonic trap . the difference results from the interaction between the light and heavy atoms . generally , for all the shapes of the trap potential , the maximal value of the momentum distribution function at @xmath42 of the light atoms is reduced with respect to the noninteracting case . for a relatively steep trap potential for the light atoms , the interaction leads to an occurrence of additional features in the momentum distribution [ see the red ( solid ) curve in the upper left panel of fig . [ fig2 ] ] . the occurrence of these additional inflection points results from a further confinement of the light atoms in an area surrounded by the heavy atoms , as can be seen in fig . [ xxx ] . ( color online ) comparison of the real - space distribution of light atoms with ( left graph , @xmath41 ) and without ( right graph ) interaction with heavy atoms ( @xmath18).,scaledwidth=50.0% ]    in order to confirm the origin of these features , we calculated the momentum distribution of the light atoms in an infinite round well with a diameter equal to the inner diameter of the ring formed by the heavy atoms . in the resulting distribution , presented in the inset in fig . [ fig2 ] , the additional peak in the center of the distribution is even more pronounced . independent of the strength of the interaction , the inhomogeneity destroys any signs of the fermi edge in the momentum distribution . when the trap becomes shallower , behavior that looks like a fermi edge appears , but this is due to the box boundary conditions at the edge of the lattice . even though there is no explicit signal in the momentum distribution function which shows the presence of different kinds of density wave ordering , the momentum distribution function does have a strong dependence on the appearance of phase separation , as can be seen in the different distribution functions in fig.[fig2 ] . for example , it is known that in a homogeneous lattice , the momentum distribution function must decrease below the noninteracting value for small momentum and increase above the noninteracting value for large momentum  @xcite when the system phase separates at low temperature . this behavior is clearly seen in all of the data , and most likely is arising from different forms of quantum confinement effects associated with the phase separation . unfortunately , it is not easy to disentangle this phase separation effect from the effect of the trap , which has a similar effect on the momentum distribution function , except in the case where the heavy particles surround the light ones and confine them with a sharp boundary . in that case , the phase separation effect causes a `` dimple '' in the momentum distribution function near @xmath43 . if the interspecies repulsion is too weak to lead to phase separation , a checkerboard configuration may be formed . [ fig3 ] shows examples of real - space configurations of the light atoms for @xmath44 . it turns out that in this case the momentum distribution is hardly affected by the interaction . [ fig4 ] presents a comparison of the momentum distribution at temperatures above and below the temperature , at which a checkerboard pattern is formed . when the temperature is lowered , the momentum distribution is enhanced for small @xmath45 , though the difference much smaller . moreover , formation of regions with density wave order does not affect the distribution in a significant way . \a )   and @xmath46 . the upper panel ( a ) corresponds to the trapping potential with @xmath18 and the lower ( b ) to @xmath47.,title="fig:",scaledwidth=30.0% ] + b )   and @xmath46 . the upper panel ( a ) corresponds to the trapping potential with @xmath18 and the lower ( b ) to @xmath47.,title="fig:",scaledwidth=30.0% ]     +   +      since the heavy atoms are localized in the proposed approach , one can not analyze their momentum distribution . instead , we determine structure factors defined by eq . [ bragg ] . in particular , we are interested in how different density - wave - ordered patterns are reflected in the structure factor . it is well known from solid state physics that different orderings produce characteristic bragg spectra . however , in the case of cold atoms the spectra are additionally affected by the confining potential that keeps the atoms inside the trap . neglecting the interspecies interaction that may lead to phase separation or ordering , at low temperature , the heavy atoms occupy the bottom of the trap forming a circular region . fourier transformation of such a configuration consists of a finite - width peak at position @xmath48 . plotted in a region @xmath49 , the peak splits into four quarters which are visible in each corner [ at @xmath50 and @xmath51 . the peak itself simply results from a finite number of heavy atoms , since according to eq . ( [ bragg ] ) , @xmath52 is equal to the average concentration of the heavy atoms . the gathering of atoms in the center of the trap broadens it , as can be seen in fig . [ circle ] .    , title="fig:",scaledwidth=23.0%],title="fig:",scaledwidth=25.0% ]    when temperature or interaction moves the atoms to more peripheral areas , the peak width shrinks , finally taking on a delta - function - like form for a random distribution of atoms . it changes the spectral weight around @xmath48 ( and equivalent points ) , affecting the magnitude of the peaks representing ordered phases . this will be discussed in more detail below . since we are mainly interested in patterns formed by the heavy atoms , we will neglect the @xmath48 peak and focus on the remainder of the spectrum . nevertheless , even in the presence of ordering , this peak is still the most pronounced feature of the spectrum . most of the results have been obtained at finite temperatures by means of the monte carlo method . these results can also be compared with zero - temperature local density approximation ( lda ) results , which we now describe . we construct the lda at @xmath53 from the homogeneous , grand canonical phase diagram , where the ground state phases are given as a function of the chemical potentials of heavy and light atoms  @xcite . the procedure is as follows . for each lattice site , we determine the local chemical potentials by subtracting the trap potential at that lattice site from a trial global chemical potential . the global chemical potential is then adjusted to produce the correct total number of heavy and light atoms in the trap . next , using the local chemical potentials , we map out the homogeneous phase diagram for each site within the trap . two sets of model parameters that give nontrivial configurations have been analyzed : @xmath54 and @xmath55 . in the former case , a checkerboard - type configuration is formed in the center of the trap , where the heavy atoms occupy , let us say , the black squares , and the light atoms are primarily on the white ones . the central part is surrounded by rings of different phases . in the latter case , the center of the trap is occupied by light atoms , while the heavy ones form a relatively thick ring composed mainly of various stripe phases . the spatial distributions have been obtained taking into account all possible periodic phases with unit cells consisting of no more than 4 lattice sites ( in all possible shapes ) . the candidate phases are presented in fig . [ configurations ] . a )   b )   c )   d )   e )   + f )   g )   h )   i )     the zero - temperature lda configurations are presented in the lower panels in figs . [ r12_9]a and [ r18_5]a , where different letters correspond to different phases marked with the same identifying letters as in fig . [ configurations ] . the bragg spectra presented in the upper panels in figs . [ r12_9]a and [ r18_5]a are calculated in the following way : first , the positions of the @xmath56-function - type peaks are determined from the fourier transforms of the configurations depicted in fig . [ configurations ] ; then the relative spectral weight proportional to the fraction of the trap occupied by the corresponding phase is assigned to the given peaks . finally , the peaks are slightly broadened to make the presentation more clear . this is necessary because the lda does not know about the finite size of the system and hence always displays perfect delta - function peaks . ( color online ) temperature dependence of the bragg spectra in the center of the trap . more precisely , this figure shows @xmath57 where the sum is calculated over the square marked by white dotted line in fig . [ r12_9 ] , scaledwidth=40.0% ]    it can be seen in fig . [ r12_9]a , that for @xmath18 and @xmath44 a checkerboard phase occupies the central part of the trap . it leads to a highly pronounced peak at @xmath58 . other phases give smaller peaks at @xmath59 , @xmath60 , @xmath61 ( and points obtained by symmetry operations ) .    for @xmath62 and @xmath41 there are no heavy atoms in the center of the trap : they are distributed in a ring of some width with two different phases that have axial stripes on its innermost part . since the stripes are invariant with respect to translations along the axis ( neglecting the finite size of the system ) , all peaks are located at @xmath63 or @xmath64 , namely at @xmath65 and @xmath66 . the outer part of the ring is densely filled with heavy atoms and its fourier transform contributes to the peak at @xmath48 .      in order to investigate how temperature affects the patterns formed by the heavy atoms , monte carlo simulations have been carried out . for a sufficiently long mc run , a number of independent heavy atom configurations are generated . the length of a single mc run depends on the temperature and the model parameters ( interaction , shape of the harmonic trap ) , but usually it is on the order of @xmath67 mc steps from which about @xmath68 independent configurations have been selected . for each configuration the fourier transform has been calculated . [ r12_9]bcd and [ r18_5]bcd present results averaged over all configurations generated for a given set of model parameters . we need to comment about the averaging procedure . in an experiment , the bragg spectra of a single configuration can be observed . however , in many cases the configuration ( and the corresponding spectra ) changes significantly between successive `` snapshots '' . in order to make our results independent of any particular distribution of the atoms , we decided to present results that can characterize the system at a given temperature and model parameters . fourier transforms of very similar configurations often have similar shape , but may have different sign , depending on the details , _ e. g. _ , on the phases of the order parameter for a checkerboard phase . therefore , in figs . [ r12_9]bcd and [ r18_5]bcd , we present averaged absolute values of the spectra , given by @xmath69 where @xmath30 and @xmath70 are the number lattice sites and the number of mc `` snapshots '' , respectively and @xmath71 is equal to 0 or 1 , depending on whether site @xmath4 is occupied or empty in the @xmath72-th `` snapshot '' . additionally , the resulting bragg spectra are self averaged using rotational and reflection symmetries of the lattice and the trap . this means that the presented spectra are calculated as averages of @xmath73 , @xmath74 , @xmath75 , @xmath76 , @xmath77 , @xmath78 , @xmath79 and @xmath80 . since we are not interested in the @xmath48 peak ( and equivalent peaks in the remaining corners ) , the false - color scales in figs . [ r12_9 ] and [ r18_5 ] are chosen in such a way that ( independent of the temperature ) only the peaks resulting from ordering are correctly represented . but of course , the scale is kept the same in panels b , c , and d.        ( color online ) temperature dependence of the bragg spectra at the edge of the trap : this figure shows @xmath81 where the sum is calculated over the rectangle marked by red dotted line in fig . [ r18_5 ] , scaledwidth=40.0% ]    in the case presented in fig . [ r12_9]a ( lda , @xmath82 ) , most of the heavy atoms form a checkerboard configuration leading to a strong peak at @xmath58 . in the peripheral areas , the concentration of heavy atoms increases and they form more dense patterns , which in turn , lead to less pronounced peaks . one can notice that the central peak in fig . [ r12_9]b is fourfold split . this results from an imperfection of the checkerboard ordering of the heavy atoms . this can be explained from an analysis of the lda : the effective chemical potential varies continuously when the distance from the center of the trap increases . hence one expects that the density of atoms ( both light and heavy ) should also vary in a continuous way . however , in some regions , the checkerboard ordering minimizes the energy , hence in that ordered region the density ( at least of the heavy atoms ) is constant . a competition between these two tendencies leads to flaws in the center of the trap . due to topological reasons , such a crack has to run across the entire ordered area . then the question is why the fourfold split is not visible in the lda results ? this is connected with the maximum size of the unit cell taken into account in the lda calculations . the @xmath83 site unit cell is too small to describe a large checkerboard area with a single line defect , but by using a larger unit cell , one should see such a configuration .    as can be expected , the ordering is reduced when the temperature increases and so are the corresponding features in the bragg spectra . figs . [ r12_9]c and [ r12_9]d illustrate how the multi - peak structure is smoothed out by destroying the real - space ordering . in order to describe this process in a more quantitative way , we calculated how the spectral weight of the central peak(s ) decreases with increasing temperature . figure [ svt_1 ] shows the spectral weight in a region @xmath84 ( marked by the white square in fig . [ r12_9]a ) as a function of temperature . the plot is shifted in such a way that the weight at @xmath85 is equal to zero . the shift is necessary since with increasing temperature less and less weight is attributed to the peak at @xmath85 , which increases the reference level more than should result from the reduction of the central peak ( this effect is especially visible in fig . [ r18_5 ] , where the color around @xmath86 becomes very bright at high temperature ) . in this paper , we have examined two simple experimental probes that can reveal information about ordered density wave phases in mixtures of light and heavy fermionic atoms ( or equivalently light fermionic and heavy `` hard core '' bosonic atoms ) on an optical lattice . namely , we examined the momentum distribution function , which comes from a time - of - flight expansion experiment and the bragg scattering signal that would come from scattered optical light that scatters off of the density wave pattern .    the momentum distribution function does not provide significant information about various ordered density wave phases , but it does provide some intuition about phase separated states . in particular , as the distribution flattens and broadens , one has an indication of phase separation setting in . furthermore , if the heavy atoms surround the lights and confine the light atoms with essentially a hard wall boundary condition , then the momentum distribution function develops a sharp dimple at low momentum which does provide a characteristic shape signalling that form of phase separation . bragg scattering is much more effective at showing the presence of ordered density wave phases , as new bragg reflection `` spots '' appear at appropriate ordering vectors for the different types of order present in the sample . the weight underneath these peaks is proportional to the strength of the ordering , and to the volume of regions which are ordered , and hence they can be used for accurate low - temperature thermometry of these systems . as @xmath40 is lowered , the weight in the peaks grows and can be calibrated via numerical calculations to produce an appropriate temperature of the system .    one caveat of this work , however , is that one must cool the mixture down to a low enough temperature that the ordering appears in the system . typically , this requires a temperature at least as low as about 1/40th of the bandwidth , and often substantially lower . this is an aggressively low temperature with current cooling technology , but hopefully can be reached as it becomes easier to manipulate entropy distributions within trapped atomic systems . finally , we also should note that direct _ in situ _ imaging via apparatus like the quantum gas microscope would provide even more convincing pictures of the ordered phases , and the fluctuations about that order , than the above proposed methods , but we are not aware of any plans to examine these kinds of mixtures with such machines in the near term . m.m.m . acknowledges support under grant no . nn  202  128736 from ministry of science and higher education ( poland ) . j. k. f. acknowledges support under aro grant number w911nf0710576 with funds from the darpa ole program and support from the mcdevitt endowment trust at georgetown university . we acknowledge m. rigol for a critical reading of the manuscript . 99 m. rigol and a. muramatsu , phys . a * 70 * , 043627 ( 2004 ) . d. a. butts and d. s. rokhsar , phys . rev . * a * 55 , 4346 ( 1997 ) . n. h. marcha , l. m. nietoc , and m. p. tosi , phys . b : condens . matter * 293 * , 308 ( 2001 ) . m. m. maka , r. lemaski , j.k . freericks , and c.j . williams , phys . * 101 * , 060404 ( 2008 ) . w. s. bakr , j.i . gillen , a. peng , s. foelling , m. greiner nature * 462 * , 74 ( 2009 ) . l.  m. falicov and j.  c. kimball , phys . * 22 * , 997 ( 1969 ) . c. ates and k. ziegler , phys . a * 71 * , 063610 ( 2005 ) . k. byczuk and d. vollhardt , phys . b * 77 * , 235106 ( 2008 ) ; ann . phys . * 18 * , 622 ( 2009 ) . m. iskin and j. k. freericks , phys . a * 80 * , 053623 ( 2009 ) . m. m. maka , k. czajka , phys . b * 74 * , 035109 ( 2006 ) ; m. m. maka , k. czajka , phys . . sol . * 242 * , 479 ( 2005 ) . daw - wei wang , phys . a * 80 * , 063620 ( 2009 ) . g. birkl , m. gatzke , i. h. deutsch , s. l. rolston , and w. d. phillips , phys . 75 * , 2823 ( 1995 ) ; matthias weidemller , axel grlitz , theodor w. hnsch , and andreas hemmerich , phys . a * 58 * , 4647 ( 1998 ) ; philipp t. ernst , sren gtze1 , jasper s. krauser , karsten pyka , dirk - sren lhmann , daniela pfannkuche , and klaus sengstock , nature physics * 6 * , 56 ( 2010 ) . t. a. corcovilos , s. k. baur , j. m. hitchcock , e. j. mueller , and r. g. hulet , phys . a * 81 * , 013415 ( 2010 ) . j. k. freericks , e. h. lieb , and d. ueltschi , commun . phys . * 227 * , 243 ( 2002 ) . g. i. watson and r. lemaski , j. phys . condens . matter * 7 * , 9521 ( 1995 ) ; r. lemaski , j. k. freericks and g. banach , j. stat . phys . * 116 * , 699 ( 2004 ) .
Abstract:

Answer: the momentum distribution is one of the most important quantities which provides information about interactions in many - body systems . at the same time it is a quantity that can easily be accessed in experiments on ultracold atoms . in this paper , we consider mixtures of light and heavy fermionic atoms in an optical lattice described effectively by the falicov - kimball model . using a monte carlo method , we study how different ordered density - wave phases can be detected by measurement of the momentum distribution of the light atoms . we also demonstrate that ordered phases can be seen in bragg scattering experiments . our results indicate that the main factor that determines the momentum distribution of the light atoms is the trap confinement . on the other hand , the pattern formed by the heavy atoms seen in the bragg scattering experiments is very sensitive to the temperature and possibly can be used in low - temperature thermometry .


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: in the functional formalism @xcite the general quantum mechanical amplitude @xmath6 is given in terms of a path integral which is simply the @xmath7 limit of the expression @xmath8 the euclidean time interval @xmath9 $ ] has been subdivided into @xmath0 equal time steps of length @xmath10 , with @xmath11 and @xmath12 . @xmath13 is the naively discretized action of the theory . in this paper we will look at theories with action of the form @xmath14 note that we use units in which @xmath15 and particle mass have been set to unity . the naively discretized action is in this case simply @xmath16 where @xmath17 , and @xmath18 . key investigations regarding numerical evaluation of path integrals were presented in the reviews of barker and henderson @xcite , kalos and whitlock @xcite , and ceperley @xcite , as well as in the papers by pollock and ceperley @xcite and barker @xcite . a modern and extensive reference on the subject of path integrals is given in the latest edition of the textbook by kleinert@xcite . as we can see , the very definition of path integrals makes it necessary to make the transition from the continuum to the discretized theory . this discretization , however , is far from unique . in fact , the details of the discretization procedure are extremely important both for analytical and numerical treatment of path integrals . this dependence on discretization procedure has been one of the principle impediments to creating a consistent mathematical theory of path integration . on the numerical side , this manifests itself in the fact that path integral simulations remain notoriously demanding of computing time  so much so that certain path integral calculations serve as benchmarks for new generations of supercomputers . for this reason let us note that we have the freedom to make two important choices that will not affect the final result , i.e. the continuum amplitude we seek to calculate . first , we have the freedom to choose the point in the interval @xmath19 $ ] in which to evaluate the potential @xmath20 . it is well known that different points correspond to different ordering prescriptions in the operator formalism . the choice of the middle point @xmath21 is the most common one . it corresponds to the symmetric or weyl ordering of operators @xmath22 and @xmath23 , and so always leads to a hermitean expression for the hamiltonian @xmath24 . two other prescriptions are also often used . the left ordering prescription evaluates the potential at @xmath25 , the left boundary of the above interval ( in the operator formalism this corresponds to taking the @xmath22 s to the left of the @xmath23 s in all the products that appear in the hamiltonian ) . similarly , one defines the right ordering prescription . although they lead to somewhat simpler looking expressions , the left and right prescriptions do not in general give hermitean hamiltonians . let us note , however , that the class of theories given by eq . ( [ action ] ) has a hamiltonian that is the sum of a @xmath22-dependent kinetic term and a @xmath23-dependant potential term and so has no ordering ambiguities . in this case different prescriptions lead to the same continuum amplitude  the discretized amplitudes do differ , but they tend to the same continuum limit .    the second , and more important freedom related to our choice of discretized action has to do with the freedom to introduce additional terms that explicitly vanish in the continuum limit . we will designate such discrete actions as effective actions . for example , the term @xmath26 where @xmath27 is regular when @xmath28 , does not change the continuum physics since it goes over into @xmath29 , i.e. it vanishes as @xmath30 . although such additional terms do not change the continuum physics they do affect the speed of convergence to that continuum limit . the aim of this paper is to give a detailed exposition of a systematic analysis that leads to the best solution in the class of all equivalent effective actions , e.g. the effective action that leads to the fastest convergence of the associated discrete amplitude @xmath31 to the continuum expression . the uncovered speedup in the path integral algorithm is a direct consequence of new analytical input that has come from the study of the relation between discretizations of differing coarseness of the same theory . we have given a brief presentation of these ideas in a recent paper @xcite . the calculations that will be presented turn out to be simplest in the mid - point prescription . before we proceed with them it is useful to spend the next section in a brief overview of known results dealing with the speed of convergence to the continuum limit . in this section we will compare the speed of convergence of several different prescriptions to the continuum limit . the naive mid - point prescription satisfies @xmath32 for all @xmath33 and @xmath34 . by naive we mean that we use the naively discretized action given in eq . ( [ actionn ] ) . on the other hand , in the naive left prescription the amplitude for @xmath35 converges much faster @xmath36 this behavior can be easily shown both analytically and numerically . we note in passing that this is strongly related to the well known result for the partition function evaluated using naive discretization in the left prescription @xcite , which follows directly from the above amplitude by integrating over @xmath33 ( to get the trace ) and writing the time of propagation @xmath37 as the inverse temperature @xmath38 , @xmath39 however , going back to the language of amplitudes , it is also easy to show that the amplitudes for different initial and final states converge slower , i.e. for @xmath40 we have @xmath41    the problems with the speed of convergence of off - diagonal amplitudes in the left prescription can be fixed very easily . we find that for all @xmath33 and @xmath34 we have @xmath42    although not related to the central investigation in this paper , let us present a proof of eq . ( [ leftright ] ) as an illustration of the use of the trotter formula @xmath43 using it we can easily show the validity of the formal expression @xmath44 . on the other hand , from eq . ( [ amplituden ] ) we find that @xmath45 . as a result we see that the left hand side of eq . ( [ leftright ] ) is simply the average of @xmath31 and @xmath46 and so has an expansion in even powers of @xmath47 . we end this section by presenting the result obtained by takahashi and imada @xcite and independently by li and broughton @xcite . in these papers the authors used a generalized form @xcite of the trotter formula to increase the speed of convergence of the discretized partition function . their final result is a derivation of a formula for the effective potential @xmath48 . the authors showed that by using this effective potential ( in the left prescription ) one gets @xmath49    a recent analysis of this method can be found in jang et al @xcite . let us note that the crucial step in the derivation of the above effective potential from the generalized trotter formula uses the cyclic property of the trace , i.e. the above increase in the speed of convergence only holds for the partition function and not the amplitudes . a direct numerical simulation shows that amplitudes calculated using this effective potential converge just as fast as the amplitudes in the naive left prescription . said another way , it is only the integral over all the diagonal amplitudes that has the @xmath50 behavior and not any individual amplitude . a recent investigation by bond et al . @xcite has uncovered a @xmath51 behavior , however , not for the case of a generic theory . at the end let us mention that several related investigations dealing with speed of convergence have focused on improvements in short time propagation @xcite or the action @xcite . the aim of this and the following section is to present a systematic exposition of the relation between different discretizations of the same path integral . throughout we will work in the mid - point prescription . we start by studying the relation between the @xmath52-fold and @xmath0-fold discretizations of a given amplitude . from eq . ( [ amplituden ] ) we see that we can write the @xmath52-fold amplitude as @xmath53 i.e. in the form of an @xmath0-fold amplitude given in terms of a new action @xmath54 determined by @xmath55 where @xmath56 is nothing but the @xmath52-fold discretization of the starting action . in the above formulas we have , for convenience , written the @xmath52-fold discretized coordinates @xmath57 in terms of @xmath58 s and @xmath59 s in the following way : @xmath60 and @xmath61 . note that we have @xmath11 , @xmath12 , while the @xmath62 remaining @xmath58 s play the role of the dynamical coordinates in the @xmath0-fold discretized theory . the @xmath59 s are the @xmath0 remaining intermediate points that we integrate over in eq . ( [ tildes ] ) . we wish to have @xmath54 belong to the same class of actions as @xmath13 . it is not difficult to show that the naively discretized action does not satisfy this requirement , i.e. the integration of eq . ( [ tildes ] ) will yield new types of terms in @xmath54 . in fact , the class of actions closed to transformation ( [ tildes ] ) is of the form @xmath63 the functions appearing in the above effective action also depend on the time step @xmath64 . we choose not to display this dependence explicitly in order to have a more compact notation . what is important is that all of these functions are regular in the @xmath65 limit . note that these effective actions are equivalent to our starting action , i.e. they all have the same continuum amplitudes as the starting theory . using eq . ( [ tildes ] ) and ( [ actions ] ) one can easily derive the following integral relation which determines the functions @xmath66 in the new action @xmath67 in terms of the related functions in the starting action :    @xmath68    where @xmath69    the above integral equation can be solved for the simple cases of a free particle and a harmonic oscillator , and gives the well known results . let us note that the integral in eq . ( [ integral ] ) is in a form that is ideal for an asymptotic expansion @xcite , whatever the potential . the time step @xmath64 plays the role of a small parameter ( in complete parallel to the role @xmath15 plays in the usual semi - classical , or loop , expansion of quantum theories ) . as is usual , the above asymptotic expansion is carried through by first taylor expanding @xmath70 around @xmath71 and then by doing the remaining gaussian integrals . we find @xmath72 where we have assumed that @xmath73 , i.e. that @xmath74 .    finally , from eq . ( [ integral ] ) and ( [ expansion ] ) we get @xmath75\ .\end{aligned}\ ] ] all that remains is to calculate the @xmath76 s using eq . ( [ f ] ) and to expand the potential and all the functions @xmath77 around the mid - point @xmath21 . for this second step we make use of the simple relations that follow from the definitions of @xmath21 and @xmath78 : @xmath79 , @xmath80 , and @xmath81 . using these relations we can expand a typical term like @xmath82 to obtain @xmath83    for example , the expansion of @xmath67 up to @xmath84 is a rather simple exercise . we find the following functional relations @xmath85+\nonumber\\ & & + \epsilon_n^2\,\left[\frac{3}{16}g_2-\frac{1}{32}v'\,^2+\frac{1}{2048}v^{(4)}+\frac{3}{128}g_1''\right]\nonumber\\ \widetilde g_1 & = & \frac{1}{4}g_1+\frac{1}{32}v''+\\ & & + \epsilon_n\left[\frac{3}{8}g_2+\frac{1}{1024}v^{(4)}-\frac{1}{64}g_1''\right]\nonumber\\ \widetilde g_2 & = & \frac{1}{16}g_2+\frac{1}{6144}v^{(4)}+\frac{1}{128}g_1''\ .\nonumber\end{aligned}\ ] ]    note that in the above relations we expanded @xmath86 up to @xmath30 , @xmath87 up to @xmath64 and @xmath88 up to @xmath89 . we also disregarded all the higher @xmath90 s . the reason for this is that the short time propagation of a generic theory satisfies @xmath91 while the @xmath77 term enters the action multiplied by @xmath92 . in general , if we expand the new action @xmath67 to @xmath93 we need to evaluate only @xmath86 ( up to @xmath94 ) and the first @xmath95 functions @xmath90 ( up to @xmath96 ) . although straight - forward , the task of calculating @xmath67 to large powers of @xmath64 is quite tedious ; using the symbolic algebra package mathematica 5.0 we have analytically solved the corresponding expressions up to @xmath5 . the memory requirements for this calculation grow exponentially with @xmath3 : the @xmath97 calculation used just under 2 gb of computer memory .    at this point it is important to comment on what has been achieved so far . evaluating @xmath67 to @xmath93 and using this new action to calculate the @xmath0-fold discretized amplitude @xmath98 we find @xmath99 so that , up to @xmath100 , this amplitude is the same as the @xmath52-fold amplitude calculated with our starting action . in this way we have halved the discretization from @xmath52 to @xmath0 . therefore , a coarser @xmath0-fold discretization using @xmath54 does the same job as the @xmath52-fold discretization of the starting theory . in the next section we will consider the iteration of this halving procedure . we will derive and solve the recursive relations that connect up the @xmath101-fold and @xmath0-fold discretizations . in particular , we will focus on the continuum limit solution when @xmath102 , i.e. the solution that connects up the continuum theory with its @xmath0-fold discretization . the iterative process of halving that starts from the @xmath101-fold discretization is governed by a recursive relation . from the @xmath103 case given in eq . ( [ halvingrelations ] ) , used in the previous section to illustrate the general discretization halving scheme , we directly get the sought - after @xmath103 level system of recursive relations    @xmath104 + \frac{\epsilon_n^2}{2^{2(s - k-1)}}\,\left[\frac{3}{16}(g_2)_{k}-\frac{1}{32}v_{k}'\,^2+\frac{1}{2048}v_{k}^{(4)}+\frac{3}{128}(g_1)_{k}''\right]\nonumber\\ ( g_1)_{k+1 } & = & \frac{1}{4}(g_1)_{k}+\frac{1}{32}v_{k}''+ \frac{\epsilon_n}{2^{s - k-1}}\left[\frac{3}{8}(g_2)_{k}+\frac{1}{1024}v_{k}^{(4)}-\frac{1}{64}(g_1)_{k}''\right]\\ ( g_2)_{k+1 } & = & \frac{1}{16}(g_2)_{k}+\frac{1}{6144}v_{k}^{(4)}+\frac{1}{128}(g_1)_{k}''\ .\nonumber\end{aligned}\ ] ]    in the above relations @xmath105 . the zeroth iterate corresponds to the starting action , the last iterate to the effective action that gives an equivalent @xmath0-fold discretization . the @xmath106 terms represent the time step of the @xmath107-th iterate in the discretization halving procedure .    although the above system of recursive relations is non - linear it is in fact quite straight - forward to solve if we remember that the system itself was derived via an expansion in @xmath64 . having this in mind we first write all the functions as expansions in powers of @xmath64 that are appropriate to the level @xmath3 we are working at . in this case we are illustrating the procedure for @xmath103 , so we have @xmath108 putting this into the @xmath103 level system of recursive relations given in eq . ( [ rec : p3 ] ) we find that @xmath109 . using the initial condition @xmath110 we find that @xmath111 for all @xmath107 . using this the remaining equations form the following set of linear recursive relations @xmath112 this system is easily solved for given initial conditions . however , what we are really interested in is the continuum limit solution which is obtained by setting @xmath113 in the above expression and taking the limit @xmath102 . by doing this we are iterating our process of discretization halving from the continuum theory down to @xmath0 . the continuum limit solution of the @xmath103 level system is simply @xmath114\nonumber\\ ( g_1)_{p=3 } & = & \frac{v''}{24 } + \epsilon_n\frac{v^{(4)}}{480}\\ ( g_2)_{p=3 } & = & \frac{v^{(4)}}{1920}\ .\nonumber\end{aligned}\ ] ] in the above expressions the label  @xmath103 " reminds us that this is the solution of the continuum limit of the recursive relations given in eq . ( [ rec : p3 ] ) describing discretization halving at the @xmath103 level . note that the continuum limit solution depends only on the initial potential @xmath20 , i.e. it is not sensitive to initial values of the @xmath77 s as these terms all vanish in the continuum limit . in this way we have obtained the effective action that gives the best @xmath0-fold discretization of the starting theory at the @xmath103 level . one can similarly obtain a set of effective actions @xmath1 , one for each value of @xmath3 . the solution for @xmath115 is given in the appendix . note that each solution contains within it all the solutions for lower levels . solutions for larger values of @xmath3 are a bit more cumbersome , however , they are just as easy to use in simulations . expressions up to @xmath97 can be found on our web site @xcite . note that one solves for the continuum limit of the level @xmath3 system of recursive relations but once for all theories , i.e. once the solution is found it works for all sufficiently smooth potentials @xmath20 . actually , the requirement for the level @xmath3 solution is that the starting potential is differentiable @xmath116 times .    the effective action satisfying the continuum limit of the discretization halving recursion relations at level @xmath3 leads to an @xmath0-fold amplitude that is equal to the continuum amplitude of the starting action up to an @xmath100 term therefore , the continuum limit solution satisfies @xmath117    expectation values can be calculated with the same precision using standard discretized estimators , provided that the discretized time step on which those observables reside is chosen appropriately . for example , the expectation value of the momentum squared @xmath118 may be calculated using the standard estimator @xmath119 if the time step from @xmath120 to @xmath121 is shortened to @xmath93 keeping all the remaining time steps unchanged . the validity of the presented analytical result will be illustrated in the following section where we present monte carlo simulations of two different models . to conclude this section : we have constructed a general procedure for calculating effective actions @xmath1 for any level @xmath3 . we have completed the procedure and found explicit values for the effective actions up to and including @xmath97 . the @xmath0-fold amplitudes of the @xmath97 effective action differ from the continuum amplitudes by a term proportional to @xmath122 . in this section we illustrate the generic results obtained in the previous sections by analyzing the speedup for different values of @xmath3 in the case of monte carlo simulations of two different models . the first model we looked at is the anharmonic oscillator with quartic coupling @xmath123 in fig . [ phi4 ] we illustrate how the discretized amplitudes @xmath124 tend to the continuum limit for levels @xmath125 ( naively discretized starting action ) , @xmath126 and @xmath127 . the top plot gives an overall view of how the discretized amplitudes for @xmath128 , @xmath129 and @xmath130 calculated using higher level effective actions systematically outperform the ones at lower level . in particular we see how they outperform the amplitude calculated using the naively discretized action . in the bottom plot we show a detail of the top plot which makes it easier to qualitatively track the differences between amplitudes calculated with effective actions at levels @xmath131 and @xmath127 . in agreement with eq . ( [ convergence ] ) , the curve fitted to the @xmath3 level data is a polynomial in @xmath47 of the form @xmath132 as derived , all the @xmath133 are ( within error bars ) equal to each other and represent the continuum amplitude @xmath134 we seek . the continuum value is represented in the plots by a dashed line . in all cases the fits were done for data with @xmath135 . the reason that the @xmath136 points were omitted is that for @xmath130 we have @xmath137 , i.e. these points do not satisfy the condition for asymptotic expansion @xmath73 . however , the @xmath136 amplitudes are quite interesting because they are algebraic expressions with no integrations . from the above plot we see that the @xmath136 amplitudes calculated with higher level effective actions give better and better approximations to the amplitude . we have seen this behavior for other potentials as well . ( color online ) ( top ) plot of discrete amplitudes @xmath138 as a function of @xmath0 for @xmath139 and @xmath127 for an anharmonic oscillator with quartic coupling @xmath140 , time of propagation @xmath130 from @xmath128 to @xmath129 , @xmath141 . ( bottom ) detail of the same plot comparing amplitudes for @xmath131 and @xmath127 . in both plots the dashed line represents the continuum limit amplitude.,title="fig:",width=321 ]   ( color online ) ( top ) plot of discrete amplitudes @xmath138 as a function of @xmath0 for @xmath139 and @xmath127 for an anharmonic oscillator with quartic coupling @xmath140 , time of propagation @xmath130 from @xmath128 to @xmath129 , @xmath141 . ( bottom ) detail of the same plot comparing amplitudes for @xmath131 and @xmath127 . in both plots the dashed line represents the continuum limit amplitude.,title="fig:",width=321 ]    for @xmath142 the @xmath136 amplitudes calculated using effective actions derived in the previous section represent very good algebraic approximations for the case of a general theory  larger levels @xmath3 give better approximations . from fig . [ phi4 ] we see that these approximations work rather well even for the marginal point @xmath130 . for @xmath143 we do need to do some integrals numerically in order to get the required amplitudes . for high level effective actions it is enough to use @xmath144 + 1 $ ] , i.e. to do only @xmath145 $ ] integrations numerically . a quantitative measure of how well the derived effective actions perform can be seen in fig . [ p1246 ] . we see explicitly that the @xmath3 level data differ from the continuum amplitudes as polynomials starting with @xmath4 . ( color online ) the deviations from the continuum limit @xmath146 as a function of @xmath0 for @xmath139 and @xmath147 ( top to bottom ) . this particular plot is for the case of an anharmonic oscillator with quartic coupling @xmath140 , time of propagation @xmath130 from @xmath128 to @xmath129 . the number of monte carlo samples used was @xmath148 for @xmath149 , @xmath150 for @xmath151 , and @xmath152 for @xmath115 . dashed lines correspond to appropriate @xmath47 polynomial fits to the data . the solid lines give the leading @xmath47 behavior . the level @xmath3 curve has an @xmath4 leading behavior.,width=321 ]    because of this , the deviations from the continuum limit @xmath146 become exceedingly small for larger values of @xmath3 making it necessary to use ever larger values of @xmath153 so that the monte carlo statistical error does not mask these extremely small deviations . for @xmath115 we see that although we used an extremely large number of monte carlo samples ( @xmath152 ) the statistical errors become of the same order as the deviations already at @xmath154 . for @xmath97 this is the case even for @xmath155 , i.e. we already get the continuum limit within a monte carlo error of around @xmath156 .    to make the deviations in fig . [ p1246 ] visible for large @xmath3 levels we needed to run a simulation with very large @xmath153 . this simulation took about a week on our 160 gflops cluster . on the other hand the simulation in fig . [ phi4 ] uses a much smaller @xmath153 and takes less than one hour to complete . in practice we see that the derived effective actions give excellent agreement with continuum limit amplitudes already for small values of @xmath0 . simulations with such values of @xmath0 take a negligible amount of time even on a single pc . the second model we consider is that of a particle moving in a modified pschl - teller potential  a well known exactly solvable model @xcite @xmath157 unlike the anharmonic oscillator the pschl - teller potential has both a continuous and discrete spectrum . the discrete eigenstates have energy @xmath158 for @xmath159 . therefore , we see that the model has critical values of coupling @xmath160 at which it acquires new bound states . ( color online ) ( top ) plot of discrete amplitudes @xmath138 as a function of @xmath0 for @xmath139 and @xmath127 for a particle in a modified pschl - teller potential with parameters @xmath161 , @xmath162 . @xmath130 , @xmath128 , @xmath129 , @xmath141 . ( bottom ) detail of the same plot comparing amplitudes for @xmath131 and @xmath127 . in both plots the dashed line represents the continuum limit amplitude . , title="fig:",width=321 ]   ( color online ) ( top ) plot of discrete amplitudes @xmath138 as a function of @xmath0 for @xmath139 and @xmath127 for a particle in a modified pschl - teller potential with parameters @xmath161 , @xmath162 . @xmath130 , @xmath128 , @xmath129 , @xmath141 . ( bottom ) detail of the same plot comparing amplitudes for @xmath131 and @xmath127 . in both plots the dashed line represents the continuum limit amplitude . , title="fig:",width=321 ]     ( color online ) same plots as in fig . [ mpt1 ] but for a modified pschl - teller potential with parameters @xmath161 and @xmath163.,title="fig:",width=321 ]   ( color online ) same plots as in fig . [ mpt1 ] but for a modified pschl - teller potential with parameters @xmath161 and @xmath163.,title="fig:",width=321 ]    fig . [ mpt1 ] displays how the discretized amplitudes @xmath124 tend to the continuum limit for levels @xmath164 for the potential with @xmath161 and @xmath162 . the same plots for the case of @xmath161 and @xmath163 ( lying on the a critical value of @xmath38 ) are given in fig . as we can see the effective actions work just as well as in the case of the anharmonic oscillator . going through a critical point like @xmath163 certainly affects the physical quantities calculated , however , the speedup algorithm is not affected in any way .    with the increase of @xmath3 level the complexity of the expressions for the effective actions grows exponentially . therefore , the increase in computation time that results from using higher @xmath3 level effective actions also grows exponentially as is shown in fig . [ time ] . ( color online ) relative increase in computation time that comes about from the increased complexity of expression for higher @xmath3 level effective actions.,width=321 ]    as we have seen , by increasing @xmath3 we drastically improve convergence to the continuum limit . an important consequence of this is that we can obtain the same precission using much smaller values of @xmath0 , i.e. much coarser discretizations . this is at the root of the speedup that we find . however , as we have seen , the exponential growth in complexity of the effective actions puts an upper bound to @xmath3 levels that can be used . from fig . [ time ] we see that @xmath97 is still far from that upper bound  the gain of eight orders of magnitude in the speed of convergence far outweighs what is roughly a ten fold increase in computation time .    at the end we briefly comment on two monte carlo algorithms developed for simulations in this section . in the first algorithm trajectories are generated by a gaussian distribution function obtained using a semi - classical expansion . the computing time of this algorithm scales as @xmath165 since it is necessary to diagonalize the quadratic form in the exponential of the distribution function . in the second algorithm we implemented the bisection method @xcite , which scales as @xmath166 . therefore , the bisection algorithm is the method of choice for large values of @xmath0 . on the other hand , our method allows us to obtain very precise results using small values of @xmath0 . in that region we have found the two algorithms to be comparable both in precision and running time .    in both algorithms we needed to use a random number generator which gives a large number of uncorrelated random numbers in a fashion suitable for parallel programming . our primary random number generator was the scalable parallel random number generator library @xcite ( sprng ) . following the good practise suggested by ferrenberg et al . @xcite we have checked all our results using a different random number generator . checks were made with the numerical recipes ran3 generator @xcite with a different seed for each mpi process . agreement was in all cases well within a 1-@xmath167 interval implying that there were no hidden systematic errors present in either the algorithms or the random number generators . we note in passing that the analytical derivations presented in this paper work equally well in both the euclidean and minkowski formalism ( with appropriate @xmath168 regularization ) , i.e. they are directly applicable to quantum systems as well as to statistical ones . however , the monte carlo simulations used to numerically document our analytical results necessarily needed to be done in the euclidean formalism . to conclude , we have presented an algorithm that leads to significant speedup of numerical procedures for the calculating of path integrals of a generic theory . the increase in speed is the result of new analytical input that has emerged from a systematic investigation of the relation between different discretizations of the same theory . we have presented an explicit procedure for obtaining a set of effective actions @xmath1 that have the same continuum limit as the starting action @xmath169 , but which approach that limit ever faster . amplitudes calculated using the @xmath0-point discretized effective action @xmath170 satisfy @xmath171 , where @xmath33 and @xmath34 are initial and final states , @xmath37 the time of propagation , and @xmath134 the sought - after amplitude of the continuum theory . we have obtained and analyzed the effective actions for @xmath5 . in this paper we quote expressions up to @xmath115 ( see the appendix ) , the rest can be found on our web site @xcite .    at the end we illustrated the obtained generic results by analyzing the speedup for different values of @xmath3 in the case of concrete monte carlo simulations of two different models : anharmonic oscillator with quartic coupling and particle in a modified pschl - teller potential . extensions of the derived algorithm to @xmath172 particles and @xmath173 dimensions , as well as to quantum field theories are both in progress . in both cases the derivation of the analogue of integral eq . ( [ integral ] ) does not seem to present a problem . the asymptotic expansion used to solve it is also directly generalizable . however , the algebraic recursive relations that determine @xmath1 will be more complex and may practically limit us to smaller values of @xmath3 . in this appendix we present the effective action at level @xmath115 . note that this solution contains within it the effective actions at all lower levels  all one needs to do is to truncate the @xmath115 solution at the appropriate order in the @xmath64 expansion . for example , the effective potential @xmath20 at the @xmath103 level is obtained from the @xmath115 level potential by disregarding term that are @xmath84 and higher . code containing @xmath1 for @xmath5 is available on our web site @xcite . @xmath174 + \epsilon_n^3\left[-\frac{v''\,^2}{360 } -\frac{v'\,v^{(3)}}{120 } +   \frac{v^{(6)}}{6720}\right]+\nonumber\\ & & + \,\epsilon_n^4\left[\frac{v'\,^2\,v''}{240 } -     \frac{23\,{v^{(3)}}^2}{40320 } -     \frac{v''\,v^{(4)}}{1680 } -     \frac{v'\,v^{(5)}}{2240 } + \frac{v^{(8)}}{241920}\right]+\nonumber\\ & & + \,\epsilon_n^5\left[\frac{v''\,^3}{5670 } + \frac{29\,v'\,v''\ ,       v^{(3)}}{20160 } + \frac{v'\,^2\,v^{(4)}}{2240 } -     \frac{47\,{v^{(4)}}^2}{1209600}\,-\right.\nonumber\\     & & \qquad\qquad-\left.\frac{19\,v^{(3)}\,v^{(5)}}{241920}-\frac{v''\,v^{(6)}}{30240}-\frac{v'\,v^{(7)}}{60480 } + \frac{v^{(10)}}{10644480}\right]\nonumber\end{aligned}\ ] ]    @xmath175+\nonumber\\ & & + \,\epsilon_n^3\left[-\frac{{v^{(3)}}^2}{4032 } - \frac{v''\,v^{(4)}}{5040 } -     \frac{v'\,v^{(5)}}{6720 } + \frac{v^{(8)}}{483840}\right]+\nonumber\\ & & + \,\epsilon_n^4\left[\frac{v''\,^3}{60480 } + \frac{v'\,v''\ ,       v^{(3)}}{3360 } + \frac{v'\,^2\,v^{(4)}}{13440 } -     \frac{13\,{v^{(4)}}^2}{806400}\,-\right.\nonumber\\     & & \qquad\qquad-\left.\frac{v^{(3)}\,v^{(5)}}{26880}-\frac{v''\,v^{(6)}}{80640 } -     \frac{v'\,v^{(7)}}{161280 } + \frac{v^{(10)}}{21288960}\right]\nonumber\\ ( g_2)_{p=6}&=&\frac{v^{(4)}}{1920}+ \epsilon_n\,\frac{v^{(6)}}{53760}+ \epsilon_n^2\left[-\frac{{v^{(3)}}^2}{32256 } -     \frac{v''\,v^{(4)}}{40320 } -     \frac{v'\,v^{(5)}}{53760 } + \frac{v^{(8)}}{1935360}\right]+\nonumber\\ & & + \,\epsilon_n^3\left[-\frac{{v^{(4)}}^2}{345600 } -     \frac{v^{(3)}\,v^{(5)}}{138240 } -     \frac{v''\,v^{(6)}}{483840 } -     \frac{v'\,v^{(7)}}{967680 } + \frac{v^{(10)}}{85155840}\right]\nonumber\\ ( g_3)_{p=6}&=&\frac{v^{(6)}}{322560}+ \epsilon_n\,\frac{v^{(8)}}{11612160}\,+\nonumber\\ & & + \,\epsilon_n^2\left[-\frac{{v^{(4)}}^2}{4147200 } -     \frac{v^{(3)}\,v^{(5)}}{1658880 } -     \frac{v''\,v^{(6)}}{5806080 } -     \frac{v'\,v^{(7)}}{11612160 } +     \frac{v^{(10)}}{510935040}\right]\nonumber\\   ( g_4)_{p=6}&=&\frac{v^{(8)}}{92897280}+ \epsilon_n\,\frac{v^{(10)}}{4087480320}\nonumber\\ ( g_5)_{p=6}&=&\frac{v^{(10)}}{40874803200}\ .\nonumber\end{aligned}\ ] ]    the investigations presented in this paper were done at the scientific computing laboratory at the institute of physics in belgrade . the authors wish to thank the ministry of science and environmental protection of the republic of serbia for financing this investigation through projects 1486 and 1899 . 23 r. p. feynman and a. r. hibbs , _ quantum mechanics and path integrals _ ( mcgraw - hill , new york , 1965 ) . r. p. feynman , _ statistical mechanics _ ( w. a. benjamin , new york , 1972 ) . j. a. barker and d. henderson , rev . phys . * 48 * , 587 ( 1976 ) . m. h. kalos and p. a. whitlock , _ monte carlo methods _ , vol . 1 : basics ( john wiley and sons , new york , 1986 ) . d. m. ceperley , rev . phys . * 67 * , 279 ( 1995 ) . e. l. pollock and d. m. ceperley , phys . b * 30 * , 2555 ( 1984 ) . j. a. barker , j. chem . phys . * 70 * , 2914 ( 1979 ) . h. kleinert , _ path integrals in quantum mechanics , statistics , polymer physics , and financial markets _ , 3rd edition ( world scientific , singapore , 2004 ) . a. bogojevi , a. bala and a. beli , phys . lett . * 94 * , 180403 ( 2005 ) f. r. krajewski and m. h. mser , phys . rev . b * 65 * , 174304 ( 2002 ) . m. takahashi and m. imada , j. phys . soc . jpn . * 53 * , 3765 ( 1984 ) . x. p. li and j. q. broughton , j. chem . phys . * 86 * , 5094 ( 1987 ) . h. de  raedt and b. de  raedt , phys . a * 28 * , 3575 ( 1983 ) . s. jang , s. jang , and g. voth , j. chem . phys . * 115 * , 7832 ( 2001 ) . s. d. bond , b. b. laird , and b. j. leimkuhler , j. comp . phys . * 185 * , 472 ( 2003 ) . n. makri and w. h. miller , chem . lett . * 151 * , 1 ( 1988 ) . n. makri and w. h. miller , j. chem . * 90 * , 904 ( 1989 ) . n. makri , chem . . lett . * 193 * , 435 ( 1992 ) . m. alford , t. r. klassen , and g. p. lepage , phys . d * 58 * , 034503 ( 1998 ) . a. erdelyi , _ asyptotic expansions _ ( dover publications , new york , 1956 ) . url http://scl.phy.bg.ac.yu / speedup/. s. flgge , _ practical quantum mechanics _ , i ( springer - verlag , berlin , 1971 ) . a. srinivasan , m. mascagni and d. ceperley , parallel computing , * 29 * , 69 ( 2003 ) . url http://sprng.cs.fsu.edu/. a.  m. ferrenberg , d.  p. landau and y.  j. wong , phys . * 69 * , 3382 ( 1992 ) . w.  h. press , s.  a. teukolsky , w.  t. vetterling and b.  p. flannery , _ numerical recipes in c : the art of scientific computing _ ( cambridge university press , 1992 ) .
Abstract:

Answer: we present and discuss a detailed derivation of a new analytical method that systematically improves the convergence of path integrals of a generic @xmath0-fold discretized theory . we develop an explicit procedure for calculating a set of effective actions @xmath1 , for @xmath2 which have the property that they lead to the same continuum amplitudes as the starting action , but that converge to that continuum limit ever faster . discretized amplitudes calculated using the @xmath3 level effective action differ from the continuum limit by a term of order @xmath4 . we obtain explicit expressions for the effective actions for levels @xmath5 . we end by analyzing the speedup of monte carlo simulations of two different models : an anharmonic oscillator with quartic coupling and a particle in a modified pschl - teller potential .