Periodic Patterns in Recurrent Networks of Varying Neuron Types

Dylan Hirsch-Shell

Abstract

We simulate a recurrent network of 100 integrate and fire neurons. Each neuron in the network forms a relatively strong inhibitory connection with every other neuron and may form a weak excitatory connection with certain other neurons. (Autapses are also allowed in both excitatory and inhibitory cases.) If all the neurons are driven with a constant current and inhibition dominates excitation, then after a brief period a periodic firing pattern emerges in the network. Next, we simulate a network with the same structural properties, except it is composed of neurons that follow the quadratic model. We want to see if we still get the periodic patterns, and if so how the characteristics of the patterns change.

Recurrent Networks

Recurrent neural networks consist of neurons connected via excitatory and inhibitory synapses. Such networks are believed to be involved in many important computations in the brain [1] as well as in serving to represent objects from the real world in the brain. [2] The topologies of the networks can vary tremendously, but when inhibition dominates overall, the networks will tend to converge upon a stable periodic pattern.[1] Example network:
Schematic of an example recurrent network
(from [1])

Integrate-and-fire Neuron

The first neuron type we looked at was the Integrate-and-Fire Neuron Model. The mathematical model for an integrate-and-fire neuron, j, is as follows:

where &tau is the membrane time constant, Vrest is the resting potential, Vj is the membrane potential, Ij is the external input, and Iexc and Iinh are the inputs from the excitatory and inhibitory synapses, respectively.
More specifically,

and

where gexc,j is the excitatory synaptic conductance and ginh,j is the inhibitory conductance, of neuron j.

The values of gexc,j and ginh,j follow the same rule. Whenever neuron i fires (i.e. Vi > V&theta) not only does Vi get reset to Vreset (where Vreset < V&theta), but also for any neuron j that i forms a connection with:

for both excitatory and inhibitory synapses.
When no connecting neurons have fired, however, the conductances show exponential decay like so:

where g0,j is the conductance after the last spike, and &tausyn is the synaptic time constant (typically much shorter than the membrane time constant). Gi,j is a constant representing the strength of the connection from i to j.

With a network of 100 of these integrate-and-fire neurons, it doesn't take long for a nice pattern to emerge:


Quadratic Neuron

We thought it would be interesting to see if similar patterns would emerge if we used a different neuron model within the same neural network structure. The model we chose was the quadratic neuron:

where U has units of potential (in our case, we set U to be near Vrest). The conductances follow the same rules as in the integrate-and-fire neuron.

Here is an example of a spike raster plot from such a network:



With the quadratic neurons, the networks behaved a bit differently. For one, usually more neurons seemed to fire at least once, but fewer remained firing after a short while, and the ones that were usually had a high firing frequency. Despite running the simulation three times longer than with the integrate-and-fire neurons, I still did not see a consistent pattern emerge. It seems to be very close, but not quite right. There are usually a few neurons that have distinctive patterns in their firing, but even those are a little inconsistent. And then there are always a few neurons randomly firing at a high rate without a discernible pattern. For instance in the above plot, neurons 3, 6, 7, 9 and 10 are never quite consistent. Perhaps making precise changes in some of the parameters of the system would result in the emergence of perfect stable patterns, however none of the changes I tried seemed to do the trick.

Conclusions

While I could not get a stable periodic pattern to appear with perfect repetition with the quadratic neural network, I still think that our results are significant. There are patterns to the firings of certain groups of neurons, to within a certain amount of error. Unfortunately, I think the system as it has been modeled is so sensitive to changes in the parameters that it is very difficult to achieve a set of parameters that behaves ideally. This is where our ability to model nature has broken down, for networks of real neurons are more robust than this. They are not hypersensitive to changes in inputs or conductances, as our model here is. It would still be good to demonstrate perfect stable periodic patterns with a myriad of different neuron models to show that our underlying assumptions about pattern-forming in recurrent networks are correct. However, again, I could not seem to do that given the limitations of my system. Perhaps my results here would be the equivalent of a fuzzy memory or a calculation gone slightly awry in a real brain.

References

[1] Dezhe Z. Jin, "Fast convergence of spike sequences to periodic patterns in recurrent networks", Physical Review Letters, 89, 208102 (2002). (PDF)

[2] Hansel D, Mato G. "Existence and stability of persistent states in large neuronal networks," Physical Review Letters, 86, 4175 (2001).


Last modified: 5/23/2003