Efficient and accurate time-stepping schemes for integrate-and-fire neuronal networks

Michael Shelley, Louis Tao

Research output: Contribution to journalArticle

Abstract

To avoid the numerical errors associated with resetting the potential following a spike in simulations of integrate-and-fire neuronal networks, Hansel et al. and Shelley independently developed a modified time-stepping method. Their particular scheme consists of second-order Runge-Kutta time-stepping, a linear interpolant to find spike times, and a recalibration of postspike potential using the spike times. Here we show analytically that such a scheme is second order, discuss the conditions under which efficient, higher-order algorithms can be constructed to treat resets, and develop a modified fourth-order scheme. To support our analysis, we simulate a system of integrate-and-fire conductance-based point neurons with all-to-all coupling. For six-digit accuracy, our modified Runge-Kutta fourth-order scheme needs a time-step of Δt = 0.5 × 10-3 seconds, whereas to achieve comparable accuracy using a recalibrated second-order or a first-order algorithm requires time-steps of 10-5 seconds or 10-9 seconds, respectively. Furthermore, since the cortico-cortical conductances in standard integrate-and-fire neuronal networks do not depend on the value of the membrane potential, we can attain fourth-order accuracy with computational costs normally associated with second-order schemes.

Original languageEnglish (US)
Pages (from-to)111-119
Number of pages9
JournalJournal of Computational Neuroscience
Volume11
Issue number2
DOIs
StatePublished - 2001

Fingerprint

Membrane Potentials
Action Potentials
Neurons
Costs and Cost Analysis

Keywords

  • Accurate time integration schemes
  • Integrate-and-fire networks
  • Numerical methods

ASJC Scopus subject areas

  • Neuroscience(all)

Cite this

Efficient and accurate time-stepping schemes for integrate-and-fire neuronal networks. / Shelley, Michael; Tao, Louis.

In: Journal of Computational Neuroscience, Vol. 11, No. 2, 2001, p. 111-119.

Research output: Contribution to journalArticle

@article{432441e7982e4e12b1d7489727d0f15b,
title = "Efficient and accurate time-stepping schemes for integrate-and-fire neuronal networks",
abstract = "To avoid the numerical errors associated with resetting the potential following a spike in simulations of integrate-and-fire neuronal networks, Hansel et al. and Shelley independently developed a modified time-stepping method. Their particular scheme consists of second-order Runge-Kutta time-stepping, a linear interpolant to find spike times, and a recalibration of postspike potential using the spike times. Here we show analytically that such a scheme is second order, discuss the conditions under which efficient, higher-order algorithms can be constructed to treat resets, and develop a modified fourth-order scheme. To support our analysis, we simulate a system of integrate-and-fire conductance-based point neurons with all-to-all coupling. For six-digit accuracy, our modified Runge-Kutta fourth-order scheme needs a time-step of Δt = 0.5 × 10-3 seconds, whereas to achieve comparable accuracy using a recalibrated second-order or a first-order algorithm requires time-steps of 10-5 seconds or 10-9 seconds, respectively. Furthermore, since the cortico-cortical conductances in standard integrate-and-fire neuronal networks do not depend on the value of the membrane potential, we can attain fourth-order accuracy with computational costs normally associated with second-order schemes.",
keywords = "Accurate time integration schemes, Integrate-and-fire networks, Numerical methods",
author = "Michael Shelley and Louis Tao",
year = "2001",
doi = "10.1023/A:1012885314187",
language = "English (US)",
volume = "11",
pages = "111--119",
journal = "Journal of Computational Neuroscience",
issn = "0929-5313",
publisher = "Springer Netherlands",
number = "2",

}

TY - JOUR

T1 - Efficient and accurate time-stepping schemes for integrate-and-fire neuronal networks

AU - Shelley, Michael

AU - Tao, Louis

PY - 2001

Y1 - 2001

N2 - To avoid the numerical errors associated with resetting the potential following a spike in simulations of integrate-and-fire neuronal networks, Hansel et al. and Shelley independently developed a modified time-stepping method. Their particular scheme consists of second-order Runge-Kutta time-stepping, a linear interpolant to find spike times, and a recalibration of postspike potential using the spike times. Here we show analytically that such a scheme is second order, discuss the conditions under which efficient, higher-order algorithms can be constructed to treat resets, and develop a modified fourth-order scheme. To support our analysis, we simulate a system of integrate-and-fire conductance-based point neurons with all-to-all coupling. For six-digit accuracy, our modified Runge-Kutta fourth-order scheme needs a time-step of Δt = 0.5 × 10-3 seconds, whereas to achieve comparable accuracy using a recalibrated second-order or a first-order algorithm requires time-steps of 10-5 seconds or 10-9 seconds, respectively. Furthermore, since the cortico-cortical conductances in standard integrate-and-fire neuronal networks do not depend on the value of the membrane potential, we can attain fourth-order accuracy with computational costs normally associated with second-order schemes.

AB - To avoid the numerical errors associated with resetting the potential following a spike in simulations of integrate-and-fire neuronal networks, Hansel et al. and Shelley independently developed a modified time-stepping method. Their particular scheme consists of second-order Runge-Kutta time-stepping, a linear interpolant to find spike times, and a recalibration of postspike potential using the spike times. Here we show analytically that such a scheme is second order, discuss the conditions under which efficient, higher-order algorithms can be constructed to treat resets, and develop a modified fourth-order scheme. To support our analysis, we simulate a system of integrate-and-fire conductance-based point neurons with all-to-all coupling. For six-digit accuracy, our modified Runge-Kutta fourth-order scheme needs a time-step of Δt = 0.5 × 10-3 seconds, whereas to achieve comparable accuracy using a recalibrated second-order or a first-order algorithm requires time-steps of 10-5 seconds or 10-9 seconds, respectively. Furthermore, since the cortico-cortical conductances in standard integrate-and-fire neuronal networks do not depend on the value of the membrane potential, we can attain fourth-order accuracy with computational costs normally associated with second-order schemes.

KW - Accurate time integration schemes

KW - Integrate-and-fire networks

KW - Numerical methods

UR - http://www.scopus.com/inward/record.url?scp=0035683779&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0035683779&partnerID=8YFLogxK

U2 - 10.1023/A:1012885314187

DO - 10.1023/A:1012885314187

M3 - Article

C2 - 11717528

AN - SCOPUS:0035683779

VL - 11

SP - 111

EP - 119

JO - Journal of Computational Neuroscience

JF - Journal of Computational Neuroscience

SN - 0929-5313

IS - 2

ER -