Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2111.09658 · Journal of Physics A Mathematical and Theoretical

Matrix product solution for a partially asymmetric 1D lattice gas with a free defect

Ivan Lobaskin, M. R. Evans, Kirone Mallick

Published Mar 21, 2022• 6 citations
arXivHTMLPDF

Abstract

Abstract A one-dimensional, driven lattice gas with a freely moving, driven defect particle is studied. Although the dynamics of the defect are simply biased diffusion, it disrupts the local density of the gas, creating nontrivial nonequilibrium steady states. The phase diagram is derived using mean field theory and comprises three phases. In two phases, the defect causes small localized perturbations in the density profile. In the third, it creates a shock, with two regions at different bulk densities. When the hopping rates satisfy a particular condition (that the products of the rates of the gas and defect are equal), it is found that the steady state can be solved exactly using a two-dimensional matrix product ansatz. This is used to derive the phase diagram for that case exactly and obtain exact asymptotic and finite size expressions for the density profiles and currents in all phases. In particular, the front width in the shock phase on a system of size L is found to scale as L 1/2, which is not predicted by mean field theory. The results are found to agree well with Monte Carlo simulations.

DOI
10.1088/1751-8121/ac5f7a

Related papers

Matrix product steady states as superposition of product shock measures in 1D driven systems

Farhad H Jafarpour, Seyedeh Raziyeh Masharian

It is known that exact travelling wave solutions exist for families of (n + 1)-states stochastic one-dimensional non-equilibrium lattice models with open boundaries provided that some constraints on the reaction rates are fulfilled. These solutions describe the diffusive motion of a product shock or a domain wall with the dynamics of a simple biased random walker. The steady state of these systems can be written in terms of linear superposition of such shocks or domain walls. These steady states can also be expressed in a matrix product form. We show that in this case the associated quadratic algebra of the system has always a twodimensional representation with a generic structure. A couple of examples for n = 1 and n = 2 cases are presented. On the macroscopic level nonlinear hydrodynamic equations, such as the Burgers equation or the Fisher equation, might exhibit shocks in some cases [1, 2]. On the microscopic level some stochastic non-equilibrium one-dimensional lattice models with open boundaries might also develop shocks provided that some constrained on the reaction rates are fulfilled [3, 4, 5]. In [6] exact travelling wave solutions have been obtained for three families of two-states drivendiffusive models in one-dimension with nearest-neighbors interactions. These three models are the Partially Asymmetric Simple Exclusion Process (PASEP) [7, 8, 9], the Branching-Coalescing Random Walk (BCRW) and the Asymmetric Kawasaki-Glauber Process (AKGP). By taking a Bernoulli shock distribution as initial distribution, the shock distribution in these models evolves in time into a linear combination of similar distributions with different shock positions. The time evolution of the shock position in these models is similar to the dynamics of a random walker moving on a one-dimensional lattice with reflecting boundaries [6]. This property of the travelling wave solution is not limited to the two-states systems with nearest-neighbor interactions. Shocks in threestates

•Nov 27, 2012•28 citations

Matrix product steady states as superposition of product shock measures in 1D driven systems

Farhad H Jafarpour, S R Masharian

It is known that exact traveling wave solutions exist for families of ( n +1)-states stochastic one-dimensional non-equilibrium lattice models with open boundaries provided that some constraints on the reaction rates are fulfilled. These solutions describe the diffusive motion of a product shock or a domain wall with the dynamics of a simple biased random walker. The steady state of these systems can be written in terms of linear superposition of such shocks or domain walls. These steady states can also be expressed in a matrix product form. We show that, in this case, the associated quadratic algebra of the system always has a two-dimensional representation with a generic structure. A couple of examples for the n = 1 and 2 cases are presented.

arXiv:0707.4341•Journal of Statistical Mechanics Theory and Experiment•Oct 24, 2007•25 citations

