Transition in the relaxation dynamics of a reversible diffusion-limited reaction

Transport and fluctuations in mass aggregation processes: Mobility-driven clustering

We calculate the bulk-diffusion coefficient and the conductivity in nonequilibrium conserved-mass aggregation processes on a ring. These processes involve chipping and fragmentation of masses, which diffuse on a lattice and aggregate with their neighboring masses on contact, and, under certain conditions, they exhibit a condensation transition. We find that, even in the absence of microscopic time reversibility, the systems satisfy an Einstein relation, which connects the ratio of the conductivity and the bulk-diffusion coefficient to mass fluctuation. Interestingly, when aggregation dominates over chipping, the conductivity or, equivalently, the mobility of masses, is greatly enhanced. The enhancement in the conductivity, in accordance with the Einstein relation, results in large mass fluctuations and can induce a mobility-driven clustering in the systems. Indeed, in a certain parameter regime, we show that the conductivity, along with the mass fluctuation, diverges beyond a critical density, thus characterizing the previously observed nonequilibrium condensation transition [Phys. Rev. Lett. 81, 3691 (1998)10.1103/PhysRevLett.81.3691] in terms of an instability in the conductivity. Notably, the bulk-diffusion coefficient remains finite in all cases. We find our analytic results in quite good agreement with simulations.

Hard core run and tumble particles on a one-dimensional lattice

We study the large scale behavior of a collection of hard core run and tumble particles on a one-dimensional lattice with periodic boundary conditions. Each particle has persistent motion in one direction decided by an associated spin variable until the direction of spin is reversed. We map the run and tumble model to a mass transfer model with fluctuating directed bonds. We calculate the steady-state single-site mass distribution in the mass model within a mean field approximation for larger spin-flip rates and by analyzing an appropriate coalescence-fragmentation model for small spin-flip rates. We also calculate the hydrodynamic coefficients of diffusivity and conductivity for both large and small spin-flip rates and show that the Einstein relation is violated in both regimes. We also show how the nongradient nature of the process can be taken into account in a systematic manner to calculate the hydrodynamic coefficients.

Aggregation-fragmentation-diffusion model for trail dynamics

We investigate statistical properties of trails formed by a random process incorporating aggregation, fragmentation, and diffusion. In this stochastic process, which takes place in one spatial dimension, two neighboring trails may combine to form a larger one, and also one trail may split into two. In addition, trails move diffusively. The model is defined by two parameters which quantify the fragmentation rate and the fragment size. In the long-time limit, the system reaches a steady state, and our focus is the limiting distribution of trail weights. We find that the density of trail weight has power-law tail P(w)∼w^{-γ} for small weight w. We obtain the exponent γ analytically and find that it varies continuously with the two model parameters. The exponent γ can be positive or negative, so that in one range of parameters small-weight trails are abundant and in the complementary range they are rare.

One-dimensional spin-anisotropic kinetic Ising model subject to quenched disorder

Large-scale Monte Carlo simulations are used to explore the effect of quenched disorder on one-dimensional, nonequilibrium kinetic Ising models with locally broken spin symmetry, at zero temperature (the symmetry is broken through spin-flip rates that differ for "+" and "-" spins). The model is found to exhibit a continuous phase transition to an absorbing state. The associated critical behavior is studied at zero branching rate of kinks, through analysis spreading of + and - spins, and of the kink density. Impurities exert a strong effect on the critical behavior only for a particular choice of parameters, corresponding to the strongly spin-anisotropic kinetic Ising model introduced by Majumdar [Phys Rev. Lett 86, 2301 (2001)]. Typically, disorder effects become evident for impurity strengths such that diffusion is nearly blocked. In this regime, the critical behavior is similar to that arising, for example, in the one-dimensional diluted contact process, with Griffiths-like behavior for the kink density. We find variable cluster exponents, which obey a hyperscaling relation, and are similar to those reported by Cafiero [Phys Rev. E 57, 5060 (1998)]. We also show that the isotropic two-component AB --> [symbol: see text] model is insensitive to reaction disorder, and that only logarithmic corrections arise, induced by strong disorder in the diffusion rate.

Static and dynamic phase transitions in multidimensional voting models on continua

A voting model (or a generalization of the Glauber model at zero temperature) on a multidimensional lattice is defined as a system composed of a lattice, each site of which is either empty or occupied by a single particle. The reactions of the system are such that two adjacent sites, one empty, the other occupied, may evolve to a state where both of these sites are either empty or occupied. The continuum version of this model in a D-dimensional region with a boundary is studied, and two general behaviors of such systems are investigated, the stationary behavior of the system, and the dominant way of relaxation of the system toward its stationary state. Based on the first behavior, a static phase transition (discontinuous changes in the stationary profiles of the system) is studied. Based on the second behavior, a dynamical phase transition (discontinuous changes in the relaxation times of the system) is studied. It is shown that the static phase transition is induced by the bulk reactions only, while the dynamical phase transition is a result of both bulk reactions and boundary conditions.

