Hodgkin Huxley neuron - The corner stone model to biophysical neurons?
The Hodgkin-Huxley model defines a four-dimensional nonlinear dynamical system describing the evolution of the membrane potential and ionic gating states. The state vector is
\[\mathbf{x}(t)=(V(t), m(t), h(t), n(t)),\]and the dynamics follow from charge conservation on the membrane combined with voltage-dependent Markov kinetics for ion channel gates.
The membrane equation is
\[C_m \frac{d V}{d t}=I_{\mathrm{ext}}(t)-\bar{g}_{\mathrm{Na}} m^3 h\left(V-E_{\mathrm{Na}}\right)-\bar{g}_{\mathrm{K}} n^4\left(V-E_{\mathrm{K}}\right)-g_L\left(V-E_L\right),\]where the exponents reflect the assumed number of independent activation gates per channel. The nonlinearities enter multiplicatively through both conductance and driving force, making the voltage equation cubic in $V$ and polynomial in the gating variables.
Each gating variable $x \in{m, h, n}$ follows first-order kinetics:
\[\frac{d x}{d t}=\alpha_x(V)(1-x)-\beta_x(V) x .\]These equations can be rewritten as
\[\frac{d x}{d t}=\frac{x_{\infty}(V)-x}{\tau_x(V)}, \quad x_{\infty}(V)=\frac{\alpha_x(V)}{\alpha_x(V)+\beta_x(V)}, \quad \tau_x(V)=\frac{1}{\alpha_x(V)+\beta_x(V)} .\]Linearization about this equilibrium yields a Jacobian whose eigenvalues determine excitability. In the classical parameter regime, the neuron behaves as an excitable system near a stable node or focus: small perturbations decay, while sufficiently large perturbations trigger a large excursion in phase space before returning to rest.
An action potential corresponds to a transient trajectory driven by a fast-slow decomposition. Sodium activation $m$ operates on the fastest timescale and can often be approximated as quasi-steady,
\[m(t) \approx m_{\infty}(V(t)),\]reducing the effective dimensionality during spike initiation. Potassium activation 𝑛 n and sodium inactivation ℎ evolve more slowly and control repolarization and refractoriness. The separation of timescales induces a relaxation-type oscillatory excursion even though the system does not possess a limit cycle under constant subthreshold input.
Threshold emerges geometrically as the intersection of voltage and gating nullclines. Rather than a fixed voltage level, it is a separatrix in state space: initial conditions on one side decay to rest, while those on the other evolve toward rapid sodium activation. This explains the variability and history dependence of spike initiation in Hodgkin–Huxley neurons.
From a stochastic perspective, the deterministic equations describe the mean-field limit of large channel populations. Channel noise can be incorporated by adding state-dependent diffusion terms to the gating equations, yielding multiplicative Langevin dynamics that perturb trajectories near the separatrix and smear threshold behavior.
Mathematically, the Hodgkin–Huxley model serves as a canonical conductance-based system: nonlinear, non-gradient, and strongly dissipative, yet structured by biophysical constraints. Reduced models—FitzHugh–Nagumo, Morris–Lecar, and integrate-and-fire variants—can be derived through systematic dimensional reduction or singular perturbation, but the full Hodgkin–Huxley equations remain the reference system in which action potentials are not abstract events, but consequences of coupled nonlinear kinetics governed by voltage-dependent ion transport.