The (Wang 1994) papers model was replicated in python by (Detorakis 2016). "The model is conductance-based and takes advantage of the
interplay between a T-type calcium current and a non-specific cation sag current and
thus, it is able to generate spindle and delta rhythms." The model also generates intermittent phase locking, non periodic firing, bursts, and tonic spike patterns.
Wang XJ (1994) Multiple dynamical modes of thalamic relay neurons: rhythmic bursting and intermittent phase-locking. Neuroscience 59(1):21-31 [PubMed]
Detorakis G (2016) [Re] Multiple dynamical modes of thalamic relay
neurons: rhythmic bursting and intermittent
phase-locking ReScience 2(1):1-11