Decay of metastable phases in a model for the catalytic oxidation of CO

TL;DR

This study uses kinetic Monte Carlo simulations to analyze the decay of metastable phases in a catalytic CO oxidation model, revealing a first-order phase transition and decay mechanisms similar to equilibrium phase transformations.

Contribution

It demonstrates that metastable decay in a nonequilibrium catalytic model follows nucleation and growth theories akin to equilibrium first-order transitions, with detailed analysis of different decay regimes.

Findings

First-order phase transition observed below a critical desorption rate.

Decay mechanisms follow Kolmogorov-Johnson-Mehl-Avrami theory.

Identification of single-droplet and multidroplet decay regimes.

Abstract

We study by kinetic Monte Carlo simulations the dynamic behavior of a Ziff-Gulari-Barshad model with CO desorption for the reaction CO + O $\to$ CO$_2$ on a catalytic surface. Finite-size scaling analysis of the fluctuations and the fourth-order order-parameter cumulant show that below a critical CO desorption rate, the model exhibits a nonequilibrium first-order phase transition between low and high CO coverage phases. We calculate several points on the coexistence curve. We also measure the metastable lifetimes associated with the transition from the low CO coverage phase to the high CO coverage phase, and {\it vice versa}. Our results indicate that the transition process follows a mechanism very similar to the decay of metastable phases associated with {\it equilibrium} first-order phase transitions and can be described by the classic Kolmogorov-Johnson-Mehl-Avrami theory of phase transformation by nucleation and growth. In the present case, the desorption parameter plays the role of temperature, and the distance to the coexistence curve plays the role of an external field or supersaturation. We identify two distinct regimes, depending on whether the system is far from or close to the coexistence curve, in which the statistical properties and the system-size dependence of the lifetimes are different, corresponding to multidroplet or single-droplet decay, respectively. The crossover between the two regimes approaches the coexistence curve logarithmically with system size, analogous to the behavior of the crossover between multidroplet and single-droplet metastable decay near an equilibrium first-order phase transition.