High-Order Unconditionally Stable IMEX BDF--FEM Schemes for the Hodgkin--Huxley Model on a Branched Axon
Abstract
We construct and analyze high-order implicit--explicit backward differentiation formula finite element (IMEX BDF--FEM) schemes for the Hodgkin--Huxley cable system on a branched metric graph. The linearly implicit schemes combine constant-step BDF methods of orders q=1,…,5 with conforming finite elements of any fixed degree p≥1. Extrapolated gate variables are smoothly masked before evaluating the conductances, preserving coercivity without projecting the numerical gates. A voltage cutoff makes the gate rates globally bounded and Lipschitz. Each timestep requires one linear voltage solve followed by independent nodal gate updates. Under exact spatial integration, we prove unconditional stability of the complete fully discrete regularized scheme in discrete history energies. Under regularity and starting-history assumptions, the voltage and reconstructed gates converge in L2(Γ) at the optimal rate O(hp+1+Δtq), without a timestep restriction or an h--Δt coupling. The same rate holds with numerical quadrature when its full residual defect is O(hp+1). Manufactured solutions recover the predicted orders through BDF5--P4. Action-potential simulations illustrate pulse propagation and splitting, while a controlled work--precision comparison within the proposed family shows that higher-order variants can reach a prescribed accuracy at much lower runtime along the tested joint-refinement paths.
Related
- Global Dynamics and Synchronization of Hodgkin-Huxley-Wilson Neural Networks
- An unconditionally stable second-order Robin partitioned method for fluid-poroelastic structure interaction
- Unconditionally L2-stable explicit--implicit--null filtered Runge--Kutta methods for general nonlinear degenerate parabolic equations
- A linear, second-order, unconditionally energy-stable scalar-auxiliary-variable scheme for the CHQ viscoelastic phase-separation model
- Stable rational cohomology splitting of moduli of branched covers of curves
BibTeX
@misc{reed2026highorder,
title = {{High-Order Unconditionally Stable IMEX BDF--FEM Schemes for the Hodgkin--Huxley Model on a Branched Axon}},
author = {Hayden Reed and Qingguang Guan},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.02423},
eprint = {2611.02423},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.02423},
}
RIS
TY - UNPB
TI - High-Order Unconditionally Stable IMEX BDF--FEM Schemes for the Hodgkin--Huxley Model on a Branched Axon
AU - Hayden Reed
AU - Qingguang Guan
PY - 2026
DA - 2026/11/02
PB - arXiv
DO - 10.48550/arxiv.2611.02423
UR - https://arxiv.org/abs/2611.02423
AB - We construct and analyze high-order implicit--explicit backward differentiation formula finite element (IMEX BDF--FEM) schemes for the Hodgkin--Huxley cable system on a branched metric graph. The linearly implicit schemes combine constant-step BDF methods of orders q=1,…,5 with conforming finite elements of any fixed degree p≥1. Extrapolated gate variables are smoothly masked before evaluating the conductances, preserving coercivity without projecting the numerical gates. A voltage cutoff makes the gate rates globally bounded and Lipschitz. Each timestep requires one linear voltage solve followed by independent nodal gate updates. Under exact spatial integration, we prove unconditional stability of the complete fully discrete regularized scheme in discrete history energies. Under regularity and starting-history assumptions, the voltage and reconstructed gates converge in L^2(Γ) at the optimal rate O(h^p+1+Δt^q), without a timestep restriction or an h--Δt coupling. The same rate holds with numerical quadrature when its full residual defect is O(h^p+1). Manufactured solutions recover the predicted orders through BDF5--P_4. Action-potential simulations illustrate pulse propagation and splitting, while a controlled work--precision comparison within the proposed family shows that higher-order variants can reach a prescribed accuracy at much lower runtime along the tested joint-refinement paths.
ER -