Nonequilibrium critical behavior in unidirectionally coupled stochastic processes

Abstract

Phase transitions from an active into an absorbing, inactive state are generically described by the critical exponents of directed percolation (DP), with upper critical dimension ${d}_{c}=4$. In the framework of single-species reaction-diffusion systems, this universality class is realized by the combined processes $\stackrel{\ensuremath{\rightarrow}}{A}A+A,$ $A+\stackrel{\ensuremath{\rightarrow}}{A}A$, and $\stackrel{\ensuremath{\rightarrow}}{A}0$. We study a hierarchy of such DP processes for particle species $A,B,\dots{}$, unidirectionally coupled via the reactions $\stackrel{\ensuremath{\rightarrow}}{A}B,\dots{}$ (with rates ${\ensuremath{\mu}}_{\mathrm{AB}},\dots{})$. When the DP critical points at all levels coincide, multicritical behavior emerges, with density exponents ${\ensuremath{\beta}}_{i}$ which are markedly reduced at each hierarchy level $i>~2$. This scenario can be understood on the basis of the mean-field rate equations, which yield ${\ensuremath{\beta}}_{i}{=1/2}^{i\ensuremath{-}1}$ at the multicritical point. Using field-theoretic renormalization-group techniques in $d=4\ensuremath{-}\ensuremath{\epsilon}$ dimensions, we identify a new crossover exponent $\ensuremath{\varphi}$, and compute $\ensuremath{\varphi}=1+O({\ensuremath{\epsilon}}^{2})$ in the multicritical regime (for small ${\ensuremath{\mu}}_{\mathrm{AB}})$ of the second hierarchy level. In the active phase, we calculate the fluctuation correction to the density exponent on the second hierarchy level, ${\ensuremath{\beta}}_{2}=1/2\ensuremath{-}\ensuremath{\epsilon}/8+O({\ensuremath{\epsilon}}^{2})$. Outside the multicritial region, we discuss the crossover to ordinary DP behavior, with the density exponent ${\ensuremath{\beta}}_{1}=1\ensuremath{-}\ensuremath{\epsilon}/6+O({\ensuremath{\epsilon}}^{2})$. Monte Carlo simulations are then employed to confirm the crossover scenario, and to determine the values for the new scaling exponents in dimensions $d<~3$, including the critical initial slip exponent. Our theory is connected to specific classes of growth processes and to certain cellular automata, and the above ideas are also applied to unidirectionally coupled pair annihilation processes. We also discuss some technical as well as conceptual problems of the loop expansion, and suggest some possible interpretations of these difficulties.