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 sing 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 nonequilibrium statistical mechanics context.