Phase description of the Huber-Braun neuron model for mammalian cold receptors

The spiking activity of mammalian cold receptors is described by the Huber-Braun neuron model. Sweeping temperature as a control parameter across a biologically relevant range this model exhibits a complex bifurcation structure seen in the sequence of interspike intervals. The model’s distinctive feature is the interaction between a fast spike generating dynamics and a slow subthreshold oscillation. Viewing the spike generation as a cycle, the dynamics may also be modeled phenomenologically by two phases, one for the spike cycle and the second for the slow subthreshold oscillation. In fact, a phase model of temperature-dependent mammalian cold receptors was already proposed by Roper et al. (2000). Here we follow their approach and investigate to what extent this model is able to reproduce the bifurcation patterns of the Huber-Braun model. Special attention is paid to the tonic firing to bursting transition observed in the low temperature range.