Volume 13 Supplement 1

Twenty First Annual Computational Neuroscience Meeting: CNS*2012

Open Access

A model for simulating Local Field Potential in the thalamus of Essential Tremor patient during deep brain stimulation

  • Ishita Basu1Email author,
  • Daniela Tuninetti1,
  • Daniel Graupe1, 2 and
  • Konstantin V Slavin3
BMC Neuroscience201213(Suppl 1):P40

DOI: 10.1186/1471-2202-13-S1-P40

Published: 16 July 2012

Local Field Potential (LFP) is simulated by adding the neuronal membrane potential of a group of model neurons surrounding a typical intra-operative micro-electrode used to localize the target for deep brain stimulation (DBS) in Essential Tremor (ET). Each neuronal membrane potential is modeled as an Ornstein Uhlenbeck Process (OUP) and the model parameters are extracted from real neuronal recording during DBS surgery following the method in [1].

Before any test stimulation through the macro ring of the microelectrode, the LFP is simulated according to:

Where, x is the membrane potential of each neuron and r is the distance from the recording tip. The simulated signal, lfp is then low pass filtered (1-40 Hz) to produce the final signal.

With stimulation on, the LFP is simulated similarly but only over time intervals in between two stimulation pulse which is free of the stimulation voltage artifact caused by amplifier saturation. The stimulation artifact time duration is estimated as:

Where, l(r) is the stimulation intensity at distance r from the recording tip and follows a similar inverse square decay as in (1), τ is the membrane time constant (an OUP model parameter) and S is the neuronal firing threshold. Thus the LFP simulated during stimulation depends on the stimulation intensity, pulse width and frequency.

The stimulation parameters used are: frequency: 160 Hz (125 samples/period), pulse width: 0.4 ms, intensity:1.5 mA

It was found that the power in the (5-12) Hz band decreases by 26.4% during stimulation. There is increased synchronization in neuronal activity of the thalamus in the theta band (4-7 Hz) for patients with ET [2]. The suppression of tremor by DBS is reflected in the reduction of power in the tremor frequency band.

Such a model , can potentially be used to find an optimal set of stimulation parameters that will produce the maximum suppression in the tremor power.
Figure 1

Spectrum of LFP before stimulation (in red) and during stimulation (in blue).

Authors’ Affiliations

Department of Electrical & Computer Engineering, University of Illinois at Chicago
Department of Bioengineering, University Illinois at Chicago
Department of Neurosurgery, University Illinois at Chicago


  1. Basu I, Graupe D, Tuninetti D, Slavin KV: Stochastic modeling of the neuronal activity in the subthalamic nucleus and model parameter identification from Parkinson patient data. Biol Cybern. 2010, 103 (4): 273-283. 10.1007/s00422-010-0397-3.View ArticlePubMedGoogle Scholar
  2. Kane A, Hutchison WD, Hodaie M, Lozano AM, Dostrovsky JO: Enhanced synchronization of thalamic theta band local field potentials in patients with essential tremor. Experimental Neurology. 2009, 217: 171-176. 10.1016/j.expneurol.2009.02.005.View ArticlePubMedGoogle Scholar


© Basu et al; licensee BioMed Central Ltd. 2012

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.