Structure-preserving numerical integrators for Hodgkin-Huxley-type systems

DOI10.1137/18M123390XMaRDI QIDQ5216792zbMATH OpenOpenAlexFDO

Authors Zhengdao Chen, Baranidharan Raman, Ari Stern

Publication date 20 February 2020

Published in SIAM Journal on Scientific Computing (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1811.00173

zbMATH Keywords

limit cyclesexponential integratorssplitting methodsHodgkin-Huxleycomputational neuroscienceVan der Pol oscillator

Mathematics Subject Classification ID

Neural biology (92C20) Computational methods for bifurcation problems in dynamical systems (37M20) Numerical bifurcation problems (65P30)

Abstract: Motivated by the Hodgkin-Huxley model of neuronal dynamics, we study explicit numerical integrators for "conditionally linear" systems of ordinary differential equations. We show that splitting and composition methods, when applied to the Van der Pol oscillator and to the Hodgkin-Huxley model, do a better job of preserving limit cycles of these systems for large time steps, compared with the "Euler-type" methods (including Euler's method, exponential Euler, and semi-implicit Euler) commonly used in computational neuroscience, with no increase in computational cost. These limit cycles are important to preserve, due to their role in neuronal spiking. Splitting methods even compare favorably to the explicit exponential midpoint method, which is twice as expensive per step. The second-order Strang splitting method is seen to perform especially well across a range of non-stiff and stiff dynamics.

Recommendations

• A comparison of six numerical methods for integrating a compartmental Hodgkin-Huxley type model
• Exponential time differencing for Hodgkin-Huxley-like ODEs
• Splitting methods
• Error analysis of a specialized numerical method for mathematical models from neuroscience
• Positive and bounded nonstandard finite difference scheme for the Hodgkin-Huxley equations

Cites work

• scientific article; zbMATH DE number 725320 (Why is no real title available?)
• Energy behaviour of the Boris method for charged-particle dynamics
• Exponential integrators
• Exponential time differencing for Hodgkin-Huxley-like ODEs
• Geometric Numerical Integration
• Implicit-explicit variational integration of highly oscillatory problems
• Invariant tori of dissipatively perturbed Hamiltonian systems under symplectic discretization
• Modified trigonometric integrators
• On the Construction and Comparison of Difference Schemes
• On the Numerical Integration of Ordinary Differential Equations by Symmetric Composition Methods
• On the Product of Semi-Groups of Operators
• On the qualitative behaviour of symplectic integrators. III: Perturbed integrable systems
• Recent progress in the theory and application of symplectic integrators
• Splitting methods
• Two and three dimensional reductions of the Hodgkin-Huxley system: Separation of time scales and bifurcation schemes

Cited in

(5)

• Exponential time differencing for Hodgkin-Huxley-like ODEs
• scientific article; zbMATH DE number 5775465 (Why is no real title available?)
• A comparison of six numerical methods for integrating a compartmental Hodgkin-Huxley type model
• Probabilistic solvers enable a straight-forward exploration of numerical uncertainty in neuroscience models
• A splitting method for SDEs with locally Lipschitz drift: illustration on the FitzHugh-Nagumo model

This page was built for publication: Structure-preserving numerical integrators for Hodgkin-Huxley-type systems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5216792)