Reaction diffusion models in one dimension with disorder

• Published in: Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics Vol. 60; no. 2 Pt A; p. 1212
• Main Authors: Doussal, P L, Monthus, C
• Format: Journal Article
• Language: English
• Published: United States 01.08.1999
• ISSN: 1063-651X
• DOI: 10.1103/physreve.60.1212

We study a large class of one-dimensional reaction diffusion models with quenched disorder using a real space renormalization group method (RSRG) which yields exact results at large time. Particles (e.g., of several species) undergo diffusion with random local bias (Sinai model) and may react upon meeting. We obtain a detailed description of the asymptotic states (i.e., attractive fixed points of the RSRG), such as the large time decay of the density of each specie, their associated universal amplitudes, and the spatial distribution of particles. We also derive the spectrum of nontrivial exponents which characterize the convergence towards the asymptotic states. For reactions which lead to several possible asymptotic states separated by unstable fixed points, we analyze the dynamical phase diagram and obtain the critical exponents characterizing the transitions. We also obtain a detailed characterization of the persistence properties for single particles as well as more complex patterns. We compute the decay exponents for the probability of no crossing of a given point by, respectively, the single particle trajectories (theta) or the thermally averaged packets (theta). The generalized persistence exponents associated to n crossings are also obtained. Specifying to the process A+A--> or A with probabilities (r,1-r), we compute exactly the exponents delta(r) and psi(r) characterizing the survival up to time t of a domain without any merging or with mergings, respectively, and the exponents deltaA(r) and psiA(r) characterizing the survival up to time t of a particle A without any coalescence or with coalescences, respectively. theta, psi, and delta obey hypergeometric equations and are numerically surprisingly close to pure system exponents (though associated to a completely different diffusion length). The effect of additional disorder in the reaction rates, as well as some open questions, are also discussed.