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: statistical mechanics , be it for classical or quantum many - body systems , relies on the reasonable assumption that all microstates of a system in equilibrium are equally likely to be realized . in order to be correctly described by a statistical ensemble , a system must therefore be either coupled to an infinite thermal bath , or ( if isolated ) effectively display some form of ergodicity as per the eigenstate thermalization hypothesis @xcite , thereby allowing it to forget its initial state by ` relaxing ' towards a well - defined equilibrium state independent of the starting conditions . the existence of these relaxation processes is not usually questioned , but their precise mechanisms occur in a generally unspecified manner . for individual systems , classical mechanics distinguishes integrable and non - integrable cases , the former possibly having stable quasiperiodic orbits protected by the kam theorem @xcite , and therefore showing no true ergodicity . for non - integrable cases , ergodicity is taken for granted as a consequence of chaotic dynamics .    in the quantum case , the question of equilibration has been the subject of many recent works , focusing in particular on the time evolution of a system after a quantum quench , at which a global parameter in the hamiltonian is suddenly changed @xcite - @xcite . this is equivalent to releasing a prepared state and letting it evolve unitarily in time according to the hamiltonian after the quench . such situations , some of which are experimentally realizable , lead to many interesting theoretical questions . under what circumstances is there a well - defined state at large times ( _ i.e. _ what , if anything , remains of the initial conditions for given quench situations ) ? can a general theory of relaxation based purely on the details of quantum dephasing be formulated ? are there different classes of systems ( _ e.g. _ integrable versus non - integrable , finite versus infinite , gapped versus gapless ) having different generic behaviour ? what remains of the usefulness of commonly used quasiparticle bases to understand the dynamics after the quench ? clearly , much progress remains to be done in order to have a full understanding of all these issues . it is our purpose here to consider a well - defined situation , which can be handled to a sufficient degree of exactness to provide some partial but reliable insights into some of these questions . the approach we will consider is based on the use of exact wavefunctions for the heisenberg magnet , and therefore exploits the integrability of this system . on the one hand , this renders some of our results non - generic . on the other hand , we are able to provide hard facts for finite as well as infinite systems , gapped versus gapless , and exclude some possible scenarios . since the theoretical understanding of the nonequilibrium dynamics of strongly - correlated systems is still in its infancy , such example cases will hopefully provide worthwhile reference points for later developments . the paper is organized as follows . in section 2 , we first define our nonequilibrium problem , and present the tools we will use to study it . we then discuss in section 3 the set of eigenstates we use to obtain quantitative results . section 4 considers the work probability distribution resulting from the quench for various system parameters . in section 5 , we study relaxation dynamics by concentrating on the loschmidt echo , which we compute . after a few words on the thermodynamic limit , we offer a discussion of our results and present our conclusions . all technical details for the computations are relegated to a series of appendices , which can safely be skipped by the reader only interested in the final results and conclusions . on the other hand , these appendices explain all calculations in sufficient detail to be reproduced by the specialist reader . we consider an isolated spin chain of @xmath2 sites , each occupied by a local spin-@xmath0 degree of freedom . for definiteness , we put the spins on a ring and impose periodic boundary conditions . the thought experiment we perform consists in preparing the quantum state of the system at @xmath3 as @xmath4 this state thus contains a magnetic domain wall between sites @xmath5 and @xmath6 , and another ( anti- ) domain wall between sites @xmath2 and @xmath7 ( in view of the periodic boundary conditions ) . this state can be prepared in different ways : we can view it as being created by an ising model with an appropriate position - dependent field , or as resulting from a sudden polarizing pulse applied on a section of an initially fully polarized chain . for times @xmath8 , we let this state evolve unitarily in time under the antiferromagnetic xxz hamiltonian @xmath9.\ ] ] note that the hamiltonian is rescaled by a factor of @xmath10 as compared to how it usually appears in the literature , in order the have a well - defined ising limit ( @xmath11 ) . for @xmath1 the spectrum of this theory is gapped , while for @xmath12 the system is in the quantum critical regime . we consider the case when the system is antiferromagnetic ( i.e. @xmath13 ) . in order to simplify the expressions and without loss of generality we put @xmath14 throughout the paper . the initial domain wall state ( [ domainwallstate ] ) is not an eigenstate of the xxz hamiltonian away from the ising limit @xmath15 . on the other hand , for any fixed value of @xmath16 , the exact eigenstates of this model form a basis in the hilbert space on which we can at least in principle decompose the initial domain - wall state to arbitrary accuracy . given this decomposition , the solution of the schrdinger equation becomes straightforward , and we obtain the exact time - dependent wavefunction after the quench as the linear decomposition @xmath17 eigenstates of @xmath18 at fixed total magnetization , and the complex amplitudes @xmath19 represent the overlaps ( vector of the quench matrix ) between the starting state and the exact eigenfunctions . since we always work with normalized wavefunctions , the coefficients @xmath19 should by definition satisfy the constraint @xmath20 this constraint will constitute an important sum rule , allowing to quantify the accuracy of our results by assessing how faithfully the resulting wavefunction is reproduced .    all the complexity of the problem is therefore hidden in two places . first , wavefunctions @xmath21 and their energy @xmath22 must be known . this is standardly handled by the bethe ansatz ( see references @xcite - @xcite and references therein ; we provide a summary of the necessary details in appendices a and b ) . second , the overlaps @xmath19 must also be known . this is a problem of much greater complexity , which in the present situation finds its solution in the framework of the algebraic bethe ansatz . the derivation of these overlaps is given in appendices c and d. note that similar overlaps ( quench matrix entries ) were also recently calculated using the algebraic bethe ansatz for the case of the interaction - quenched richardson model @xcite . here , we also obtain these overlaps using integrability , but using a different method based on the explicit structure of the solutions to the bethe equations .    despite being in possession of these two fundamental building blocks , one major difficulty remains . since the hilbert space is exponentially large in system size @xmath2 , the summation in ( [ timeevolphi ] ) is difficult to handle , and must in practice be truncated in order to reach sizes sufficiently large to allow the discrimination between finite- and infinite - size behaviour . we will show in the next section that this truncation is possible in the situation we consider . in particular , this puts us in position to study the long - time average of observables according to any prescription desired , a common one being @xmath23 and to make reliable observations on the relaxation / thermalization of the prepared state . the time evolution of this domain wall state has been studied in similar settings . for gapless spin chains this type of quench is studied in @xcite . an exact analysis has been carried out for the xx - chain in @xcite . the short time regime for the xxz chain is studied using tdmrg in @xcite , and an exact diagonalization was performed for the gapped xxz chain with open boundary conditions @xcite . we here offer a complement to these studies , consisting in results from integrability . the space of eigenstates of the xxz chain is spanned by bethe wavefunctions , each described by a set of rapidity parameters @xmath24 obtained as solution to the bethe equations . bethe wavefunctions are generically quite complicated objects , and can contain both unbound and bound magnons . we refer the reader to the standard literature @xcite as well as to our appendices for the basic detail . in summary , solutions of the bethe ansatz equations can be classified using strings ( see appendix [ a : bae ] for a discussion ) . an @xmath25string in a set of @xmath26 complex rapidities sharing the same real part and invariant under complex conjugation . these strings can be interpreted as bound states of ` mass ' @xmath26 . in the ising - limit it can be shown that an @xmath26-string corresponds to @xmath26 adjacent down spins . from a perturbative analysis , we expect at large anisotropy @xmath16 a well - defined hierarchy of states in terms of the importance of their overlap with the domain wall state . this hierarchy is presented in table [ tab : hierarchy ] . the most important states are those consisting of a single @xmath27string , which occupy a dispersionless line in the thermodynamic limit . they are followed by successively more complicated partitionings of the @xmath5 rapidities into more and more individual bound or unbound states , each partitioning representing a whole continuum of excitations . the number of particles in the state is defined here to be the number of elements in the partitioning . from a perturbative analysis it can be easily shown that for sufficiently large @xmath5 the system is insensitive to small changes of @xmath5 ( _ i.e. _ , of how distant the domain and anti - domain walls are from each other ) . this allows us to perform the calculation with an @xmath5 slightly below @xmath28 without loss of generality . the reason for this choice is that for these values the string hypothesis is better satisfied for the majority of string states . .the hierarchy of states , in order of importance of contributions to the normalization sum rule ( [ sumrule ] ) . each partitioning of @xmath5 rapidities into strings leads to an independent excitation class , only the few shown here being of relevance to the current setup . [ cols="^,^,^,^",options="header " , ] in this work we studied the relaxation dynamics of the gapped xxz chain after a quench from an initial domain wall state . the techniques that were used are based on the algebraic bethe ansatz for a finite chain . some of the results were extended to the thermodynamic limit . the key ingredient for studying this quench was deriving a numerically efficient expression for the overlap between the initial state and eigenstates of the xxz hamiltonian . overlaps for different ising states as initial state could be computed following the same logic , however in the general case of a finite density of domain walls this will become an intractable combinatorial problem . since the energies of the relevant part of the spectrum are intensive a good comparison with the thermodynamic limit is possible . one important result is that the bandwidth of the @xmath27strings is nonzero for finite systems allowing the long time average of the loschmidt echo going to zero , something that does not happen in the thermodynamic limit . hence we conclude that the limits @xmath29 and @xmath30 do not commute , which was also pointed out in @xcite .    in the thermodynamic limit we derived a lower bound for the infinite time average of the loschmidt echo and argue that this bound is the exact value . this led us to the most important conclusion of this paper , namely that the system for this quench will not thermalize . the reason why the system does not thermalize can be explained from its spectrum , and does not directly rely on the fact that the system is integrable . imagine for instance a gapped system with lowest ( or highest ) band having zero bandwidth in the thermodynamic limit , such as a system with macroscopically - degenerate ground states . consider a quench such that @xmath31 lies right between the lowest ( or highest ) energy band and the next one . a lower bound for the overlap @xmath32 of the lowest band can then be obtained by simple reasoning . let @xmath33 be the energy of the lowest energy band of zero bandwidth . we denote the band gap by @xmath34 . the expectation value of the initial hamiltonian is @xmath35 . we then choose @xmath36 . in this case we can obtain a lower bound for the total overlap of the @xmath33 states . @xmath37 from which follows the lower bound @xmath38 in other words , such a quench must occupy the lowest ( resp . highest ) dispersionless band with @xmath39 amplitude . since these states are dispersionless and separated from other states by a band gap , they can not dephase , and therefore can not effectively relax . the large - time asymptotic state will therefore always maintain memory of the ground state occupation amplitude of this band , which is a quench - dependent quantity . note that such situations trivially can not show eigenstate thermalization according to the arguments of @xcite , in which such circumstances were argued to be non - generic and therefore not treated . the condition used for deriving this lower bound is quite strong . in figure [ fig : spectrum_overlaps ] one can see that for @xmath40 that @xmath31 does not lie in the gap anymore , but the overlaps of the @xmath27string are still substantially big . one can also question whether the condition of an isolated peak is essential . from the work probability distribution it is clear that the contributions of states made up of different strings can be distinguished . simply looking at energy levels is therefore insufficient to assess the presence or absence of relaxation , and one must also take into account the values of the wavefunction overlaps . if we consider for example @xmath26 strings with a diverging length in the thermodynamic limit , they all will become coherent . although they are embedded in a continuum , we can reasonably expect a non - zero height peak of zero width in @xmath41 , which again will prevent thermalization .    in the case of evolution under a gapless hamiltonian , our method loses its efficiency , since all overlaps scale to zero and the number of states to take into account grows accordingly , preventing an efficient truncation of the sum in [ timeevolphi ] . the absence of a gap thus opens the door to dephasing involving arbitrarily complex excitation continua , in correspondence with anderson s orthogonality catastrophe principle . therefore it is expected that @xmath42 in this case , and the possibility of thermalization remains , although it can not be quantified here . the presence of a gap is thus determinantal to the long - time asymptotics of dynamics after the quench release of the domain wall state , but is not very sensitive to system size for large enough systems . to summarize we considered a quench starting from a domain wall state under the time evolution of the gapped xxz chain . we observed that the long time behavior can not be described by a statistical ensemble . this can be understood by the fact that the spectrum is gapped , and does not directly rely on the fact that the system is integrable . furthermore we were able to assess finite size effects by considering both results for the finite and infinite systems . general conclusions for the long time behavior of the xxz model after quench can not however be drawn from the results presented here . for instance , the two quasi - degenerate ground states of the xxz chain we are dealing with are separated from all other states by the gap , so one might expect a similar effect if one would perform a quench starting from the nel state . however , the low energies are extensive in contrast with the high energy excitations which are intensive . the effect of the gap is much less pronounced and one might at least expect that the staggered magnetization will vanish for all @xmath1 , in correspondence with the results of @xcite . in future work , we will consider this kind of initial state as well as the time evolution of more elaborate observables . in the algebraic bethe ansatz the eigenstates are completely characterized by a set of rapidities @xmath43 , solutions to the bethe equations . these solutions can be either real or complex . the complex solutions typically come in groups , invariant under complex conjugation , which represent bound states . at a small density of down spins or large anisotropy @xmath16 , most solutions take the form of so - called strings , in which a number of rapidities share a real centre while the imaginary parts are equally spaced according to the so - called string hypothesis . an @xmath25string takes the form :    @xmath44    where @xmath45 , and @xmath46 is a deviation . when the string hypothesis holds ( i.e. when all deviations @xmath46 are sufficiently small ) the bethe equations can be rewritten as equations for the real parts of the rapidities only , yielding the bethe - gaudin - takahashi equations @xcite . the logarithmic version of the bethe - gaudin - takahashi equations for @xmath1 is given by @xmath47 with the dispersion and scattering kernels @xmath48 @xmath49 solutions of the bethe - gaudin - takahashi equations can be classified by the set of quantum numbers @xmath50 , which are integers ( half - odd integers ) if @xmath51 is odd ( even ) . the partitioning into strings satisfies the constraint @xmath52 . the energy of an @xmath25string @xmath53 is @xmath54 the total energy and momentum of an eigenstate containing strings are given by @xmath55 @xmath56 for finite @xmath57 and large @xmath26 the energy @xmath58 becomes dispersionless . an @xmath25string can be interpreted as a magnon bound state with a mass proportional to @xmath26 .      in the thermodynamic limit @xmath59 the band structure is easily obtained from . for a given string configuration @xmath60 one can define : @xmath61 by putting all @xmath53 to @xmath62 . similarly one can define @xmath63 by putting all @xmath53 to @xmath64 . the @xmath27string is a special case since it has a vanishing bandwidth in the limit @xmath65 . in the ising limit ( @xmath66 ) the energy of a state made of @xmath67 different strings is @xmath68 . this corresponds in the ising language to a state made of @xmath67 distinct blocks of down spins and clarifies the bound states interpretation in this limit . in order to evaluate the spectrum numerically a correspondence between quantum numbers and rapidities is needed . for the isotropic chain @xmath69 every set of quantum numbers corresponds to a unique solution of the bethe - gaudin - takahashi equations . a complete classification of solutions in terms of quantum numbers was given in @xcite . for the case of the xxz chain with @xmath70 this classification in not known in general . for solutions with only real solutions this was done in @xcite . here we generalize the result to string solutions . note that a genuine proof that the bethe - gaudin - takahashi equations produce the right number of eigenstates for @xmath1 remains to be found . in @xcite it was argued that @xmath72 if @xmath73 . equal rapidities do not yield proper bethe wavefunctions , and we therefore only need to consider sets of distinct quantum numbers . in the parameterization chosen in the bethe - gaudin - takahashi equations all rapidities are restricted to an interval of width @xmath74 , i.e. for all strings of given length @xmath26 we have : @xmath75 . from this requirement and assuming that the bethe - gaudin - takahashi equations are monotonic in @xmath76 it follows that : @xmath77 with @xmath78 . we call @xmath71 the bandwidth of the quantum numbers @xmath79 .      since @xmath81 as far as the wavefunctions are concerned , there is no one - to - one mapping of quantum numbers and wavefunctions . namely , the shift @xmath82 leaves the wavefunction invariant although the quantum numbers in the bethe - gaudin - takahashi equations change . for every string length @xmath26 we can define a transformation : @xmath83 with : @xmath84 from this transformation we see that @xmath85 . for every string length @xmath26 there are @xmath86 bands of quantum numbers which can not be transformed into each other by @xmath87 , whereas in the xxx case there is only one such band for every string length . if we restrict to only one string length it follows that the number of choices for the quantum numbers @xmath88 in a single band is : @xmath89 . in case of @xmath80 bands the number of solutions is : @xmath90      to find the minimal set of quantum numbers one should first determine using @xmath87 the bandwidth @xmath71 and the number of bands @xmath80 for every string length @xmath26 . however , the transformation @xmath87 not only affects the quantum numbers of length @xmath26 but also all the others . to take this effect into account one should also consider transformations of the form @xmath91 , leading to additional state exclusions . in general these are very difficult to determine . the explicit classification for two string types is given below . we consider @xmath67 @xmath25strings and @xmath92 @xmath93strings such that : @xmath94 . the bandwidths of the two sets of quantum numbers are : @xmath95 the number of bands are @xmath96 and @xmath97 . to take care of the degeneracies coming from different string sectors we consider the transformation : @xmath98 this turns out to be the only extra transformation one needs to consider . define : @xmath99 . from the conditions : @xmath100 and @xmath101 follows that we have to exclude : @xmath102 and @xmath103 . the total number of unique solutions for this case is @xmath104 the algebraic bethe ansatz deals with the problem of diagonalizing simultaneously the transfer matrix @xmath105 for all values of @xmath106 , from which all conserved quantities can be obtained @xcite . all operators can be expressed in terms of four non - local hilbert space operators @xmath107 and @xmath108 . the algebraic bethe ansatz requires a pseudo - vacuum @xmath109 which in case of the xxz chain is the state with all spins up , such that : @xmath110 the transfer matrix is constructed as : @xmath111 . the vacuum eigenvalues of the @xmath112 and @xmath108 operators are : @xmath113 @xmath114 the set @xmath115 are inhomogeneity parameters , which should be set to @xmath116 in order to obtain the xxz chain . states , respectively dual states can be constructed as : @xmath117 in order to represent an eigenstate the rapidities @xmath118 should satisfy the bethe equations . these states are not automatically normalized , we therefore write : @xmath119 . the norm of a bethe state in case of strings is given by @xcite : @xmath120 where @xmath121 is called the reduced gaudin matrix : @xmath122 an extremely useful formula is slavnov s expression for the scalar product between an eigenstate @xmath123 ( represented by a set of rapidities @xmath124 that satisfy the bethe equations ) and an arbitrary state @xmath125 ( without restrictions on @xmath126 ) @xcite : @xmath127 with @xmath128      the initial state can be written in the algebraic bethe ansatz language using the inverse mapping @xcite : @xmath129 making use of @xmath130 and @xmath131 we can write the state as : @xmath132 which is a normalized state . the overlap between an eigenstate for finite @xmath16 and the initial state @xmath133 can be computed using slavnov s theorem for the scalar product . in the case when @xmath134 the matrix @xmath135 simplifies : @xmath136 for the xxz spin chain all inhomogeneity parameters should be set to @xmath137 . the limit should be taken carefully , there are @xmath138 poles coming from : @xmath139 . since the columns of the matrix @xmath135 are also @xmath138-fold degenerate we can apply lhpital s rule . first we notice that we can easily extract factors @xmath140 from the determinant , since applying lhpital s rule to those terms does not remove any degeneracies . @xmath142 note that this determinant of @xmath143 is the same as that in the partition function of the 6-vertex model with domain wall boundary conditions @xcite . however in this case both the sets @xmath144 and @xmath145 play the role of boundary conditions and satisfy no bethe equations . to be able to apply lhpital we are interested in @xmath146 . first we rewrite @xmath147 : @xmath148 next , consider @xmath149 using this result it is easy to derive for @xmath26 even : @xmath150 and for @xmath26 odd we get a similar result . in general the coefficients @xmath151 are complicated expressions . however it is straightforward to derive that @xmath152 for @xmath26 even . since in this case there is no mixing between the left and right terms in . it will turn out that theses coefficients are the only ones we need , since all terms with coefficient @xmath153 for @xmath154 can be removed using row manipulations while leaving the determinant invariant . differentiating the denominator : @xmath139 is much simpler . first we put @xmath155 . if we now consider @xmath156 we get a single zero in the limit @xmath157 . we can continue this logic and see that @xmath115 gives rises to @xmath158 zero s , hence we get a factor @xmath159 . we also notice that the factors @xmath160 from the numerator cancel against the factors from the denominator so we are left with : @xmath161 the matrix @xmath162 has a similar structure to a vandermonde determinant . this makes the matrix ill - conditioned and is therefore not suitable for a numerical evaluation . however , for a given string structure the determinant can be written in terms of sums of vandermonde determinants . in the following sections this is done explicitly for the most important cases . the normalized overlap is :    @xmath163    these expressions hold for all three parameterizations in .      in the case where the solution is an @xmath27string : @xmath164 the determinant of @xmath162 is exactly a vandermonde determinant . a vandermonde matrix is defined as : @xmath165 with all @xmath166 distinct . its determinant is given by : @xmath167 introducing the notation @xmath168 the determinant of @xmath162 can be written as an @xmath169 vandermonde matrix determinant : @xmath170 the reduced gaudin determinant for an m - string takes the simple form : @xmath171 this results in the normalized overlap for an m - string : @xmath172 the absolute value squared of the overlap is @xmath173      we consider states constructed from two strings with lengths @xmath67 and @xmath92 such that @xmath175 , and @xmath176 : @xmath177 and @xmath178 . some convenient notations are : @xmath179 now the overlap matrix takes the form : @xmath180 expanding this determinant in the first column results in a sum over @xmath181 vandermonde determinants . after factoring out common factors we obtain : @xmath182      consider a state constructed of @xmath26 strings with lengths : @xmath183 . introducing the notation @xmath184 the overlaps can be written as a sum over vandermonde determinants : @xmath185 the bethe - gaudin - takahashi equations for the @xmath27strings reduces to the simple equation : @xmath186 for even @xmath5 and @xmath70 the @xmath2 possible quantum numbers @xmath187 are @xmath188 which follows from . in the thermodynamic limit we can define a density function @xmath189 for an m - string solution centered around @xmath106 : @xmath190 such that @xmath191 . this function @xmath192 should be interpreted as the number of solutions there are with a rapidity in the domain @xmath193 . the total contribution of all @xmath2 @xmath27string solutions is : @xmath194 where @xmath195 is the hypergeometric function : @xmath196 in the limit when both @xmath5 and @xmath2 are sent to infinity a simpler expression can be obtained : @xmath197      the isotropic limit @xmath198 should be taken with care , since some m - strings solutions correspond to different string configurations . for the xxx case we have the following restriction for the quantum number of an m - string : @xmath199 . from the bethe - gaudin - takahashi equation : @xmath200 we see that in the limit @xmath201 and finite @xmath5 allowed values for @xmath106 go from minus to plus infinity . in this limit we can derive the weight function : @xmath202 so the contribution of the @xmath5-string for finite @xmath5 becomes : @xmath203 both authors gratefully acknowledge support from the stichting voor fundamenteel onderzoek der materie ( fom ) in the netherlands . we thank g. palacios for useful discussions .
Abstract:

