## Brazilian Journal of Physics

##
*Print version* ISSN 0103-9733*On-line version* ISSN 1678-4448

### Braz. J. Phys. vol.30 no.4 São Paulo Dec. 2000

#### http://dx.doi.org/10.1590/S0103-97332000000400004

**Exact solution of asymmetric diffusion with N classes of particles of arbitrary size and hierarchical order **

F. C. Alcaraz

*Departamento de Física, Universidade Federal de São Carlos, 13565-905, São Carlos, SP Brazil*

R. Z. Bariev

*Departamento de Física, Universidade Federal de São Carlos, 13565-905, São Carlos, SP Brazil*

*The Kazan Physico-Technical Institute of the Russian Academy of Sciences, Kazan 420029, Russia*

**Received on 5 August, 2000**

The exact solution of the asymmetric exclusion problem withNdistinct classes of particles (c= 1, 2, ...,N), with hierarchical order is presented. In this model the particles (size 1) are located at lattice points, and diffuse with equal asymmetric rates, but particles in a classcdo not distinguish those in the classesc'>cfrom holes (empty sites). We generalize and solve exactly this model by considering the molecules in each distinct classc= 1, 2, ...,Nwith sizess(_{c}s= 0, 1, 2, ...), in units of the lattice spacing. The solution is derived via a Bethe ansatz of nested type._{c}

** I Introduction **

The similarity between the master equation describing time fluctuations in nonequilibrium problems and the Schrödinger equation describing the quantum fluctuations of quantum spin chains turns out to be fruitful for both areas of research [1]-[15]. Since many quantum chains are known to be exactly integrable through the Bethe ansatz, this provides exact information on the related stochastic model. At the same time classical physical intuition and probabilistic methods successfully applied to nonequilibrium systems give new insights into the physical and algebraic properties of quantum chains.

An example of this fruitful interchange is the problem of asymmetric diffusion of hard-core particles on the one dimensional lattice ( see [16, 17, 18] for reviews). This model is related to the exactly integrable anisotropic Heisenberg chain in its ferromagnetic regime [19] (XXZ model). However if we demand this quantum chain to be invariant under a quantum group symmetry *U*_{q}(*SU*(2)), we have to introduce, for the equilibrium statistical system, unusual surface terms, which on the other hand have a nice and simple interpretation for the related stochastic system [3, 4].

In the area of exactly integrable models it is well known that one of the possible extensions of the spin- XXZ chain to higher spins is the anisotropic spin-S Sutherland model (grading _{1} = _{2} = ... = _{2s+1} = 1) [20]. On the other hand in the area of diffusion limited reactions a simple extension of the asymmetric diffusion problem is the problem of diffusion with particles belonging to *N* distinct classes (*c* = 1, 2, ..., *N*) with hierarchical order [22]-[24] . In this problem a mixture of hard-core particles diffuses on the lattice. Particles belonging to a class *c* (*c* = 1, ...,* N*) ignore the presence of those in classes *c*¢ > *c*, i.e., they see them in the same way as they see the holes (empty sites). In [3] it was shown that for open boundary conditions the anisotropic spin-1 Sutherland model and this last stochastic model, in the case *N* = 2, are exactly related. The Hamiltonian governing the quantum or time fluctuations of both models is given in terms of generators of a Hecke algebra, invariant under the quantum group *U*_{q}*SU*(3). In fact this relation can be extended to arbitrary values of *N*, and the quantum chain associated to the stochastic model is invariant under the quantum *U*_{q}(*SU*(*N* + 1)) group. In this paper we derive through the Bethe ansatz the exact solution of the associated quantum chain, on a closed lattice. Recently [15] (see also [14]) we have shown that without losing its exact integrability, we can consider the problem of asymmetric diffusion with an arbitrary mixture of molecules with different sizes (even zero), as long they do not interchange positions, that is, there is no reactions. In this paper we extend the asymmetric diffusion problem with *N* types of particles with hierarchical order, to the case where the particles in each class have an arbitrary size, in units of the lattice spacing. Unlike the case of asymmetric diffusion problem, we have in this case a nested Bethe ansatz [25]. A pedagogical presentation for the simplest case *N* = 2 was presented in [26].

