skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Generalized rate-code model for neuron ensembles with finite populations

Abstract

We have proposed a generalized Langevin-type rate-code model subjected to multiplicative noise, in order to study stationary and dynamical properties of an ensemble containing a finite number N of neurons. Calculations using the Fokker-Planck equation have shown that, owing to the multiplicative noise, our rate model yields various kinds of stationary non-Gaussian distributions such as {gamma}, inverse-Gaussian-like, and log-normal-like distributions, which have been experimentally observed. The dynamical properties of the rate model have been studied with the use of the augmented moment method (AMM), which was previously proposed by the author from a macroscopic point of view for finite-unit stochastic systems. In the AMM, the original N-dimensional stochastic differential equations (DEs) are transformed into three-dimensional deterministic DEs for the means and fluctuations of local and global variables. The dynamical responses of the neuron ensemble to pulse and sinusoidal inputs calculated by the AMM are in good agreement with those obtained by direct simulation. The synchronization in the neuronal ensemble is discussed. The variabilities of the firing rate and of the interspike interval are shown to increase with increasing magnitude of multiplicative noise, which may be a conceivable origin of the observed large variability in cortical neurons.

Authors:
 [1]
  1. Department of Physics, Tokyo Gakugei University, Koganei, Tokyo 184-8501 (Japan)
Publication Date:
OSTI Identifier:
21072434
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physical Review. E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics; Journal Volume: 75; Journal Issue: 5; Other Information: DOI: 10.1103/PhysRevE.75.051904; (c) 2007 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND TECHNOLOGY; BRAIN; COMPUTERIZED SIMULATION; DISTRIBUTION; FLUCTUATIONS; FOKKER-PLANCK EQUATION; GAUSS FUNCTION; MOMENTS METHOD; NERVE CELLS; NEURAL NETWORKS; NOISE; STOCHASTIC PROCESSES; SYNCHRONIZATION

Citation Formats

