arXiv:cond-mat/0305691 · Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics
Persistence properties of a system of coagulating and annihilating random walkers
Supriya Krishnamurthy, R. Rajesh, Oleg V. Zaboronski
Abstract
We study a d-dimensional system of diffusing particles that on contact either annihilate with probability $1/(q\ensuremath{-}1)$ or coagulate with probability $(q\ensuremath{-}2)/(q\ensuremath{-}1).$ In one dimension, the system models the zero-temperature Glauber dynamics of domain walls in the q-state Potts model. We calculate $\overline{P}(m,t),$ the probability that a randomly chosen lattice site contains a particle whose ancestors have undergone exactly $(m\ensuremath{-}1)$ coagulations. Using perturbative renormalization group analysis for $d<2,$ we show that, if the number of coagulations m is much less than the typical number $M(t),$ then $\overline{P}(m,t)\ensuremath{\sim}{m}^{\ensuremath{\zeta}/d}{t}^{\ensuremath{-}\ensuremath{\theta}},$ with $\ensuremath{\theta}=dQ+Q(Q\ensuremath{-}1/2)\ensuremath{\epsilon}+O({\ensuremath{\epsilon}}^{2}),$ $\ensuremath{\zeta}=(2Q\ensuremath{-}1)\ensuremath{\epsilon}+(2Q\ensuremath{-}1)(Q\ensuremath{-}1)(1/2+AQ){\ensuremath{\epsilon}}^{2}+O({\ensuremath{\epsilon}}^{3}),$ where $Q=(q\ensuremath{-}1)/q,$ $\ensuremath{\epsilon}=2\ensuremath{-}d$ and $A=\ensuremath{-}0.006\dots{}.$ $M(t)$ is shown to scale as $M(t)\ensuremath{\sim}{t}^{d/2\ensuremath{-}\ensuremath{\delta}},$ where $\ensuremath{\delta}=d(1\ensuremath{-}Q)+(Q\ensuremath{-}1)(Q\ensuremath{-}1/2)\ensuremath{\epsilon}+O({\ensuremath{\epsilon}}^{2}).$ In two dimensions, we show that $\overline{P}(m,t)\ensuremath{\sim}\mathrm{ln}{(t)}^{Q(3\ensuremath{-}2Q)}\mathrm{ln}{(m)}^{(2Q\ensuremath{-}{1)}^{2}}{t}^{\ensuremath{-}2Q}$ for $m\ensuremath{\ll}{t}^{2Q\ensuremath{-}1}.$ We also derive an exact nonperturbative relation between the exponents: namely $\ensuremath{\delta}(Q)=\ensuremath{\theta}(1\ensuremath{-}Q).$ The one-dimensional results corresponding to $\ensuremath{\epsilon}=1$ are compared with results from Monte Carlo simulations.