Skip to content


BMC Neuroscience

Open Access

Is the Langevin phase equation an efficient model for stochastic limit cycle oscillators in real neurons?

  • Keisuke Ota1Email author,
  • Toshiaki Omori2, 3,
  • Shigeo Watanabe4,
  • Hiroyoshi Miyakawa4,
  • Masato Okada2, 3 and
  • Toru Aonishi1, 3
BMC Neuroscience200910(Suppl 1):P236

Published: 13 July 2009


The Langevin phase equation /dt = 1 + Z(ϕ)(G(t) + σs(t)), where ϕ is the phase, which is disturbed by a perturbation G and Langevin force s of intensity σ, and Z is the phase response curve (PRC), has been deemed to be a good model for stochastic limit cycle oscillators [1], and it has been extensively used in theoretical neuroscience as a model neural oscillator [2]. Inspired by the theoretical research, experimental researchers have measured PRCs, but none of them have identified the Langevin phase equation for real neurons directly. In fact, biological experiments have yet to show whether this equation is a good model for neural oscillators.


Here, we demonstrate that the Langevin phase equation is a good model for neural oscillators in rat hippocampal CA1 pyramidal neurons, through two steps in the same neuron. (A) Estimation step: We estimated the parameters of the Langevin phase equation, i.e., a PRC and the intensity of the Langevin force, from physiological noisy data (phase shifts disturbed by one-shot rectangle pulses 20 pA, 5 msec) by using the MAP estimation algorithm, from our previous study [3]. Figure 1 shows the estimated PRC (a) and the hyper parameters (b). The estimated PRC was fully positive, which indicates that hippocampal CA1 pyramidal neurons could be classified as Type I excitability neurons, while the estimated Langevin force intensity suggests that the inherent noise in the recording neuron is 9.0 pA. (B) Prediction step: We injected two different periodic perturbations, consisting successive rectangular pulses. One consisted of 20 pA, 5 msec pulses (the same values used in the PRC estimation), the other 10 pA, 20 msec pulses. We obtained histograms of phase differences between spikes and pulses for each perturbation. Then we ascertained whether the Fokker Planck equation (FPE) derived from the Langevin phase equation with the estimated PRC and Langevin force intensity could predict these histograms.
Figure 1

(a) Estimated PRC (solid line) from the phase shifts (crosses) disturbed by one-shot perturbations (20 pA, 5 msec) in a pyramidal neuron of rat hippocampal CA1. (b) The landscape of the marginal likelihood, where the maximum represents the estimated values of hyper parameters.

Results and discussion

The distribution derived from FPE was in good agreement with the experimental histogram (Fig. 2), even when the shape of the perturbation was not the same as the one used in the PRC measurement (Fig. 2b). This result suggests that the Langevin phase equation can describe the stochastic dynamics of neural oscillators, no matter what perturbation disturbs them.
Figure 2

Comparison of distributions predicted by the Langevin phase equation (green line) and phase difference histograms calculated from physiological data (red line). The upper right of each graph displays the shape of the periodic perturbation, (a) 20 pA, 5 msec and (b) 10 pA, 20 msec.

Authors’ Affiliations

Dept. of Computational Intelligence and System Science, Tokyo Tech., Yokohama, Japan
Dept. of Complexity Science and Engineering, Univ. of Tokyo, Kashiwa, Japan
Brain Science Institute, RIKEN, Wako, Japan
Dept. of Life Science, Tokyo Univ. of Pharmacy and Life Science, Hachioji, Japan


  1. Nakao H, Teramae J, Ermentrout GB: Comment on "phase reduction of stochastic limit cycle oscillators". []
  2. Ermentrout GB, Galan RF, Urban NN: Relating neural dynamics to neural coding. Phys Rev Lett. 2007, 99: 248103-10.1103/PhysRevLett.99.248103.PubMed CentralPubMedView ArticleGoogle Scholar
  3. Ota K, Omori T, Aonishi T: MAP estimation algorithm for phase reponse curves based on analysis of the observation process. J Comput Neurosci.Google Scholar


© Ota et al; licensee BioMed Central Ltd. 2009

This article is published under license to BioMed Central Ltd.