Hasegawa, Hideo. Generalized rate-code model for neuron ensembles with finite populations. United States: N. p., 2007. Web. doi:10.1103/PHYSREVE.75.051904.
Hasegawa, Hideo. Generalized rate-code model for neuron ensembles with finite populations. United States. doi:10.1103/PHYSREVE.75.051904.
Hasegawa, Hideo. Tue . "Generalized rate-code model for neuron ensembles with finite populations". United States. doi:10.1103/PHYSREVE.75.051904.
@article{osti_21072434,
title = {Generalized rate-code model for neuron ensembles with finite populations},
author = {Hasegawa, Hideo},
abstractNote = {We have proposed a generalized Langevin-type rate-code model subjected to multiplicative noise, in order to study stationary and dynamical properties of an ensemble containing a finite number N of neurons. Calculations using the Fokker-Planck equation have shown that, owing to the multiplicative noise, our rate model yields various kinds of stationary non-Gaussian distributions such as {gamma}, inverse-Gaussian-like, and log-normal-like distributions, which have been experimentally observed. The dynamical properties of the rate model have been studied with the use of the augmented moment method (AMM), which was previously proposed by the author from a macroscopic point of view for finite-unit stochastic systems. In the AMM, the original N-dimensional stochastic differential equations (DEs) are transformed into three-dimensional deterministic DEs for the means and fluctuations of local and global variables. The dynamical responses of the neuron ensemble to pulse and sinusoidal inputs calculated by the AMM are in good agreement with those obtained by direct simulation. The synchronization in the neuronal ensemble is discussed. The variabilities of the firing rate and of the interspike interval are shown to increase with increasing magnitude of multiplicative noise, which may be a conceivable origin of the observed large variability in cortical neurons.},
doi = {10.1103/PHYSREVE.75.051904},
journal = {Physical Review. E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics},
number = 5,
volume = 75,
place = {United States},
year = {Tue May 15 00:00:00 EDT 2007},
month = {Tue May 15 00:00:00 EDT 2007}
}
  • The authors propose a new class of cluster growth models where growth sites have a finite lifetime tau, which contains as special cases the Eden model (tau = infinity) and the kinetic growth walk (tau = 1). For finite but large tau values the growth process can be characterized by a crossover time t/sub x/; for times below t/sub x/ an Eden-type cluster is formed, while for times above t/sub x/ the growth process belongs to the universality class of the self-avoiding random walk. The crossover time increases monotonically with tau. They develop a scaling theory for the time evolutionmore » of the mean end-to-end distance between the seed and the last-added site, and for the average number of growth sites by which the kinetics of the growth process can be characterized. They test this scaling theory by extensive Monte Carlo simulations. They also extend their results to inhomogeneous media (percolation systems).« less
  • The effects of millimeter waves (mm-waves, 75 GHz) and temperature elevation on the firing rate of the BP-4 pacemaker neuron of the pond snail Lymnaea stagnalis were studied by using microelectrode techniques. The open end of a rectangular waveguide covered with a thin Teflon film served as a radiator. Specific absorption rates (SARs), measured in physiological solution at the radiator outlet, ranged from 600 to 4,200 W/kg, causing temperature rises from 0.3 to 2.2 C, respectively. Irradiation at an SAR of 4,200 W/kg caused a biphasic change in the firing rate, i.e., a transient decrease in the firing rate followedmore » by a gradual increase to a new level that was 68 {+-} 21% above control. The biphasic changes in the firing rate were reproduced by heating under the condition that the magnitude (2 C) and the rate of temperature rise were equal to those produced by the irradiation. The addition of 0.05 mM of ouabain caused the disappearance of transient responses of the neuron to the irradiation. It was shown that the rate of temperature rise played an important role in the development of a transient neuronal response. The threshold stimulus for a transient response of the BP-4 neutron found in warming experiments was a temperature rise of 0.0025 C/s.« less
  • To test the hypothesis that the functional state of hypothalamic LHRH neurons and pituitary gonadotrophs might alter their radiosensitivity, we determined the experimental conditions under which the gonadotropin response to castration could be impaired by a single dose of cranial irradiation. Single doses of cranial irradiation greater than 2000 rads were lethal to unshielded rats. Shielding of the oropharynx and esophagus allowed the animals to survive doses up to 5000 rads. Doses between 2000 and 5000 rads had no effect on basal gonadotropin levels for as long as 3 months after irradiation. Irradiation caused a dose- and time-dependent impairment, however,more » in the gonadotropin response to castration. Impairment of the gonadotropin levels of castrate animals occurred in animals that were irradiated either before or after castration. However, rats irradiated in the castrate state showed a decreased susceptibility to irradiation damage. Additionally, stimulation of the pituitary by LHRH agonist (LHRHa) 3 h before irradiation significantly reduced the impairment of gonadotropin secretion 12-20 weeks after irradiation (P less than 0.05). Thus, increased functional activity of the rat hypothalamus or pituitary at the time of irradiation, induced by either castration or acute LHRHa administration, was associated with some protection against the gonadotropin-lowering effect of irradiation. Based upon these data, we hypothesize that stimulation of gonadotropin secretion at the time of therapeutic cranial irradiation in humans might protect against subsequent impairment of gonadotropin secretion.« less
  • A Monte Carlo charge breeding code (MCBC) is being developed by FAR-TECH, Inc. to model the capture and charge breeding of 1+ ion beam in an electron cyclotron resonance ion source (ECRIS) device. The ECRIS plasma is simulated using the generalized ECRIS model which has two choices of boundary settings, free boundary condition and Bohm condition. The charge state distribution of the extracted beam ions is calculated by solving the steady state ion continuity equations where the profiles of the captured ions are used as source terms. MCBC simulations of the charge breeding of Rb+ showed good agreement with recentmore » charge breeding experiments at Argonne National Laboratory (ANL). MCBC correctly predicted the peak of highly charged ion state outputs under free boundary condition and similar charge state distribution width but a lower peak charge state under the Bohm condition. The comparisons between the simulation results and ANL experimental measurements are presented and discussed.« less