The paper is organized as follows. In the next section we introduce the generalized asymmetric model with *N* types of particles with hierarchical order and derive the associated quantum chain. In section 3 the Bethe ansatz solution of the model is presented. Finally in section 4 we present our conclusions, with some possible generalizations of the stochastic problem considered in this paper, and some perspectives on future work.

** II The generalized asymmetric diffusion model with N classes of particles with hierarchical order **

A simple extension of the asymmetric exclusion model, in which hard-core particles diffuse on the lattice, is the problem where a mixture of particles belonging to different classes (*c* = 1, 2, ..., *N*) diffuses on the lattice.. This problem in the case where we have only *N* = 2 classes was used to describe shocks [22]-[24] in nonequilibrium and also has a stationary probability distribution that can be expressed via the matrix-product ansatz [27]. In [28] it was also shown that the stationary state of the case *N* = 3 can also be expressed by the matrix-product ansatz. In this model we have *n*_{1}, *n*_{2}, ... *n*_{N} molecules belonging to the classes *c* = 1, 2, ..., *N*, respectively. All classes of molecules diffuse asymmetrically, but with the same asymmetrical rates, whenever they encounter empty sites (holes) at nearest-neighbor sites. However, when molecules of different classes, *c* and *c*¢ (*c* < *c*¢), are at their minimum separation, the molecules of class *c* exchange position with the same rate as they diffuse, and consequently the molecules in the class *c* see no difference between molecules belonging to the classes *c*¢ > *c* and holes.

We now introduce a generalization of the above model, where instead of having unit size, the molecules in each distinct class *c* = 1, 2, ..., *N* have in general distinct sizes *s*_{1}, *s*_{2}, ..., *s*_{N} (*s*_{1}, ..., *s*_{N} = 1, 2, ¼), respectively, in units of lattice spacing. In Fig. 1 we show some examples of molecules of different sizes. We may think of a molecule of size *s* as formed by *s* monomers (size 1), and for simplicity, we define the position of the molecule as the center of its leftmost monomer. The molecules have a hard-core repulsion: the minimum distance *d*_{ab}, in units of the lattice spacing, between molecules a and b, with a in the left, is given by *d*_{ab} = *s*_{a}. In order to describe the occupancy of a given configuration of molecules we define at each lattice site *i* (*i* = 1, 2, ..., *L*) a variable b_{i} (*i* = 1, 2, ¼, *L*), taking the values b_{i} = 0, 1, ..., *N*. The values b = 1, 2, ¼, *N* represent sites occupied by molecules of class *c* = 1, 2, ¼,* N*, respectively. On the other hand the value b = 0 represents an empty site or an excluded one, due to the finite size of the molecules. As an example, in a chain of *L* = 8 sites, the configuration in which a particle of class 1, with size *s*_{1} = 2 is at site 1, and another particle, of class 2, with size *s*_{2} = 3 is at site 3, is represented by {b} = {1, 0, 0, 2, 0, 0, 0, 0}. Thus the allowed configurations are given by the set {b_{i}} (*i* = 1, ¼, *L*), where for each pair (b_{i}, b_{j}) ¹ 0 with *j* > *i* we should have *j* -*i* ³ *s*_{bi}.

Figura.1 Example of configurations of molecules with distinct sizes s in a lattice of size L = 6. The coordinates of the molecules are denoted by the black squares. |

The time evolution of the probability distribution *P*({b}, *t*), of a given configuration {b} is given by the master equation

where G({b}® {b¢}) is the transition rate for configuration {b} to change to {b¢}. In the present model we only allow, whenever the constraint of excluded volume is satisfied, the particles to diffuse to nearest-neighbor sites, or to exchange positions. The possible motions are diffusion to the right

diffusion to the left

and interchange of particles

As we see from (4), particles belonging to a given class *c* interchange positions with those of class *c*¢ > *c* with the same rate as they interchange positions with the empty sites (diffusion). We should remark however that unless the particles in class *c*¢ have unit size (*s*_{c} = 1), the net effect of these particles in those of class *c* is distinct from the effect produced by the holes, since as the result of the exchange the particles in class *c* will move by *s*_{c¢} lattice size units, accelerating its diffusion.