Answer: we study the dynamics of a quench - prepared domain wall state released into a system whose unitary time evolution is dictated by the hamiltonian of the heisenberg spin-@xmath0 gapped antiferromagnetic chain . using exact wavefunctions and their overlaps with the domain wall state allows us to describe the release dynamics to high accuracy , up to the long - time limit , for finite as well as infinite systems . the results for the infinite system allow us to rigorously prove that the system in the gapped regime ( @xmath1 ) can not thermalize in the strict sense .


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 concept of chiral symmetry breaking in strong interactions goes back to the pioneering work by nambu @xcite , and it has been widely explored and discussed since that time . restoration of chiral symmetry at high densities was proposed by lee and wick @xcite already in 1974 . theoretical procedures for studying chiral symmetry restoration at high temperatures were proposed on the same year by dolan and jackiw @xcite and by weinberg @xcite . qcd at finite temperature and density has attracted much interest . for latest studies of the qcd phase diagram see for instance @xcite . physical applications are to high energy heavy ion collisions @xcite and to the physics of the early universe @xcite . a review of lattice calculations can be found for instance in ref . the nambu - jona - lasinio model has been used under different approximations to obtain indications on the chiral transition ( see for instance @xcite ) . a number of peculiar features may emerge due to the chiral phase transition and recently much attention has been attracted by the possibility of disoriented chiral condensate formation ( see for instance @xcite ) . it is strongly suspected that a single transition occurs rather than separate transitions for deconfinement and chiral - symmetry . we had suggested a heuristic argument indicating that , at least for zero density , the critical temperature for chiral transition , @xmath3 , coincides with that for deconfinement , @xmath4 @xcite . the order parameters usually considered cover extreme and opposite ranges . the thermally averaged polyakov loop is suitable in the limit of infinite quark masses to describe the transition from the confined to the non - confined phase . the other extreme is the limit of vanishing quark masses , where the thermally averaged quark - antiquark bilinears are the typical order parameters for chiral symmetry transition . we shall deal here with the chiral transition , concentrating on the thermally averaged quark bilinears at finite temperatures as order parameters for chirality . for light massive quarks , such as @xmath5 and @xmath6 , rigorously , chiral symmetry is already broken in the lagrangian , but we can still retain the notion of phase transition , by looking at the region of @xmath7 where the condensate has a rapid variation . the current quark mass plays a role analogous to that of an external magnetic field in the ferromagnetic transition , as it explicitly violates the chiral symmetry whose restoration characterizes the phase transition .    the analysis will be based on a composite operator formalism at finite temperature . the formalism makes use of an effective action for composite operators @xcite . actually the effective potential admits of a landau expansion around the critical point and thus the behaviour of the condensate is well reconstructed by knowing the coefficients ( which are infrared safe ) and the critical exponents . we shall compare with the results following from a landau expansion and show the general agreement .    within the composite operator formalism we had already discussed the @xmath7 dependence of @xmath8 @xcite , in the whole range of temperatures up to @xmath3 . we shall here mainly deal with the ratio of the scalar to pseudoscalar mass @xmath9 at varying temperature . this ratio has a peculiar theoretical interest as a sensible indicator of the degree of chiral symmetry breaking . let us for the moment neglect the small @xmath5 and @xmath6 quark masses . when chiral symmetry is restored the two mesons are degenerate in mass . the ratio is then equal to unity . at zero temperature instead in the broken chiral phase the pion is a goldstone and has vanishing mass , while the sigma has a finite mass from the chiral condensation . thus one expects the mass ratio to decrease from @xmath10 to one at the chiral transition . quark masses will change this picture quantitatively leading to a decrease from a finite value to a value close to one . beyond some temperature , before approaching the critical value for the transition , the sigma will not have phase space left to decay into two pions . its instability will decrease when increasing the temperature , and the channel into two pions will finally be suppressed at some temperature in the vicinity but lower than the critical temperature . the decay channel into two photons will still be available and become the dominant channel @xcite . it is too early to say whether a possible experimental signature for the transition may be inferred from such behaviour . an accurate determination of the mass ratio versus temperature is in itself of theoretical interest as it constitutes a significant parameter for the degree of chiral symmetry breaking at finite temperature . this we shall do in this note in the composite operator formalism by making use of parameter values obtained from fits to zero temperature qcd . in section 2 we shall summarize the main results of the application of the composite operator formalism in qcd . in section 3 we determine the relevant qcd parameters and discuss the hadronic masses at non zero temperature . the temperature dependence of the observables is discussed in section 4 , where we also calculate the temperature at which the sigma can no longer decay via strong interaction . following ref .  @xcite the zero temperature euclidean effective action for an @xmath11 qcd - like gauge theory is @xmath12 - { \rm tr}\left [ { \bf \sigma}{\bf s}\right ] -{\bf \gamma}_2({\bf s } ) + counterterms \label{effact}\ ] ] where @xmath13 , * m * is the bare quark mass matrix and @xmath14 is the sum of all the two - particle irreducible vacuum diagrams with fermionic propagator @xmath15 and @xmath16 . at two - loop level @xmath17 , where @xmath18 is the gauge boson propagator , so that @xmath19 , @xmath20 = 2{\bf \gamma}_2 $ ] and one can rewrite eq . ( [ effact ] ) in terms of @xmath21 @xmath22 + { \bf \gamma}_2({\bf \sigma } ) + counterterms\nonumber\\ & = & - { \rm tr}\ln\left[{\bf s}_0^{-1 } -{\bf \sigma } \right ] + { 1\over 2}{\rm tr}\left({\bf \sigma}\delta^{-1}{\bf \sigma } \right ) + counterterms \label{effact1}\end{aligned}\ ] ]    here the variable @xmath21 plays the role of a dynamical variable . at the minimum of the functional action , that is when the schwinger - dyson equation is satisfied , @xmath21 is nothing but the fermion self - energy . a parametrization for @xmath21 , employed in @xcite was @xmath23 with a suitable ansatz for @xmath24 , and with * s * and * p * scalar and pseudoscalar constant fields respectively to be taken as the variational parameters . the effective potential one obtains from eq . ( [ effact1 ] ) ( see ref . @xcite ) is @xmath25-\nonumber\\ & & - n~{\rm tr}\ln\left[i\hat{k } - \left({\bf m}+{\bf \sigma}_s \right)-i\gamma_5{\bf \sigma}_p\right ] + \delta z~{\rm tr}({\bf m}~{\bf s } ) \label{effpot}\end{aligned}\ ] ] where @xmath26 is the quadratic casimir of the fermion representation ( for @xmath27 @xmath28 ) . furthermore @xmath29 , @xmath30 , @xmath31 ( @xmath32 , @xmath33 , @xmath34 gell - mann matrices , @xmath35 ) , @xmath36 is the gauge coupling constant and @xmath37 is the four - dimensional volume . in eq . ( [ effpot ] ) @xmath38 has a divergent part to compensate the leading divergence proportional to @xmath39 in the logarithmic term . we remark that both the current mass * m * and the self - energy @xmath21 are , in general , matrices in flavour space . however , as discussed in refs . @xcite if we neglect the mixing between different flavours originating , for instance , from terms such as the t hooft determinant @xcite , it follows that only the flavour diagonal elements of the fermion self - energy and mass can be different from zero at the minimum . with vanishing off - diagonal terms , the effective potential decomposes into the sum of @xmath40 contributions ( @xmath40  =  number of flavours ) , one for each flavour . therefore , to study the minima of the effective potential , it is formally sufficient to consider a single contribution . of course , the choice of a given flavour number will reflect in the particular parameters assumed . in the present paper , as in ref.@xcite , we will take @xmath41 and a number of colour @xmath42 . the value of the parameters will be specified later on .    as far as the choice for the function @xmath24 is concerned , in qcd , the operator product expansion suggests ( neglecting logarithmic corrections ) to take for @xmath24 a momentum behaviour as @xmath43 for large @xmath44 . we have chosen as a variational ansatz @xcite @xmath45 where @xmath46 is a momentum scale which is expected to be of the order of @xmath47 .    to extend the zero - temperature theory to finite temperature we can still work with continuous energies by substituting for the sum over discrete energies @xmath48 ( where @xmath49 ) a sum of integrals over continuous energies by means of the poisson s formula @xcite @xmath50    this substitution corresponds to the imaginary - time formalism after using the poisson s formula and allows for an equivalent version af the dolan - jackiw finite temperature feynman rules @xcite .    in conclusion the final form for the effective potential , for a quark of mass @xmath51 , is ( see ref . @xcite ) @xmath52-\nonumber\\ & & - 2n~\sum_{n=-\infty}^{n=+\infty}~(-)^{n}~\int{d^{4}k\over ( 2\pi)^{4}}~ \ln\left[k^2 + \left(m + \sigma_s\right)^2+\sigma^{2}_{p}\right ] + \delta z~ms \label{effpot1}\end{aligned}\ ] ] where @xmath53 and @xmath54 . let us now comment on the choice for the gauge coupling constant . as suggested by asymptotic freedom and renormalization - group considerations , we expect the strong forces to weaken at high temperature @xcite . we shall then assume that in the uv region the coupling constant @xmath55 depends logarithmically on the temperature @xmath7 . we take into account this assumption by writing @xmath56 where @xmath57 and we will discuss the parameters @xmath58 , @xmath46 and @xmath59 later on .    at @xmath60 the effective potential does not acquire any extra divergence with respect to the @xmath61 case . the renormalization at @xmath61 can be performed by adding a counterterm and requiring that the derivative of the effective potential with respect to the term which breaks explicitly the chiral symmetry , evaluated at the minimum , satisfies for each flavour the renormalization condition @xcite @xmath62 the fermion condensate @xmath63 is related to the minimum @xmath64 of the effective potential , renormalized at the scale @xmath46 through the relation ( see ref . @xcite ) @xmath65 and @xmath66 . at finite @xmath7 we do not have any additional divergence with respect to the @xmath61 case ; nevertheless , in order to satisfy the generalization of the normalization condition ( [ normcond ] ) at finite @xmath7 we have to add a finite counterterm to that determined at @xmath61 @xcite . finally we recall that with the appropriate normalization for the pion field ( see next section ) , eq . ( [ normcond ] ) and its generalization at finite @xmath7 is also equivalent to the adler - dashen formula @xmath67 to compute the masses of the scalar and pseudoscalar mesons ( @xmath68 and @xmath69 respectively ) one has to take the second derivative of the effective potential ( [ effpot1 ] ) with respect to the scalar field @xmath70 and the pseudoscalar field @xmath71 , evaluated at the minimum ( for the sake of simplicity we are assuming that the up and down quarks are degenerate in mass and @xmath72 in eq . ( [ effpot1 ] ) ) . the actual values of the masses will be obtained by multiplying the second derivative by the appropriate factor @xmath73 that relates the physical fields @xmath74 ( @xmath75 ) to @xmath71 ( @xmath70 ) according to @xmath76 . this factor can be obtained in terms of the decay constant @xmath77 through standard arguments of current algebra @xcite . one gets @xmath78 where @xmath79 is the extremum of the effective potential in the presence of a bare mass .    to derive a more physically transparent expression for the masses of the mesons as expressed by eq . ( [ masses ] ) through the second derivative of the effective potential , we will use the gap equation and the generalization at finite @xmath7 of the normalization condition ( [ normcond ] ) . the extremum condition is @xmath80 } ~e^{\displaystyle{in\beta k_{0}}}\right]+~m~\delta z=0 \label{extrcond}\\ { \partial v\over \partial p}=0 \rightarrow   { \bar p}=0 \label{extrcond1}\end{aligned}\ ] ] where @xmath81 . as far as the normalization condition is concerned , by using the relation ( [ condenst ] ) between the scalar field at the minimum and the scalar condensate we shall write in general @xmath82 \label{deltaz}\ ] ] with @xmath83 and where @xmath84 is the minimum of the effective potential in the massless case . now , by using the gap equations ( [ extrcond])-([extrcond1 ] ) and the normalization condition ( [ deltaz ] ) in the explicit expressions ( [ masses ] ) for the meson masses , we can eliminate @xmath38 and @xmath85 to finally obtain , for a quark of mass @xmath51 : @xmath86 ~e^{\displaystyle{in\beta k_{0 } } } \label{emmepai2}\ ] ] @xmath87^{2 } } ~e^{\displaystyle{in\beta k_{0 } } } \label{emmesig2}\ ] ]    we notice that eq . ( [ emmepai2 ] ) reproduces , in the small mass limit , the adler - dashen formula ( [ adash ] ) and differs from this formula only by terms of order @xmath88 .    in order to fix the values of the parameters @xmath89 and @xmath90 we follow the same procedure of ref .  @xcite . in the massive case the effective potential is uv divergent and from the normalization condition ( in the small - mass limit ) and gap equation one is able to fix , through a self - consistency relation , the values @xmath91 . by inserting this value in the gap equation for massless quarks , one finds that chiral symmetry is spontaneously broken ( at @xmath61 ) via a minimum of @xmath92 located at @xmath93 and the point @xmath94 is a local maximum . to determine the mass scale @xmath46 and the quark masses from the experimental data , one has to derive the explicit expressions for the masses and for the decay constants of the pseudoscalar octet mesons which are the pseudo - goldstone bosons of chiral symmetry breaking . these expressions constitute a system of coupled equations which we have solved by an approximation method . the experimental inputs are the decay constant and mass of the charged pion , the charged kaon mass and the electromagnetic mass difference between the neutral and the charged kaon . the outputs of the numerical fit for the octet meson masses ( agreement within 3% ) are the masses of the @xmath5 , @xmath6 and @xmath70 quarks and the mass parameter @xmath46 . the values we get are the following @xcite @xmath95 finally , the parameter @xmath59 , which appears in the expression ( [ coupling ] ) for the running coupling constant , has to be determined on phenomenological grounds . if we compare our model in the low-@xmath7 regime with the results of ref . @xcite we find @xmath96 for @xmath40 flavours , which gives @xmath97 for @xmath41 @xcite .    with these values for the relevant parameters it is easy to determine the masses of the pseudoscalar and scalar mesons at @xmath61 ( we are using the value @xmath98 for the pion decay constant ) @xmath99 and ( again at @xmath61 ) the fermionic condensates at @xmath46 for massless and massive quarks @xmath100 the value we have obtained for the mass of the scalar meson , @xmath101 , accords very well with those obtained by other authors within different approaches to the study of dynamical chiral symmetry breaking in qcd @xcite . it is well known that as the temperature @xmath7 increases we expect some changes in the hadronic properties of matter as a consequence of a possible phase transition . for instance , a physical quantity as @xmath8 will decrease as @xmath7 increases @xcite . this means that the pion tends to decouple from matter ( quarks and leptons ) as @xmath7 increases . as far as the pion mass is concerned , @xmath103 has a weak dependence on the temperature . in fact the value of @xmath103 is dominated by the current quark mass for temperatures below the critical value , whereas it becomes independent of @xmath51 at the critical temperature or above it . this should be expected because , for @xmath104 , the pion becomes an ordinary resonance and its mass is dominated by @xmath47 and not by @xmath51 . this is a clear sign that the pion loses its goldstone nature once the approximate chiral symmetry is restored @xcite . the scalar meson mass @xmath105 will decrease in association with partial chiral symmetry restoration in hot matter . in fact the scalar @xmath68 meson can be described as a @xmath106 quasi bound state and its mass is about twice the constituent quark mass @xcite . thus the @xmath68-mass will decrease substantially with increasing temperature and this behaviour is mainly determined by the @xmath7-dependence of the constituent quark mass . the decreasing of the @xmath68-mass is expected from the restoration of chiral symmetry at @xmath107 . in fact @xmath108 for @xmath109 at @xmath107 . as a consequence @xmath105 will become smaller than @xmath110 at some temperature @xmath111 and the @xmath68-meson , which has a large width at @xmath61 due to the decay @xmath112 , would appear as a sharp resonance at high temperature since the amplitude for @xmath112 would vanish above @xmath0 ( near the critical point of the chiral phase transition ) . thus one can expect to have a better chance of seeing the @xmath68-meson emerge as a clear , narrow state at non zero temperature . in fact , since for @xmath113 the decay width of @xmath68 by strong interaction vanishes , only the electromagnetic width coming from @xmath114 is left @xcite .    within our model we are able to make some rough estimate of the temperature @xmath0 at which this phenomenon takes place by simply studying the behaviour of the meson masses ( [ emmepai2 ] ) , ( [ emmesig2 ] ) as a function of the temperature . however , since we are , as already said , mainly interested in determining the temperature @xmath0 at which @xmath1 , instead of evaluating separately @xmath105 and @xmath103 it will be sufficient to study the ratio @xmath115 @xmath116^{2}}~{\bar \sigma}_{s}^{2 } ~e^{\displaystyle{in\beta k_{0 } } } } \over\displaystyle { \langle{\bar\psi}\psi\rangle_t~-12~ { \sum_{n}(-)^{n}\int { d^{4}k\over ( 2\pi)^{4 } } \left[{{\bar \sigma}_{s}\over k^{2}+\left(m+{\bar \sigma}_{s}\right)^{2 } } ~-~{{\bar s}\over   s_{0}}~ { \sigma_{0s}\over k^{2}+\sigma^2_{0s}}\right]}~ e^{\displaystyle{in\beta k_{0 } } } } } \label{ratiosq}\ ] ] by straightforward calculations at the leading order in @xmath117 we obtain @xmath118^{2 } } \label{ratiosq1}\ ] ] furthermore , as we expect that the temperature @xmath0 will be not far from the critical temperature of the chiral phase transition , we can start our investigation by performing a mean - field expansion _ a - la _ landau of the effective potential in the small quark - mass limit around the critical point . following ref . @xcite we have @xmath119 \label{landau}\ ] ] where @xmath120 . if we now compute the masses of the scalar and pseudoscalar mesons from the expression ( [ landau ] ) of the effective potential we get @xmath121 by noticing that the generalization at finite @xmath7 of the normalization condition ( [ normcond ] ) on the effective potential reads now @xmath122 and that the coefficient @xmath123 is given by the expression ( see ref . @xcite ) @xmath124^{2 } } ~e^{\displaystyle{in\beta m y_{0 } } } \label{aquattro}\ ] ] it is easy to check that eq . ( [ landmasratio ] ) is nothing but eq . ( [ ratiosq1 ] ) . just to have a rough idea of the value of the temperature @xmath0 we can push our approximation further by evaluating @xmath123 and @xmath125 at @xmath126 and using for @xmath127 the expression for the massless minimum around @xmath126 @xcite @xmath128 in eq . ( [ landmasratio ] ) . the numerical values for the coefficients relevant for the evaluation of the mass ratio are @xmath129 and @xmath130 @xcite .    by now imposing the condition @xmath131 we have a constraint on the temperature ratio @xmath132 @xmath133 and this constraint fixes @xmath134 . the approximations we have made up to now ( we have evaluated @xmath123 and @xmath125 at @xmath126 , and we have chosen for @xmath127 the behaviour as given by eq . ( [ essezerot ] ) ) allow us to study @xmath9 only for @xmath7 very close to @xmath126 . to study the mass ratio @xmath9 as a function of @xmath7 for all the values of the temperature from @xmath61 up to @xmath7 very close to @xmath126 we have to go back to eq.([ratiosq1 ] ) which is valid for any @xmath7 at the leading order in the quark mass @xmath51 . to evaluate @xmath135 we can perform the integration in @xmath136 by using the cauchy theorem and summing up the series ; then , through a numerical calculation we can finally perform the integration left . in fig . [ fig1 ] we plot the mass ratio @xmath9 as a function of @xmath7 for a mass @xmath137 ; from the numerical analysis we can evaluate again the temperature @xmath0 @xmath138 and this constraint now fixes @xmath139 . we remark that this result for @xmath0 is in complete agreement with that previously derived by studying the landau expansion of the effective potential for @xmath7 very close to @xmath126 . finally , if one wants to study @xmath140 including all the orders in the quark mass @xmath51 and for any temperature it is necessary to go back to the general formula ( [ ratiosq ] ) . in fig . [ fig2 ] it is plotted the behaviour of the mass ratio as a function of @xmath7 as derived by the explicit evaluation of eq . ( [ ratiosq ] ) . in this case the temperature @xmath0 we are looking for is determined again by the condition @xmath141 which fixes @xmath142 . curiously this result exactly coincides with that obtained through the landau expansion and the evaluation of the coefficients @xmath123 , @xmath125 at @xmath126 and by putting for @xmath84 the expression ( [ essezerot ] ) as suggested by the mean - field theory . we have applied the composite operator formalism to qcd at finite temperature . after carrying out consistently the renormalization in the effective potential , we have calculated the masses of the scalar and pseudoscalar mesons . the relevant qcd parameters have been derived from the expressions for the masses and decay constants of the pseudoscalar multiplet , by inserting as inputs the mass and decay constant of the pion , the charged kaon mass and the electromagnetic mass difference between neutral and charged kaon .    a main purpose of the calculation was to determine the temperature @xmath0 at which @xmath1 . such temperature will be reached before the critical point , since the scalar meson mass @xmath105 is expected to decrease along with partial chiral symmetry restoration in hot matter , and at the same time the pion will lose its goldstone nature.the temperature @xmath0 is also expected not to lie far from the critical temperature of the chiral transition . in principle the @xmath68-meson would appear as a narrow state at high enough temperature . for @xmath113 its partial width by strong interaction has to vanish . the result of our calculation gives @xmath2 , in an approximation including current quark masses and with the full effective potential given by the formalism . we find that this result almost coincides with a much simpler treatment using a mean - field expansion of the effective potential _ a - la _ landau around the critical point in the limit of small quark mass . it is perhaps premature to discuss whether a possible experimental signature for the qcd phase transition may be inferred from the behaviour of the sigma width . the mass ratio itself is in any case a significant parameter to describe the degree of chiral symmetry breaking at different temperatures and therefore it has to be considered as a quantity of intrinsic theoretical interest . stephanov in the proceedings of workshop on qcd at finite baryon density : a complex system with a complex action , bielefeld , germany , 27 - 30 apr 1998 , hep - ph/9806509 ; j. berges . mit - ctp-2758 , talk given at workshop on qcd at finite baryon density : a complex system with a complex action , bielefeld , germany , 27 - 30 apr 1998 , hep - ph/9807362 k. rajagopal , talk given at 3rd workshop on continuous advances in qcd ( qcd 98 ) , minneapolis , mn , 16 - 19 apr , 1998 , hep - ph/9808348 , m.harada and a. shibata , ifp-762-unc , jul 1998 hep - ph/9807408 .    for general reviews and references see : h. satz , h. j. specht and r. stock , proceedings of the international conference `` quark matter 87 '' , z. phys . * c38 * , 1 ( 1988 ) ; m. jacob , proceedings of the international conference `` quark matter 88 '' , nucl . 498a * , 1 ( 1989 ) ; p. carruthers and j. rafelski , proceedings of the international conference `` hadronic matter in collision '' tucson arizona , world scientific , singapore ( 1989 ) ; e. v. shuryak , `` the qcd vacuum '' , world scientific , singapore ( 1988 ) ; d. gross , r. pisarski , and l. yaffe , rev . mod . phys . * 53 * , 43 ( 1981 ) ; j. cleymans , r. v. gavai , e. suhonen , phys . rep . * 130 * , 217 ( 1986 ) ; l. mclerran , rev . * 58 * , 1021 ( 1986 ) ; f. karsch , z. phys . * c38 * , 147 ( 1988 ) .    for reviews and references see : a. d. linde , rep . . phys . * 42 * , 389 ( 1979 ) ; b. sinha and s. raha , proceedings of the international conference on physics and astrophysics of quark gluon plasma , world scientific , singapore ( 1988 ) . j. c. collins and m. j. perry , phys . lett . * 34 * , 1353 ( 1975 ) ; p. d. morley and m. b. kislinger , phys . rep . * 51 * , 63 ( 1979 ) ; h. matsumoto , y. nakano and h. umezawa , phys . rev . * d38 * , 3211 ( 1984 ) ; h. nakkagawa , physica , * 158a * , 235 ( 1989 ) ; o. k. kalashnikov , z. phys . * 39c * , 427 ( 1988 ) . r. casalbuoni , in _ proceedings of the international symposium on composite models of quarks and leptons _ , tokyo , japan , 1985 , edited by h. terazawa and m. yasu ( institute for nuclear study , university of tokyo , tokyo , 1985 ) .
Abstract:

