University of Cambridge > Talks.cam > Isaac Newton Institute Seminar Series > Fast Reactive Brownian Dynamics

Fast Reactive Brownian Dynamics

Add to your list(s) Download to your calendar using vCal

If you have a question about this talk, please contact INI IT.

SDBW04 - Spatially distributed stochastic dynamical systems in biology

I will describe a particle-based algorithm for reaction-diffusion problems that combines Brownian dynamics with a Markov reaction process. The microscopic model simulated by our Split Reactive Brownian Dynamics (SRBD) algorithm is based on the Doi or volume-reactivity model. This model applies only to reactions with at most two reactants, which is physically realistic. Let us consider the simple reaction A+B->product. In the Doi model, particles are independent spherical Brownian walkers (this can be relaxed to account for hydrodynamic interactions), and while an A and a B particle overlap, there is a Poisson process with a given microscopic reaction rate for the two particles to react and give a product. Our goal is to simulate this complex Markov process in dense systems of many particles, in the presence of multiple reaction channels.

Our algorithm is inspired by the Isotropic Direct Simulation Monte Carlo (I-DSMC) method and the next subvolume method. Strang splitting is used to separate diffusion and reaction; this is the only approximation made in our method. In order to process reactions without approximations, with the particles frozen in place, we use an event-driven algorithm. We divide the system into a grid of cells such that only particles in neighboring cells can react. Each cell schedules the next potential reaction to happen involving a particle in that cell and a particle in one of the neighboring cells, and an event queue is used to select the next cell in which a reaction may happen.

Note that, while a grid of cells is used to make the algorithm efficient, the results obtained by the SRBD method are grid-independent and thus free of grid artifacts, such as loss of Galilean invariance and sensitivity of the results to the grid spacing. I will compare our SRBD method with grid-based methods, such as©RDME and a variant of RDME that we call Split Brownian Dynamics with Reaction Master Equation (S-BD-RME), on a problem involving the spontaneous f

This talk is part of the Isaac Newton Institute Seminar Series series.

Tell a friend about this talk:

This talk is included in these lists:

Note that ex-directory lists are not shown.

 

© 2006-2021 Talks.cam, University of Cambridge. Contact Us | Help and Documentation | Privacy and Publicity