The master equation (1) can be written as a Schrödinger equation in Euclidean time (see Ref. [3] for general application for two body processes)

if we interpret | *P* > º *P*({b}, *t*) as the associated wave function. If we represent b_{i} as |b > _{i} the vector |b > _{1 }Ä |b > _{2} Ä ¼ Ä |b > _{L} will give us the associated Hilbert space. The process (2)-(4) gives us the Hamiltonian (see Ref. [3] for general applications)

with

and periodic boundary conditions. The matrices *E*^{a,b} are (*N* + 1) × (*N *+ 1) matrices with a single nonzero element (*E*^{a, b})_{i,j} = d_{a, i}d_{b, j} (a, b, *i*,* j* = 0, ..., *N*). The projector in (6), projects out from the associated Hilbert space the vectors |{b} > which represent forbidden positions of the molecules due to their finite size, which mathematically means that for all *i*,*j* with b_{i}, b_{j} ¹ 0, |*i* - *j*| ³ *s*_{bi} (*j* > *i*). The constant *D* in (6) fixes the time scale; for simplicity we chose *D* = 1. A particular simplification of (6) occurs when the molecules in all classes have the same size *s*_{1} = *s*_{2} = ... = *s*_{N} = *s*. In this case the Hamiltonian can be expressed as an anisotropic nearest-neighbor interaction spin-*N*/2 *SU*(*N* + 1) chain. Moreover in the case where their sizes are unity (*s* = 1) the model can be related to the anisotropic version [21] of the *SU*(*N* + 1) Sutherland model [20] with twisted boundary conditions.

**III The Bethe ansatz equations **

We present in this section the exact solution of the general quantum chain (6). A pedagogical presentation for the particular case where *N* = 2 was presented in [26].

Due to the conservation of particles in the diffusion and interchange processes the total number of particles *n*_{1}, *n*_{2}, ..., *n*_{N} in each class are good quantum numbers and consequently we can split the associated Hilbert space into block disjoint sectors labeled by the numbers *n*_{1}, *n*_{2}... *n*_{N} (*n*_{i} = 0, 1, ¼; *i* = 1, ..., *N*). We therefore consider the eigenvalue equation

where

and *n* = *n*_{i} is the total number of particles. In (10) |*x*_{1},* Q*_{1}; ¼;* x*_{n}, *Q*_{n} > means the configuration where a particle of class *Q*_{i} (*Q*_{i} = 1, 2, ..., *N*) is at position *x*_{i} (*x*_{i} = 1, ¼, *L*). The summation {*Q*} = {*Q*_{1}, ¼, *Q*_{n}} extends over all permutations of the *n* integer numbers {1, 2, ..., *N*} in which *n*_{i} terms have the value *i* (*i* = 1, 2, ..., *N*) , while the summation {*x*} = {*x*_{1}, ¼, *x*_{n}} runs, for each permutation {*Q*}, in the set of the *n* nondecreasing integers satisfying

Before getting the results for general values of *n* let us consider initially the cases where we have 1 or 2 particles.

**n = 1.** For one particle on the chain, in any class *c* = 1, 2, ..., *N*, as a consequence of the translational invariance of (6) it is simple to verify directly that the eigenfunctions are the momentum-*k* eigenfunctions

with

and energy given by

**n =2.** For two particles of classes *Q*_{1} and *Q*_{2} (*Q*_{1},* Q*_{2} = 1, 2, ...,* N*) on the lattice, the eigenvalue equation (9) gives us two distinct relations depending on the relative location of the particles. The first relation applies to the case in which a particle of class *Q*_{1} (size *s*_{Q1}) is at position *x*_{1} and a particle *Q*_{2} (size *s*_{Q2}) is at position *x*_{2}, where *x*_{2} > *x*_{1} + *s*_{Q1}. We obtain in this case the relation

where we have used the relation _{+} + _{-} = 1. This last equation can be solved promptly by the ansatz

