Simulations for trapping reactions with subdiffusive traps and subdiffusive particles

While there are many well known and extensively tested results involving diffusion-limited binary reactions, reactions involving subdiffusive reactant species are far less well understood. Subdiffusive motion is characterized by a mean square displacement ∼ tγ with 0<γ<1. Recently, we calculated the asymptotic survival probability P(t) of a (sub)diffusive particle (γ') surrounded by (sub)diffusive traps (γ) in one dimension. These are among the few known results for reactions involving species characterized by different anomalous exponents. Our results were obtained by bounding, above and below, the exact survival probability by two other probabilities that are asymptotically identical (except when γ' = 1 and 0<γ<2/3). Using this approach, we were not able to estimate the time of validity of the asymptotic result, nor the way in which the survival probability approaches this regime. Toward this goal, here we present a detailed comparison of the asymptotic results with numerical simulations. In some parameter ranges the asymptotic theory describes the simulation results very well even for relatively short times. However, in other regimes more time is required for the simulation results to approach asymptotic behaviour, and we arrive at situations where we are not able to reach asymptotes within our computational means. This is regrettably the case for γ' = 1 and 0<γ<2/3, where we are therefore not able to prove or disprove even conjectures about the asymptotic survival probability of the particle