Trapping with biased diffusion species

We analyze a trapping reaction with a single penetrable trap, in a one dimensional lattice, where both species (particles and trap) are mobile and have a drift velocity. We obtain the density as seen from a reference system attached to the trap and from the laboratory frame. In addition we study the nearest neighbor distance to the trap. We exploit a stochastic model previously developed, and compare the results with numerical simulations, resulting in an excellent agreement.Comment: 6 pages, 7 Postscript figure

Decay Process for Three - Species Reaction - Diffusion System

We propose the deterministic rate equation of three-species in the reaction - diffusion system. For this case, our purpose is to carry out the decay process in our three-species reaction-diffusion model of the form $A+B+C\to D$. The particle density and the global reaction rate are also shown analytically and numerically on a two-dimensional square lattice with the periodic boundary conditions. Especially, the crossover of the global reaction rate is discussed in both early-time and long-time regimes.Comment: 6 pages, 3 figures, Late

Finite-Size Scaling Studies of Reaction-Diffusion Systems Part III: Numerical Methods

The scaling exponent and scaling function for the 1D single species coagulation model $(A+A\rightarrow A)$ are shown to be universal, i.e. they are not influenced by the value of the coagulation rate. They are independent of the initial conditions as well. Two different numerical methods are used to compute the scaling properties: Monte Carlo simulations and extrapolations of exact finite lattice data. These methods are tested in a case where analytical results are available. It is shown that Monte Carlo simulations can be used to compute even the correction terms. To obtain reliable results from finite-size extrapolations exact numerical data for lattices up to ten sites are sufficient.Comment: 19 pages, LaTeX, 5 figures uuencoded, BONN HE-94-0

Reaction Front in an A+B -> C Reaction-Subdiffusion Process

We study the reaction front for the process A+B -> C in which the reagents move subdiffusively. Our theoretical description is based on a fractional reaction-subdiffusion equation in which both the motion and the reaction terms are affected by the subdiffusive character of the process. We design numerical simulations to check our theoretical results, describing the simulations in some detail because the rules necessarily differ in important respects from those used in diffusive processes. Comparisons between theory and simulations are on the whole favorable, with the most difficult quantities to capture being those that involve very small numbers of particles. In particular, we analyze the total number of product particles, the width of the depletion zone, the production profile of product and its width, as well as the reactant concentrations at the center of the reaction zone, all as a function of time. We also analyze the shape of the product profile as a function of time, in particular its unusual behavior at the center of the reaction zone

Localization-delocalization transition of a reaction-diffusion front near a semipermeable wall

The A+B --> C reaction-diffusion process is studied in a system where the reagents are separated by a semipermeable wall. We use reaction-diffusion equations to describe the process and to derive a scaling description for the long-time behavior of the reaction front. Furthermore, we show that a critical localization-delocalization transition takes place as a control parameter which depends on the initial densities and on the diffusion constants is varied. The transition is between a reaction front of finite width that is localized at the wall and a front which is detached and moves away from the wall. At the critical point, the reaction front remains at the wall but its width diverges with time [as t^(1/6) in mean-field approximation].Comment: 7 pages, PS fil

Behavior of the reaction front between initially segregated species in a two-stage reaction

The large-time asymptotic behavior of a two-stage reaction (A+B→R, B+R→S) with initially segregated reactants is described. The concentration of the reactants is found to be significantly less than the initial concentrations in a depletion zone of width proportional to t[sup. 1/2], where t is time; the reaction takes place in a thinner zone of width proportional to t[sup. 1/6]. Similarity solutions for the chemical concentration profiles in the reaction zone are calculated, and are compared with numerical simulations of the full partial differential reaction-diffusion equations. The large-time asymptotic scalings reported here are the same as in the absence of the secondary reaction, but we find that the location of the reaction zone is significantly shifted due to the secondary reaction. The reaction zone may behave in an exotic fashion at large time, moving first one way, then reversing its direction.Stephen M. Cox and Matthew D. Fin