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: collective phenomena in systems far from thermal equilibrium have been of considerable interest in recent years @xcite . unlike systems in thermal equilibrium where the gibbs picture provides a theoretical framework within which such phenomena can be studied , here no such framework exists and one has to resort to studies of specific models in order to gain some understanding of the phenomena involved .    one class of such models is driven diffusive systems ( dds ) @xcite . driven by an external field these systems do not generically obey detailed balance so that the steady state has non - vanishing currents . theoretical studies of dds have revealed basic differences between systems in thermal equilibrium and systems far from thermal equilibrium . for example , it is well known that one dimensional ( @xmath0 ) systems in thermal equilibrium with short - range interactions do not exhibit phenomena such as phase transitions , spontaneous symmetry breaking ( ssb ) and phase separation ( except in the limit of zero temperature or in the context of long - range interactions ) @xcite . in contrast , some examples of noisy @xmath0 dds with local dynamics have been found to exhibit such phenomena . one example of a noisy system which exhibits ssb in @xmath0 is the asymmetric exclusion model of two types of charge studied in @xcite . in this model , two types of charge are biased to move in opposite directions on a @xmath0 lattice with open ends . the charges interact via a hard - core interaction , and are injected at one end of the lattice and ejected at the other end . this model is symmetric under the combined operations of charge conjugation and parity ( pc symmetry ) . however , this symmetry is broken in the steady state , where the currents of the two charges are not equal . the reason for symmetry breaking in this model lies to some extent in the open boundaries . other examples of models in which there is ssb in @xmath0 have also been found in the context of cellular automata @xcite and surface growth @xcite . in the latter , ssb was due to the fact that one of the rates for a local dynamical move in the models is zero . once this zero rate changes to a non - zero rate ssb disappears . a closely related problem to spontaneous symmetry breaking , is that of _ phase separation _ in @xmath0 noisy systems . this has been observed in driven diffusive models with inhomogeneities , such as defect sites @xcite or particles @xcite . in these models it has been found that macroscopic regions of high densities are formed near the defect , much like a high density of cars behind a slow car in a traffic jam @xcite . here the phase separation is triggered by the defects . it is of interest to study whether phase separation can occur in @xmath0 noisy homogeneous systems such as on a ring geometry with no defects , where _ all _ possible local transition rates which are consistent with the symmetry and conservation laws of the model are non vanishing . recently , lahiri and ramaswamy have introduced a lattice model in the context of sedimenting colloidal crystals , where phase separation is found to take place without any inhomogeneities @xcite . in this model , there are two rings coupled to each other and particles on each ring undergo an asymmetric exclusion process . the hopping rate between sites @xmath1 and @xmath2 on each ring depends on the occupation at the @xmath1th site on the other ring . however , this model is studied mainly using monte carlo simulations and no analytical results are available so far .    in a recent letter @xcite we introduced a simple three - species driven diffusive model exhibiting phase separation and spontaneous breaking of the translational symmetry on a ring . in the model nearest - neighbor particles exchange with given rates and the numbers of each species are conserved under the dynamics . the rates of all local dynamical moves that obey the conservation laws are non zero . an argument indicating that generically the system phase separates , thus breaking the translational invariance , was given for the case when none of the species of particle has zero density . in the special case of equal number of particles of each type , it was shown that the local dynamics obeys detailed balance with respect to a long - range asymmetric ( chiral ) hamiltonian . in this special case , using the hamiltonian , we have found the steady state of the model exactly and have been able to prove the existence of phase separation analytically . the existence of a hamiltonian for this special case is of interest in the light of speculation that non - equilibrium systems exhibiting generic long - range correlations might be described by effective hamiltonians containing long - range interactions @xcite . here we explicitly demonstrate that for the special case where the three densities are equal the model is exactly described by a long - range asymmetric hamiltonian . the model not only has long - range correlations but has generic long - range _ order_. the mechanism found in this study suggests that systems with dynamical rules defined completely locally and _ a priori _ without respect to any hamiltonian , may have a steady state where the configuration space is sampled according to a measure that is intrinsically global . the hamiltonian also allows us to identify the analog of a temperature in the microscopic dynamics as related to the drive of the system ; for zero drive , that is symmetric diffusion of the particles , the effective temperature is infinite and phase separation is lost . we note that a related but distinct three - species model has recently been introduced by arndt _ et al_. this model also exhibits phase separation@xcite .    in the present work we analyze in detail the @xmath3 species model which was introduced in @xcite and then generalize it to larger @xmath4 . we provide the complete proof of phase separation which follows from the exact calculation of the partition sum in the thermodynamic limit . we also provide numerical evidence of phase separation in the general case where the densities of the three particles are not equal .    in order to study the coarsening process monte - carlo simulations are performed . however , simulation of the microscopic model is hampered by slow dynamics which makes it difficult to access the scaling regime . the system becomes trapped in metastable states comprising several domains of each type of particle . the number of metastable states is exponentially large in the system size . the lifetimes of the metastable states increase exponentially with the average domain size as the fully phase separated state is approached . thus the model provides an example of slow dynamics in a system without any quenched disorder @xcite .    to ameliorate the difficulty of numerically studying such slow dynamics we employ a toy model wherein it is the domains that are updated rather than the individual particles this allows the long - time scaling behavior of the domain size to be investigated and to confirm a logarithmic growth of the average domain size with time . the toy model also affords a mean - field solution for the long - time dynamical behavior , that again confirms the scaling behavior .    returning to the case of equal numbers of particles of different species it is of interest to investigate the steady - state behavior in finite - size systems . we have found it convenient to do this by employing a matrix product technique previously used to solve the steady state of asymmetric exclusion processes @xcite . however , in the case of three species the simplest form of this technique @xcite is applicable only to a limited class of systems @xcite . for the present model we generalize the matrix product to a product of rank 6 tensors and write the steady state by taking an appropriate contraction . the partition sum and steady - state correlation functions can be conveniently computed numerically using this tensor product ansatz . the paper is organized as follows : in section ii we define the model introduced in @xcite and we present an argument which indicates that the system should phase separate as long as none of the species of particles has zero density . in section iii we study the special case where the model satisfies detailed balance and explicitly write down the steady - state weight for the three - species model . the existence of phase separation in the model for any non - infinite temperature is proved analytically by calculating some bounds on the two - point correlation functions . section iv contains numerical evidence for phase separation in the general case where the densities of the three species of particles are not equal . the toy model , which facilitates efficient monte carlo simulations , is used to study the dynamics of phase separation . a mean - field analysis of this toy model is presented , the details being left to appendix a. in section v we present results for finite systems obtained via the tensor product ansatz . using these results we study finite - size scaling in the system . in section vi we address phase separation in systems with more than three species of particle and a proof of detailed balance for special cases is given in appendix b. we conclude in section vii and discuss some open questions . we start by defining a three - species model which exhibits phase separation in @xmath0 . consider a one - dimensional , ring - like ( periodic ) lattice of length @xmath5 where each site is occupied by one of the three types of particles , @xmath6 , @xmath7 , or @xmath8 . the model evolves under a random sequential update procedure which is defined as follows : at each time step two neighboring sites are chosen randomly and the particles at these sites are exchanged according to the following rates @xmath9 the particles thus diffuse asymmetrically around the ring . the dynamics conserves the number of particles , @xmath10 and @xmath11 of the three species . the @xmath12 case is special . here the diffusion is symmetric and every local exchange of particles takes place with the same rate as the reverse move . the system thus obeys detailed balance reaching a steady state in which all microscopic configurations ( compatible with the number of particles @xmath10 and @xmath11 ) are equally probable . this state is homogeneous , and no phase separation takes place . we now present a simple argument suggesting that for @xmath13 the steady state of the system is not homogeneous in the thermodynamic limit . for simplicity the case @xmath14 is examined . as a result of the bias in the exchange rates an @xmath6 particle prefers to move to the left inside a @xmath7 domain and to the right inside a @xmath8 domain . similarly the motion of @xmath7 and @xmath8 particles in foreign domains is biased . consider the dynamics starting from a random initial configuration . the configuration is composed of a random sequence of domains of @xmath6 , @xmath7 , and @xmath8 particles . due to the bias a local configuration in which an @xmath6 domain is placed to the right of a @xmath7 domain is unstable and the two domains exchange places on a relatively short time scale which is linear in the domain size . similarly , @xmath15 and @xmath16 domains are unstable too . on the other hand @xmath17 , @xmath18 and @xmath19 configurations are stable and long lived . thus after a relatively short time the system reaches a state of the type @xmath20 in which @xmath21 and @xmath8 domains are located to the right of @xmath22 and @xmath7 domains , respectively . the evolution of this state takes place via a slow diffusion process in which , for example , the time scale for an @xmath6 particle to cross an adjacent @xmath7 domain is @xmath23 , where @xmath24 is the size of the @xmath7 domain . the system therefore coarsens and the average domain size increases with time as @xmath25 @xcite . eventually the system phase separates into three domains of the three species of the form @xmath26 .    in a finite system the phase separated state may further evolve and become disordered due to fluctuations . however , the time scale for this to happen grows exponentially with the system size . for example it would take a time of order of @xmath27 for the @xmath6 domain in the totally phase separated state to break up into smaller domains . hence in the thermodynamic limit , this time scale diverges and the phase separated state remains stable provided the density of each species is non - zero . note that there are always small fluctuations about a totally phase separated state . however , these fluctuations affect the densities only near the domain boundaries . they result in a finite width for the domain walls . the fact that any phase separated state is stable for a time exponentially long in the system size amounts to a breaking of the translational symmetry . since the exchange rates are asymmetric , the system generically supports a particle current in the steady state . to see this , consider the @xmath6 domain in the phase separated state . an @xmath6 particle near the @xmath28 boundary can traverse the entire @xmath7 domain to the right with an effective rate proportional to @xmath29 . once it crosses the @xmath7 domain it will move through the @xmath8 domain with rate @xmath30 . similarly an @xmath6 particle near the @xmath31 boundary can traverse the entire @xmath8 domain to the left with a rate proportional to @xmath32 . once the domain is crossed it moves through the @xmath7 domain with rate @xmath30 . hence the net @xmath6 particle current is of the order of @xmath33 . since this current is exponentially small in system size , it vanishes in the thermodynamic limit . for the case of @xmath34 , this argument suggests that the current is strictly zero for any @xmath5 . in sections iii and v we study this case in detail . the arguments presented above suggesting phase separation for @xmath35 may be easily extended to @xmath36 . in this case , however , the phase separated state is @xmath37 rather than @xmath38 . this may be seen by noting that the dynamical rules are invariant under the transformation @xmath39 together with @xmath40 . in this section we show that the dynamics ( [ eq : dynamics ] ) , for the special case @xmath34 , satisfies detailed balance . the corresponding hamiltonian , which determines the steady - state distribution , is found to have long - range asymmetric interactions . using this hamiltonian , we calculate analytically the partition sum and bounds on the correlation functions in the thermodynamic limit . these are then used to prove the existence of phase separation in the model . later , in section [ matrixans ] we study finite systems for this case and the approach to the thermodynamic limit . the general argument presented in the previous section suggests that for the special case @xmath34 , the steady state carries no current for any system size . we demonstrate this explicitly by showing that the _ local dynamics _ of the model satisfies detailed balance with respect to a _ long - range asymmetric _ hamiltonian @xmath41 . we define the occupation variables @xmath42 and @xmath43 as follows : @xmath44 the variables @xmath45 and @xmath43 are defined similarly . clearly the relation @xmath46 is satisfied . a microscopic configuration is thus described by a set @xmath47 . using these variables , we will show that the hamiltonian @xmath41 and the steady - state distribution @xmath48 corresponding to the dynamics ( [ eq : dynamics ] ) for the case @xmath49 are given by @xmath50 \ ; , \ ] ] @xmath51 here @xmath52 is the partition sum given by @xmath53 , where the sum is over all configurations in which @xmath34 . note that the hamiltonian @xmath41 does not determine the dynamics of the system , it just governs the steady - state distribution as given in ( [ eq : hamilton ] ) and ( [ eq : weight ] ) . ( [ eq : weight ] ) suggests that @xmath54 serves as a temperature variable with @xmath55 . thus , @xmath56 is the infinite - temperature limit . the hamiltonian ( [ eq : hamilton ] ) is written in a form which is not manifestly translationally invariant . however , careful examination reveals that when the relation @xmath34 is taken into account , the hamiltonian as given by ( [ eq : hamilton ] ) is indeed translationally invariant ( see appendix b ) . therefore site @xmath57 may be chosen arbitrarily . an expression for @xmath41 which is manifestly translationally invariant will be derived at the end of this section . note that @xmath58 since the lhs yields the number of @xmath19 ( and @xmath15 ) pairs in the system . using this relation the hamiltonian may also be written in a form where the cyclic symmetry is more apparent : @xmath59-(n/3)^2 \ ; .\ ] ]    the proof of eqs . ( [ eq : hamilton],[eq : weight ] ) is straightforward . this is done by considering a nearest - neighbor particle exchange and verifying that detailed balance is satisfied with respect to ( [ eq : weight ] ) . we start by considering nearest - neighbor sites in the interior of the lattice , namely pairs other than @xmath60 . for example consider the exchange @xmath61 taking place at two adjacent sites @xmath62 and @xmath63 , where @xmath64 . this exchange results in the contribution of one more @xmath65 term in @xmath66 and hence the energy of the resulting configuration is higher by @xmath57 . it is easy to see using eq . ( [ eq : weight ] ) that @xmath67 , as required by detailed balance . similar relations are easily derived for exchange of @xmath18 and @xmath19 pairs . now consider an exchange taking place between sites @xmath5 and @xmath57 , say @xmath68 . according to ( [ eq : hamilton ] ) this exchange costs an energy of @xmath69 . therefore the exchange satisfies the detailed balance condition @xmath70 only when @xmath71 . similarly , by considering the exchanges @xmath72 and @xmath73 , one deduces that the detailed balance condition is satisfied for any exchange at sites @xmath5 and @xmath57 as long as @xmath34 . in appendix b we consider the most general nearest - neighbor exchange rates for @xmath4 species and arbitrary densities and derive conditions ( [ condition ] ) for exchange rates which satisfy detailed balance .    to write @xmath41 in a manifestly translationally invariant form we define @xmath74 as the hamiltonian in which site @xmath75 is the origin . namely , @xmath76 \ ; , \ ] ] where the summation over @xmath1 and @xmath77 is modulo @xmath5 . summing ( [ hi0 ] ) over all @xmath75 and dividing by @xmath5 , one obtains , @xmath78 where in the summation the value of the site index @xmath79 is modulo @xmath5 . in the hamiltonian ( [ eq : h_tinvar ] ) the interaction is linear in the distance between the particles , and thus is long - ranged . the distance is measured in a preferred direction from site @xmath1 to site @xmath80 . moreover it is asymmetric in the sense that @xmath66 is not invariant under the parity operation . it is also related to chiral hamiltonians @xcite .      before proceeding further to evaluate the partition sum and some correlation functions associated with the hamiltonian ( [ eq : hamilton ] ) let us make a few observations . the ground state of the hamiltonian is given by the fully separated state @xmath81 and its translationally related states . the degeneracy of the ground state is thus @xmath5 and its energy is zero . a simple way of evaluating the energy of an arbitrary configuration is obtained by noting that nearest - neighbor ( nn ) exchanges @xmath82 and @xmath83 cost one unit of energy each while the reverse exchanges result in an energy gain of one unit . the energy of an arbitrary configuration may thus be evaluated by starting with the ground state and performing nn exchanges until the configuration is reached , keeping track of the energy changes at each step of the way . the highest energy is @xmath84 and it corresponds to the totally phase separated configuration @xmath85 and its @xmath5 translations .    in considering the excited states of the hamiltonian ( [ eq : hamilton ] ) we note that the model exhibits a set of metastable states which correspond to local minima of the energy : any exchange of nn particles results in an increase of the energy . in these states no @xmath86 and @xmath15 nn pairs exist ; only @xmath87 and @xmath19 nn pairs may be found in addition to @xmath88 , @xmath89 , and @xmath90 . any metastable state is thus composed of a sequence of domains in which @xmath21 and @xmath8 domains follow @xmath22 and @xmath7 domains , respectively . therefore each metastable state has an equal number of domains , @xmath91 , of each type with @xmath92 . the @xmath93 case corresponds to the ground state while @xmath94 corresponds to the @xmath95 state , composed of a total of @xmath5 domains each of length @xmath57 . ( the total number of domains in an @xmath91-state is @xmath96 . )    for calculating the free energy and some correlation functions corresponding to the hamiltonian ( [ eq : hamilton ] ) we find it useful to first derive some bounds for the number @xmath97 of @xmath91-states and their energies . in the following such bounds are presented . they are then used , in the next section , to evaluate the free energy and correlation functions of the model .    to obtain a bound for @xmath97 we note that the number of ways of dividing @xmath98 @xmath6 particles into @xmath91 domains is @xmath99 . the number of ways of combining @xmath91 divisions of each of the three types of particles is clearly @xmath100 ^ 3 $ ] . there are at most @xmath5 ways of placing this string of domains on a lattice to obtain a metastable state ( the number of ways need not be equal to @xmath5 since the string may possess some translational symmetry ) . one therefore has @xmath101 ^ 3 \le{\cal n}(s ) \le n \left[{n/3 - 1}\choose { s-1}\right]^3\ ; . \label{nmss}\ ] ] thus , the total number of metastable states is exponential in @xmath5 . we now consider the energy of the metastable states . it is easy to convince oneself that among all @xmath91-states , none has energy lower than the following configuration , @xmath102 where the @xmath103 rightmost domains are of size @xmath57 and the three leftmost domains are of size @xmath104 each . the energy of this state , @xmath105 satisfies the following recursion relation @xmath106 with @xmath107 . to see this one notes that the @xmath91-state may be created from the @xmath108-state by first moving a @xmath7 particle from the leftmost @xmath7 domain across @xmath109 @xmath8 particles to the right and then moving an @xmath6 particle from the leftmost @xmath6 domain to the right across the adjacent @xmath7 and @xmath8 domains . the energy cost of these moves is @xmath109 , yielding ( [ rr ] ) . the recursion relation ( [ rr ] ) , together with @xmath110 , is then readily solved to give @xmath111 the energy of all metastable @xmath91-states is larger or equal to @xmath105 as given by eq.([emss ] ) . in the following section we use the bounds ( [ nmss ] ) and ( [ emss ] ) to calculate the partition sum and some correlation functions corresponding to the hamiltonian ( [ eq : hamilton ] ) . in this section we prove that , in the large @xmath5 limit and for all @xmath35 , the partition sum is given by @xmath112 ^ 3 \ ; , \label{psum}\ ] ] where @xmath113 the partition sum for @xmath36 may be obtained by replacing @xmath54 by @xmath114 in ( [ psum ] ) . note that the partition sum is linear and not exponential in @xmath5 , meaning that the free energy is not extensive . this is a result of the long - range interaction in the hamiltonian and the fact that the energy excitations are localized near the domain boundaries , as will be shown in the following . for @xmath54 close to @xmath57 , @xmath115 has an essential singularity , @xmath116}\;.\ ] ] this suggests that extensivity of the free energy could be restored in the double limit @xmath117 and @xmath118 with @xmath119 finite . this scaling behavior is only suggestive since expression ( [ psum ] ) may not be valid in this limit . for @xmath12 all configurations with @xmath34 are equally probable , so that the partition sum is given by @xmath120 which goes like @xmath121 for large @xmath5 . it is instructive to first present the proof of ( [ psum ] ) for small values of @xmath54 . this proof will then serve as the basis for the proof for any @xmath35 . we start by noting that in calculating the partition sum ( [ psum ] ) , configurations with energy larger then @xmath122 ( @xmath123 ) may be neglected in the thermodynamic limit . for simplicity we first demonstrate this for @xmath124 , although later we show it for any @xmath35 . the contribution to the partition sum from these energy states , @xmath125 , is given by @xmath126 where @xmath127 is the number of configurations of energy @xmath128 . clearly , the number of possible configurations in the system is bounded crudely from above by @xmath121 . this bound implies @xmath129 thus , for @xmath124 , the contribution to the partition sum arising from energies larger than @xmath122 is exponentially small in @xmath5 and may be neglected in the thermodynamic limit . the calculation of @xmath52 is thus reduced to calculating a truncated partition sum in which only energies up to @xmath122 are summed over . to proceed we consider @xmath130 and take into account configurations with energy less than @xmath131 . this simplifies the calculations considerably since all configurations with energy @xmath132 may be decomposed into @xmath5 disjoint sets of states , each corresponding to a unique ground state ( see previous section for a discussion of ground states ) . we label the sets by @xmath133 , the position of the rightmost @xmath6 particle in the @xmath6 domain of the ground state belonging to the set ( see fig . [ ground ] for an example of an @xmath134 ground state ) . each state in a specific set can be obtained from the corresponding ground state by exchanging nearest neighbors so that the energy always increases along the intermediate states . note that this is correct only if excitations of energy less than @xmath131 are considered . this is because not all higher energy states can be reached by uphill steps from a ground state . thus , using translational invariance , the partition sum can be written as @xmath135 where @xmath136 is the truncated partition sum of one of the @xmath5 sets of configurations . we proceed to calculate @xmath136 . this is done by considering all the possible energy excitations with energy less than @xmath131 above one ground state . consider a specific domain boundary , say @xmath17 . excitations of energy @xmath128 at this boundary can be created by moving one or more @xmath6 particles into the @xmath7 domain ( this is equivalent to moving @xmath7 particles into the @xmath6 domain ) . an @xmath6 particle moving into the @xmath7 domain is considered as a walker . the excitation energy increases linearly with the distance the walker has moved . thus , in this picture an excitation of energy @xmath128 is created by @xmath137 walkers , traveling a total distance @xmath128 . the number of excitations of energy @xmath128 is then given by the number of ways , @xmath138 , of partitioning an integer @xmath128 into the sum of a sequence of non - increasing positive integers . taking into account excitations at all three boundaries , an excitation of energy @xmath128 in the system is created by three independent excitations of energy @xmath139 and @xmath140 at the different domain boundaries such that @xmath141 . the number of excitations of this form is just given by @xmath142 . thus , @xmath136 is given by @xmath143 taking the thermodynamic limit we obtain @xmath144 using a well known result from number theory , attributed to euler , for the generating function of @xmath138 @xcite , @xmath145 and using ( [ thelim ] ) and ( [ exp ] ) , eq.([psum ] ) is obtained . so far we have proved that for @xmath146 , eq.([psum ] ) is exact in the thermodynamic limit . we now extend these results for any @xmath35 . first we have to show that the states ignored in the previous calculation for @xmath146 may be ignored for all @xmath35 . to do this we calculate upper and lower bounds on @xmath52 and show they converge for large enough @xmath5 .    for this we have to consider the entire energy spectrum of the hamiltonian . any configuration of the system which is neither a ground state nor a metastable state can be obtained from at least one ground state ( @xmath93 ) or a metastable state ( @xmath147 ) as follows : starting from this @xmath91-state exchange nearest neighbors such that the energy always increases along the path until the configuration is reached . in what follows it is demonstrated that none of the configurations which can be obtained from @xmath91-states , with @xmath147 , by the above procedure of particle exchange , contributes to the partition sum in the thermodynamic limit . an upper bound on the partition sum may be calculated as follows : using the same steps of derivation used for computing @xmath136 , it is straightforward to show that the contribution to the partition sum from an @xmath91-state and associated configurations is _ at most _ @xmath148^{-3s}$ ] . the prefactor @xmath149 arises from the minimum energy ( [ emss ] ) of this metastable state . therefore by considering the contributions from all the @xmath91-states and using ( [ nmeta ] ) the following bound is found @xmath150 ^ 3 + \sum_{s=2}^{n/3}n { n/3 - 1\choose s-1 } q^{(s-1)n/3 -s(s-1)/2 } [ ( q_{\infty})]^{-3s}\;. \label{zupperb1}\ ] ] the second term on the rhs represents the contribution from excitations around the metastable states . replacing @xmath151 by an upper bound @xmath152 one can sum the binomial series . the resulting expression is exponentially small in @xmath5 for any @xmath35 . a lower bound on @xmath52 can be calculated by neglecting configurations with energy greater than @xmath131 as follows , @xmath153 ^ 3 - n \sum_{m = n/3 - 1}^{\infty}q^m\sum_{m_i=0}^mp(m_1)p(m_2)p(m_3)\delta_{m_1+m_2+m_3,m } \\ & > & n/[(q)_{\infty}]^3 - n \sum_{m = n/3 - 1}^{\infty}q^m ( mp(m))^3 \ ; . \label{lowbound}\end{aligned}\ ] ] the asymptotic behavior of @xmath138 @xcite is given by @xmath154 thus , for large @xmath5 the lower bound ( [ lowbound ] ) converges to ( [ psum ] ) as does the upper bound ( [ zupperb1 ] ) . whether or not a system has long - range order in the steady state can be found by studying the decay of two - point density correlation functions . for example the probability of finding an @xmath6 particle at site @xmath1 and a @xmath7 particle at site @xmath77 is , @xmath155 where the summation is over all configurations @xmath156 in which @xmath34 . due to symmetry many of the correlation functions will be the same , for example @xmath157 . a sufficient condition for the existence of phase separation is @xmath158 since @xmath159 we wish to show that @xmath160 . in fact we will show below that for any given @xmath161 and for sufficiently large @xmath5 , @xmath162 this result not only demonstrates that there is phase separation , but also that each of the domains is pure . namely the probability of finding a particle a large distance inside a domain of particles of another type is vanishingly small in the thermodynamic limit .    to prove eq . ( [ corrf3 ] ) , we use the relation @xmath163 , and show that the correlation function @xmath164 is of @xmath165 and @xmath166 is of @xmath167 . here we show only the proof for @xmath164 , since the proof of @xmath166 is similar . we also restrict ourselves to @xmath168 , which is sufficient for proving eq . ( [ corrf3 ] ) . we have already seen that the contribution to the partition sum from the metastable states and excitations above them are exponentially small in the system size and hence may be neglected . therefore , for calculating the correlation function it is sufficient to consider the @xmath5 ground states and excitations above them which may be reached by moves which only increase the energy . as we have seen , these states form @xmath5 disjoint sets of states , each associated with one of the ground states . using this we now show that @xmath169 . for this purpose we use a restricted partition sum @xmath170 , which is defined as the partition sum @xmath136 calculated with the constraint that one of the walkers , say of type @xmath6 , has traveled at least distance @xmath91 . it is given as @xmath171 by @xmath172 here @xmath173 is the number of partitions of integer @xmath128 with the constraint that in all the partitions the integer @xmath91 occurs at least once . noting that @xmath174 it is easy to show that @xmath175 where @xmath176 . we now proceed to derive a bound for @xmath177 . recall that @xmath24 is the position of the rightmost @xmath6 particle in the @xmath6 domain in the ground state labeled @xmath24 . if we define @xmath178 as the correlation function calculated within the set of states labeled @xmath24 , we can write @xmath179 up to exponentially small corrections in the system size . for convenience we break the summation over @xmath24 into @xmath180 sums according to the values of @xmath181 and @xmath182 in the ground state . these @xmath180 parts correspond to ( i ) ground states where @xmath183 ; ( ii ) ground states where @xmath184 ; ( iii ) ground states where @xmath185 and ( iv ) ground states where @xmath186 . we now consider each of them in detail and give an upper bound for @xmath187 in each case . \(i ) _ ground states where @xmath183 _ : in this case the site @xmath57 is inside the @xmath6 domain and site @xmath161 is inside the @xmath7 domain . since we consider only @xmath168 , these states correspond to the @xmath24 ground states with @xmath188 . using the fact that @xmath189 one finds @xmath190 ( ii ) _ ground states where @xmath184 _ : in principle site @xmath161 might be either inside the @xmath6 domain or inside the @xmath8 domain . however , since site @xmath57 is in the @xmath6 domain and we consider only @xmath191 , site @xmath161 must be in the @xmath6 domain . the ground states @xmath24 for which this takes place satisfy @xmath192 . clearly , only the excited states where @xmath193 contribute to @xmath187 . in such excited states one of the @xmath7 walkers travels at least a distance @xmath194 into the @xmath6 domain ( see fig . [ a10bk1 ] ) . for this case we can give the upper bound @xmath195 . from eq . ( [ psums1 ] ) , @xmath196 . hence , @xmath197 \(iii ) _ ground states where @xmath185 _ : again , since @xmath168 site @xmath57 has to be inside the @xmath7 domain . the values of @xmath24 satisfying this condition are in the range @xmath198 . in this case only excited states in which one of the @xmath6 walkers travels at least a distance @xmath199 into the @xmath7 domain will contribute to @xmath200 . hence we can use the upper bound @xmath201 , in this case . therefore , @xmath202 \(iv ) _ ground states where @xmath186 _ : there are three possibilities here . @xmath203 site @xmath57 is inside the @xmath8 domain and site @xmath161 is inside the @xmath6 domain ( @xmath204 ) , @xmath205 both the sites @xmath57 and @xmath161 are inside the @xmath8 domain ( @xmath206 ) , @xmath207 site @xmath57 inside the @xmath7 domain and site @xmath161 inside the @xmath8 domain ( @xmath208 ) . since all these are consistent with @xmath209 , all these cases can occur . it can be shown that the minimal energy needed to create an excited state where @xmath210 and @xmath193 is @xmath211 for the case @xmath203 , @xmath212 for the case @xmath205 and @xmath213 for the case @xmath207 . the resulting expression for the bound is @xmath214 the summations on the rhs of eqs . ( [ ubound2]-[ubound4 ] ) can be carried out explicitly . to leading order , the summations gives @xmath215 for each of eqs . ( [ ubound2 ] ) and ( [ ubound3 ] ) . the summation on the rhs of eq . ( [ ubound4 ] ) vanishes exponentially in the thermodynamic limit . using eqs . ( [ a1br]-[ubound3 ] ) , we get the following expression for the upper bound on @xmath216 @xmath217 \;.\ ] ] therefore @xmath218 . similarly one can show that @xmath219 . thus for all @xmath35 , @xmath220 , proving the existence of a complete phase separation . we have demonstrated that in the thermodynamic limit the system is phase separated when @xmath34 . the general arguments given in section ii indicate that when the global densities of the three species are non - vanishing and @xmath221 , the system phase separates , even when the three densities are not equal . the argument suggests that the typical time , @xmath222 , in which the system leaves a specific phase separated configuration increases exponentially with the system size . thus , a phase separated state is stable in the thermodynamic limit . in the following we use monte carlo simulations to support these arguments . the time , @xmath222 , can be measured using the auto - correlation function defined as , @xmath223 where @xmath224 and @xmath225 are the values of the occupation variables @xmath42 and @xmath43 at time @xmath226 , and @xmath227 denotes an average over histories of evolution . clearly , @xmath228 while @xmath229 , the value of the autocorrelation between two independent configurations . thus , @xmath222 may be defined as the decay time of @xmath230 to @xmath231 when at @xmath232 the system is totally phase separated . we have measured the time scale @xmath222 using monte carlo simulations for different system sizes for @xmath34 and for @xmath233 for several @xmath54 values . an example of such measurements for @xmath234 , @xmath235 and @xmath236 is presented in fig . [ fig_corrt ] . in the plot @xmath222 is plotted versus system size for several values of @xmath54 . it agrees with the exponential growth of @xmath222 with the system size suggested by the simple argument of section ii . the same behavior seems to occur for all @xmath237 and different choices of @xmath238 , @xmath239 and @xmath240 . therefore we conclude that the monte carlo simulations support the claim that for any @xmath241 the system will phase separate into three domains in the thermodynamic limit , even when the number of particles of each species is not equal . in the thermodynamic limit the translational symmetry is spontaneously broken in this state . due to the slow dynamics , which reflects escape from metastable states , monte carlo simulations could be performed only for relatively small system size ( @xmath242 ) . in order to study the coarsening process for larger systems we employ , in the following , a toy model which mimics the dynamics of the model ( [ eq : dynamics ] ) . the toy model can be conveniently simulated for systems larger by about two orders of magnitude . we now construct a simple toy model which captures the essential physics of the coarsening process in the model at large times and enables us to simulate systems much larger than those accessible by monte carlo simulation . using the toy model we examine another characteristic scale of the system . namely , the average domain size @xmath243 as a function of the time @xmath226 . the results support the simple argument leading to a domain growth law @xmath244 . a mean - field version of the toy model is then solved analytically . we consider a system at time @xmath226 such that the average domain size , @xmath243 , is much larger than the domain wall width . at these time scales , the domain walls can be taken as sharp and we may consider only events which modify the size of domains . this means that the dynamics of the system can be approximated by considering only the movement of particles between neighboring domains of the same species . using this we represent a configuration by a sequence of domains of the form @xmath245 , where the @xmath1th domain of , say @xmath6 , particles is represented by @xmath246 , as shown in fig . [ fig_toy1 ] . the exchange of particles between domains , say @xmath247 and @xmath248 , takes place at a rate dictated by the size of the domains @xmath249 and @xmath250 which separate them . since intermediate configurations of the form @xmath251 rearrange on short time scales compared with the evolution between metastable states , only metastable configurations are considered in the toy model . events in which a domain splits into two are ignored .    using these ideas we define the dynamics of the toy model as follows : at each time step two neighboring domains of the same species of particle are chosen randomly , say @xmath246 and @xmath252 . let @xmath253 , @xmath254 and @xmath255 denote the lengths of the domains @xmath256 and @xmath257 respectively . the length of the domain chosen is then modified by carrying out one of the following processes : @xmath258{ll }      \left . \begin{array}[]{ll }      1 ) & a_i \rightarrow a_i - 1    \\        & a_{i+1}\rightarrow a_{i+1 } + 1   \\ \end{array } \right\ } & { \rm with~rate } \quad q^{b_i } \\ & \\      \left . \begin{array}[]{ll }      2 ) & a_i \rightarrow a_i + 1     \\        & a_{i+1}\rightarrow a_{i+1 } - 1   \\ \end{array } \right\ } & { \rm with~rate } \quad q^{c_i } \\ \end{array}\ ] ] where , as before , @xmath14 is considered . if @xmath253 becomes zero , then delete the domain @xmath246 from the list of domains , and merge @xmath259 and @xmath257 with @xmath260 and @xmath261 , respectively . then for @xmath262 , shift the indices of the domains from @xmath77 to @xmath263 , so that @xmath264 becomes @xmath265 . the rules for updating @xmath266 and @xmath267 domains can be obtained from ( [ toy_rule ] ) using cyclic permutations and a slight change of indices . note that the toy model is only relevant to the description of the coarsening dynamics . this is because here , once the system is left with three domains , it remains in that state .    to simulate the toy model efficiently , an algorithm suitable for rare event dynamics must be used due to the small rate of events @xcite . we use an algorithm which is performed by repeating the following steps :    1 .   list all possible events @xmath268 and assign to them rates @xmath269 according to the rules of the model . 2 .   choose an event @xmath128 with probability @xmath270 where @xmath271 . advance time by @xmath272 , where @xmath273 . the algorithm would be equivalent to a usual monte carlo simulation , where @xmath57 time step is equivalent to one monte carlo sweep , if in step @xmath274 , @xmath275 would be drawn from a poisson distribution @xmath276 $ ] . however , here we make an approximation by using the deterministic choice @xmath277 . we have simulated the dynamics for lattices of size up to @xmath278 . for simplicity we consider the case @xmath34 . an example of a typical behavior of the average domain size as a function of @xmath226 is shown in fig . [ fig_domsize ] . one can see that after an initial transient growth time the data fits very well with a @xmath279 behavior . simulations for different @xmath54 values indicate that , @xmath280 with @xmath281 . the toy model enables one to verify the scaling behavior ( [ loga ] ) and estimate the constant @xmath282 . this would be very difficult to do by simulation of the full model ( [ eq : dynamics ] ) .      here we present the solution of a mean - field version of the toy model based on ideas presented by rutenberg and bray @xcite and derrida _ et al._@xcite in the study of the ordering dynamics in a one - dimensional scalar model . to construct the mean - field model we notice that since all steps in the toy model which involve exchange of particles between domains occur at a rate exponentially small in the size of domains , one can consider a model where dynamics occurs only in the vicinity of the smallest domains . the mean field approximation assumes that different domains are uncorrelated and does not distinguish between domains of different species . the second assumption relaxes the conservation of particles of each species . thus , in contrast to the systems studied by @xcite and @xcite we do not expect the mean field to become exact in the scaling limit . we define the mean - field model as follows :    1 . pick one of the smallest domains @xmath283 . 2 .   pick @xmath284 domains @xmath285 and @xmath286 randomly and treat them as the neighbors of @xmath283 . 3 .   pick @xmath274 more domains randomly say , @xmath287 and @xmath288 . 4 .   eliminate @xmath289 and @xmath290 from the system and change the length of the domains @xmath287 and @xmath288 by @xmath291    steps @xmath57 to @xmath180 are performed simultaneously for all of the smallest domains in the system . here @xmath292 is the length of domain @xmath1 , and @xmath293 is the length of the smallest domain . steps @xmath57 to @xmath274 choose the smallest domain @xmath283 and its nearest neighbors ( see fig . [ fig_toy2 ] ) . step @xmath180 uses the fact that the dynamics occur in the model only in the vicinity of the smallest domain and eliminates the three domains @xmath283 , @xmath294 and @xmath295 joining them appropriately with the other domains @xmath296 , @xmath297 and @xmath298 . note that this mean - field dynamics does not take into account the time taken for these events to happen . this will be done later when we derive the growth law of domains .    to solve the mean - field model we follow the method used in @xcite . let @xmath299 be the number of domains of size @xmath24 , irrespective of the type of particle it consists of . let @xmath300 be the length of the smallest domain and @xmath301 be the total number of domains at time @xmath226 . we will assume that @xmath299 has the following scaling form in the large - time limit . @xmath302 a solution of the model given in appendix [ meanfield ] yields @xmath303 where @xmath304 is given by @xmath305    the growth law for @xmath306 can be derived following @xcite : after the elimination of the smallest domain , @xmath307 increases by @xmath57 . this happens at a rate @xmath308 , namely , the inverse time required by a typical domain to cross a distance @xmath309 ( thus , causing the annihilation of @xmath283 ) . using @xmath310 we write @xmath311 from this one can obtain the scaling form of the average domain size , @xmath312    note that in this eq . @xmath304 does not depend , according to the mean - field solution , on either @xmath226 or @xmath54 . one can see from ( [ sform ] ) that for large @xmath24 , @xmath313 which was confirmed by the simulations of the toy model where @xmath314 ( see ( [ loga ] ) ) . a numerical evaluation of ( [ avexpp ] ) yields @xmath315 as compared with @xmath316 obtained from the toy model simulations . in section iii the partition function , @xmath52 , and correlation function , @xmath317 for finite @xmath161 , have been calculated in the thermodynamic limit . it is also of interest to obtain results for finite systems for the study of finite - size effects and the approach to the thermodynamic limit . recently a matrix ansatz method has been introduced to study one - dimensional non - equilibrium systems @xcite . it has been shown that in certain three - species models the steady - state weight and correlation functions can be represented as a product of matrices @xcite . in the ansatz a specific matrix is associated with each type of particle . then the unnormalized probability of a certain configuration is obtained from a matrix product . the matrices corresponding to the different species of particles satisfy an algebra derived from the dynamics of the model . a scalar , _ i.e. _ the weight of a configuration or some correlation function , is usually obtained by performing a trace over the product of matrices , or by multiplying both sides of the product of matrices by vectors . generalizing this method to replace matrices by tensors @xcite , we have been able to obtain recursion relations for the partition function and correlation function for finite systems for the special case @xmath34 . the recursion relations are then used to obtain the partition function and correlation function @xmath318 for any @xmath161 in small systems . the results are used to study the scaling of the correlation function near the critical point @xmath12 ( infinite temperature ) , where the typical domain wall width diverges . it is convenient to consider the unnormalized weights , @xmath319 , defined through @xmath320 where @xmath321 is the probability of being in configuration @xmath322 . the partition sum @xmath52 is given by @xmath323 where the sum is over all configurations with @xmath34 . we generalize the matrix ansatz and construct the steady - state weight , @xmath324 , from a product of tensors each corresponding to a particle located in a specific place on the lattice . the contraction of the tensors yields a tensor which is then contracted with ` left ' and ` right ' tensors to generate a scalar . the three tensors which represent the different type of particle are defined as rank @xmath325 tensors through the following tensor products of square matrices @xmath326 here @xmath327 is a unit matrix . the matrices @xmath328 and @xmath329 will be chosen in what follows to satisfy a commutation relation which will be dictated by the detailed balance condition .    to define @xmath324 we introduce the following notation : the contraction of two rank @xmath325 tensors @xmath330 and @xmath331 , where @xmath332 and @xmath333 are square matrices , according to the rule @xmath334 is denoted by @xmath335 . the contraction of a rank 6 tensor @xmath336 with a ` left ' rank @xmath274 tensor @xmath337 , where @xmath338 are transposed vectors , and a ` right ' rank @xmath274 tensor , @xmath339 , where @xmath340 are vectors , defined through @xmath341 is denoted by @xmath342 .    using these definitions we write the steady - state weight of the system as @xmath343 \vert { \cal v } \rangle \ ; , \ ] ] where @xmath344 , @xmath45 and @xmath43 are the occupation variables defined in ( [ ocup ] ) . the expression states that a tensor @xmath345 is present at place @xmath1 in the tensor product if site @xmath1 is occupied by an @xmath6 particle , a tensor @xmath346 is present if site @xmath1 is occupied by a @xmath7 particle , and a tensor @xmath347 is present if site @xmath1 is occupied by a @xmath8 particle . the action of the tensor product on @xmath348 and @xmath349 ( to be specified later ) produces the scalar @xmath324 . it is straightforward to show using detailed balance that a necessary condition for eq . ( [ eq : ansatz ] ) to be the steady - state weight is that the following commutation relations are satisfied between @xmath350 , @xmath351 , and @xmath352 : @xmath353 using ( [ eq : tensors ] ) one can verify that these commutation relations are satisfied provided that @xmath354 . this deformed commutator is of relevance in other stochastic systems @xcite . a representation of the matrices @xmath328 and @xmath329 which satisfies this commutation relation can be obtained as follows : let @xmath355 denote a basis set ( @xmath356 ) forming a vector space . in this basis we choose the matrices so that , @xmath357 while for @xmath358 , @xmath359 . an explicit form for @xmath329 and @xmath328 is given by the following @xmath360 square matrices @xmath361 to obtain @xmath324 , @xmath362 and @xmath349 have to be specified . this should obviously be done so that @xmath324 is non - zero if the ansatz is to give a non - trivial result . we consider a general tensor product which corresponds to some configuration . the product has @xmath98 tensors of each type @xmath363 and @xmath352 , which results in a tensor product of three matrix products . using ( [ eq : tensors ] ) it can be seen that each matrix product contains @xmath98 matrices of each type @xmath328 , @xmath329 and @xmath327 . since @xmath329 and @xmath327 are diagonal , while @xmath328 acts to its left as a lowering matrix ( see ( [ eq : matrix ] ) ) , choosing @xmath364 , and @xmath365 will give a non - zero @xmath366 . this makes clear that the minimal size choice for the vector - spaces is @xmath367 . under this choice it is easy to see that in the ground states one has @xmath368 which corresponds to a ground state energy @xmath84 . however , the choice of @xmath369 and @xmath370 is determined only up to some multiplicative factors . these factors may be used to shift the ground state energy of the system . for example choosing @xmath369 as before with @xmath371 will shift the ground state energy to @xmath372 . in the following the factors are taken to be @xmath57 . finally we would like to remark that usually when using the matrix ansatz for systems with periodic boundary conditions it is often convenient to use a trace of the matrix product @xcite . in this case this is not possible since the trace of our tensor product is always zero .      in terms of the ansatz ( [ eq : ansatz ] ) the partition function @xmath52 is given by @xmath373 note that we are using the canonical ensemble since all tensor products with unequal number of particles do not contribute to @xmath52 this is easily seen since in these cases there are always more than @xmath98 @xmath328 matrices acting on one of the vectors @xmath374 .    to obtain the partition function we derive a recursion relation for @xmath375 one can see that @xmath376 . rewriting @xmath377 as @xmath378 and using relations ( [ eq : tensors ] ) and ( [ eq : matrix ] ) the following recursion relation can be derived @xmath379 the boundary conditions for this recursion relation is given by the no particle partition function @xmath380 if @xmath381 and is zero otherwise .    for small systems ( up to @xmath382 ) for which the recursion relation is tractable analytically on mathematica , we obtained the partition function @xmath383 as a polynomial in @xmath54 . as expected the first @xmath384 terms of the polynomial match the first @xmath384 terms of the expansion of ( [ psum ] ) up to a factor of @xmath385 due to the energy shift in the ground state . for larger @xmath5 we solve the recursion relation numerically . we note that ( [ eq : reccursion1 ] ) could have been derived directly from the definition of the partition function without recourse to the tensor ansatz . however , we believe the utility of the ansatz lies in the ease with which correlation functions can be manipulated and relations such as that of the next subsection derived . the correlation function @xmath317 is given in terms of the ansatz by @xmath386 using relations ( [ eq : tensors]),([eq : commutation ] ) and ( [ eq : matrix ] ) we obtain @xmath387 where we define the object @xmath388 through , @xmath389 a recursion relation for @xmath388 can be obtained , similarly to the recursion relation for the partition function ( [ eq : reccursion1 ] ) . using ( [ eq : tensors ] ) and ( [ eq : matrix ] ) gives @xmath390 the boundary conditions for the recursion relation are obtained by noting that @xmath391 , _ i.e. _ , @xmath392 . using the same methods one can obtain recursion relations for all other correlation functions . the recursion relations are solved numerically for finite systems . this is done by first solving numerically for @xmath393 , and then using the result as boundary conditions for the recursion relation ( [ eq : u recurr ] ) . owing to ( [ a1ar ] ) we are ultimately interested in @xmath394 . using these recursion relations we have calculated the correlation function @xmath395 , which is a measure of the phase separation in the system , for various system sizes . when @xmath396 is close to zero the system is phase separated . in the disordered case , @xmath12 , the value of the correlation function is @xmath397 , approaching @xmath398 in the thermodynamic limit . the results are shown in the inset of fig . [ order ] one can see that the system is phase separated for small values of @xmath54 while for @xmath54 close to @xmath57 the system is disordered . the range of @xmath54 values for which the system is phase separated increases as the system size increases . the natural scaling variable near the critical point @xmath12 is @xmath399 , the ratio between the domain size @xmath98 and the domain wall width @xmath400 . in fig.([order ] ) the correlation function is plotted as a function of the scaling variable . one can see that the data collapse improves as the system size increases . this scaling variable was also suggested by the form of the partition function ( [ psum ] ) in section iii c. in the following we discuss possible generalization of the model to @xmath401 species . to demonstrate how this might be done we first discuss the case @xmath402 . we then comment briefly on @xmath403 . we now define a @xmath180-species model and argue that it phase separates . consider a ring where each site is occupied by either an @xmath6 , @xmath7 , @xmath8 or @xmath404 particle . the model evolves according to the following procedure : at each time step two nearest neighbors are chosen randomly and exchanged according to the rates ,    @xmath405    as before the model conserves the number of particles of each species . note that several other generalizations of the model to four species are possible . however , for simplicity , we discuss only the model defined by ( [ eq : four ] ) with @xmath35 . we now argue that the system phase separates into a configuration of the form @xmath406 ( where each letter now indicates a domain ) as long as the densities of particles of each species are non - zero . note that in the model @xmath17 , @xmath18 , @xmath407 , @xmath408 , @xmath15 and @xmath409 boundaries are stable while reverse boundaries @xmath410 , @xmath16 , @xmath411 , @xmath412 , @xmath19 and @xmath413 are unstable . as in the case of the three species model the system , starting from a random initial condition , evolves on a short time scale ( _ i.e. _ , which is not determined by the size of the system ) into a metastable configuration where only stable domains are present . this configuration then slowly coarsens by slow diffusion of particles through neighboring domains . the system will finally reach the most stable state where the number of domains is minimal . one can easily check that this configuration is given by @xmath406 . note that the system may exhibit other metastable states . for example , a state composed of @xmath414 is also stable under the choice of rates ( [ eq : four ] ) . however , since this state is composed of more domains than the @xmath180-domain state , some of the domains are necessarily smaller . according to the argument presented in section ii the relaxation time of this sequence ( proportional to @xmath415 where @xmath128 is the typical domain size ) is much shorter than the relaxation time of the @xmath180-domain state . therefore the @xmath180-domain state is more stable so that the system will finally evolve into it .    in considering @xmath403 models one finds that for some choices of transition rates several states with the minimum number of domains may become metastable . for example , for @xmath416 it is possible to choose transition rates for which both @xmath417 and @xmath418 are locally stable . the relative stability ( and thus the resulting phase separated state ) may be found by determining the relaxation time of these states using simple considerations such as those presented in section ii .    as is the case of @xmath3 detailed balance is found to be satisfied for certain densities and transition rates for @xmath419 . the condition for this is derived in appendix [ dbappen ] . in this case the relative stability of metastable states could be determined by comparing free energies . in this paper a model of three species of particles diffusing on a ring previously introduced in @xcite has been studied . the model is governed by local dynamics in which all moves compatible with the conservation of the three densities are allowed . we argue that phase separation should occur as long as all densities are non - zero . in the special case of equal densities we find that the steady state generated by the local stochastic dynamics is exactly given by a long - range asymmetric hamiltonian . phase separation for this case is explicitly demonstrated . the model provides an explicit example for the mechanism leading to phase separation or breaking of ergodicity in systems with _ local _ stochastic dynamics . although we did not succeed in solving the steady state in the case of non - equal densities of particles , there is a strong evidence that phase separation still occurs . in order to investigate further the case of equal densities we employed a generalized matrix ansatz to calculate correlation function for finite - size systems . the novel structure of the ansatz may give some clue as to handle other @xmath0 models which have so far resisted solution . the dynamics of phase separation reduces to a coarsening problem where the typical domain size grows logarithmically in time . this results from the elimination of the domains at a rate exponentially small in their size . the slow dynamics poses a problem of how to access the scaling regime numerically . with direct numerical simulations only small systems can be studied ( see fig . however , by employing a toy model in which domains rather than individual sites are updated one can simulate much larger systems and probe the scaling regime ( see fig .  5 ) . such ideas of updating domains have been used before in the study of coarsening @xcite . with the aid of the toy model it should be possible to study other aspects of the scaling regime associated with the slow dynamics and escape from metastable states . generalizations are possible to models with @xmath419 species . we have discussed some possibilities and have shown that phase separation may take place , although the structure of the set of metastable states is more complicated . as was the case for @xmath3 , conditions for detailed balance with respect to a long - range hamiltonian may be determined . the problem of phase separation and coarsening is of interest also in the broader context of phase transitions in one dimensional systems . here the existence of conserved quantities results in certain local transition rates being zero . it would be interesting to generalize this study to models in which no conserved quantity exists and all local rates are non - vanishing . also , another open problem is to calculate the steady state of the present model in the case of non - equal densities . * acknowledgments * : mre is a royal society university research fellow and thanks the weizmann institute and the einstein center for warm hospitality during several visits . the support of the minerva foundation , munich germany , the israeli science foundation , and the israel ministry of science is gratefully acknowledged . computations were performed on the sp2 at the inter - university high performance computer center , tel aviv . we thank s. wiseman for helpful advice on programming and m. j. e. richardson for careful reading of the manuscript . to solve the mean - field model we follow the method used in @xcite . let @xmath299 be the number of domains of size @xmath24 , irrespective of the type of particle it consists of . let @xmath300 be the length of the smallest domain and @xmath301 be the total number of domains at time @xmath226 . we will assume that @xmath299 has the following scaling form in the large - time limit . @xmath420 after the elimination of the smallest domain as given by ( [ toy2 ] ) , @xmath421 and @xmath4 change according to @xmath422 using the scaling form ( [ nlform ] ) we have , @xmath423\ ] ] where @xmath424 and we have expanded in @xmath425 . now substituting these in ( [ modif ] ) , it is straightforward to show that @xmath426 using the laplace transform , @xmath427~f(x)\ ] ] one can show that @xmath428 satisfies the differential equation , @xmath429[2 \phi(p ) + e^{-p}]\ ] ] since @xmath430 , by expanding ( [ phip ] ) to order @xmath431 one obtains @xmath432 . the solution of ( [ phip ] ) with boundary conditions @xmath433 and @xmath434 for @xmath435 is @xmath436 where @xmath437/t = -\log(x ) - \gamma - \sum_{n=1}^{\infty } ( -x)^n/(n!\ n)$ ] and @xmath438 is the euler constant . inverse laplace transform of ( [ solu ] ) gives the domain size distribution . from ( [ solu ] ) and using the expansion @xmath439 , where the average is with respect to @xmath440 , we get @xmath441 we now define the most general @xmath4 species model , where @xmath442 . let @xmath443 denote a variable at site @xmath1 of a ring of size @xmath5 , which takes values @xmath444 . @xmath445 means that site @xmath1 is occupied by a particle of type @xmath128 . the system evolves by a random sequential , nearest - neighbor exchange dynamics , with the following rates :          consider a particle exchange between sites @xmath62 and @xmath63 , where @xmath451 and @xmath452 ( _ i.e. _ in the bulk , note that site @xmath57 is chosen arbitrarily ) . expanding the product in ( [ eq : sweight ] ) , it is easy to verify that @xmath453 since this hold for any @xmath454 , and is irrespective of the number of particles of each species , the steady - state weight ( [ eq : sweight ] ) satisfies detailed balance for all nearest - neighbor exchanges in the bulk . if the weights ( [ eq : sweight ] ) are translationally invariant then detailed balance will also hold for exchanges between sites @xmath57 and @xmath5 .    thus , to complete the proof of detailed balance it is sufficient to demand that ( [ eq : sweight ] ) is translationally invariant . to do this we relabel sites @xmath455 . the weight then becomes @xmath456 where @xmath457 is identical to @xmath458 . rewriting this equation by relabeling the indices we obtain , @xmath459 \ \prod_{k=1}^{n-1}\frac{q(x_k , x_n)}{q(x_n , x_k ) } \ ; .\ ] ] comparing ( [ eq : sweight3 ] ) with ( [ eq : sweight ] ) and noting for example that , @xmath460^{n_l } \ ; , \ ] ] one can see that ( [ eq : sweight ] ) is translational invariant if @xmath461^{n_{l } } = 1 \ ; , \ ] ] for every @xmath462 . thus , detailed balance holds if ( [ condition ] ) is satisfied . we note that for given densities @xmath463 the manifold of solutions for the rates is of @xmath464 dimensions .
Abstract:

Answer: a driven system of three species of particle diffusing on a ring is studied in detail . the dynamics is _ local _ and conserves the three densities . a simple argument suggesting that the model should phase separate and break the translational symmetry is given . we show that for the special case where the three densities are equal the model obeys detailed balance and the steady - state distribution is governed by a hamiltonian with asymmetric long - range interactions . this provides an explicit demonstration of a simple mechanism for breaking of ergodicity in one dimension . the steady state of finite - size systems is studied using a generalized matrix product ansatz . the coarsening process leading to phase separation is studied numerically and in a mean - field model . the system exhibits slow dynamics due to trapping in metastable states whose number is exponentially large in the system size . the typical domain size is shown to grow logarithmically in time . generalizations to a larger number of species are discussed .