# American Institute of Mathematical Sciences

June  2018, 15(3): 807-825. doi: 10.3934/mbe.2018036

## A multi-base harmonic balance method applied to Hodgkin-Huxley model

 Normandie Univ, UNIHAVRE, LMAH, FR-CNRS-3335, ISCN, 76600 Le Havre, France

* Corresponding author: V. Lanza.

Received  November 23, 2016 Revised  October 09, 2017 Published  December 2017

Our aim is to propose a new robust and manageable technique, called multi-base harmonic balance method, to detect and characterize the periodic solutions of a nonlinear dynamical system. Our case test is the Hodgkin-Huxley model, one of the most realistic neuronal models in literature. This system, depending on the value of the external stimuli current, exhibits periodic solutions, both stable and unstable.

Citation: Aymen Balti, Valentina Lanza, Moulay Aziz-Alaoui. A multi-base harmonic balance method applied to Hodgkin-Huxley model. Mathematical Biosciences & Engineering, 2018, 15 (3) : 807-825. doi: 10.3934/mbe.2018036
##### References:

show all references

##### References:
Bifurcation diagram of the HH model, showing the stable (solid line) and unstable (dotted line) branches of the periodic solutions of HH model. For each periodic solution the minimum and the maximum values of the potential $V$ over one period are represented. Depending on the values of $I$, two regions with different dynamical behaviors can be identified.
Zoom for $I \in~[0, ~I_{2}]$ of Fig 1. HH model exhibits one equilibrium point, one stable limit cycle (solid line) and up to 3 unstable ones (dotted lines).
(a) The stable periodic solution detected by the HB method for $I =6.25$ exhibits a sort of Gibbs phenomenon. (b) Zoom showing the small oscillations, sign of a non accurate approximation of the limit cycle, despite the exploitation of 50 harmonics.
(a) Stable and (b) unstable periodic solutions for different values of $I$, in a neighborhood of (a) $I_1$ and (b) $I_2$, respectively.
Time series of (a) the stable periodic solution for $I = 152.2500$ and (b) the unstable periodic solution for $I = 9.71889$.
Stable (solid line) and unstable (dashed line) limit cycles near the first saddle-node of cycles bifurcation, for (a) $I =6.2649$ both solutions are almost coincident, and for (b) $I =6.2716$.
Projection of two unstable limit cycles on the $(V, n)$ plane for (a) $I = 7.92198548\lesssim I_3$ and (b) $I = I_3 = 7.92198549$.
Projection of two unstable limit cycles on the $(V, n)$ plane for (a) $I = I_4 = 7.84654752$ and (b) $I = 7.84654876\lesssim I_4$.
(a)-(b) Floquet multipliers for the stable limit cycles and unstable limit cycles, respectively, associated to the first saddle node of cycles bifurcation for $I\in [6.2792, 6.7872]$. As $I$ increases, in (a) the multiplier $\mu_4$ starts from the value +1 and then enters in the unit circle, while in (b) the multiplier $\mu_4$ starts to the value +1 and becomes bigger and bigger. (c)-(d) Floquet multipliers for the two unstable limit cycles associated to the second saddle node of cycles bifurcation for $I\in [7.921985465, 7.921985491]$. Here, in both cases, the third multiplier is outside the unit circle (this makes the limit cycle unstable) and is not shown, since it takes very high values with respect to the others. As in the previous case, as $I$ decreases, the multiplier $\mu_4$ starts from the value +1 and either (c) enters in the unit circle, or (d) takes higher and higher values. (e)-(f) Floquet multipliers for the two unstable limit cycles associated to the third saddle node of cycles bifurcation for $I\in [7.846557778, 7.846616827]$. Also in this case, for both limit cycles, the third multiplier is not represented. As $I$ increases, the multiplier $\mu_4$ starts from the value +1 and either (e) escapes from, or (f) enters in the unit circle.
Floquet multipliers near the period-doubling bifurcation for different values of $I\in [7.92197743, 7.92197799]$. By decreasing $I$, the multiplier $\mu_4$ crosses the unit cycle through $-1$.
By decreasing the value of I, the multipliers $\mu_4$ decreases, crosses the value -1 for $I =7.92197768$ and enters into the unit circle.
 I $\mu_1$ $\mu_2$ $\mu_3$ $\mu_4$ 7.92197799 1.000 0.000 -2940.687 -1.041 7.92197793 1.000 -0.000 -2964.042 -1.033 7.92197787 1.000 0.000 -2987.386 -1.025 7.92197781 1.000 0.000 -3010.719 -1.017 7.92197775 1.000 -0.000 -3034.042 -1.009 7.92197768 1.000 0.000 -3057.354 -1.001 7.92197762 1.000 -0.000 -3080.655 -0.993 7.92197756 1.000 0.000 -3103.946 -0.986 7.92197750 1.000 0.000 -3127.225 -0.978 7.92197743 1.000 0.000 -3150.494 -0.9713
 I $\mu_1$ $\mu_2$ $\mu_3$ $\mu_4$ 7.92197799 1.000 0.000 -2940.687 -1.041 7.92197793 1.000 -0.000 -2964.042 -1.033 7.92197787 1.000 0.000 -2987.386 -1.025 7.92197781 1.000 0.000 -3010.719 -1.017 7.92197775 1.000 -0.000 -3034.042 -1.009 7.92197768 1.000 0.000 -3057.354 -1.001 7.92197762 1.000 -0.000 -3080.655 -0.993 7.92197756 1.000 0.000 -3103.946 -0.986 7.92197750 1.000 0.000 -3127.225 -0.978 7.92197743 1.000 0.000 -3150.494 -0.9713
 [1] Li-Bin Liu, Ying Liang, Jian Zhang, Xiaobing Bao. A robust adaptive grid method for singularly perturbed Burger-Huxley equations. Electronic Research Archive, 2020, 28 (4) : 1439-1457. doi: 10.3934/era.2020076 [2] Hongyu Cheng, Shimin Wang. Response solutions to harmonic oscillators beyond multi–dimensional brjuno frequency. Communications on Pure & Applied Analysis, 2021, 20 (2) : 467-494. doi: 10.3934/cpaa.2020222 [3] Mengyu Cheng, Zhenxin Liu. Periodic, almost periodic and almost automorphic solutions for SPDEs with monotone coefficients. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021026 [4] Zhihua Liu, Yayun Wu, Xiangming Zhang. Existence of periodic wave trains for an age-structured model with diffusion. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2021009 [5] Rong Chen, Shihang Pan, Baoshuai Zhang. Global conservative solutions for a modified periodic coupled Camassa-Holm system. Electronic Research Archive, 2021, 29 (1) : 1691-1708. doi: 10.3934/era.2020087 [6] Dong-Ho Tsai, Chia-Hsing Nien. On space-time periodic solutions of the one-dimensional heat equation. Discrete & Continuous Dynamical Systems - A, 2020, 40 (6) : 3997-4017. doi: 10.3934/dcds.2020037 [7] Yi Guan, Michal Fečkan, Jinrong Wang. Periodic solutions and Hyers-Ulam stability of atmospheric Ekman flows. Discrete & Continuous Dynamical Systems - A, 2021, 41 (3) : 1157-1176. doi: 10.3934/dcds.2020313 [8] Jiangtao Yang. Permanence, extinction and periodic solution of a stochastic single-species model with Lévy noises. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020371 [9] Thazin Aye, Guanyu Shang, Ying Su. On a stage-structured population model in discrete periodic habitat: III. unimodal growth and delay effect. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2021005 [10] Chao Wang, Qihuai Liu, Zhiguo Wang. Periodic bouncing solutions for Hill's type sub-linear oscillators with obstacles. Communications on Pure & Applied Analysis, 2021, 20 (1) : 281-300. doi: 10.3934/cpaa.2020266 [11] Sishu Shankar Muni, Robert I. McLachlan, David J. W. Simpson. Homoclinic tangencies with infinitely many asymptotically stable single-round periodic solutions. Discrete & Continuous Dynamical Systems - A, 2021  doi: 10.3934/dcds.2021010 [12] Michal Fečkan, Kui Liu, JinRong Wang. $(\omega,\mathbb{T})$-periodic solutions of impulsive evolution equations. Evolution Equations & Control Theory, 2021  doi: 10.3934/eect.2021006 [13] Liupeng Wang, Yunqing Huang. Error estimates for second-order SAV finite element method to phase field crystal model. Electronic Research Archive, 2021, 29 (1) : 1735-1752. doi: 10.3934/era.2020089 [14] Martin Kalousek, Joshua Kortum, Anja Schlömerkemper. Mathematical analysis of weak and strong solutions to an evolutionary model for magnetoviscoelasticity. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 17-39. doi: 10.3934/dcdss.2020331 [15] Helmut Abels, Johannes Kampmann. Existence of weak solutions for a sharp interface model for phase separation on biological membranes. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 331-351. doi: 10.3934/dcdss.2020325 [16] Philippe Laurençot, Christoph Walker. Variational solutions to an evolution model for MEMS with heterogeneous dielectric properties. Discrete & Continuous Dynamical Systems - S, 2021, 14 (2) : 677-694. doi: 10.3934/dcdss.2020360 [17] Yi-Long Luo, Yangjun Ma. Low Mach number limit for the compressible inertial Qian-Sheng model of liquid crystals: Convergence for classical solutions. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 921-966. doi: 10.3934/dcds.2020304 [18] Lu Xu, Chunlai Mu, Qiao Xin. Global boundedness of solutions to the two-dimensional forager-exploiter model with logistic source. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020396 [19] Guoyuan Chen, Yong Liu, Juncheng Wei. Nondegeneracy of harmonic maps from ${{\mathbb{R}}^{2}}$ to ${{\mathbb{S}}^{2}}$. Discrete & Continuous Dynamical Systems - A, 2020, 40 (6) : 3215-3233. doi: 10.3934/dcds.2019228 [20] Mehdi Bastani, Davod Khojasteh Salkuyeh. On the GSOR iteration method for image restoration. Numerical Algebra, Control & Optimization, 2021, 11 (1) : 27-43. doi: 10.3934/naco.2020013

2018 Impact Factor: 1.313