Perturbative calculation of one-point functions of one-dimensional single-species reaction-diffusion systems

Perturbations around autonomous one-dimensional single-species reaction-diffusion systems are investigated. It is shown that the parameter space corresponding to the autonomous systems is divided into two parts: In one part, the system is stable against the perturbations, in the sense that the largest relaxation time of the one-point functions changes continuously with perturbations. In the other part, however, the system is unstable against perturbations, so that any small perturbation drastically modifies the large-time behavior of the one-point functions.

Reaction front in an A+B-->C reaction-subdiffusion process

We study the reaction front for the process A+B-->C in which the reagents move subdiffusively. Our theoretical description is based on a fractional reaction-subdiffusion equation in which both the motion and the reaction terms are affected by the subdiffusive character of the process. We design numerical simulations to check our theoretical results, describing the simulations in some detail because the rules necessarily differ in important respects from those used in diffusive processes. Comparisons between theory and simulations are on the whole favorable, with the most difficult quantities to capture being those that involve very small numbers of particles. In particular, we analyze the total number of product particles, the width of the depletion zone, the production profile of product and its width, as well as the reactant concentrations at the center of the reaction zone, all as a function of time. We also analyze the shape of the product profile as a function of time, in particular, its unusual behavior at the center of the reaction zone.

Adsorption of reactive particles on a random catalytic chain: an exact solution

We study equilibrium properties of a catalytically activated annihilation A+A-->0 reaction taking place on a one-dimensional chain of length N (N--> infinity ) in which some segments (placed at random, with mean concentration p) possess special, catalytic properties. Annihilation reaction takes place as soon as any two A particles land onto two vacant sites at the extremities of the catalytic segment, or when any A particle lands onto a vacant site on a catalytic segment while the site at the other extremity of this segment is already occupied by another A particle. Noncatalytic segments are inert with respect to reaction and here two adsorbed A particles harmlessly coexist. For both "annealed" and "quenched" disorder in placement of the catalytic segments, we calculate exactly the disorder-averaged pressure per site. Explicit asymptotic formulas for the particle mean density and the compressibility are also presented.

Cluster approximation solution of a two-species annihilation model

A two-species reaction-diffusion model, in which particles diffuse on a one-dimensional lattice and annihilate when meeting each other, has been investigated. Mean-field equations for general choice of reaction rates have been solved exactly. The cluster mean-field approximation of the model is also studied. It is shown that the general form of large time behavior of one- and two-point functions of the number operators are determined by the diffusion rates of the two type of species, and is independent of annihilation rates.

Generalized empty-interval method applied to a class of one-dimensional stochastic models

In this work we study, on a finite and periodic lattice, a class of one-dimensional (bimolecular and single-species) reaction-diffusion models that cannot be mapped onto free-fermion models. We extend the conventional empty-interval method, also called interparticle distribution function (IPDF) method, by introducing a string function, which is simply related to relevant physical quantities. As an illustration, we specifically consider a model that cannot be solved directly by the conventional IPDF method and that can be viewed as a generalization of the voter model and/or as an epidemic model. We also consider the reversible diffusion-coagulation model with input of particles and determine other reaction-diffusion models that can be mapped onto the latter via suitable similarity transformations. Finally we study the problem of the propagation of a wave front from an inhomogeneous initial configuration and note that the mean-field scenario predicted by Fisher's equation is not valid for the one-dimensional (microscopic) models under consideration.

Exactly solvable models through the empty-interval method

The most general one dimensional reaction-diffusion model with nearest-neighbor interactions that can be solved exactly through empty-interval method has been introduced. Assuming translationally invariant initial conditions, the probability that n consecutive sites are empty, E(n), has been exactly obtained. Here, however, we do not consider reactions changing two empty neighboring sites. In the thermodynamic limit, the large-time behavior of the system has also been investigated. Releasing translationally invariance, the evolution equation for the probability that n consecutive sites, starting from the site k, are empty, E(k,n), is obtained. In the thermodynamic limit, the large time behavior of the system is also considered. Finally, the continuum limit of the model is considered and the empty-interval probability function is obtained.

Subdiffusion-limited A+A reactions

We consider the coagulation dynamics A+A-->A and A+A <==> A and the annihilation dynamics A+A-->0 for particles moving subdiffusively in one dimension. This scenario combines the "anomalous kinetics" and "anomalous diffusion" problems, each of which leads to interesting dynamics separately and to even more interesting dynamics in combination. Our analysis is based on the fractional diffusion equation.

Method of intervals for the study of diffusion-limited annihilation, A+A-->0

We introduce a method of intervals for the analysis of diffusion-limited annihilation, A+A-->0, on the line. The method leads to manageable diffusion equations whose interpretation is intuitively clear. As an example, we treat the following cases: (a) annihilation in the infinite line and in infinite (discrete) chains; (b) annihilation with input of single particles, adjacent particle pairs, and particle pairs separated by a given distance; (c) annihilation, A+A-->0, along with the birth reaction A-->3A, on finite rings, with and without diffusion.