11Brain neurons exhibit complex electroresponsive properties -including intrinsic subthreshold 12 oscillations and pacemaking, resonance and phase-reset -which are thought to play a critical role 13 in controlling neural network dynamics. Although these properties emerge from detailed 14representations of molecular-level mechanisms in "realistic" models, they cannot usually be 15 generated by simplified neuronal models (although these may show spike-frequency adaptation 16 and bursting). We report here that this whole set of properties can be generated by the extended 17 generalized leaky integrate-and-fire (E-GLIF) neuron model. E-GLIF derives from the GLIF 18 model family and is therefore mono-compartmental, keeps the limited computational load typical 19 of a linear low-dimensional system, admits analytical solutions and can be tuned through gradient-20 descent algorithms. Importantly, E-GLIF is designed to maintain a correspondence between model 21 parameters and neuronal membrane mechanisms through a minimum set of equations. In order to 22 test its potential, E-GLIF was used to model a specific neuron showing rich and complex 23 electroresponsiveness, the cerebellar Golgi cell, and was validated against experimental 24 electrophysiological data recorded from Golgi cells in acute cerebellar slices. During simulations, 25 E-GLIF was activated by stimulus patterns, including current steps and synaptic inputs, identical 26 to those used for the experiments. The results demonstrate that E-GLIF can reproduce the whole 27 set of complex neuronal dynamics typical of these neurons -including intensity-frequency curves, 28 spike-frequency adaptation, depolarization-induced and post-inhibitory rebound bursting, 29 spontaneous subthreshold oscillations, resonance and phase-reset, -providing a new effective tool 30to investigate brain dynamics in large-scale simulations. 31Word count (body of the text): 8047. Figures: 9. Tables: 5 . 32 33 65 Gerstner, 2005). However, the nonlinearity entailed more difficulties in optimizing model 66 parameters and in computational efficiency. Therefore, recently, new linear adaptive point models 67have been developed (Generalized LIF, GLIF), with spike-triggered currents and moving threshold 68 as the source of adaptation (Mihalaş and Niebur, 2009) and with stochastic processes in firing 69 emission (Pozzorini et al., 2015; Rössert et al., 2016). The possibility to use a linear and 70 analytically solvable neuron model is fundamental when simulating large-scale SNNs, since 71 computational efficiency can be enhanced without severe loss in spike time accuracy and realism 72 (Hanuschkin et al., 2010). However, GLIF can hardly generate phenomena like subthreshold 73 oscillations, resonance and phase-reset, which are critical for large-scale network entrainment and 74 communication (Buzsáki, 2004; Buzsáki and Draghun, 2004). 75We propose here an extended GLIF (E-GLIF) model, which achieves a sound compromise 76 between model complexity, biological plausibility and computational efficienc...