Abstract:Recent advances in computational neuroscience have demonstrated the usefulness and importance of stochastic, spatial reaction-diffusion simulations. However, ever increasing model complexity renders traditional serial solvers, as well as naive parallel implementations, inadequate. This paper introduces a new generation of the STochastic Engine for Pathway Simulation (STEPS) project (http://steps.sourceforge.net/), denominated STEPS 4.0, and its core components which have been designed for improved scalability, performance, and memory efficiency. STEPS 4.0 aims to enable novel scientific studies of macroscopic systems such as whole cells while capturing their nanoscale details. This class of models is out of reach for serial solvers due to the vast quantity of computation in such detailed models, and also out of reach for naive parallel solvers due to the large memory footprint. Based on a distributed mesh solution, we introduce a new parallel stochastic reaction-diffusion solver and a deterministic membrane potential solver in STEPS 4.0. The distributed mesh, together with improved data layout and algorithm designs, significantly reduces the memory footprint of parallel simulations in STEPS 4.0. This enables massively parallel simulations on modern HPC clusters and overcomes the limitations of the previous parallel STEPS implementation. Current and future improvements to the solver are not sustainable without following proper software engineering principles. For this reason, we also give an overview of how the STEPS codebase and the development environment have been updated to follow modern software development practices. We benchmark performance improvement and memory footprint on three published models with different complexities, from a simple spatial stochastic reaction-diffusion model, to a more complex one that is coupled to a deterministic membrane potential solver to simulate the calcium burst activity of a Purkinje neuron. Simulation results of these models suggest that the new solution dramatically reduces the per-core memory consumption by more than a factor of 30, while maintaining similar or better performance and scalability.

On a variable step size modification of Hines' method in computational neuroscience

Design Fast Algorithms For Hodgkin-Huxley Neuronal Networks

Library-based Numerical Reduction of the Hodgkin–Huxley Neuron for Network Simulation

Refactorization of a variable step, unconditionally stable method of Dahlquist, Liniger and Nevanlinna

Multiscale modeling via split-step methods in neural firing

Stability Analysis of Time Varying Delayed Stochastic Hopfield Neural Networks in Numerical Simulation.

Exponential Time Differencing Algorithm for Pulse-Coupled Hodgkin-Huxley Neural Networks.

The Seven-League Scheme: Deep Learning for Large Time Step Monte Carlo Simulations of Stochastic Differential Equations

On the convergence of adaptive approximations for stochastic differential equations

Library-based Fast Algorithm for Simulating the Hodgkin-Huxley Neuronal Networks

STEPS 4.0: Fast and memory-efficient molecular simulations of neurons at the nanoscale

A Self-Organizing State-Space-Model Approach for Parameter Estimation in Hodgkin-Huxley-Type Models of Single Neurons

Asymptotic Behaviour of Time Stepping Methods for Phase Field Models

Euler Methods with Different Step-Sizes for Different Equations for Stochastic Differential Equations

Adaptive stepsize algorithms for Langevin dynamics

Adaptive Step Size for Hybrid Monte Carlo Algorithm

Mean Square Stability of Two Classes of Theta Methods for Numerical Computation and Simulation of Delayed Stochastic Hopfield Neural Networks.

Efficient optimization of ODE neuron models using gradient descent

Adaptive step-size selection for state-space probabilistic differential equation solvers

Laws of Large Numbers and Langevin Approximations for Stochastic Neural Field Equations

An Antithetic Multilevel Monte Carlo-Milstein Scheme for Stochastic Partial Differential Equations