# 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] Cecilia Cavaterra, Denis Enăchescu, Gabriela Marinoschi. Sliding mode control of the Hodgkin–Huxley mathematical model. Evolution Equations & Control Theory, 2019, 8 (4) : 883-902. doi: 10.3934/eect.2019043 [2] Qingyun Wang, Xia Shi, Guanrong Chen. Delay-induced synchronization transition in small-world Hodgkin-Huxley neuronal networks with channel blocking. Discrete & Continuous Dynamical Systems - B, 2011, 16 (2) : 607-621. doi: 10.3934/dcdsb.2011.16.607 [3] Jian Wu, Jiansheng Geng. Almost periodic solutions for a class of semilinear quantum harmonic oscillators. Discrete & Continuous Dynamical Systems - A, 2011, 31 (3) : 997-1015. doi: 10.3934/dcds.2011.31.997 [4] Dag Lukkassen, Annette Meidell, Peter Wall. On the conjugate of periodic piecewise harmonic functions. Networks & Heterogeneous Media, 2008, 3 (3) : 633-646. doi: 10.3934/nhm.2008.3.633 [5] Daniel Vasiliu, Jianjun Paul Tian. Periodic solutions of a model for tumor virotherapy. Discrete & Continuous Dynamical Systems - S, 2011, 4 (6) : 1587-1597. doi: 10.3934/dcdss.2011.4.1587 [6] Arturo Hidalgo, Lourdes Tello. On a climatological energy balance model with continents distribution. Discrete & Continuous Dynamical Systems - A, 2015, 35 (4) : 1503-1519. doi: 10.3934/dcds.2015.35.1503 [7] Graziano Crasta, Benedetto Piccoli. Viscosity solutions and uniqueness for systems of inhomogeneous balance laws. Discrete & Continuous Dynamical Systems - A, 1997, 3 (4) : 477-502. doi: 10.3934/dcds.1997.3.477 [8] Laura Caravenna. Regularity estimates for continuous solutions of α-convex balance laws. Communications on Pure & Applied Analysis, 2017, 16 (2) : 629-644. doi: 10.3934/cpaa.2017031 [9] Piotr Gwiazda, Piotr Orlinski, Agnieszka Ulikowska. Finite range method of approximation for balance laws in measure spaces. Kinetic & Related Models, 2017, 10 (3) : 669-688. doi: 10.3934/krm.2017027 [10] Dinh-Liem Nguyen. The factorization method for the Drude-Born-Fedorov model for periodic chiral structures. Inverse Problems & Imaging, 2016, 10 (2) : 519-547. doi: 10.3934/ipi.2016010 [11] Ernest Fontich, Rafael de la Llave, Yannick Sire. A method for the study of whiskered quasi-periodic and almost-periodic solutions in finite and infinite dimensional Hamiltonian systems. Electronic Research Announcements, 2009, 16: 9-22. doi: 10.3934/era.2009.16.9 [12] Pierluigi Colli, Antonio Segatti. Uniform attractors for a phase transition model coupling momentum balance and phase dynamics. Discrete & Continuous Dynamical Systems - A, 2008, 22 (4) : 909-932. doi: 10.3934/dcds.2008.22.909 [13] James Walsh, Christopher Rackauckas. On the Budyko-Sellers energy balance climate model with ice line coupling. Discrete & Continuous Dynamical Systems - B, 2015, 20 (7) : 2187-2216. doi: 10.3934/dcdsb.2015.20.2187 [14] Nataliia V. Gorban, Olha V. Khomenko, Liliia S. Paliichuk, Alla M. Tkachuk. Long-time behavior of state functions for climate energy balance model. Discrete & Continuous Dynamical Systems - B, 2017, 22 (5) : 1887-1897. doi: 10.3934/dcdsb.2017112 [15] Daniela Calvetti, Jenni Heino, Erkki Somersalo, Knarik Tunyan. Bayesian stationary state flux balance analysis for a skeletal muscle metabolic model. Inverse Problems & Imaging, 2007, 1 (2) : 247-263. doi: 10.3934/ipi.2007.1.247 [16] Shangbing Ai. Multiple positive periodic solutions for a delay host macroparasite model. Communications on Pure & Applied Analysis, 2004, 3 (2) : 175-182. doi: 10.3934/cpaa.2004.3.175 [17] Rodrigo Donizete Euzébio, Jaume Llibre. Periodic solutions of El Niño model through the Vallis differential system. Discrete & Continuous Dynamical Systems - A, 2014, 34 (9) : 3455-3469. doi: 10.3934/dcds.2014.34.3455 [18] Klemens Fellner, Wolfang Prager, Bao Q. Tang. The entropy method for reaction-diffusion systems without detailed balance: First order chemical reaction networks. Kinetic & Related Models, 2017, 10 (4) : 1055-1087. doi: 10.3934/krm.2017042 [19] Jiawei Dou, Lan-sun Chen, Kaitai Li. A monotone-iterative method for finding periodic solutions of an impulsive competition system on tumor-normal cell interaction. Discrete & Continuous Dynamical Systems - B, 2004, 4 (3) : 555-562. doi: 10.3934/dcdsb.2004.4.555 [20] Marco Spadini. Branches of harmonic solutions to periodically perturbed coupled differential equations on manifolds. Discrete & Continuous Dynamical Systems - A, 2006, 15 (3) : 951-964. doi: 10.3934/dcds.2006.15.951

2018 Impact Factor: 1.313