# A Diffusion Based Embedding of the Stochastic Simulation Algorithm in Continuous Space

@article{Thomas2020ADB, title={A Diffusion Based Embedding of the Stochastic Simulation Algorithm in Continuous Space}, author={Marcus Thomas and Russell Schwartz}, journal={Biophysical Journal}, year={2020}, volume={118} }

A variety of simulation methodologies have been used for modeling reaction-diffusion dynamics -- including approaches based on Differential Equations (DE), the Stochastic Simulation Algorithm (SSA), Brownian Dynamics (BD), Green's Function Reaction Dynamics (GFRD), and variations thereon -- each offering trade-offs with respect to the ranges of phenomena they can model, their computational tractability, and the difficulty of fitting them to experimental measurements. Here, we develop a… Expand

#### Figures and Tables from this paper

#### One Citation

Self-assembly of finite-sized colloidal aggregates.

- Materials Science, Physics
- Soft matter
- 2020

A more versatile design allowing us to fabricate different types of finite-sized aggregates, including colloidal molecules or core-shell clusters, starting from finite density suspensions of isotropic colloids in bulk is validated. Expand

#### References

SHOWING 1-10 OF 57 REFERENCES

Efficient reactive Brownian dynamics.

- Physics, Medicine
- The Journal of chemical physics
- 2018

We develop a Split Reactive Brownian Dynamics (SRBD) algorithm for particle simulations of reaction-diffusion systems based on the Doi or volume reactivity model, in which pairs of particles react… Expand

Fast Monte Carlo Simulation Methods for Biological Reaction-Diffusion Systems in Solution and on Surfaces

- Computer Science, Medicine
- SIAM J. Sci. Comput.
- 2008

It is shown how spatially realistic Monte Carlo simulations of biological systems can be far more cost-effective than often is assumed, and provide a level of accuracy and insight beyond that of continuum methods. Expand

Efficient stochastic simulation of chemical kinetics networks using a weighted ensemble of trajectories.

- Computer Science, Biology
- The Journal of chemical physics
- 2013

The "weighted ensemble" (THE AUTHORS) simulation strategy is applied to a series of systems-biology models that range in complexity from a one-dimensional system to a system with 354 species and 3680 reactions, and results significantly enhance the precision with which rare or slow processes can be characterized. Expand

Perspective: Stochastic algorithms for chemical kinetics.

- Medicine, Chemistry
- The Journal of chemical physics
- 2013

The perspective on stochastic chemical kinetics is outlined, paying particular attention to numerical simulation algorithms, and some of the more promising strategies for dealing stochastically with stiff systems, rare events, and sensitivity analysis are taken. Expand

Green's-function reaction dynamics: a particle-based approach for simulating biochemical networks in time and space.

- Medicine
- The Journal of chemical physics
- 2005

The main idea of GFRD is to exploit the exact solution of the Smoluchoswki equation to set up an event-driven algorithm, which combines in one step the propagation of the particles in space with the reactions between them. Expand

Efficient stochastic sampling of first-passage times with applications to self-assembly simulations.

- Biology, Medicine
- The Journal of chemical physics
- 2008

Two methods for accelerating sampling in SSA models are developed: an exact method and a scheme allowing for sampling accuracy up to any arbitrary error bound, based on the analysis of the eigenvalues of continuous time Markov models that define the behavior of the SSA. Expand

The small-voxel tracking algorithm for simulating chemical reactions among diffusing molecules.

- Chemistry, Medicine
- The Journal of chemical physics
- 2014

A molecule-tracking algorithm that has the same theoretical underpinnings and strategic aims as eGFRD, but a different implementation procedure, called the small-voxel tracking algorithm (SVTA), which combines the well known voxel-hopping method for simulating molecular diffusion with a novel procedure for rectifying the unphysical predictions of the diffusion equation on the small spatiotemporal scale of molecular collisions. Expand

Approximate accelerated stochastic simulation of chemically reacting systems

- Computer Science
- 2001

With further refinement, the τ-leap method should provide a viable way of segueing from the exact SSA to the approximate chemical Langevin equation, and thence to the conventional deterministic reaction rate equation, as the system size becomes larger. Expand

A diffusional bimolecular propensity function.

- Medicine, Chemistry
- The Journal of chemical physics
- 2009

We derive an explicit formula for the propensity function (stochastic reaction rate) of a generic bimolecular chemical reaction in which the reactant molecules move about by diffusion, as solute… Expand

A Brownian dynamics algorithm for colloids in curved manifolds.

- Physics, Medicine
- The Journal of chemical physics
- 2014

The many-particle Langevin equation, written in local coordinates, is used to derive a Brownian dynamics simulation algorithm to study the dynamics of colloids moving on curved manifolds and is applied to strong correlated systems, namely, paramagnetic colloids along a circle and soft colloids on a sphere. Expand