Reaction-controlled diffusion: Monte Carlo simulations

TL;DR

This study uses Monte Carlo simulations to analyze a reaction-controlled diffusion model, revealing subdiffusive behavior of particles near phase transitions and highlighting the impact of fluctuations on particle displacement distributions.

Contribution

It provides the first detailed simulation analysis of the reaction-controlled diffusion model, confirming mean-field predictions and uncovering non-Gaussian displacement distributions.

Findings

A particles exhibit subdiffusive behavior with alpha_A = alpha_B.

Displacement distributions are non-Gaussian and scale dynamically.

Fluctuations enhance particle localization and large-distance traversal.

Abstract

We study the coupled two-species non-equilibrium reaction-controlled diffusion model introduced by Trimper et al. [Phys. Rev. E 62, 6071 (2000)] by means of detailed Monte Carlo simulations in one and two dimensions. Particles of type A may independently hop to an adjacent lattice site provided it is occupied by at least one B particle. The B particle species undergoes diffusion-limited reactions. In an active state with nonzero, essentially homogeneous B particle saturation density, the A species displays normal diffusion. In an inactive, absorbing phase with exponentially decaying B density, the A particles become localized. In situations with algebraic decay rho_B(t) ~ t^{-alpha_B}, as occuring either at a non-equilibrium continuous phase transition separating active and absorbing states, or in a power-law inactive phase, the A particles propagate subdiffusively with mean-square displacement <x(t)^2_A> ~ t^{1-alpha_A}. We find that within the accuracy of our simulation data, \alpha_A = \alpha_B as predicted by a simple mean-field approach. This remains true even in the presence of strong spatio-temporal fluctuations of the B density. However, in contrast with the mean-field results, our data yield a distinctly non-Gaussian A particle displacement distribution n_A(x,t) that obeys dynamic scaling and looks remarkably similar for the different processes investigated here. Fluctuations of effective diffusion rates cause a marked enhancement of n_A(x,t) at low displacements |x|, indicating a considerable fraction of practically localized A particles, as well as at large traversed distances.