The steady-state bimolecular annihilation reaction A + B -*· 0 on two-dimensional surfaces is studied via computer simulations.
In the simulations A and B are randomly adsorbed on vacant sites, and reaction takes place whenever A and B reach
nearest-neighbor sites, either directly following adsorption or through diffusion. It is found that both with and without diffusion
the reactants segregate into separate islands of As and Bs. The islands vary in size and exhibit highly ramified shapes.
Moreover, the islands are self-similar with a fractal dimension D =
1.89 (similar to percolation, but also other clusters).
D is found to be independent of the diffusion rate K. Other fractal dimensions, e.g., of the hull" differ from those of percolating
clusters. The steady-state coverage 0* =
0*A + 0*B decreases with K, as expected (0*A =
0*B, corresponding to equal fluxes
of A and B is the only physical solution). For systems with immobile particles (K =
0) we find 0* s 0.59 and 0* s 0.49
for the square and the triangular lattices, respectively, similar to the percolation thresholds on these lattices. The long-time
characteristics of the system (D, 0*. etc.) are independent of the initial conditions of the simulation, indicating that the system
reaches a stable steady state. Furthermore, for the large systems simulated (typically 500 X 500 lattice sites) it was found
that the long-time behavior is independent of the input mode. Namely, the same results are obtained for conserved (i.e.,
exactly balanced) and nonconserved (statistically balanced) A,B input mechanisms, indicating that on the time scale of the
simulations (104 Monte Carlo steps) the apparent steady state (for nonconserved input) is essentially identical with the
true steady state (for the conserved input).