with energy

where *k*_{1},* k*_{2}, and are free parameters to be fixed. In (16) the summation is over the permutations *P* = *P*_{1},* P*_{2} of (1, 2). The second relation applies when *x*_{2} = *x*_{1} + *s*_{Q1}. In this case instead of (15) we have

If we now substitute the ansatz (16) with the energy (17), the constants and , initially arbitrary, should now satisfy

where

At this point it is convenient to consider separately the case where *Q*_{1} = *Q*_{2} from those where *Q*_{1} ¹ *Q*_{2}. If *Q*_{1} = *Q*_{2} = *Q* (*Q* = 1, ..., *N*) eq. (19) gives

and the cases *Q*_{1} ¹ *Q*_{2} give us the equations

Performing the above summation we obtain, after lengthy but straightforward algebra, the following relation among the amplitudes

where

Equations (21) and (22) can be written in a compact form

with

where we have introduced the *S* matrix. From (21) and (23) this S matrix has only *N*(2*N* - 1) non zero elements, namely

Equations (23) do not fix the "wave numbers" *k*_{1} and *k*_{2}. In general, these numbers are complex, and are fixed due to the cyclic boundary condition

which from (16) give the relations

This last equation, when solved by exploiting (23)-(25), gives us the possible values of *k*_{1} and *k*_{2}, and from (17) the eigenenergies in the sector with 2 particles. Instead of solving these equations for the particular case *n* = 2 let us now consider the case of general *n*.

**General n.** The above calculation can be generalized for arbitrary occupation {*n*_{1},* n*_{2}, ..., *n*_{N}} of particles in classes 1, 2, ..., *N*, respectively. The ansatz for the wave function (10) becomes

where the sum extends over all permutations *P* of the integers 1, 2, ¼,* n*, and *n* = *n*_{i} is the total number of particles.

Application of the translation operator to the above wave functions implies that (10) are also eigenfunctions of the momentum operator with eigenvalues

For the components |*x*_{1},* Q*_{1}; ¼; *x*_{n}, *Q*_{n} > where *x*_{i+1} -*x*_{i} > *s*_{Qi} for *i* = 1, 2, ¼,* n*, it is simple to see that the eigenvalue equation (9) is satisfied by the ansatz (28) with energy

On the other hand if a pair of particles of class *Q*_{i}, *Q*_{i+1} is at positions *x*_{i}, *x*_{i+1}, where *x*_{i+1} = *x*_{i} +* s*_{Qi}, equation (9) with the ansatz (28) and the relation (30) give us the generalization of relation (23), namely

with *S* given by eq. (25). Inserting the ansatz (28) in the boundary condition

we obtain the additional relation

which together with (31) should give us the energies.

Successive applications of (31) give us in general distinct relations between the amplitudes. For example relate to by performing the permutations abg ® bag ® bga ® gba or abg ® agb ® gab ® gba, and consequently the *S*-matrix should satisfy the Yang-Baxter [19, 29] equation

for a, a¢, a¢¢, b, b¢, b¢¢ = 1, 2, ..., *N* and *S* given by (25). Actually the relation (34) is a necessary and sufficient condition [19, 29] to obtain a non-trivial solution for the amplitudes in Eq. (31).

We can verify by a long and straightforward calculation that for arbitrary number of classes *N* and values of the sizes *s*_{1},* s*_{2}, ..., *s*_{N}, the *S* matrix (25), satisfies the Yang-Baxter equation (34), and consequently we may use relations (31) and (33) to obtain the eigenenergies of the Hamiltonian (6). Applying relation (31) *n* times on the right of equation (33) we obtain a relation between the amplitudes with the same ordering in the lower indices:

where we have introduced the harmless extra sum

(see [26] for illustrations of the above equations). In order to fix the values of {*k*_{j}} we should solve (35), i.e., we should find the eigenvalues L(*k*) of the matrix

with periodic boundary condition

The Bethe-ansatz equations which fix the set {*k*_{l}} will be given from (35) by

