After noticing a typo (hello missing negative sign!) in the paper, and debugging two typos in my implementation, I have what seem to be two working models (VCN bushy cells Types I-c and II). Here I’m applying a constant 100 pA stimulus to them. Already I see significant differences between the two cell types: