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

DOI10.1137/18M123390XzbMATH Open1431.65238arXiv1811.00173OpenAlexW3100501652MaRDI QIDQ5216792FDOQ5216792

Authors: Zhengdao Chen, Baranidharan Raman, Ari Stern 

Publication date: 20 February 2020

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

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.

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

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

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)

Cites Work

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

Cited In (5)

• Exponential time differencing for Hodgkin-Huxley-like ODEs
• Title not available (Why is that?)
• 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)