The matrix *T*(*k*) has dimension *N*^{n} ×* N*^{n} and can be interpreted as the transfer matrix of an inhomogeneous *N*(2*N *- 1)-vertex model in a two dimensional lattice with periodic boundary conditions in the horizontal direction (*n* sites). Due to the special form of the S matrix (25) the eigenvalues of (37) are invariant under a local gauge transformation where for each factor *S*(*k*_{Pl},* k*) in (37):

where (*l* = 1, ..., *L*; a = 1, ..., *N*) are arbitrary functions. If we perform the transformation (40) with the special choice

the equivalent transfer matrix to be diagonalized is given by

where

with the twisted boundary condition

with twisted phase

The matrix in (43) and (44) is obtained from those in (25) by taking the size of all particles equal to unity. In this way the problem is transformed into the evaluation of the eigenvalues of a regular (all particles with size 1) inhomogeneous transfer matrix *T*_{0} with *n*(2*N* - 1) non-zero vertex and twisted boundary condition.

**Diagonalization of T_{0}(k)**

The simplest way to diagonalize *T*_{0} is through the introduction of the monodromy matrix (*k*) [25], which is a transfer matrix of the inhomogeneous vertex model under consideration, where the first and last link in the horizontal direction are fixed to the values m_{1} and m_{N + 1} (m_{1}, m_{N + 1} = 1, 2, ¼, *N*), that is

The monodromy matrix (*k*) has coordinates {*Q*},{*Q*¢} in the vertical space (*N*^{n} dimensions) and coordinates m_{1}, m_{N + 1} in the horizontal space (*N*^{2} dimensions). This matrix satisfies the following important relations

for m_{1}, n_{1}, m_{N + 1}, n_{N + 1} = 1, 2, ¼, *N*.This relation follows directly from successive applications of the Yang-Baxter equations (34) (see [26], for a graphical representation of these equations).

In order to exploit relation (47) let us denote the components of the monodromy matrix in the horizontal space by

where a, b = 1, 2, ¼, *N *- 1. Clearly the transfer matrix *T*_{0}(*k*) of the inhomogeneous lattice with twisted boundary conditions, we want to diagonalize, is given by

As a consequence of (47) the matrices , *B*^{a}, *C*^{a} and *D* in (48) obey some algebraic relations. By setting (n_{1}, m_{1}, n_{N + 1}, m_{N + 1}) = (*N*, a, g, b) in (47) we obtain

with (a, b = 1, ¼, *N *- 1). By setting (n_{1}, m_{1}, n_{N + 1}, m_{N + 1}) = (*N*, *N*, *N*, a) we obtain

where (a = 1, ¼,* N* - 1). The diagonalization of *T*_{0}(*k*) in (49) will be done by exploiting the above relations. This procedure is known in the literature as the algebraic Bethe ansatz [25]. The first step in this method follows from the identification of a reference state |W > , which should be an eigenstate of (*k*) and *D*(*k*), and hence *T*_{0}(*k*), but not of *B*^{a}(*k*). In the present case a suitable reference state is |W > = |{a_{l} = *N*} > _{l = 1, ¼, n}, which corresponds to a state with N-class particles only. It is simple to calculate

where

and | > = |{a_{l ¹ i} = *N*}, a_{i} = a > . The matrices *B*^{a}(*k*) act as creation operators in the reference ("vacuum") state, by creating particles of class a (1, 2, ¼, *N*) in a sea of particles of Nth class |W > . We then expect that the eigenvectors of *T*_{0}(*k*) corresponding to *m*_{1} (1, 2, ¼, *n*) particles, belonging to classes distinct from *N*, can be expressed as

where {*k*_{l}^{(1)}, *l* = 1, ¼, *m*_{1}} and *F*_{b1, ¼, bm1} are variables to be fixed by the eigenvalue equation

Using (50) successively, and (52), (53) we obtain

where the "unwanted terms" are those ones which are not expressed in the "Bethe basis" produced by the *B*^{a}(*k*_{j}^{(1)}) operators. Similarly, using (51) successively and (52)-(53) we obtain

The relations (56) and (57) when used in (54)-(55) give us

where