Bethe Ansatz solution for a defect particle in the asymmetric exclusion process

Bernard Derrida, M. R. Evans

Abstract: The asymmetric exclusion process on a ring in one-dimension is considered with a single defect particle. The steady state has previously been solved by a matrix product method. Here we use the Bethe ansatz to solve exactly for the long time limit behaviour of the generating function of the distance travelled by the defect particle. This allows us to recover steady state properties known from the matrix approach such as the velocity, and obtain new results such as the diffusion constant of the defect particle. In the case where the defect particle is a second class particle we determine the large deviation function and show that in a certain range the distribution of the distance travelled about the mean is Gaussian. Moreover the variance (diffusion constant) grows as L 1/2 where L is the system size. This behaviour can be related to the superdiffusive spreading of excess mass fluctuations on an infinite system. In the case where the defect particle produces a shock, our expressions for the velocity and the diffusion constant coincide with those calculated previously for an infinite system by Ferrari and Fontes.

•Jan 1, 1999•46 citations

Exact solutions for driven-diffusive systems: Excess-mass formation and phase transitions

Marko Wölki

A totally asymmetric exclusion process on a ring is investigated in which particles can move one or two sites. Special attention is spent to the high-speed case where particles are not allowed to move a single site if they could move two sites. The stationary state is calculated exactly in the framework of the matrix-product ansatz. Independently of the update this process evolves into subspace of the configuration space and can lead to the formation of an excess hole. One observes two phases where its velocity takes different values which are calculated exactly from the normalization-generating function. Numerically computed density profiles show an interesting algebraic form as a limit of a shock profile. For continuous time the process turns out to be related to the ASEP with a single defect particle and for synchronous update it leads to a natural defect dynamics. For the general definition of the process with maximum velocity two, from an exact analysis of the two-particle sector we find that the distribution of inter-particle distance can oscillate. This is underlined in the thermodynamic limit by an improved mean-field theory which shows a remarkably good agreement with simulations. The main focus of this work is the exact solution of lattice models by matrix-product ansatz. We finally obtain stationary states also for related processes. Especially considerable is a new formulation of steady states for parallel dynamics like the ASEP with open boundaries as a product of a pair-factorized and a matrix-product state.

•Jul 27, 2009•0 citations

Bethe ansatz solution for a defect particle in the asymmetric exclusion process

Bernard Derrida, M. R. Evans

The asymmetric exclusion process on a ring in one dimension is considered with a single defect particle. The steady state has previously been solved by a matrix product method. Here we use the Bethe ansatz to solve exactly for the long time limit behaviour of the generating function of the distance travelled by the defect particle. This allows us to recover steady state properties known from the matrix approach such as the velocity, and obtain new results such as the diffusion constant of the defect particle. In the case where the defect particle is a second-class particle we determine the large deviation function and show that in a certain range the distribution of the distance travelled about the mean is Gaussian. Moreover, the variance (diffusion constant) grows as L 1/2 where L is the system size. This behaviour can be related to the superdiffusive spreading of excess mass fluctuations on an infinite system. In the case where the defect particle produces a shock, our expressions for the velocity and the diffusion constant coincide with those calculated previously for an infinite system by Ferrari and Fontes.

arXiv:cond-mat/9902133•Journal of Physics A Mathematical and General•Jan 1, 1999•72 citations

Towards a matrix product Ansatz in two dimensions

Chandraniva Guha Ray, Aikya Banerjee, P. K. Mohanty

