TY - JOUR
T1 - Analytic modeling of neural tissue
T2 - II. Nonlinear membrane dynamics
AU - Schwartz, B. L.
AU - Brown, S. M.
AU - Muthuswamy, J.
AU - Sadleir, R. J.
N1 - Publisher Copyright:
© 2022 Author(s).
PY - 2022/11/1
Y1 - 2022/11/1
N2 - Computational modeling of neuroactivity plays a central role in our effort to understand brain dynamics in the advancements of neural engineering such as deep brain stimulation, neuroprosthetics, and magnetic resonance electrical impedance tomography. However, analytic solutions do not capture the fundamental nonlinear behavior of an action potential. What is needed is a method that is not constrained to only linearized models of neural tissue. Therefore, the objective of this study is to establish a robust, straightforward process for modeling neurodynamic phenomena, which preserves their nonlinear features. To address this, we turn to decomposition methods from homotopy analysis, which have emerged in recent decades as powerful tools for solving nonlinear differential equations. We solve the nonlinear ordinary differential equations of three landmark models of neural conduction - Ermentrout-Kopell, FitzHugh-Nagumo, and Hindmarsh-Rose models - using George Adomian's decomposition method. For each variable, we construct a power series solution equivalent to a generalized Taylor series expanded about a function. The first term of the decomposition series comes from the models' initial conditions. All subsequent terms are recursively determined from the first. We show rapid convergence, achieving a maximal error of <10-12 with only eight terms. We extend the region of convergence with one-step analytic continuation so that our complete solutions are decomposition splines. We show that this process can yield solutions for single- and multi-variable models and can characterize a single action potential or complex bursting patterns. Finally, we show that the accuracy of this decomposition approach favorably compares to an established polynomial method, B-spline collocation. The strength of this method, besides its stability and ease of computation, is that, unlike perturbation, we make no changes to the models' equations; thus, our solutions are to the problems at hand, not simplified versions. This work validates decomposition as a viable technique for advanced neural engineering studies.
AB - Computational modeling of neuroactivity plays a central role in our effort to understand brain dynamics in the advancements of neural engineering such as deep brain stimulation, neuroprosthetics, and magnetic resonance electrical impedance tomography. However, analytic solutions do not capture the fundamental nonlinear behavior of an action potential. What is needed is a method that is not constrained to only linearized models of neural tissue. Therefore, the objective of this study is to establish a robust, straightforward process for modeling neurodynamic phenomena, which preserves their nonlinear features. To address this, we turn to decomposition methods from homotopy analysis, which have emerged in recent decades as powerful tools for solving nonlinear differential equations. We solve the nonlinear ordinary differential equations of three landmark models of neural conduction - Ermentrout-Kopell, FitzHugh-Nagumo, and Hindmarsh-Rose models - using George Adomian's decomposition method. For each variable, we construct a power series solution equivalent to a generalized Taylor series expanded about a function. The first term of the decomposition series comes from the models' initial conditions. All subsequent terms are recursively determined from the first. We show rapid convergence, achieving a maximal error of <10-12 with only eight terms. We extend the region of convergence with one-step analytic continuation so that our complete solutions are decomposition splines. We show that this process can yield solutions for single- and multi-variable models and can characterize a single action potential or complex bursting patterns. Finally, we show that the accuracy of this decomposition approach favorably compares to an established polynomial method, B-spline collocation. The strength of this method, besides its stability and ease of computation, is that, unlike perturbation, we make no changes to the models' equations; thus, our solutions are to the problems at hand, not simplified versions. This work validates decomposition as a viable technique for advanced neural engineering studies.
UR - http://www.scopus.com/inward/record.url?scp=85144821712&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85144821712&partnerID=8YFLogxK
U2 - 10.1063/5.0124414
DO - 10.1063/5.0124414
M3 - Article
AN - SCOPUS:85144821712
SN - 2158-3226
VL - 12
JO - AIP Advances
JF - AIP Advances
IS - 11
M1 - 115019
ER -