Metastability in Schloegl's second model for autocatalysis: Lattice-gas realization with particle diffusion

The role of fluctuations in bistability and oscillations during the H2 + O2 reaction on nanosized rhodium crystals

A combined experimental and theoretical study is presented of fluctuations observed by field ion microscopy in the catalytic reaction of water production on a rhodium tip. A stochastic approach is developed to provide a comprehensive understanding of the different phenomena observed in the experiment, including burst noise manifesting itself in a bistability regime, noisy oscillations, and nanopatterns with a cross-like oxidized zone separating the surface into four quadrants centered on the {111} facets. The study is based on a stochastic model numerically simulating the processes of adsorption, desorption, reaction, and transport. The surface diffusion of hydrogen is described as a percolation process dominated by large clusters corresponding to the four quadrants. The model reproduces the observed phenomena in the ranges of temperature, pressures, and electric field of the experiment.

Schloegl's second model for autocatalysis on hypercubic lattices: Dimension dependence of generic two-phase coexistence

Schloegl's second model on a (d ≥ 2)-dimensional hypercubic lattice involves: (i) spontaneous annihilation of particles with rate p and (ii) autocatalytic creation of particles at vacant sites at a rate proportional to the number of suitable pairs of neighboring particles. This model provides a prototype for nonequilibrium discontinuous phase transitions. However, it also exhibits nontrivial generic two-phase coexistence: Stable populated and vacuum states coexist for a finite range, pf(d)<p<pe(d), spanned by the orientation-dependent stationary points for planar interfaces separating these states. Analysis of interface dynamics from kinetic Monte Carlo simulation and from discrete reaction-diffusion equations (dRDEs) obtained from truncation of the exact master equation, reveals that pe(f)∼0.2113765+ce(f)/d as d→∞, where Δc=ce-cf≈0.014. A metastable populated state persists above pe(d) up to a spinodal p=ps(d), which has a well-defined limit ps(d→∞)=1/4. The dRDEs display artificial propagation failure, absent in the stochastic model due to fluctuations. This feature is amplified for increasing d, thus complicating our analysis.

Catalytic conversion reactions mediated by single-file diffusion in linear nanopores: hydrodynamic versus stochastic behavior

We analyze the spatiotemporal behavior of species concentrations in a diffusion-mediated conversion reaction which occurs at catalytic sites within linear pores of nanometer diameter. Diffusion within the pores is subject to a strict single-file (no passing) constraint. Both transient and steady-state behavior is precisely characterized by kinetic Monte Carlo simulations of a spatially discrete lattice-gas model for this reaction-diffusion process considering various distributions of catalytic sites. Exact hierarchical master equations can also be developed for this model. Their analysis, after application of mean-field type truncation approximations, produces discrete reaction-diffusion type equations (mf-RDE). For slowly varying concentrations, we further develop coarse-grained continuum hydrodynamic reaction-diffusion equations (h-RDE) incorporating a precise treatment of single-file diffusion in this multispecies system. The h-RDE successfully describe nontrivial aspects of transient behavior, in contrast to the mf-RDE, and also correctly capture unreactive steady-state behavior in the pore interior. However, steady-state reactivity, which is localized near the pore ends when those regions are catalytic, is controlled by fluctuations not incorporated into the hydrodynamic treatment. The mf-RDE partly capture these fluctuation effects, but cannot describe scaling behavior of the reactivity.