is a (*N*-1)^{m1}-dimensional transfer matrix of a inhomogeneous vertex model, with inhomogeneties (notice the reverse order of the inohomogeneties, when compared with (43)) and twisted boundary conditions (boundary phases F_{a}, a = 1, ¼, *N* - 1).

In order to proceed we need now to diagonalize the new transfer matrix *T*_{1}(*k*), that is we must solve

and then (58) give us

where, using the fact that (*k*_{l}, *k*) = 1,

In order to prove that L^{(0)} and |,* F* > are the eigenvalues and eigenvectors of *T*_{0}(*k*), we should fix {, ¼, } by requiring that the "unwanted terms" in (61) vanish. Although for *N* = 2 this calculation is not complicated [26] for arbitrary *N* it is not simple. Since the expression (62) for the eigenvalues should be valid for arbitrary values of *k* we can obtain L^{(1)}(*k*_{j}^{(1)}) in an alternative way from the following trick [31]. At *k* = *k*_{j}^{(1)} (*j* = 1, ¼,* m*_{1}) the denominators of the factors in (62) vanish ((*k*_{j}^{(1)},* k*_{j}^{(1)}) = 0 , *l* ¹ *N*), and since we should have a finite result, we have the conditions

Notice that our result in (63) does not depend on the particular ordering of the additional variables *k*_{j}^{(1)} (*j* = 1, ¼, *m*_{1}). This means that if instead of the ordering chosen in (54), we chose the reverse order, namely,

we would obtain the same results (61)-(63), but now *T*_{1} is the transfer matrix, with boundary condition specified by the phase F_{a}, of a problem with (*N*-1) species and inhomogeneities , ¼, (notice we have now the same order of the inhomogeneties as in (43)). This means that the eigenvalue L(*k*) = L^{(0)}(*k*) of the transfer matrix of the problem with *N* classes and inhomogeneities () º (*k*_{1}, *k*_{2}, ¼, *k*_{n}) is related to the eigenvalue L^{(1)}(*k*) of the problem with (*N* - 1) classes and inhomogeneities , , ¼, . Iterating these calculations we obtain the generalization of the relation (62) and the condition (63)

which connects the eigenvalues of the inhomogeneous transfer matrix *T*_{l}(*k*) and *T*_{l+1}(*k*), with inhomogeneities {*k*_{j}^{(l)}} and , related with the problem with (*N *- *l*) and (*N *- *l* - 1) classes of particles, respectively.

However from (39) and (42)-(43), in order to obtain the Bethe-ansatz equations for our original problem we need the eigenvalues of the transfer matrices evaluated at *k*_{j} (*j* = 1, ¼, *n*), i.e., L^{(0)}(*k*_{j}), which are given by

since (*k*_{j},* k*_{Pj}) = 0. The conditions that fix the variables (*k*_{j}^{(1)}, *j* = 1, ¼, *m*_{1}) are given by (63). In the left side of this equation we have L^{(1)}(*k*_{j}^{(1)}), which are the eigenvalues of the transfer matrix *T*_{1} of the model with (*N* - 1) classes of particles and inhomogeneities {*k*_{j}^{(1)}, *j* = 1, ¼, *m*_{1}}, evaluated at the partcular point *k*_{j}^{(1)}. This value can be obtained from (65) which gives a generalization of (67)

The condition (63) is then replaced by

where now we need to find the relations that fix {*k*_{j}^{(2)}}. Iterating this process we find the generalization of (69)

Equations (67) and (70) give us the eigenvalues of the transfer matrix *T*_{0}(*k*) evaluated at the points {*k*_{j}}, i. e. L^{(0)}(*k*_{j}). Inserting the above results in (42) and then in (39) we obtain the Bethe-ansatz equations of our original problem.

The eigenenergies of the Hamiltonian (6) in the sector containing *n*_{i} particles in class *i* (*i* = 1, 2, ..., *N*) ( *n* = * n*_{j}) and total momentum *p* = (*l* = 0, 1, ..., *L* - 1) are given by

where {*k*_{j}^{(0)} = *k*_{j}, *j* = 1, ..., *n*} are obtained from the solutions {*k*_{j}^{(l)}, *l* = 0, ..., *N* - 1; *j* = 1, ..., *m*_{l}} of the Bethe ansatz equations

