RESUMEN
We study by kinetic Monte Carlo simulations the dynamic behavior of a Ziff-Gulari-Barshad model with CO desorption for the reaction CO + O-->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 vice versa. Our results indicate that the transition process follows a mechanism very similar to the decay of metastable phases associated with 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.
RESUMEN
We present a kinetic Monte Carlo study of the dynamical response of a Ziff-Gulari-Barshad model for CO oxidation with CO desorption to periodic variation of the CO pressure. We use a square-wave periodic pressure variation with parameters that can be tuned to enhance the catalytic activity. We produce evidence that, below a critical value of the desorption rate, the driven system undergoes a dynamic phase transition between a CO2 productive phase and a nonproductive one at a critical value of the period and waveform of the pressure oscillation. At the dynamic phase transition the period-averaged CO2 production rate is significantly increased and can be used as a dynamic order parameter. We perform a finite-size scaling analysis that indicates the existence of power-law singularities for the order parameter and its fluctuations, yielding estimated critical exponent ratios beta/nu approximately 0.12 and gamma/nu approximately 1.77. These exponent ratios, together with theoretical symmetry arguments and numerical data for the fourth-order cumulant associated with the transition, give reasonable support for the hypothesis that the observed nonequilibrium dynamic phase transition is in the same universality class as the two-dimensional equilibrium Ising model.
RESUMEN
Using both analytical and simulational methods, we study low-temperature nucleation rates in kinetic Ising lattice-gas models that evolve under two different Arrhenius dynamics that interpose between the Ising states a transition state representing a local energy barrier. The two dynamics are the transition-state approximation [T. Ala-Nissila, J. Kjoll, and S. C. Ying, Phys. Rev. B 46, 846 (1992)] and the one-step dynamic [H. C. Kang and W. H. Weinberg, J. Chem. Phys. 90, 2824 (1989)]. Even though they both obey detailed balance and are here applied to a situation that does not conserve the order parameter, we find significant differences between the nucleation rates observed with the two dynamics, and between them and the standard Glauber dynamic [R. J. Glauber, J. Math. Phys. 4, 294 (1963)], which does not contain transition states. Our results show that great care must be exercised when devising kinetic Monte Carlo transition rates for specific physical or chemical systems.