A comb-shaped chaotic region has been simulated in multiple two-dimensional parameter spaces using the Hindmarsh–Rose(HR) neuron model in many recent studies, which can interpret almost all of the previously simulated bifurcation processes with chaos in neural firing patterns. In the present paper, a comb-shaped chaotic region in a twodimensional parameter space was reproduced, which presented different processes of period-adding bifurcations with chaos with changing one parameter and fixed the other parameter at different levels. In the biological experiments, different period-adding bifurcation scenarios with chaos by decreasing the extra-cellular calcium concentration were observed from some neural pacemakers at different levels of extra-cellular 4-aminopyridine concentration and from other pacemakers at different levels of extra-cellular caesium concentration. By using the nonlinear time series analysis method, the deterministic dynamics of the experimental chaotic firings were investigated. The period-adding bifurcations with chaos observed in the experiments resembled those simulated in the comb-shaped chaotic region using the HR model. The experimental results show that period-adding bifurcations with chaos are preserved in different two-dimensional parameter spaces, which provides evidence of the existence of the comb-shaped chaotic region and a demonstration of the simulation results in different two-dimensional parameter spaces in the HR neuron model. The results also present relationships between different firing patterns in two-dimensional parameter spaces.
The coexistence of a resting condition and period-1 firing near a subcritical Hopf bifurcation point, lying between the monostable resting condition and period-1 firing, is often observed in neurons of the central nervous systems. Near such a bifurcation point in the Morris–Lecar(ML) model, the attraction domain of the resting condition decreases while that of the coexisting period-1 firing increases as the bifurcation parameter value increases. With the increase of the coupling strength, and parameter and initial value dependent synchronization transition processes from non-synchronization to compete synchronization are simulated in two coupled ML neurons with coexisting behaviors: one neuron chosen as the resting condition and the other the coexisting period-1 firing. The complete synchronization is either a resting condition or period-1 firing dependent on the initial values of period-1 firing when the bifurcation parameter value is small or middle and is period-1 firing when the parameter value is large. As the bifurcation parameter value increases, the probability of the initial values of a period-1 firing neuron that lead to complete synchronization of period-1 firing increases, while that leading to complete synchronization of the resting condition decreases. It shows that the attraction domain of a coexisting behavior is larger, the probability of initial values leading to complete synchronization of this behavior is higher. The bifurcations of the coupled system are investigated and discussed. The results reveal the complex dynamics of synchronization behaviors of the coupled system composed of neurons with the coexisting resting condition and period-1 firing, and are helpful to further identify the dynamics of the spatiotemporal behaviors of the central nervous system.