First order phase transition in a modified Ziff-Gulari-Barshad model with self-oscillating reactant coverages

Using kinetic Monte Carlo simulations, we study the effect of oscillatory kinetics due to surface reconstructions on Ziff-Gulari-Barshad (ZGB) model discontinuous phase transition. To investigate the transition, we do extensive finite size scaling analysis. It is found that the discontinuous transition still exists. On inclusion of desorption in the model, the order-parameter probability distribution broadens but remains bimodal. That is, the first-order phase transition becomes weaker with increase in desorption rate.