We investigate the quantum reaction-diffusion dynamics of fermionic particles which coherently hop in a one-dimensional lattice and undergo annihilation reactions. The latter are modelled as dissipative processes which involve losses of pairs 2A→∅, triplets 3A→∅, and quadruplets 4A→∅ of neighboring particles. When considering classical particles, the corresponding decay of their density in time follows an asymptotic power-law behavior. The associated exponent in one dimension is different from the mean-field prediction whenever diffusive mixing is not too strong and spatial correlations are relevant. This specifically applies to 2A→∅, while the mean-field power-law prediction just acquires a logarithmic correction for 3A→∅ and is exact for 4A→∅. A mean-field approach is also valid, for all the three processes, when the diffusive mixing is strong, i.e., in the so-called reaction-limited regime. Here we show that the picture is different for quantum systems. We consider the quantum reaction-limited regime and we show that for all the three processes power-law behavior beyond mean field is present as a consequence of quantum coherences, which are not related to space dimensionality. The decay in 3A→∅ is further, highly intricate, since the power-law behavior therein only appears within an intermediate time window, while at long times the density decay is not power law. Our results show that emergent critical behavior in quantum dynamics has a markedly different origin, based on quantum coherences, to that applying to classical critical phenomena, which is, instead, solely determined by the relevance of spatial correlations.
Quantum reaction-transport systems consist of coherently propagating particles that irreversibly react upon encounter. Their relaxation is commonly classified as reaction-limited or transport-limited, depending on the relative timescales of reaction and particle transport. We show that this expectation fails in low dimensions by studying the quantum binary annihilation,
A+A→∅, where reaction processes acquire singular fluctuation corrections below the upper critical dimension
Dc=2. Consequently, fluctuations dominate the asymptotic kinetics for
D<2. In one dimension, they render mean-field relaxation transient and drive the system toward a transport-limited regime with
n∼t−1/2, governed by a quantum-Zeno scale even for arbitrarily weak loss. At
D=2, the kinetics acquires logarithmic corrections, whereas above two dimensions mean-field scaling is asymptotically restored. At and below the upper critical dimension, mean-field behavior can nevertheless persist over parametrically long crossover times before the asymptotic fluctuation-dominated regime emerges. We also show that the same fluctuations generate effective elastic collisions despite the absence of microscopic coherent interactions. In
D>1, these collisions redistribute momentum populations and are parametrically faster than losses in the transport-limited regime, allowing the system to approach a quasi-stationary thermal state.
Abstract We investigate quantum reaction–diffusion (RD) systems in one-dimension with bosonic particles that coherently hop in a lattice, and when brought in range react dissipatively. Such reactions involve binary annihilation ( A + A → ∅ ) and coagulation ( A + A → A ) of particles at distance d. We consider the reaction-limited regime, where dissipative reactions take place at a rate that is small compared to that of coherent hopping. In classical RD systems, this regime is correctly captured by the mean-field approximation. In quantum RD systems, for noninteracting fermionic systems, the reaction-limited regime recently attracted considerable attention because it has been shown to give universal power law decay beyond mean field for the density of particles as a function of time. Here, we address the question whether such universal behavior is present also in the case of the noninteracting Bose gas. We show that beyond mean-field density decay for bosons is possible only for reactions that allow for destructive interference of different decay channels. Furthermore, we study an absorbing-state phase transition induced by the competition between branching A → A + A , decay A → ∅ and coagulation A + A → A . We find a stationary phase-diagram, where a first and a second-order transition line meet at a bicritical point which is described by tricritical directed percolation. These results show that quantum statistics significantly impact on both the stationary and the dynamical universal behavior of quantum RD systems.
We consider the quantum reaction-diffusion dynamics in
d spatial dimensions of a Fermi gas subject to binary annihilation reactions $A+A\ensuremath{\rightarrow}\ensuremath{\emptyset}$. These systems display collective nonequilibrium long-time behavior, which is signalled by an algebraic decay of the particle density. Building on the Keldysh formalism, we devise a field theoretical approach for the reaction-limited regime, where annihilation reactions are scarce. Combining a perturbative expansion of the dissipative interaction with Euler-hydrodynamic scaling limit, we derive a description in terms of a large-scale universal kinetic equation. Our approach shows how the time-dependent generalized Gibbs ensemble assumption, which is often employed for treating low-dimensional nonequilibrium dissipative systems, emerges from systematic diagrammatics. It also allows us to exactly compute---for arbitrary spatial dimension---the decay exponent of the particle density. The latter is based on the large-scale description of the quantum dynamics and it differs from the mean-field prediction even in dimension larger than one. We moreover consider spatially inhomogeneous setups involving an external potential. In confined systems the density decay is accelerated towards the mean-field algebraic behavior, while for deconfined scenarios the power-law decay is replaced by a slower nonalgebraic decay.
We generalize nonperturbative renormalization group methods to nonequilibrium critical phenomena. Within this formalism, reaction-diffusion processes are described by a scale-dependent effective action, the flow of which is derived. We investigate branching and annihilating random walks with an odd number of offspring. Along with recovering their universal physics (described by the directed percolation universality class), we determine their phase diagrams and predict that a transition occurs even in three dimensions, contrarily to what perturbation theory suggests.
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
dc=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,…, 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.