On Matrix Product States for Periodic Boundary Conditions

The possibility of a matrix product representation for eigenstates with energy and momentum zero of a general m-state quantum spin Hamiltonian with nearest neighbour interaction and periodic boundary condition is considered. The quadratic algebra used for this representation is generated by 2m operators which fulfil m^2 quadratic relations and is endowed with a trace. It is shown that {\em not} every eigenstate with energy and momentum zero can be written as matrix product state. An explicit counter-example is given. This is in contrast to the case of open boundary conditions where every zero energy eigenstate can be written as a matrix product state using a Fock-like representation of the same quadratic algebra.Comment: 7 pages, late

The Study of Shocks in Three-States Driven-Diffusive Systems: A Matrix Product Approach

We study the shock structures in three-states one-dimensional driven-diffusive systems with nearest neighbors interactions using a matrix product formalism. We consider the cases in which the stationary probability distribution function of the system can be written in terms of superposition of product shock measures. We show that only three families of three-states systems have this property. In each case the shock performs a random walk provided that some constraints are fulfilled. We calculate the diffusion coefficient and drift velocity of shock for each family.Comment: 15 pages, Accepted for publication in Journal of Statistical Mechanics: Theory and Experiment (JSTAT

Universality properties of the stationary states in the one-dimensional coagulation-diffusion model with external particle input

We investigate with the help of analytical and numerical methods the reaction A+A->A on a one-dimensional lattice opened at one end and with an input of particles at the other end. We show that if the diffusion rates to the left and to the right are equal, for large x, the particle concentration c(x) behaves like As/x (x measures the distance to the input end). If the diffusion rate in the direction pointing away from the source is larger than the one corresponding to the opposite direction the particle concentration behaves like Aa/sqrt(x). The constants As and Aa are independent of the input and the two coagulation rates. The universality of Aa comes as a surprise since in the asymmetric case the system has a massive spectrum.Comment: 27 pages, LaTeX, including three postscript figures, to appear in J. Stat. Phy

Finite-dimensional representation of the quadratic algebra of a generalized coagulation-decoagulation model

The steady-state of a generalized coagulation-decoagulation model on a one-dimensional lattice with reflecting boundaries is studied using a matrix-product approach. It is shown that the quadratic algebra of the model has a four-dimensional representation provided that some constraints on the microscopic reaction rates are fulfilled. The dynamics of a product shock measure with two shock fronts, generated by the Hamiltonian of this model, is also studied. It turns out that the shock fronts move on the lattice as two simple random walkers which repel each other provided that the same constraints on the microscopic reaction rates are satisfied.Comment: Minor revision

Stochastic Models on a Ring and Quadratic Algebras. The Three Species Diffusion Problem

The stationary state of a stochastic process on a ring can be expressed using traces of monomials of an associative algebra defined by quadratic relations. If one considers only exclusion processes one can restrict the type of algebras and obtain recurrence relations for the traces. This is possible only if the rates satisfy certain compatibility conditions. These conditions are derived and the recurrence relations solved giving representations of the algebras.Comment: 12 pages, LaTeX, Sec. 3 extended, submitted to J.Phys.

Construction of a matrix product stationary state from solutions of finite size system

Stationary states of stochastic models, which have $N$ states per site, in matrix product form are considered. First we give a necessary condition for the existence of a finite $M$-dimensional matrix product state for any ${N,M}$. Second, we give a method to construct the matrices from the stationary states of small size system when the above condition and $N\le M$ are satisfied. Third, the method by which one can check that the obtained matrices are valid for any system size is presented for the case where $M=N$ is satisfied. The application of our methods is explained using three examples: the asymmetric exclusion process, a model studied in [F. H. Jafarpour: J. Phys. A: Math. Gen. 36 (2003) 7497] and a hybrid of both of the models.Comment: 22 pages, no figure. Major changes: sec.3 was shortened; the list of references were changed. This is the final version, which will appear in J.Phys.

Discontinuous Phase Transition in an Exactly Solvable One-Dimensional Creation-Annihilation System

An exactly solvable reaction-diffusion model consisting of first-class particles in the presence of a single second-class particle is introduced on a one-dimensional lattice with periodic boundary condition. The number of first-class particles can be changed due to creation and annihilation reactions. It is shown that the system undergoes a discontinuous phase transition in contrast to the case where the density of the second-class particles is finite and the phase transition is continuous.Comment: Revised, 8 pages, 1 EPS figure. Accepted for publication in Journal of Statistical Mechanics: theory and experimen

First Order Phase Transition in a Reaction-Diffusion Model With Open Boundary: The Yang-Lee Theory Approach

A coagulation-decoagulation model is introduced on a chain of length L with open boundary. The model consists of one species of particles which diffuse, coagulate and decoagulate preferentially in the leftward direction. They are also injected and extracted from the left boundary with different rates. We will show that on a specific plane in the space of parameters, the steady state weights can be calculated exactly using a matrix product method. The model exhibits a first-order phase transition between a low-density and a high-density phase. The density profile of the particles in each phase is obtained both analytically and using the Monte Carlo Simulation. The two-point density-density correlation function in each phase has also been calculated. By applying the Yang-Lee theory we can predict the same phase diagram for the model. This model is further evidence for the applicability of the Yang-Lee theory in the non-equilibrium statistical mechanics context.Comment: 10 Pages, 3 Figures, To appear in Journal of Physics A: Mathematical and Genera