Abstract Matrix product ansatz (MPA) is a powerful framework for constructing exact steady state weights of one dimensional non-equilibrium stochastic processes; but its generalization to higher dimensions is limited. Here, we introduce the MPA formalism for two dimensions (2D). As a concrete application, we introduce and exactly solve a non-conserved assisted exclusion model (NAEM) in one and two dimensions with constrained hopping and local birth–death dynamics: a particle can hop to a neighbouring site only when exactly one of its neighbouring sites is vacant, while creation and annihilation occur exclusively at sites whose neighbours are all occupied. The MPA yields exact steady-state weights and provides a systematic method to compute observables such as density moments and particle currents. In the particle-conserving limit, the system undergoes an absorbing phase transition at the critical density ρ c =1/2 with order-parameter exponent β=3. We further show that the steady state of the NAEM maps exactly onto the well-studied hard-square lattice gas with nearest-neighbour exclusion, thereby providing a nonequilibrium dynamical route to realizing equilibrium states of constrained lattice gases. Our work generalizes matrix-product methods beyond one dimension, establishing a systematic approach to exact solutions of interacting stochastic systems in 2D.

arXiv:2605.20160•Journal of Physics A Mathematical and Theoretical•Oct 23, 2026•0 citations
Abstract A one-dimensional, driven lattice gas with a freely moving, driven defect particle is studied. Although the dynamics of the defect are simply biased diffusion, it disrupts the local density of the gas, creating nontrivial nonequilibrium steady states. The phase diagram is derived using mean field theory and comprises three phases. In two phases, the defect causes small localized perturbations in the density profile. In the third, it creates a shock, with two regions at different bulk densities. When the hopping rates satisfy a particular condition (that the products of the rates of the gas and defect are equal), it is found that the steady state can be solved exactly using a two-dimensional matrix product ansatz. This is used to derive the phase diagram for that case exactly and obtain exact asymptotic and finite size expressions for the density profiles and currents in all phases. In particular, the front width in the shock phase on a system of size L is found to scale as L 1/2, which is not predicted by mean field theory. The results are found to agree well with Monte Carlo simulations.
It is known that exact travelling wave solutions exist for families of (n + 1)-states stochastic one-dimensional non-equilibrium lattice models with open boundaries provided that some constraints on the reaction rates are fulfilled. These solutions describe the diffusive motion of a product shock or a domain wall with the dynamics of a simple biased random walker. The steady state of these systems can be written in terms of linear superposition of such shocks or domain walls. These steady states can also be expressed in a matrix product form. We show that in this case the associated quadratic algebra of the system has always a twodimensional representation with a generic structure. A couple of examples for n = 1 and n = 2 cases are presented. On the macroscopic level nonlinear hydrodynamic equations, such as the Burgers equation or the Fisher equation, might exhibit shocks in some cases [1, 2]. On the microscopic level some stochastic non-equilibrium one-dimensional lattice models with open boundaries might also develop shocks provided that some constrained on the reaction rates are fulfilled [3, 4, 5]. In [6] exact travelling wave solutions have been obtained for three families of two-states drivendiffusive models in one-dimension with nearest-neighbors interactions. These three models are the Partially Asymmetric Simple Exclusion Process (PASEP) [7, 8, 9], the Branching-Coalescing Random Walk (BCRW) and the Asymmetric Kawasaki-Glauber Process (AKGP). By taking a Bernoulli shock distribution as initial distribution, the shock distribution in these models evolves in time into a linear combination of similar distributions with different shock positions. The time evolution of the shock position in these models is similar to the dynamics of a random walker moving on a one-dimensional lattice with reflecting boundaries [6]. This property of the travelling wave solution is not limited to the two-states systems with nearest-neighbor interactions. Shocks in threestates
It is known that exact traveling wave solutions exist for families of ( n +1)-states stochastic one-dimensional non-equilibrium lattice models with open boundaries provided that some constraints on the reaction rates are fulfilled. These solutions describe the diffusive motion of a product shock or a domain wall with the dynamics of a simple biased random walker. The steady state of these systems can be written in terms of linear superposition of such shocks or domain walls. These steady states can also be expressed in a matrix product form. We show that, in this case, the associated quadratic algebra of the system always has a two-dimensional representation with a generic structure. A couple of examples for the n = 1 and 2 cases are presented.
Abstract: The asymmetric exclusion process on a ring in one-dimension is considered with a single defect particle. The steady state has previously been solved by a matrix product method. Here we use the Bethe ansatz to solve exactly for the long time limit behaviour of the generating function of the distance travelled by the defect particle. This allows us to recover steady state properties known from the matrix approach such as the velocity, and obtain new results such as the diffusion constant of the defect particle. In the case where the defect particle is a second class particle we determine the large deviation function and show that in a certain range the distribution of the distance travelled about the mean is Gaussian. Moreover the variance (diffusion constant) grows as L 1/2 where L is the system size. This behaviour can be related to the superdiffusive spreading of excess mass fluctuations on an infinite system. In the case where the defect particle produces a shock, our expressions for the velocity and the diffusion constant coincide with those calculated previously for an infinite system by Ferrari and Fontes.
A totally asymmetric exclusion process on a ring is investigated in which particles can move one or two sites. Special attention is spent to the high-speed case where particles are not allowed to move a single site if they could move two sites. The stationary state is calculated exactly in the framework of the matrix-product ansatz. Independently of the update this process evolves into subspace of the configuration space and can lead to the formation of an excess hole. One observes two phases where its velocity takes different values which are calculated exactly from the normalization-generating function. Numerically computed density profiles show an interesting algebraic form as a limit of a shock profile. For continuous time the process turns out to be related to the ASEP with a single defect particle and for synchronous update it leads to a natural defect dynamics. For the general definition of the process with maximum velocity two, from an exact analysis of the two-particle sector we find that the distribution of inter-particle distance can oscillate. This is underlined in the thermodynamic limit by an improved mean-field theory which shows a remarkably good agreement with simulations. The main focus of this work is the exact solution of lattice models by matrix-product ansatz. We finally obtain stationary states also for related processes. Especially considerable is a new formulation of steady states for parallel dynamics like the ASEP with open boundaries as a product of a pair-factorized and a matrix-product state.
The asymmetric exclusion process on a ring in one dimension is considered with a single defect particle. The steady state has previously been solved by a matrix product method. Here we use the Bethe ansatz to solve exactly for the long time limit behaviour of the generating function of the distance travelled by the defect particle. This allows us to recover steady state properties known from the matrix approach such as the velocity, and obtain new results such as the diffusion constant of the defect particle. In the case where the defect particle is a second-class particle we determine the large deviation function and show that in a certain range the distribution of the distance travelled about the mean is Gaussian. Moreover, the variance (diffusion constant) grows as L 1/2 where L is the system size. This behaviour can be related to the superdiffusive spreading of excess mass fluctuations on an infinite system. In the case where the defect particle produces a shock, our expressions for the velocity and the diffusion constant coincide with those calculated previously for an infinite system by Ferrari and Fontes.
Abstract Matrix product ansatz (MPA) is a powerful framework for constructing exact steady state weights of one dimensional non-equilibrium stochastic processes; but its generalization to higher dimensions is limited. Here, we introduce the MPA formalism for two dimensions (2D). As a concrete application, we introduce and exactly solve a non-conserved assisted exclusion model (NAEM) in one and two dimensions with constrained hopping and local birth–death dynamics: a particle can hop to a neighbouring site only when exactly one of its neighbouring sites is vacant, while creation and annihilation occur exclusively at sites whose neighbours are all occupied. The MPA yields exact steady-state weights and provides a systematic method to compute observables such as density moments and particle currents. In the particle-conserving limit, the system undergoes an absorbing phase transition at the critical density ρ c =1/2 with order-parameter exponent β=3. We further show that the steady state of the NAEM maps exactly onto the well-studied hard-square lattice gas with nearest-neighbour exclusion, thereby providing a nonequilibrium dynamical route to realizing equilibrium states of constrained lattice gases. Our work generalizes matrix-product methods beyond one dimension, establishing a systematic approach to exact solutions of interacting stochastic systems in 2D.