Neurons are poised near the edge of chaos

This paper shows the action potential (spikes) generated from the Hodgkin-Huxley equations emerges near the edge of chaos consisting of a tiny subset of the locally active regime of the HH equations. The main result proves that the eigenvalues of the 4 × 4 Jacobian matrix associated with the mathematically intractable system of four nonlinear differential equations are identical to the zeros of a scalar complexity function from complexity theory. Moreover, we show the loci of a pair of complex-conjugate zeros migrate continuously as a function of an externally applied DC current excitation emulating the net synaptic excitation current input to the neuron. In particular, the pair of complex-conjugate zeros move from a subcritical Hopf bifurcation point at low excitation current to a super-critical Hopf bifurcation point at high excitation current. The spikes are generated as the excitation current approaches the vicinity of the edge of chaos, which leads to the onset of the subcritical Hopf bifurcation regime. It follows from this in-depth qualitative analysis that local activity is the origin of spikes.