Abstract. Recent technological advances including brain imaging (higher resolution in space and time), miniaturization of integrated circuits (nanotechnologies), and acceleration of computation speed (Moore's Law), combined with interpenetration between neuroscience, mathematics, and physics have led to the development of more biologically plausible computational models and novel therapeutic strategies. Today, mathematical models of irreversible medical conditions such as Parkinson's disease (PD) are developed and parameterised based on clinical data. How do these evolutions have a bearing on deep brain stimulation (DBS) of patients with PD? We review how the idea of DBS, a standard therapeutic strategy used to attenuate neurological symptoms (motor, psychiatric), has emerged from past experimental and clinical observations, and present how, over the last decade, computational models based on different approaches (phase oscillator models, spiking neuron network models, population-based models) have started to shed light onto DBS mechanisms. Finally, we explore a new mathematical modelling approach based on neural field equations to optimize mechanisms of brain stimulation and achieve finer control of targeted neuronal populations. We conclude that neuroscience and mathematics are crucial partners in exploring brain stimulation and this partnership should also include other domains such as signal processing, control theory and ethics.