Phase diagram and critical behavior of the pair annihilation model

TL;DR

This study investigates the critical behavior of the pair annihilation model with diffusion across multiple dimensions, combining analytical approximations and simulations to understand how diffusion influences survival thresholds and phase transitions.

Contribution

It provides new simulation data confirming the critical diffusion thresholds and the universality class of the model, clarifying the dependence of the critical creation rate on diffusion.

Findings

Critical diffusion threshold D* in 3D is approximately 0.333.

In 2D, lambda_c decreases exponentially as D approaches 1.

The model belongs to the directed percolation universality class.

Abstract

We study the critical behavior of the pair annihilation model (PAM) with diffusion in one, two and three dimensions, using the pair approximation (PA) and Monte Carlo simulation. Of principal interest is the dependence of the critical creation rate, lambda_c, on the diffusion probability D, in particular, whether survival is possible at arbitrarily small creation rates, for sufficiently rapid diffusion. Whilst the PA predicts that in any spatial dimension d \geq 1, lambda_c -> 0 at some diffusion probability D* < 1, Katori and Konno [Physica A {\bf 186}, 578 (1992)] showed rigorously that for d \leq 2, one has lambda_c > 0 for any D<1. Our simulation results are consistent with this theorem. In two dimensions, the extinction region becomes narrow as D approaches unity, following lambda_c \propto \exp[- \rm{const.}/(1-D)^\gamma], with gamma = 1.41(2). In three dimensions we find D* = 0.333(3). The simulation results confirm that the PAM belongs to the directed percolation universality class.