Answer: the composite operator formalism is applied to qcd at finite temperature to calculate the masses of scalar and pseudoscalar mesons . in particular the ratio of the sigma mass to the pion mass is an interesting measure of the degree of chiral symmetry breaking at different temperatures . we calculate the temperature @xmath0 at which @xmath1 , above which the sigma partial width into two pions vanishes . we find @xmath2 ( where @xmath3 is the critical temperature for the chiral phase transition ) , within the full effective potential given by the formalism . we find that an expansion _ a - la _ landau of the effective potential around the critical point in the limit of small quark mass provides for a very good determination of @xmath0 . ugva - dpt 1998/09 - 1015    firenze preprint - dff - 320/09/1998


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: error correcting codes are an important device for communicating information through noisy channels . we are living in what some social commentators have called the `` information age '' , with apparently ever increasing appetite for broadband internet , mobile phone data transmission , satellite television , computational power etc . hence , it is no surprise that such codes have are of significant practical use . a family of error - correcting codes based on an identification with the ever popular ising model from statistical mechanics were proposed by sourlas @xcite , and novel coding methods were examined @xcite . a further important advance , in terms of practical significance , was the move to finite rate , finite length sourlas type codes @xcite - which used interactions between @xmath0 bits and @xmath1 interactions per bit . while these codes showed good theoretical performance , the number of successfully transmitted bits was significantly higher ( for certain channel noises ) when @xmath2 . however , in this case the basin of attraction ( boa ) for a local decoding algorithm is very small meaning that decoding is computationally extremely challenging for reasonable code lengths . a more directly useful version of the above theory involves cascading sourlas codes @xcite . the basic idea is to use multiple values of @xmath0 in the interactions . for @xmath3 ( i.e. 2-body interactions ) there is an energy landscape dominated by the state which overlaps the sent message . once an algorithm has found this state with good overlap , higher body interactions can be invoked which increase this overlap substantially provided that the state is within the boa of the high overlap state . belief propagation ( bp ) @xcite , which is closely related to the bethe approximation and the cavity ( replica ) method in finite connectivity disordered systems @xcite , has been the principle algorithm used for decoding the finite rate sourlas codes including the case of the cascade version . sometimes the quenched disorder in this problem , in the form of corrupted bits which cause frustration in the graph of bit interactions , brings about a certain disruption to the bp dynamics , which can lead to deterioration of error correction performance . recently , a finite temperature survey propagation algorithm ( sp ) @xcite was applied to a biased 3-body sourlas code @xcite . the latter found that this more advanced algorithm did lead to a shift in the critical channel noise for computationally feasible decoding when bp did not converge . they also noted that averaging the beliefs in time , or time averaged belief propagation ( tabp ) @xcite could improve decoding in such cases as well . however , the reason for this improvement has not been fully clarified yet . the objective of this article is to pinpoint the major contributing factors for such improvement in practical decoding . for this , we will investigate how the overlap between the original and estimated messages in the first state of the cascading codes , for which only @xmath4 interactions are considered , are influenced by the employed decoding algorithms when the almeida - thouless ( at ) stability is broken @xcite . for simplicity of numerical experiments , we will mainly focus on the zero temperature cases . although the model that we will deal with does not directly accord with those in the preceding work @xcite , the two systems have a common feature that the at stability is broken and , therefore , insights gained from our simpler system is probably useful for understanding what occurs in the biased 3-body interaction system at finite temperature . by applying sp to the first stage of a cascading sourlas code , we will find that we can increase the overlap with the original message below the at - transition , which in turn improves the chance of being within the boa for good decoding once we introduce higher body interactions . it has also recently been shown that the first step of replica symmetry breaking ( 1rsb ) approach leads to higher values of the magnetisation for disordered ferromagnetically biased @xmath5 2-spin ising systems on finite connectivity random graphs @xcite . thus , we strongly speculate that the higher first stage overlap gained by introduction of the rsb ansatz in the cascading scheme is the main factor for the improvement of the error correcting performance . this speculation also suggests that algorithms based on increasing levels of replica symmetry breaking ( rsb ) would improve matters further , albeit at high computational cost . we make an initial investigation of this by a brief analysis of a 2rsb survey propagation ( 2sp : although this could rightly be called a survey of surveys propagation ) . we will also argue that tabp corresponds to survey propagation of a specific type for which the parisi parameter , which specifies the size of subgroups of replicas in the 1rsb ansatz , is set to zero and is never bettered by the simple bp algorithm , as is observed in cdma multiuser detection problems @xcite . this article is organized as follows . in the next section , we introduce the model that we will study . in sections 3 and 4 , two decoding schemes based on bp and sp are presented in the case of vanishing temperature . in sections 5 and 6 , we will show numerically that the sp based algorithm with zero parisi parameter ( sp0 ) is , in practice , equivalent to tabp and provides a better performance than bp when the at stability is broken . this relationship of sp0 to other algorithms is provided analytically in the vanishing code rate limit in appendix a. the general proof of this relationship , unfortunately , not been managed yet . we also show numerically that 2sp can lead to an improvement of the decoding performance in the situation where the at stability broken . this implies that a further improvement can be expected by introducing a higher level rsb ansatz although the computational cost for decoding will be sacrificed . the final section is devoted to a summary . we assume that we have some message source vector @xmath6 ( for simplicity of exposition we work within the ising formalism ) . we encode this source vector via a transformation to an @xmath7-bit vector @xmath8 where @xmath9 and where we have @xmath10 . the rate @xmath11 of such a code is given by @xmath12 . for each set @xmath13 , given its size , we populate it with indices at random from @xmath14 with each index appearing exactly once in any given set @xmath13 , and each index appearing @xmath15 times in each @xmath0 body interaction . we take the channel to be a binary symmetric channel so that when the encoded message vector @xmath8 is passed through the channel each bit is flipped independently with probability @xmath16 and we denote the received message by @xmath17 . the object of any decoding algorithm is to estimate @xmath18 with some @xmath19 using information from @xmath17 . a natural performance measure is the overlap between our estimate and the true message @xmath20 . one way to proceed with decoding is to look for the ground state of the following hamiltonian : @xmath21 since this hamiltonian is invariant under the transformation @xmath22 and @xmath23 we can gauge away the original message and are left with a multispin ferromagnetically biased @xmath5 spin glass model .    for the cascading approach we define a new auxiliary hamiltonian @xmath24 where @xmath25 if @xmath26 and @xmath27 otherwise . the initial step is to find minima of @xmath28 which is relatively easier than finding minima of @xmath29 . we then use this minima as an initial condition for local search algorithms looking for minima of the full hamiltonian ( [ eq : hamiltonian ] ) . an approximate method of finding the groundstate is available via belief propagation , which is equivalent to the cavity method at the level of replica symmetry for locally tree like graphs , which the present code is in the long code length limit . the same results can be found in the ensemble of graphs in the thermodynamic ( long code length ) limit by use of the replica method , however , since we are interested in decoding particular examples the former method is preferable . the method has been discussed a great deal in recent literature @xcite so we will be relatively brief . we define cavity fields @xmath30 to be the field on site @xmath31 ( decoded estimator @xmath32 ) in the absence of link @xmath13 ( received bit @xmath33 ) and the message @xmath34 to be the message from received bit @xmath33 to decoded bit @xmath32 . we are working with a discrete energy model and a self - consistent solution to these equations is given by restricting the values of @xmath35 to integers . then , the belief propagation ( bp ) update equations describing these variables are : @xmath36 where @xmath37 denotes the set @xmath38 and we take the convention @xmath39 . other than the trivial paramagnetic state ( @xmath40 ) , fixed points of the message passing equations ( [ eq : cavityupdates_rs ] ) correspond to local minima of the energy of the hamiltonian ( [ eq : hamiltonian ] ) . however , it is certainly not guaranteed that the equations ( [ eq : cavityupdates_rs ] ) will converge at all . if the messages converge , or after some specified number of updates if they do not , the estimator for the original message bits is given by @xmath41 consider a specific example , that of a code with @xmath42 and @xmath43 so each original message bit @xmath44 is contained within three 2-body interactions and three 3-body interactions . due to the 3-body interactions we find that without using the cascading approach we never find the state with high overlap - even at zero channel noise - due to the higher body interactions , which is similar to the result found for purely ferromagnetic 3-body interactions in @xcite , where even for these purely ferromagnetic interactions the structure of the interactions led to glassiness . however , if we first allow the system to equilibrate using the energy landscape given by only the 2-body interactions ( so perform bp on the auxiliary hamiltonian ) we find that for a certain range of @xmath45 we obtain non - zero values of @xmath7 . since for purely 2-body interactions there is a reflection symmetry in the hamiltonian there is no bias at this stage between @xmath46 . however , when the 3-body interactions are added to the hamiltonian it is possible to infer which of these two states gives a better overlap overall , namely the state with lower energy . bp is then performed on the hamiltonian with the 3-body interactions added . in practice , we add the the 3-body interactions in three batches , at each batch we add one 3-body interaction to the hamiltonian for each spin ( so @xmath47 interactions ) . since the initial state for the combined system is one with finite @xmath7 , it is closer to the ground state with high overlap than a random initial condition so bp has a better chance of reaching the ground state . past some critical channel noise we expect there to be an ergodicity breaking transition corresponding to replica symmetry breaking , at which point the update equations ( [ eq : cavityupdates_rs ] ) will no longer be sufficient to correctly characterise the system . for example , if we consider the initial stage of the algorithm where messages are passed purely between sites and 2-body links then it is known that for @xmath42 ( three 2-body interactions ) the at - transition occurs at @xmath48 @xcite . below this point a survey propagation algorithm @xcite should improve at least the initial overlap @xmath7 ( the at - point for the full system will not be the same as that for the initial system ) - at most until the transition from the mixed phase ( glassy ferromagnetic pure states ) to the spin glass phase . by considering the finite temperature replica symmetric transition lines away from the paramagnetic phase , a bifurcation analysis of the belief propagation equations yields critical temperatures for the @xmath42 regular sourlas code : @xmath49 which implies the triple point is at @xmath50 . according to the parisi - toulouse hypothesis of no re - entrance this gives us the zero temperature transition point between the mixed phase and the spin glass phase .    at the first level of replica symmetry breaking , survey propagation with parisi s replica symmetry breaking parameter equal to @xmath51 ( spx ) rather than solving for a single message from sites to links ( and links to sites ) we have to solve for a distribution of messages this distribution is over the messages in different ergodic sectors . while true ergodicity breaking will not occur in a finite length message , the breaking effects are sufficiently strong to mean that it does occur to all intents and purposes on any practical timescale for a message of reasonable length . the key idea behind the spx is that when performing an iteration of our messages , the energies in different ergodic sectors will vary , and as they vary the probability of finding such a given state also varies . the spx equations keep track of these changes assuming that the true nature of the states is 1rsb ( which is not necessarily the case ) . if we again restrict ourselves to integer values for the messages ( so in fact the integrals are really sums over the appropriate domain ) then the iteration equations are : @xmath52 \delta[h - \sum_b u_b ] \rme^{x    \left|\sum_b u_b\right|}\\ q_{a { { \rightarrow}}i}^t(u ) \propto \int [ \prod_{j \in a \setminus i } \rmd h_j p_{j    { { \rightarrow}}a}^t(h_j ) ] \delta[u - { { \mbox{sign}}}(j_{a } \prod_{j } h_j ) ] \rme^{x    \delta_{u,0 } } \label{eq : spx}\end{aligned}\ ] ] the above pair of equations are iterated until convergence ( or for some sufficient number of iterations ) . at that point we can construct our estimate for the decoded message via : @xmath53 \delta[h - \sum_b u_b ] \rme^{x    \left|\sum_b u_b\right|}\end{aligned}\ ] ] the approach thus far has still left the choice of @xmath51 . in physical problems , the general approach is to extremeise the generalized free energy wrt @xmath51 @xcite as these states are the ground states of the system . however , despite using a physical approach , we are interested in the error correcting abilities of this system and it does not necessarily hold that extremising the zero temperature free energy corresponds to maximising the overlap @xmath7 . in fact , we have found that minimizing the energy @xmath54 gives a good guide to choice of @xmath51 which is reasonable since our initial aim was to decode via minimizing the energy of the hamiltonian ( [ eq : hamiltonian ] ) . if the computational effort required to fit the model is prohibitive , choosing a particular value of @xmath51 for spx ( we found about 0.15 was reasonable ) still gives a marked improvement over bp in decoding across a range of channel noises @xmath45 . in figure [ fig : overlapvsx ] we plot variation in the overlap from its mean value to that of the energy and the free energy ( we take 10 times the variation of the free energy so that the scales are similar ) . the energy ( which is observable ) seems to give a reasonable indication of the overlap ( which is not observable in practical situations ) . ( 300,85 ) ( -45,15)=60 ( 0,6)(0,0)@xmath51 ( 55,15)=60 ( 95,6)(0,0)@xmath51    it is to be expected that to correctly describe state space up to the transition from a mixed ( glassy ferromagnetic ) phase to a spin glass phase full rsb needs to be taken into account . however , this is not feasible computationally . this does make the next obvious step a 2rsb algorithm , 2spxy , where the messages will be surveys of surveys . for a more detailed discussion on applications to @xmath55-sat see @xcite . we restrict our discussion for simplicity to the first stage of the decoding procedure ( where all interactions are 2-body ) . then we have update equations ; @xmath56 \propto \int \prod_{b \in \mathcal{n}(i )    \setminus a } \rmd \mathcal{q}_{b { { \rightarrow}}i}^{t}[q_{b { { \rightarrow}}i } ] z[\{q\ } ; y]^{x / y }    \delta[p - p[\{q\ } ; y]]\\ \mathcal{q}_{a { { \rightarrow}}i}^t[q ] = \int \rmd \mathcal{p}_{j { { \rightarrow}}a}^t[p ] \delta[q    - q[p;a]]\\ p[\{q\};y]= \frac{1}{z[\{q\};y ] } \int \prod_{b \in \mathcal{n}(i ) \setminus a } \rmd u_b    q_b(u_b ) \delta[h - \sum_b u_b ] \rme^{y ( |\sum_b u_b| - \sum_b |u_b|)}\\ q[p;a](u ) = \int \rmd h p(h ) \delta[u - { { \mbox{sign}}}(j_a h)]\end{aligned}\ ] ] where @xmath57 $ ] is the normalising constant for @xmath58 $ ] and @xmath51 and @xmath59 are the zero temperature parisi parameters . the integrals over @xmath60 are sums over the appropriate integers while the integrals over measures are over distributions on the integers . we will compare these different approaches numerically in the following . in general there is no guarantee that bp will converge on a graph with loops . however , experimentally it appears that for large disordered random graphs the transition to non - convergence is closely linked to the at - transition  @xcite . one can also rationalise this in terms of correlation lengths . if the correlation length is less than log(@xmath61 ) then locally the graph appears to be a tree ( apart from at a finite number of sites ) , so bp should converge provided the graph is sufficiently large ( of course this is very large as the length scale is log(@xmath61 ) ) . below the at transition due to the proliferation of pure states there are necessarily correlations on all length scales throughout the graph , so one would not expect bp to converge as it is affected by the loops . one simple , robust approach to deal with estimation when bp has not converged is to use time averaged bp ( tabp ) , which corresponds to taking averages in time of the messages , and performing inference with these average messages @xcite . ( 200,105 ) ( 0,10)=80 ( 57,6)(0,0)@xmath51 ( -10,50)(0,0)@xmath62    another algorithm which improves on bp but is simpler than full spx is sp0 , as used for the cdma problem in @xcite . sp0 does not resort to any reweighing of messages as the messages are passed but gains its improvement from a `` majority rules '' approach . i.e. when decoding in a region where bp does not converge , inference is performed at each bit by using the majority decision of the different instances . for the fully connected cdma problem , the local fields are gaussian and it is possible to show analytically that sp0 improves inference compared to bp analytically  @xcite . in the present situation due to the finite connectivity nature of the problem the effective fields are not gaussian distributed , thus we have not been able to show that sp0 improves on bp analytically ( basically we do not have sufficiently general control of @xmath63 ) . however , in the high @xmath1 limit our present model reduces to the sk - model @xcite and in this case we again find that sp0 is strictly better in inference terms than bp below the at - line , which is shown in appendix a.    in fact , tabp and sp0 are intuitively quite similar . the first averages over bp states in time , while the latter averages over different bp states in parallel from different initial conditions and with mixing of different solutions at each update step . provided the dynamics of bp are sufficiently chaotic that it explores all states one could expect both averages to give similar ( or even identical ) results . to test this hypothesis we compare the distribution of local fields given by sp0 to the distribution given by tabp . more specifically we examine the order parameter @xmath64\end{aligned}\ ] ] with the hellinger distance @xmath65 . in , figure [ fig : distances ] we look at @xmath62 vs @xmath51 and see that we find very good agreement for the hypothesis that sp0 is equivalent to tabp which gives an alternative explanation to its improved performance over bp and provides an new method for implementing sp0 . ( 200,105 ) ( 0,10)=80 ( 57,6)(0,0)@xmath45 ( -10,50)(0,0)@xmath7    we have used four different algorithms for the initial decoding , bp , tabp ( or sp0 ) , spx and 2spxy , corresponding to different levels of replica symmetry breaking and model fitting . for each algorithm there is essentially only one pertinent question which is whether the overlap reached after initial decoding is sufficient to lead to good full decoding . while it is tempting to think of some critical @xmath7 - in reality this is not sufficient to describe the state in full detail and other factors such as @xmath45 and the decoding algorithms themselves will all contribute to changes in the @xmath7 required ( in addition to finite size fluctuations ) to reach the high overlap state . however , due to issues of computational complexity we have tested all four algorithms on the initial decoding only . the 2spxy algorithm has only been investigated very briefly there , for a single value of @xmath51 and @xmath59 rather than optimising over them for each codeword received . we compare the results of these initial overlaps in figure [ fig : initial_mag ] . we see that each improvement in the algorithm does lead to an increased overlap after the initial stage , between the beginning of the mixed phase and the at point . the 2spxy does appear ( as far as we can tell from the data ) to give an improvement on the spx algorithm for a limited range of channel noise . next we look in figure [ fig : final_mag ] at the final overlap reached after full decoding with bp , sp0 and spx . we see that improving the quality of the algorithm leads to a marked improvement in the critical channel noise enlarging the space for which the cascading code approach will work . ( 200,105 ) ( 0,10)=80 ( 57,6)(0,0)@xmath45 ( -10,50)(0,0)@xmath7 we have examined how taking replica symmetry breaking into account can improve decoding of cascading sourlas codes . the basis behind our idea is that below the at transition , the true overlap at the minima of the hamiltonian ( [ eq : hamiltonian ] ) is above its replica symmetric estimate so increasing steps of replica symmetry breaking will serve to improve the overlap between the estimate of the message and the true message . we have also shown that tabp is closely related ( if not identical ) to sp0 - this throws up two obvious areas for further investigation , one to find under what conditions these algorithms are identical and secondly to show that both of them are strict improvements on bp below the at transition generally . another area to improve is in terms of the 2sp , where further optimisation of the parisi parameters should lead to further improvements , albeit at some considerable computational cost . the graph structure used here is not expected to be optimal - after the initial stage using higher body interactions ( i.e. @xmath66 ) should give improvements in the final overlap , optimising over possible architectures should also improve matters , but the principles developed here will help in any case where the channel noise is sufficiently high that decoding is taking place below the at transition . jplh would like to thank tommaso castellani and bastian wemmenhove for helpful discussion of their work . yk acknowledges support by grants - in - aid nos . 14084205 and 17340116 from mext / jsps , japan . here , we show analytically that sp0 is never bettered by bp in terms of the error correction performance when the number of connections tends to infinity . such a property is likely to hold in the case of sparse connectivity as well , analytical proof of which has not yet been managed .      for simplicity and generality , we focus on regular @xmath4 systems at general _ temperature @xmath67 . extension to other @xmath0 values and the case of vanishing temperature is straightforward . spx of this system can be expressed as @xmath68 \delta [ h - \sum_b u_b ] ( 2 \cosh ( \beta \sum_b u_b))^{\frac{x}{\beta } } , \label{spx_sk1 } \\ q_{a { { \rightarrow}}i}^t(u ) & \propto & \int [ \prod_{j \in a \backslash i } d h_j p_{j { { \rightarrow}}a}^t(h_j ) ] \delta [ u - u(h_j , j_a;\beta ) ] , \label{spx_sk2}\end{aligned}\ ] ] where @xmath69 . in order to obtain a simpler expression for @xmath70 , we introduce appropriate scalings @xmath71 and @xmath72 , where @xmath73 and @xmath74 . this corresponds to situations of low code rate communication that copes with strong noise despite restricted signal power . let us denote @xmath75 , @xmath76 and characterize the system utilizing two macroscopic variables @xmath77 and @xmath78 based on the law of large numbers . these , in conjunction with the central limit theorem , indicate that eq . ( [ spx_sk1 ] ) is reduced to @xmath79 where @xmath80 in turn , eq . ( [ spx_reduced1 ] ) provides expressions of microscopic variables @xmath81 eqs . ( [ spx_sk_simple1])-([spx_reduced2_2 ] ) constitute the simpler expression of spx in the dense connectivity limit . utilizing @xmath82 and @xmath83 , the estimate for the decoded message at the @xmath84th update is provided as @xmath85 where @xmath86 . the performance of the spx decoding algorithm ( [ spx_sk_simple1])-([spx_reduced2_2 ] ) can be examined by following time evolution of several relevant macroscopic variables  @xcite . the maximum entropy principle , in conjunction of the law of large numbers , makes it possible to regard the microscopic variable @xmath87 as gaussian random numbers , the average and variance of which are provided independently of site and check indices @xmath31 , @xmath13 as @xmath88 and @xmath89 , respectively . we assume that the gauge transformation @xmath22 and @xmath90 is already carried out , which implies that the original message is specified by @xmath91 . these indicate that @xmath92 , @xmath93 and the macroscopic overlap to the original message @xmath94 are represented by @xmath95 respectively , where @xmath96 and @xmath97 . on the other hand , the self - averaging property  @xcite and eqs . ( [ 1rsb_q1 ] ) , ( [ 1rsb_q0 ] ) and ( [ 1rsb_m ] ) imply that @xmath98 and @xmath99 are provided by @xmath100 respectively , where @xmath101 denotes the configuration average with respect to the channel noise and lattice configuration . the update rule for @xmath102 is already provided by eq . ( [ spx_sk_simple2 ] ) . using these macroscopic variables , the error probability per bit for spx , @xmath103 , is calculated from @xmath104 where @xmath105 and @xmath106 for @xmath107 and @xmath108 , respectively . it may be remarked that eqs . ( [ spx_sk_simple2 ] ) and ( [ 1rsb_q1])-([1rsb_e ] ) accord with the macroscopic description of the spx dynamics for the sk model  @xcite . now , we focus on two specific parisi parameters @xmath110 and @xmath111 . a distinctive property of these two parameter choices is that only four out of the six variables are relevant in eqs . ( [ spx_sk_simple2 ] ) and ( [ 1rsb_q1])-([1rsb_e ] ) , which can be reduced to an identical set of four equations @xmath112 where @xmath113 , @xmath114 for @xmath110 while @xmath115 , @xmath116 for @xmath111 . as @xmath83 is a non - negative variable , this implies that @xmath117 , which indicates the variance of the local fields , for @xmath111 can not be smaller than that for @xmath110 while @xmath118 , which represents the signal strength of the original message in the local fields , is identical in both cases . applying this result to eq . ( [ bit_error_rate ] ) leads to a conclusion that sp0 does not provide worse performance than sp@xmath109 .    in terms of the decoding problem , sp@xmath109 and bp provide identical performance . therefore , the above argument indicates that sp0 is never overcome by bp with respect to the error correction performance . in order to show this , we characterise the auxiliary distributions @xmath119 and @xmath120 utilising the first moments of @xmath121 and @xmath122 , respectively , as @xmath123 although higher moments are necessary for completely describing general spx dynamics , eqs . ( [ spx_sk1 ] ) and ( [ spx_sk2 ] ) indicate that for the specific choice of @xmath111 , @xmath124 and @xmath125 constitute a closed update rule as @xmath126 which is nothing but the bp dynamics of finite temperature @xmath67 and are reduced to eqs . ( [ eq : cavityupdates_rs0 ] ) and ( [ eq : cavityupdates_rs ] ) in the zero temperature limit @xmath127 . this implies that as long as only the first moment is concerned , which is the case for the current decoding problem , sp@xmath109 and bp are equivalent .
Abstract:

Answer: we investigate how insights from statistical physics , namely survey propagation , can improve decoding of a particular class of sparse error correcting codes . we show that a recently proposed algorithm , time averaged belief propagation , is in fact intimately linked to a specific survey propagation for which parisi s replica symmetry breaking parameter is set to zero , and that the latter is always superior to belief propagation in the high connectivity limit . we briefly look at further improvements available by going to the second level of replica symmetry breaking .