and

and *m*_{l} = * n*_{j}, *l* = 0, ¼, *N* (*m*_{0} = *N*,* m*_{N} = 0). It is interesting to observe that in the particular case where *n*_{2} = *n*_{3} = ... = *n*_{N} = 0 we obtain the Bethe-ansatz equations, recently derived [15] (see also [14]), for the asymmetric diffusion problem with particles of size *s*_{1}. Also the case *s*_{1} = *s*_{2} = ... = *s*_{N} = 1 gives us the corresponding Bethe ansatz equations for the standard problem of *N* types of particles in hierarchical order. The Bethe-ansatz solution in the particular case of N=2 with a a single particle of class 2 (*n*_{1} = *n* - 1, *n*_{2} = 1) was derived recently [32]. The Bethe-ansatz equations for the fully asymmetric problem are obtained by setting in (72)-(73) _{+} = 1 and _{-} = 0.

** IV Conclusions and generalizations **

We obtained through the Bethe ansatz the exact solution of the problem in which particles belonging to *N* distinct classes with hierarchical order diffuse as well interchange positions with rates depending on their relative hierarchy. We show that the exact solution can also be derived in the general case where the particles have arbitrary sizes.

Some extensions of our results can be made. A first and quite interesting generalization of our model happens when we allow molecules in any class to have size *s* = 0. Molecules of size zero do not occupy space on the lattice, having no hard-core exclusion effect. Consequently we may have, at a given lattice point, an arbitrary number of them. The Bethe-ansatz solution presented in the previous section is extended directly in this case (the equations are the same) and the eigenenergies are given by fixing in (71)-(73) the appropriate sizes of the molecules. It is interesting to remark that particles of a given class *c*¢ (2, 3, ..., *N*), with size *s*_{c¢} = 0, contrary to the case *s*_{c¢} > 1, where they "accelerate" the diffusion of particles in classes *c* < *c*¢, they now "retard" the diffusive motion of these particles. The quantum Hamiltonian in the cases where the particles have size zero is obviously not given by (6) but can be written in terms of spin *S* = ¥ quantum chains. Another further extension of our model is obtained by considering an arbitrary mixture of molecules, where molecules in the same hierarchy may have distinct sizes. The results presented in [15] correspond to the particular case of this generalization where *N* = 1 (simple diffusion). For general *N* the *S* matrix we obtain in (25) is also a solution of the Yang-Baxter equation (34), but the diagonalization of the transfer matrix of the associated inhomogeneous vertex model is more complicated. The Bethe-ansatz equations in the case of asymmetric diffusion, with particles of unit size [10, 11], or with arbitrary size [15], were used to obtain the finite-size corrections of the mass gap *G*_{N} of the associated quantum chain. The real part of these finite-size corrections is governed by the dynamical critical exponent *z*, i. e.,

The calculation of the exponent *z* for the model presented in this paper, with particles of arbitrary sizes, is presently in progress [30].

**Acknowledgements**

This work was supported in part by Conselho Nacional de Desenvolvimento Científico e Tecnológico - CNPq - Brazil and by the Russian Foundation of Fundamental Investigation ( Grant 99-02-17646).

**References **

[1] A. A. Lushnikov, Zh. Éksp. Teor. Fiz. **91**, 1376 (1986) [ Links ][Sov. Phys. JETP **64**, 811 (1986)], [ Links ]Phys. Lett. A **120**, 135 (1987). [ Links ]

[2] G. Schütz, J. Stat. Phys. **71**, 471 (1993). [ Links ]

[3] F. C. Alcaraz, M. Droz, M. Henkel, and V. Rittenberg, Ann. Phys. (N.Y.) **230**, 250 (1994). [ Links ]

[4] F. C. Alcaraz and V. Rittenberg, Phys. Lett. B **324**, 377 (1993). [ Links ]

[5] F. C. Alcaraz, Int. J. Mod. Phys. B **8**, 3449 (1994). [ Links ]

