Abstract:Space is a very important aspect in the simulation of biochemical systems; recently, the need for simulation algorithms able to cope with space is becoming more and more compelling. Complex and detailed models of biochemical systems need to deal with the movement of single molecules and particles, taking into consideration localized fluctuations, transportation phenomena, and diffusion. A common drawback of spatial models lies in their complexity: models can become very large, and their simulation could be time consuming, especially if we want to capture the systems behavior in a reliable way using stochastic methods in conjunction with a high spatial resolution. In order to deliver the promise done by systems biology to be able to understand a system as whole, we need to scale up the size of models we are able to simulate, moving from sequential to parallel simulation algorithms. In this paper, we analyze Smoldyn, a widely diffused algorithm for stochastic simulation of chemical reactions with spatial resolution and single molecule detail, and we propose an alternative, innovative implementation that exploits the parallelism of Graphics Processing Units (GPUs). The implementation executes the most computational demanding steps (computation of diffusion, unimolecular, and bimolecular reaction, as well as the most common cases of molecule-surface interaction) on the GPU, computing them in parallel on each molecule of the system. The implementation offers good speed-ups and real time, high quality graphics output.

Boltzmann Generators and the New Frontier of Computational Sampling in Many-Body Systems

A Boltzmann generator for the isobaric-isothermal ensemble

Transferable Boltzmann Generators

Efficient and Unbiased Sampling of Boltzmann Distributions via Consistency Models

Flow Perturbation to Accelerate Unbiased Sampling of Boltzmann distribution

Active learning of Boltzmann samplers and potential energies with quantum mechanical accuracy

Scalable Lattice Sampling using Factorized Generative Models

Stein Boltzmann Sampling: A Variational Approach for Global Optimization

Stable generative modeling using Schrödinger bridges

Transition Path Sampling with Boltzmann Generator-based MCMC Moves

Scalable Normalizing Flows Enable Boltzmann Generators for Macromolecules

Thermodynamic Interpolation: A generative approach to molecular thermodynamics and kinetics

Markov-chain sampling for long-range systems without evaluating the energy

Energy based diffusion generator for efficient sampling of Boltzmann distributions

Generative methods for sampling transition paths in molecular dynamics

Smoldyn on Graphics Processing Units: Massively Parallel Brownian Dynamics Simulations

Fast and Functional Structured Data Generators Rooted in Out-of-Equilibrium Physics

Enhanced Sampling of Configuration and Path Space in a Generalized Ensemble by Shooting Point Exchange

How to generate an object under an ordinary Boltzmann distribution via an exponential Boltzmann sampler