This title appears in the Scientific Report :
2018
Please use the identifier:
http://dx.doi.org/10.1016/j.jpdc.2017.11.019 in citations.
A scalable algorithm for simulating the structural plasticity of the brain
A scalable algorithm for simulating the structural plasticity of the brain
The neural network in the brain is not hard-wired. Even in the mature brain, new connections between neurons are formed and existing ones are deleted, which is called structural plasticity. The dynamics of the connectome is key to understanding how learning, memory, and healing after lesions such as...
Saved in:
Personal Name(s): | Rinke, Sebastian (Corresponding author) |
---|---|
Butz-Ostendorf, Markus / Hermanns, Marc-André / Naveau, Mikael / Wolf, Felix | |
Contributing Institute: |
Jülich Supercomputing Center; JSC JARA - HPC; JARA-HPC |
Published in: | Journal of parallel and distributed computing, 120 (2017) S. 251-266 |
Imprint: |
Amsterdam [u.a.]
Elsevier
2017
|
DOI: |
10.1016/j.jpdc.2017.11.019 |
Document Type: |
Journal Article |
Research Program: |
SimLab Neuroscience Scalable Performance Analysis of Large-Scale Parallel Applications Supercomputing and Modelling for the Human Brain Computational Science and Mathematical Methods |
Publikationsportal JuSER |
The neural network in the brain is not hard-wired. Even in the mature brain, new connections between neurons are formed and existing ones are deleted, which is called structural plasticity. The dynamics of the connectome is key to understanding how learning, memory, and healing after lesions such as stroke work. However, with current experimental techniques even the creation of an exact static connectivity map, which is required for various brain simulations, is very difficult. One alternative is to use network models to simulate the evolution of synapses between neurons based on their specified activity targets. This is particularly useful as experimental measurements of the spiking frequency of neurons are more easily accessible and reliable than biological connectivity data. The Model of Structural Plasticity (MSP) by Butz and van Ooyen is an example of this approach. However, to predict which neurons connect to each other, the current MSP model computes probabilities for all pairs of neurons, resulting in a complexity . To enable large-scale simulations with millions of neurons and beyond, this quadratic term is prohibitive. Inspired by hierarchical methods for solving -body problems in particle physics, we propose a scalable approximation algorithm for MSP that reduces the complexity to without any notable impact on the quality of the results. We show that an MPI-based parallel implementation of our scalable algorithm can simulate the structural plasticity of up to neurons—four orders of magnitude more than the naïve version. |