“…In the past three decades, numerous simplified neuron models had been fantastically extended from the classical Hodgkin-Huxley model [1] to reconstruct the main dynamical characteristics of neuronal electrical activities [2][3][4][5][6][7][8], among which the two-and three-dimensional Hindmarsh-Rose (HR) neuron models are effective and available for dynamical analysis in electrical activities of biological neurons [9,10]. In the last few years, a wide variety of the HR neuron models, such as original three-dimensional HR models [10][11][12][13][14][15][16], extended or nonlinear feedback coupled HR models [17][18][19][20], time delayed HR models [20][21][22], fractional-order HR models [23,24], and memristor based HR models under electromagnetic radiations [9,[25][26][27], have been proposed and further studied by bifurcation analysis methods for understanding the dynamics of electrical activities among neurons [8]. For this reason, bifurcation analysis theory plays an essential role in describing mode transitions between spiking and bursting in the neuronal electrical activities [9][10][11][12][13][14][15][16][17][18]…”