A renumbering method to decrease matrix banding in equations describing branched neuron-like structures.

R M Eichler West; G L Wilcox

Home/Archive/PubMed

PubMed · 8884608

A renumbering method to decrease matrix banding in equations describing branched neuron-like structures.

R M Eichler West·G L Wilcox

View original source ↗Export metadata ↓

Abstract

The solution to matrix equations which describe branched neuron-like structures can be made more efficient by minimizing matrix banding. This can be accomplished through the reordering of the compartmental numbering system. The renumbering method presented here extends upon the numbering method of Hines ((1984) Int. J. Biomed. Comput., 15: 69-76). A demonstration of efficient numbering will be presented for several general cases of branching structures. Theoretical computational savings can be estimated for the above structures. An algorithm to renumber a matrix already in Hines form will be described. Branched nerve equations, electrical networks and chemical reaction models are examples of systems which can benefit from this application.

Explore related subjects

Models, Neurological Neurons Journal Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

R M Eichler West, G L Wilcox. 1996. A renumbering method to decrease matrix banding in equations describing branched neuron-like structures.. https://pubmed.ncbi.nlm.nih.gov/8884608/

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2025

Balanced state of networks of winner-take-all units.

Irregularly timed action potentials, or spikes, are pervasively observed in the brain activity of awake mammals. However, the role of this temporal irregularity in neural computation is still not well understood. In canonical network models irregular spiking emerges via balanced, fluctuating input currents, leading to collective responses that track inputs linearly. How networks characterized by irregular spiking could support flexible nonlinear dynamics needed for general-purpose computation remains under ongoing debate. Here we characterize the dynamics of networks whose elementary unit is not a single neuron but a small group of neurons, with distinct tunings, that compete at each timestep via a winner-take-all (WTA) interaction. While WTA has long been proposed as an elementary functional motif in the brain and represents a powerful computational primitive, how large networks of such units behave has received less investigation. We show that these networks, like classic excitatory-inhibitory balanced networks, exhibit a chaotic fluctuation-driven regime characterized by sustained irregular activity resembling realistic cortical spiking, which we interpret as a multidimensional balance spread over several competing neural populations with different tunings. We develop a mean-field theory for the network, which shows how irregular spiking sustained by time-varying input fluctuations can support flexible nonlinear collective dynamics. Using the theory we predict and verify network regimes in which input fluctuations alone yield multistability, stable sequence generation, or complex heterogeneous firing rate dynamics-three core dynamical primitives thought to underlie memory-dependent neural computation-via consistent Poisson-like spiking produced through chaos. Thus, networks of WTA units support a chaotic fluctuation-driven regime characterized by irregular spiking that can power complex nonlinear collective dynamics. This represents a new model of brain activity capable of simultaneously reproducing realistic spike trains and diverse nonlinear firing rate patterns well posed for flexible computation, and which can be trained or fit to data.

Models, Neurological↗

PubMed1996

Chaos and synchrony in a model of a hypercolumn in visual cortex.

Neurons in cortical slices emit spikes or bursts of spikes regularly in response to a suprathreshold current injection. This behavior is in marked contrast to the behavior of cortical neurons in vivo, whose response to electrical or sensory input displays a strong degree of irregularity. Correlation measurements show a significant degree of synchrony in the temporal fluctuations of neuronal activities in cortex. We explore the hypothesis that these phenomena are the result of the synchronized chaos generated by the deterministic dynamics of local cortical networks. A model of a "hypercolumn" in the visual cortex is studied. It consists of two populations of neurons, one inhibitory and one excitatory. The dynamics of the neurons is based on a Hodgkin-Huxley type model of excitable voltage-clamped cells with several cellular and synaptic conductances. A slow potassium current is included in the dynamics of the excitatory population to reproduce the observed adaptation of the spike trains emitted by these neurons. The pattern of connectivity has a spatial structure which is correlated with the internal organization of hypercolumns in orientation columns. Numerical simulations of the model show that in an appropriate parameter range, the network settles in a synchronous chaotic state, characterized by a strong temporal variability of the neural activity which is correlated across the hypercolumn. Strong inhibitory feedback is essential for the stabilization of this state. These results show that the cooperative dynamics of large neuronal networks are capable of generating variability and synchrony similar to those observed in cortex. Auto-correlation and cross-correlation functions of neuronal spike trains are computed, and their temporal and spatial features are analyzed. In other parameter regimes, the network exhibits two additional states: synchronized oscillations and an asynchronous state. We use our model to study cortical mechanisms for orientation selectivity. It is shown that in a suitable parameter regime, when the input is not oriented, the network has a continuum of states, each representing an inhomogeneous population activity which is peaked at one of the orientation columns. As a result, when a weakly oriented input stimulates the network, it yields a sharp orientation tuning. The properties of the network in this regime, including the appearance of virtual rotations and broad stimulus-dependent cross-correlations, are investigated. The results agree with the predictions of the mean field theory which was previously derived for a simplified model of stochastic, two-state neurons. The relation between the results of the model and experiments in visual cortex are discussed.

Models, Neurological↗

PubMed1996

Coding of odor intensity in a steady-state deterministic model of an olfactory receptor neuron.

The coding of odor intensity by an olfactory receptor neuron model was studied under steady-state stimulation. Our model neuron is an elongated cylinder consisting of the following three components: a sensory dendritic region bearing odorant receptors, a passive region consisting of proximal dendrite and cell body, and an axon. First, analytical solutions are given for the three main physiological responses: (1) odorant-dependent conductance change at the sensory dendrite based on the Michaelis-Menten model, (2) generation and spreading of the receptor potential based on a new solution of the cable equation, and (3) firing frequency based on a Lapicque model. Second, the magnitudes of these responses are analyzed as a function of odorant concentration. Their dependence on chemical, electrical, and geometrical parameters is examined. The only evident gain in magnitude results from the activation-to-conductance conversion. An optimal encoder neuron is presented that suggests that increasing the length of the sensory dendrite beyond about 0.3 space constant does not increase the magnitude of the receptor potential. Third, the sensitivities of the responses are examined as functions of (1) the concentration at half-maximum response, (2) the lower and upper concentrations actually discriminated, and (3) the width of the dynamic range. The overall gain in sensitivity results entirely from the conductance-to-voltage conversion. The maximum conductance at the sensory dendrite appears to be the main tuning constant of the neuron because it determines the shift toward low concentrations and the increase in dynamic range. The dynamic range of the model cannot exceed 5.7 log units, for a sensitivity increase at low odor concentration is compensated by a sensitivity decrease at high odor concentration.

Models, Neurological↗