[6] M. D. Grynberg and R. B. Stinchcombe, Phys. Rev. Lett. **74**, 1242 (1995). [ Links ]

[7] K. Krebs, M. P. Pfannmüller, B. Wehefritz, and H. Henrichsen, J. Stat. Phys. **78**, 1429 (1995). [ Links ]

[8] J. E. Santos, G. M. Schütz, and R. B. Stinchcombe, J. Chem. Phys. **105**, 2399 (1996). [ Links ]

[9] M. J. de Oliveira, T. Tomé, and R. Dickman, Phys. Rev. A **46**, 6294 (1992). [ Links ]

[10] L. H. Gwa and H. Spohn, Phys. Rev. Lett. **68**, 725 (1992), [ Links ]Phys. Rev. A **46**, 844 (1992). [ Links ]

[11] D. Kim, Phys. Rev. E **52**, 3512 (1995). [ Links ]

[12] D. Kim, J. Phys. A **30**, 3817 (1997). [ Links ]

[13] F. C. Alcaraz, S. Dasmahapatra and V. Rittenberg, J. Phys. A. **31**, 845 (1998). [ Links ]

[14] T. Sasamoto and M. Wadati, J. Phys. A **31**, 6057 (1998). [ Links ]

[15] F. C. Alcaraz and R. Z. Bariev, Phys. Rev. E **60**, 79 (1999). [ Links ]

[16] B. Derrida, Physics Reports **301**, 65 (1998). [ Links ]

[17] T. M. Liggett, *Stochastic Interacting Systems: Contact, Voter and Exclusion Process* (Springer Verlag, 1999) [ Links ]

[18] G. M. Schütz, *Exactly solvable models for many-body systems far from equilibrium* in "Phase Transition and Critical Phenomena". Eds. C. Domb and J.L. Lebowitz Vol 19 (to appear). [ Links ]

[19] C. N. Yang, Phys. Rev. Lett. **19**, 1312 (1967). [ Links ]

[20] B. Sutherland, Phys. Rev. B **12**, 3795 (1975). [ Links ]

[21] S. V. Pokrovskii and A. M. Tsvelick, Zh. Eksp. Teor. Fiz. **93**, 2232 (1987) [ Links ](Sov. Phys.-JETP **66**, 6 (1987) [ Links ])

[22] C. Boldrighini, G. Cosimi, S. Frigio, and G. Nuñes, J. Stat. Phys. **55**, 611 (1989). [ Links ]

[23] P. A. Ferrari, C. Kipnis and E. Saada, Ann. Prob. **19**, 226 (1991). [ Links ]

[24] P. A. Ferrari, Prob. Theory Relat. Fields **91**, 81 (1992). [ Links ]

[25] L. A. Takhtajan and L. D. Faddeev, Russ. Math. Surv. **34**,11 (1979); [ Links ]V. E. Korepin, I. G. Izergin and N. M. Bogoliubov, *Quantum Inverse Scattering Method, Correlation Functions and Algebraic Bethe Ansatz* (Cambridge University Press, Cambridge, 1992). [ Links ]

[26] F. C. Alcaraz and R. Z. Bariev, Braz. J. Phys. **30**,13 (2000) [ Links ]

[27] B. Derrida, M. R. Evans, V. Hakim, and V. Pasquier, J. Phys. A **26**, 1493 (1993). [ Links ]

[28] K. Mallick, S. Mallick and N. Rajewsky, J. Phys. A **32**, 48 (1999) [ Links ]

[29] R. J. Baxter, *Exactly Solved Models in Statistical Mechanics* (Academic Press, New York, 1982). [ Links ]

[30] F. C. Alcaraz and R. Z. Bariev, (to be published).

[31] P. P. Kulish and N. Yu. Reshetikhin, Zh. Eksp. Teor. Fiz. **80**, 214 (1981); [ Links ]C. L. Schultz, Physica A **122**, 71 (1983); [ Links ]M. Gaudin *La fonction d'onde Bethe* (Paris, Masson, 1983). [ Links ]

[32] B. Derrida and M. R. Evans, J. Phys. A **32**, 4833 (1999). [ Links ]