Regensburg 2013 – scientific programme

Parts | Days | Selection | Search | Updates | Downloads | Help

DY: Fachverband Dynamik und Statistische Physik

DY 5: Reaction-Diffusion Systems

DY 5.3: Talk

Monday, March 11, 2013, 15:30–15:45, H48

Efficient kinetic Monte Carlo method for reaction-diffusion problems with spatially varying annihilation rates — •Karsten Schwarz and Heiko Rieger — Theoretische Physik, Universität des Saarlandes, Saarbrücken

We present an efficient Monte Carlo method to simulate reaction-diffusion processes with spatially varying particle annihilation or transformation rates as it occurs for instance in the context of motor-driven intracellular transport. Like Green's function reaction dynamics and first-passage time methods, our algorithm avoids small diffusive hops by propagating sufficiently distant particles in large hops to the boundaries of protective domains. Since for spatially varying annihilation or transformation rates the single particle diffusion propagator is not known analytically, we present an algorithm that generates efficiently either particle displacements or annihilations with the correct statistics, as we prove rigorously. The numerical efficiency of the algorithm is demonstrated with an illustrative example.