Power laws from linear neuronal cable theory: power spectral densities of the soma potential, soma membrane current and single-neuron contribution to the EEG.
Power laws, that is, power spectral densities (PSDs) exhibiting 1/f(α) behavior for large frequencies f, have been observed both in microscopic (neural membrane potentials and currents) and macroscopic (electroencephalography; EEG) recordings. While complex network behavior has been suggested to be at the root of this phenomenon, we here demonstrate a possible origin of such power laws in the biophysical properties of single neurons described by the standard cable equation. Taking advantage of the analytical tractability of the so called ball and stick neuron model, we derive general expressions for the PSD transfer functions for a set of measures of neuronal activity: the soma membrane current, the current-dipole moment (corresponding to the single-neuron EEG contribution), and the soma membrane potential. These PSD transfer functions relate the PSDs of the respective measurements to the PSDs of the noisy input currents. With homogeneously distributed input currents across the neuronal membrane we find that all PSD transfer functions express asymptotic high-frequency 1/f(α) power laws with power-law exponents analytically identified as α∞(I) = 1/2 for the soma membrane current, α∞(p) = 3/2 for the current-dipole moment, and α∞(V) = 2 for the soma membrane potential. Comparison with available data suggests that the apparent power laws observed in the high-frequency end of the PSD spectra may stem from uncorrelated current sources which are homogeneously distributed across the neural membranes and themselves exhibit pink (1/f) noise distributions. While the PSD noise spectra at low frequencies may be dominated by synaptic noise, our findings suggest that the high-frequency power laws may originate in noise from intrinsic ion channels. The significance of this finding goes beyond neuroscience as it demonstrates how 1/f(α) power laws with a wide range of values for the power-law exponent α may arise from a simple, linear partial differential equation.
Dept. of Mathematical Sciences and Technology, Norwegian University of Life Sciences, Ås, Norway; Letten Centre and GliaLab, Institute of Basic Medical Sciences, University of Oslo, Oslo, Norway; Centre for Molecular Medicine Norway, University of Oslo, Oslo, Norway.
Klas H
KH
Pettersen
Dept. of Mathematical Sciences and Technology, Norwegian University of Life Sciences, Ås, Norway; Department of Neuroscience and Pharmacology, University of Copenhagen, Copenhagen, Denmark; Dept. of Computational Biology, School of Computer Science and Communication, Royal Institute of Technology (KTH), Stockholm, Sweden.
Henrik
H
Lindén
Dept. of Mathematical Sciences and Technology, Norwegian University of Life Sciences, Ås, Norway; Inst. of Neuroscience and Medicine (INM-6), Computational and Systems Neuroscience & Inst. for Advanced Simulation (IAS-6), Theoretical Neuroscience, Jülich Research Centre and JARA, Jülich, Germany.
Tom
T
Tetzlaff
Dept. of Mathematical Sciences and Technology, Norwegian University of Life Sciences, Ås, Norway; Dept. of Physics, University of Oslo, Oslo, Norway.
Gaute T
GT
Einevoll
D019295:Computational Biology; D003198:Computer Simulation; D004569:Electroencephalography; D008564:Membrane Potentials / Q000502:physiology*; D008959:Models, Neurological; D009474:Neurons / Q000502:physiology*