it has been hypothesized that , in many areas of the brain , having different brain functionality , repeatable precise spatiotemporal patterns of spikes play a crucial role in coding and storage of information .
temporally structured replay of spatiotemporal patterns have been observed to occur during sleep , both in the cortex and hippocampus ( nadasdy et al .
2009 ) , and it has been hypothesized that this replay may subserve memory consolidation . the sequential reactivation of hippocampal place cells , corresponding to previously experienced behavioral trajectories , has been observed also in the awake state ( awake replay ) ( diba and buzsaki 2007 ; davidson et al .
awake replay may reflect trajectories through either the current environment or previously , spatially remote , visited environments .
a possible interpretation is that spatiotemporal patterns , stored in the plastic synaptic connections of hippocampus , are retrieved when a cue activates the emergence of a stored pattern , allowing these patterns to be replayed and then consolidated in distributed circuits beyond the hippocampus ( carr et al .
cross - correlogram analysis revealed that in prefrontal cortex the time scale of reactivation of firing patterns during post - behavioral sleep was compressed five- to eightfold relative to waking state ( euston et al .
2007 ; schwindel and mcnaughton 2011 ) , a similar compression effect may also be seen in primary visual cortex ( ji and wilson 2007 ) .
internally generated spatiotemporal patterns have also been observed in the rat hippocampus during the delay period of a memory task , showing that the emergence of consistent pattern of activity may be a way to maintain important information during a delay in a task ( pastalkova et al .
2008 ) , i.e. patterns with precise relative phases of the spikes of neurons participating to a collective oscillation , or precise phases of spikes relatively to the ongoing oscillation .
first experimental evidence of the importance of spike phases in neural coding was observed in experiments on theta phase precession in rat s place cells ( okeefe and recce 1993 ; okeefe and burgess 2005 ) , showing that spike phase is correlated with rat s position .
recently , the functional role of oscillations in the hippocampal - entorinal cortex circuit for path - integration has been deeply investigated ( okeefe and recce 1993 ; lengyel et  al .
2006 ) , showing that place cells and grid cells form a map in which precise phase relationship among units plays a central role . in particular it has been shown ( burgess et al . 2007 ; blair et al . 2008 ; welday et al .
2011 ) that both spatial tuning and phase - precession properties of place cells can arise when one has interference among oscillatory cells with precise phase relationship and velocity - modulated frequency . further evidence of phase coding comes from the experiments on spike - phase coding of natural stimuli in auditory and visual primary cortex ( montemurro et al . 2008 ; kayser et al .
2009 ) , and from experiments on short - term memory of multiple objects in prefrontal cortices of monkeys ( siegel et al .
these experimental works support the hypothesis that collective oscillations may underlie a phase dependent neural coding and an associative memory behavior which is able to recognize the phase coded patterns .
the importance of precise timing relationships among neurons , which may carry information to be stored , is supported also by the evidence that precise timing of few milliseconds is able to change the sign of synaptic plasticity ( markram et al .
the dependence of synaptic modification on the precise timing and order of pre- and post - synaptic spiking has been demonstrated in a variety of neural circuits of different species .
many experiments show that a synapse can be potentiated or depressed depending on the precise relative timing of the pre- and post - synaptic spikes .
this timing dependence of magnitude and sign of plasticity , observed in several types of cortical ( markram et al .
2001 ) and hippocampal ( bi and poo 1998 ; magee and johnston 1997 ; sjostrom et al .
1998 ; bi and poo 2001 ) neurons , is usually termed spike timing dependent plasticity ( stdp ) .
the role of stdp has been investigated both in supervised learning framework ( legenstein et al . 2005 ) , in unsupervised framework in which repeating patterns are detected by downstream neurons ( masquelier et al .
2011 ) , cortical development ( song and abbot 2001 ) , generation of sequences ( fiete et al . 2010 ; verduzco- flores et al . 2011 ) and polychronous activity ( izhikevich 2006 ) , and in an associative memory framework with binary units ( scarpetta et al .
however , this is the first time that this learning rule has been used to make a if network to work as associative memory for phase - coded patterns of spike , each of which becomes a dynamic attractor of the network .
notably , in a phase coded pattern not only the order of activation matters , but the precise spike timing intervals between units .
we therefore present a possibility to build a circuit with stable phase relationships between the spikes of a population of if neurons , in a robust way with respect to noise and changes of frequency .
the first important result of the paper is the measurement of the storage capacity of the model , i.e. the maximum number of distinct spatiotemporal patterns that can be stored and selectively retrieved , since it has never been computed in a spiking model for spatiotemporal patterns .
. 2001 , and references therein ) have focused on storage capacity of binary model with static binary patterns ( hopfield 1982 ) , and much efforts have been done to use more biophysical models and patterns ( gerstner et al .
1996 , 1993 ; hopfield 1995 ; de almeida et al . 2007 ; amit and treves 1989 ; battaglia and treves 1998 ; anishchenko and treves 2006 ; borisyuk and hoppensteadt 1998 ; leibold and kempter 2006 ; memmesheimer and timme 2006 ; olmi et al .
2010 ; scarpetta et al . 2002 , 2011b ) , but , up to our knowledge , without any calculation of the storage capacity of spatiotemporal patterns in if spiking models .
notably , by introducing an order - parameter which measures the overlap between phase coded spike trains , we are able quantitatively measure of the overlap between the stored pattern and the replay activity , and to compute the storage capacity as a function of the model parameters . another important result is the study of the different regimes observed by changing the excitability parameters of the network . in particular , we find that near the region of the parameter space where the network tends to become unresponsive and silent there is a regime in which the network responds selectively to cue presentation with a short transient replay of the phase - coded pattern .
differently , in the region of higher excitability , the patterns are replayed persistently and selectively , and eventually with more then one spike per cycle .
the paper is organized as follows : section 2 introduces the leaky - integrate - and - fire ( if ) neuronal model ; section 3 describes the stdp learning rule used to design the connections ; in section 4 we study the emergence of collective dynamics and introduce an order parameter to measure the overlap between the collective dynamics and the stored phase coded patterns ; section 5 reports on the storage capacity of the network , i.e. the maximum number of patterns that can be stored and selectively retrieved in the network ; the parameter space and the different working regions are also investigated in section 5 ; in section 6 we study the robustness of the retrieval dynamics wrt noise and heterogeneity ; section 7 reports on the implication of this model in the framework of oscillatory interference model of path - integration ; summary and discussion are outlined in section 8 .
we consider a recurrent neural network with n(n    1 ) possible connections jij , where n is the number of neural units .
the connections jij are designed during the learning mode , when the connections change their efficacy according to a learning rule inspired to the stdp .
after the learning stage , the connections values are frozen , and the collective dynamics is studied .
this distinction in two stages , plastic connection in the learning mode and frozen connections in the dynamics mode , is a useful framework to simplify the analysis .
it also finds some neurophysiological motivations in the effects of neuromodulators , such as dopamine and acetylcholine ( hasselmo 1993 , 1999 ) , which regulate excitability and plasticity .
the single neuron model is a leaky integrate - and - fire ( if ) ( gerstner and kistler 2002 ) .
this simple choice , with few parameters for each neuron , is suitable to study the emergence of collective dynamics and the diverse regimes of the dynamics , instead of focusing on the complexity of the neuronal internal structure .
we use the spike response model ( srm ) formulation ( gerstner and kistler 2002 ; gerstner et al .
1993 ) of the if model , which allows us to use an event - driven programming and makes the numerical simulations faster with respect to a differential equation formulation . in this picture
, the postsynaptic membrane potential is given by : 
 1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ h_i(t)=\sum\limits_{j } j_{ij } \sum\limits _ { { \hat t}_j > { \hat t}_i } \epsilon\left(t-{\hat
t}_j\right ) , \label{if } $ $ \end{document}where jij are the synaptic connections , (t ) describes the response kernel to incoming spikes on neuron i , and the sum over \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\hat t}_j$\end{document } runs over all presynaptic firing times following the last spike of neuron i. namely , each presynaptic spike j , with arrival time \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\hat{t}_j$\end{document } , is supposed to add to the membrane potential a postsynaptic potential of the form \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$j_{ij } \epsilon(t-{\hat t}_j)$\end{document } , where 
 2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{lll } \epsilon\left(t-{\hat t}_j\right)\nonumber\\ \quad = k \left[\ ! \exp\left(-\frac{t-{\hat t}_j}{\tau_m}\right ) - \exp\left(-\frac{t-{\hat
t}_j}{\tau_s}\right ) \!\right ] \theta\left(t-{\hat t}_j\right ) \label{tre } \end{array}$$\end{document } where m is the membrane time constant ( here 10 ms ) , s is the synapse time constant ( here 5 ms ) ,  is the heaviside step function , and k is a multiplicative constant chosen so that the maximum value of the kernel is 1 .
the sign of the synaptic connection jij sets the sign of the postsynaptic potential s change , so there s inhibition for negative jij and excitation for positive jij . when the membrane potential hi(t ) exceeds the spiking threshold \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^i$\end{document } , a spike is scheduled , and the membrane potential is reset to the resting value zero .
we use the same threshold th for all the units , except in section 6 where different values \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^i$\end{document } are used and the robustness w.r.t .
clearly the spiking threshold th of the neurons is related to the excitability of the network , an increase of the value of th is also equivalent to a decrease of k , the size of the unitary postsynaptic potential , or , equivalently to a global decrease in the scaling factor of synaptic connections jij .
numerical simulations of this dynamics are performed for a network with p stored patterns , where connections jij are determined via a learning rule described in the next paragraph .
we found that a few number of spikes , given a in proper time order , are able to selectively induce the emergence of a persistent collective spatiotemporal pattern , which replays one of the stored pattern ( see section 4 ) .
( 2007 ) , the average change in the connection jij , occurring in the time interval [    tlearn,0 ] due to periodic spike trains of period t , with tlearn  >  >  t , was formulated as follows : 
 3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \delta j_{ij } = \frac{t}{t_{\rm learn } } \int\limits_{-t_{\rm learn}}^{0}d t \int\limits_{-t_{\rm learn}}^{0}d t^\prime \ , x_i(t ) a(t - t^\prime ) x_j(t^\prime ) \label{lr } $ $ \end{document}where t / tlearn is a normalization factor , xj(t ) is the activity of the pre - synaptic neuron at time t , and xi(t ) the activity of the post - synaptic one .
it means that the probability for unit i to have a spike in the interval ( t , t  +  t ) is proportional to xi(t)t in the limit t0 .
the learning window a(  =  t    t ) is the measure of the strength of synaptic change when a time delay  occurs between pre and post - synaptic activity . to model the experimental results of stdp in hippocampal neurons , the learning window a( ) should be an asymmetric function of  , mainly positive ( ltp ) for   >  0 and mainly negative ( ltd ) for   <  0 .
equation ( 3 ) holds for activity pattern x(t ) which represents instantaneous firing rate , and is suitable to use in analog rate models ( scarpetta et al .
2001 , 2002 , 2008 ; yoshioka et al . 2007 ; scarpetta and marinaro 2005 ) and spin network models ( scarpetta et al . 2010 , 2011b ) .
differently , here , being interested in spiking neurons , the patterns to be stored are defined as precise periodic sequence of spikes , i.e. spike - phase coded patterns .
namely , activity of the neuron j is a spike train at times \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t^\mu_j$\end{document } , 
 4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_j^\mu(t)=\sum\limits_n \delta(t- ( t^\mu_j+ n t^\mu ) ) , \label{spiketr1 } $ $ \end{document}where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t^\mu_j+n t^\mu$\end{document } is the set of spikes times of unit j in the pattern  with period t , and frequency   =  1/t .
( 3 ) , the change in the connections jij due to the learning of the pattern  when the time duration of the learning process tlearn is longer then a single period t , is simply given by 
 5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ j_{ij}^\mu=\sum\limits_{n=-\infty}^{\infty } a(t^\mu_j - t^\mu_i+ n t^\mu ) .
\label{lr2 } $ $ \end{document}the window a( ) , shown in fig .  1 ,
is given by 
 6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ a(\tau ) = \left \ { \begin{array}{ll } a_p e^{-\tau / t_p } - a_{\tiny d } e^{-\eta \tau / t_p } & \quad \textrm{for}\quad \tau>0 \\ a_p
/ t_{\tiny d } } & \quad \textrm{for}\quad \tau<0 , \end{array}\right . \label{at } $ $ \end{document}with the same parameters used in abarbanel et al .
( 2002 ) to fit the experimental data of bi and poo ( 1998 ) , namely ap = [1/tp + \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\eta / t_d]^{-1}$\end{document } , ad = [/tp + \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$1/t_d]^{-1}$\end{document } , with tp  =  10.2 ms , td  =  28.6 ms ,
this function satisfies the balance condition \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\int_{-\infty}^\infty a(\tau ) d\tau = 0$\end{document}. notably , when a( ) is used in eq .
( 5 ) to learn phase - coded patterns with uniformly distributed phases , then the property \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\int a(\tau ) dt = 0$\end{document } assures that in the connection matrix the summed excitation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(1/n)\sum_{i , j_{ij}>0 } j _
{ ij}$\end{document } and the summed inhibition \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(1/n)\sum_{i , j_{ij}<0 } j _ { ij}$\end{document } are equal in the thermodynamic limit , and therefore it assures a balance between excitation and inhibition . 
 fig . 
1(a ) plot of the learning window a( ) used in the learning rule ( see eqs .
( 3 ) , ( 5 ) and ( 6 ) ) to model stdp effects .
the parameters of the function a( ) ( eq . ( 6 ) ) are determined by fitting the experimental data reported in bi and poo ( 1998 ) ( a ) plot of the learning window a( ) used in the learning rule ( see eqs .
( 3 ) , ( 5 ) and ( 6 ) ) to model stdp effects .
the parameters of the function a( ) ( eq . ( 6 ) ) are determined by fitting the experimental data reported in bi and poo ( 1998 ) writing eqs . 
( 35 ) , implicitly we have assumed that , with periodic phase - coded spike trains used to induce plasticity , the effects of all separate spike pairs sum linearly , each weighted by the same stdp window reported in fig .  1 .
1 are indeed typically measured by giving an order of 100 pairs of spikes repeatly , with fixed phase relationship , and fixed frequency in a proper range .
however , in different situations , for instance if the frequency is too low or to high ( sjostrom et al .
2001 ) , or in case of few spike pairs ( wittenberg and wang 2006 ) , the timing dependence of plasticity is not well described by the bidirectional window used here , and a more detailed model is needed to account for integration of spike pairs when arbitrary trains are used ( see shouval et al .
the spikes patterns used in this work are periodic spatiotemporal sequences , made up of one spike per cycle and each of which has a phase \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi^\mu_j$\end{document } randomly chosen from a uniform distribution in [ 0,2 ) . in each pattern
, information is coded in the precise time delay between spikes of unit i and unit j , which corresponds to a precise phase relationship among units i and j. a spatiotemporal pattern represented in this way is often called phase coded pattern . pattern s information is coded in the spiking phases which , in turn , shape the synaptic connectivity responsible of the emerging dynamics and the memory formation .
the set of timing of spikes of unit j can be defined as \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t^\mu_j + n t^\mu = ( \phi^\mu_j)/(2\pi \nu^\mu ) + n/\nu^\mu$\end{document } , where  is the oscillation frequency of the neurons .
thus , each pattern  is represented through the frequency  and the specific phases of spike \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi^\mu_j$\end{document } of the neurons j  =  1, .. ,n .
the change in the connection jij provided by the learning of pattern  is given by 
 7\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{rll } j_{ij}^\mu&=&\displaystyle\sum\limits_{n } a(t^\mu_j - t^\mu_i+n t^\mu)\nonumber\\ & = & \sum\limits_n a\left(\frac{\phi^\mu_j } { 2\pi \nu^\mu}-\frac{\phi^\mu_i}{2\pi\nu^\mu } + n/\nu^\mu \right ) .
\label{conn } \end{array}$$\end{document } when multiple phase coded patterns are stored , the learned connections are simply the sum of the contributions from individual patterns , namely 
 8\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ j_{ij}=\sum\limits_{\mu=1}^p j^{\mu}_{ij}. \label{connp } $ $ \end{document}note that ring - like topology with strong unidirectional connections is formed only in the case p=1 , when a single pattern is stored . when multiple patterns are stored in the same connectivity , with phases of one pattern uncorrelated with the others , bidirectional connections are possible , and the more the stored patterns , the less the ring - like is the connectivity . even in the cases
when the connectivity is not ring - like the network is still able to retrieve each of the p stored patterns in a proper range of threshold values ( see storage capacity in section 5 ) .
we study a recurrent network with n leaky integrate and fire units , with connections fixed to the values calculated in eqs .
the results show that , within a well specified range of parameters , our if network is able to work as an associative memory for spike - phase patterns . in order to check if the network is able to retrieve selectively each of the stored patterns , we give an initial signal , made up of m    n spikes , taken from the stored pattern  , and we check if this initial short cue is able to selectively trigger a collective sustained activity that is the replay of the same stored pattern  , i.e. checking if the sustained activity has spikes aligned to the phases \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi_i^\mu$\end{document } of pattern . an example of successful selective retrieval process is shown in fig . 
2 where , depending on the partial cue presented to the network , a different collective activity emerges with the phases of the firing neurons which resemble one or another of the stored patterns . 

2examples of selective successful retrieval ( ( c ) , ( d ) ) of two stored patterns ( ( a ) , ( b ) ) .
the raster plot of 50 units ( randomly chosen ) are shown sorted on the vertical axis according to increasing values of phase \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi_i^1$\end{document } of the first stored pattern   =  1 .
the network has n  =  3000 if neurons , th  =  70 and connections given by eqs .
( 7 ) and ( 8) with p  =  5 stored patterns at   =  3  hz .
two of the stored patterns used during the learning mode are shown in ( a ) , ( b ) .
the dynamics emerging after a short train of m  =  n/10 spikes with phases similar to the pattern shown in ( a ) and ( b ) , is shown in ( c ) and ( d ) respectively .
the dynamics of the network , after a transient , is periodic of period t. the spikes which belong to the trigger are shown in pink in ( c ) and ( d ) , the other different colors represent the value of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t_i / t\!\mod 4$\end{document } , where ti is the time of the spike of the unit i during the emerging spontaneous dynamics .
figure ( c ) shows that when the network dynamic is stimulated by a partial cue of pattern   =  1 , the neurons oscillate with phase alignments resembling pattern   =  1 , but at different frequency .
otherwise , in ( d ) , when the partial cue is taken from pattern   =  2 , the neurons phase relationships , even if periodic , are uncorrelated with pattern   =  1 , and recall the phase of pattern   =  2 examples of selective successful retrieval ( ( c ) , ( d ) ) of two stored patterns ( ( a ) , ( b ) ) .
the raster plot of 50 units ( randomly chosen ) are shown sorted on the vertical axis according to increasing values of phase \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi_i^1$\end{document } of the first stored pattern   =  1 .
the network has n  =  3000 if neurons , th  =  70 and connections given by eqs .
( 7 ) and ( 8) with p  =  5 stored patterns at   =  3  hz .
two of the stored patterns used during the learning mode are shown in ( a ) , ( b ) .
the dynamics emerging after a short train of m  =  n/10 spikes with phases similar to the pattern shown in ( a ) and ( b ) , is shown in ( c ) and ( d ) respectively .
the dynamics of the network , after a transient , is periodic of period t. the spikes which belong to the trigger are shown in pink in ( c ) and ( d ) , the other different colors represent the value of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t_i / t\!\mod 4$\end{document } , where ti is the time of the spike of the unit i during the emerging spontaneous dynamics .
figure ( c ) shows that when the network dynamic is stimulated by a partial cue of pattern   =  1 , the neurons oscillate with phase alignments resembling pattern   =  1 , but at different frequency .
otherwise , in ( d ) , when the partial cue is taken from pattern   =  2 , the neurons phase relationships , even if periodic , are uncorrelated with pattern   =  1 , and recall the phase of pattern   =  2 in this work the cue is a stimulation with m spikes , with m  =  n/10 , at times \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t^\mu_i= t_{\rm stim } \phi_i^\mu$\end{document } , 0  < 
i  <  m , with tstim  =  50 ms . in the example shown in fig . 
2 the short stimulation ( which lasts less then 5 ms , shown in pink in all the figures ) has the effect to selectively trigger the sustained replay of pattern . note that the retrieval dynamics has the same phase relationship among units than the stored pattern , but the replay may happen on a time scale different from the scale used to store the pattern , and the collective spontaneous dynamics is a time compressed ( or dilated ) replay of the stored pattern . indeed
, the period of the collective periodic pattern which emerges during retrieval stage may be different then the period of the periodic pattern used in the learning stage . in the example of fig .  2 the time scale of the retrieval dynamics ( fig . 
2(c ) , ( d ) ) is faster then the time scale used to learn the patterns ( fig . 
in the following we will study the factors affecting the time scale during retrieval , given the time scale of the pattern used during learning .
clearly , regions of the parameter space in which the network is unable to retrieve selectively the patterns also exist . in these regions
the retrieval dynamics may correspond to a mixture of patterns or a spurious state , i.e. a state which is not correlated with any of the stored patterns because the number of stored patterns exceeded the storage capacity of the network . as discussed below , the storage capacity , defined as the maximum number of encoded and successfully retrieved patterns , depends on the frequency used during the learning stage ( which affects connectivity ) , and on the spiking threshold of the units ( which affects excitability and network dynamics ) .
3(b ) the network has too low excitability and the response is not persistent , while in fig . 
3(a ) the emerging dynamics is not correlated with any of the stored patterns . 
 fig . 
3example of neural response in two case of failure of retrieval . a spurious state emerge in ( a ) , while a short transient response emerges in the case shown in ( b ) .
n  =  3000 and   =  3  hz as in previous picture , while the values th and p are th  =  10 , p  =  5 in ( a ) and th  =  95 , p  =  5 in ( b ) .
the dynamics emerging after a short train of m  =  n/10 spikes with phases equal to the stored pattern ( pattern shown in fig . 
for clarity , the raster plot of only 50 ( randomly chosen ) units are shown , sorted according to increasing value of phase i of the stored pattern example of neural response in two case of failure of retrieval . a spurious state emerge in ( a ) , while a short transient response emerges in the case shown in ( b ) .
n  =  3000 and   =  3  hz as in previous picture , while the values th and p are th  =  10 , p  =  5 in ( a ) and th  =  95 , p  =  5 in ( b ) .
the dynamics emerging after a short train of m  =  n/10 spikes with phases equal to the stored pattern ( pattern shown in fig . 
2(a ) ) , is not a self - sustained retrieval of the pattern . for clarity ,
the raster plot of only 50 ( randomly chosen ) units are shown , sorted according to increasing value of phase i of the stored pattern to measure quantitatively the success of the retrieval , in analogy with the hopfield model , we introduce an order parameter , which estimates the overlap between the network collective activity during the spontaneous dynamics and the stored phase - coded pattern .
this quantity is 1 when the phases j of neurons j coincides with the stored phases \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi_j^\mu$\end{document } , and is close to zero when the phases are uncorrelated with the stored ones .
therefore , we consider the following dot product |m(t)|  =  <  (t)|  > where  is the vector having components \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$ e^{i \phi_j^\mu}$\end{document } , namely : 
 9\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ |m^\mu(t)|= \left|\frac{1}{n}\sum\limits_{\stackrel{j=1,\ldots , n}{\scriptscriptstyle t - t^*<t_j^*<t } } e^{-i 2 \pi t_j^*/t^ * } e^{i \phi_j^\mu}\right| \label{nn } $ $ \end{document}where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t_j^*$\end{document } is the spike timing of neuron j during the spontaneous dynamics , and t is an estimation of the period of the collective spontaneous periodic dynamics .
( 9 ) is equal to 1 when the phase - coded pattern is perfectly retrieved ( i.e. same sequence and phase relationships among spikes , even though on a different time scale ) , while is of order \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\simeq 1/\sqrt{n}$\end{document } when phases of spikes are uncorrelated to the stored phases .
the order parameter |m| allow us to measure the network storage capacity in the space of parameters th and . note that the value of m(t ) between two periodic spike trains measures the similarity in the sequence of spiking neurons and in the phase lag between the spikes , being invariant by a simple change in time scale .
this is a suitable choice especially when the replay of a spatio temporal pattern has to be detected independently from the compression of the time scale .
note that if we have a spike train that is not periodic , we can not define the period , however we can define the order parameter eq . 
. this can be useful in the case when one looks for a short replay hidden in a not - periodic spike train , such in many experimental situations .
the value of m(t ) after a transient converges to a stable value which is close to one when pattern  is retrieved ( for example in fig . 
2(c ) at large times m  =  1 , and m  =  0.01 ) while m(t ) is of order  \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\simeq 1/\sqrt{n}$\end{document } for all  in the case of failure of retrieval .
3(a ) m(t ) after the transient has values in the range 0.010.02 for all  because the emerging dynamics is a spurious state not correlated with any of the stored phase - patterns , while m(t ) is zero in fig . 
, the storage capacity of the model is analyzed considering the maximum number of patterns that the network is able to store and selectively recall .
in particular we investigate the role of two model parameters : the frequency of the stored patterns  , and the spiking threshold th affecting the excitability of the network .
numerical simulations of the if network with n  =  3000 neurons were performed by systematically changing the value of the spiking threshold th , the connections jij , and for different number of patterns p and frequency . here we propose results for a unique value of the spiking threshold th for all neurons , however the behavior is also robust with respect to a variability in the threshold values among neurons , as reported in the next section .
network storage capacity is defined as c  =  pmax / n , where n is the number of neurons and pmax is the maximum number of patterns that can be stored and successfully retrieved with an overlap |m| larger than a certain value , which measures the degree of similarity .
given that in our simulations the overlap |m(t)| at large times has mostly two possible values , close to one ( success ) or close to \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$1/\sqrt n$\end{document } ( failure ) , we fixed the desired similarity value to 0.5 , since the whole storage capacity analysis is very robust with respect to this parameter ( since the transition between low values and high values of |m(t)| as a function of p is sharp ) .
patterns with random phases were extracted and used to define the network connections jij with the rule eq .
after the stimulation with a short train of m  =  n/10 spikes taken at times ti from the first pattern , the dynamics is simulated and the overlap defined in eq .
if the overlap |m(t)| , averaged over 50 runs , is greater than 0.5 at time \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t>\bar t$\end{document } ( where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\bar t= 600$\end{document } ms for all the simulations ) , then we consider the retrieval successful for that pattern .
the maximum value of p , for which the network is able to successfully replay each of the stored patterns , defines the storage capacity of the network . the storage capacity as a function of the spiking threshold th and storing frequency  is reported in fig . 
the largest capacity is achieved when the frequency of the stored patterns during learning is     8 hz and the spiking threshold of the units during retrieval is th    130 , which provides a capacity max  =  pmax / n  =  0.016 . 

4(a ) storage capacity in a network of n  =  3000 units , as a function of the spiking threshold th and oscillation frequency  of stored patterns .
the maximum number of patterns successful retrievable pmax is shown in color - coded legend , the value grows from pmax  =  0 ( dark blue ) to pmax  =  50 ( strong red ) .
( b ) the storage capacity pmax as a function of the frequency of stored patterns , once fixed the threshold th to the optimal value for each frequency ( a ) storage capacity in a network of n  =  3000 units , as a function of the spiking threshold th and oscillation frequency  of stored patterns .
the maximum number of patterns successful retrievable pmax is shown in color - coded legend , the value grows from pmax  =  0 ( dark blue ) to pmax  =  50 ( strong red ) .
( b ) the storage capacity pmax as a function of the frequency of stored patterns , once fixed the threshold th to the optimal value for each frequency in fig .
4(b ) we show the storage capacity pmax as a function of the frequency  once fixed the value of the threshold th to the optimal value , corresponding to highest capacity for each frequency .
the optimal storing frequencies and threshold values depend on the time constants of the model , such as the s,m of the if units and the temporal shape of the learning kernel a( ) , whereas different shapes of a( ) may subserve to different storing frequency ranges . in this work s,m and a( )
are set to the values described in section 2 , and the emergent collective dynamics is studied as a function of the other network parameters .
4(b ) shows that for the learning kernel a( ) used here , there is peak in the storage capacity around 8 hz , in the range 220 hz .
figure 4(a ) also proves that , for each stored frequency , a large interval of spiking threshold values th exists for which the network is still able to work properly as associative memory for phase - coded patterns .
the associative memory properties as a function of the spiking threshold are reported in fig .
5 , when the oscillation frequency of the patterns stored during learning is   =  3 hz .
5 corresponds to cases in which the retrieval is successful and the cue is able to selectively activate the self - sustained replay of the stored pattern ( with an order parameter m larger than 0.5 ) .
when spiking threshold changes in the range 10  <  th  <  90 the storage capacity changes between pmax / n  =  1/3000 and pmax / n  =  29/3000 . 
 fig . 
5(a ) storage capacity at   =  3 hz : the region of successful retrieval as a function of spiking threshold and number of patterns is marked in green .
the region with persistent activity not correlated with any of the stored pattern is marked in red ( spurious states ) , and the region in which the network responds with only a short transient and then becomes silent is marked in blue ( see examples in fig . 
( b ) the probability that the size stot of the network response , measured as the number of the spikes that follow the cue stimulation , is larger than n , with n  =  0,n/2,n , is shown as a function of spiking threshold th , in a network with   =  3  hz and p  =  1 . as always in this paper the number of units is n  =  3000 .
the figure shows that near \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^{\rm crit } \simeq 90$\end{document } there s a transition from a region of persistent replay to a region of silence ( a ) storage capacity at   =  3 hz : the region of successful retrieval as a function of spiking threshold and number of patterns is marked in green .
the region with persistent activity not correlated with any of the stored pattern is marked in red ( spurious states ) , and the region in which the network responds with only a short transient and then becomes silent is marked in blue ( see examples in fig . 
( b ) the probability that the size stot of the network response , measured as the number of the spikes that follow the cue stimulation , is larger than n , with n  =  0,n/2,n , is shown as a function of spiking threshold th , in a network with   =  3  hz and p  =  1 .
as always in this paper the number of units is n  =  3000 .
the figure shows that near \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^{\rm crit } \simeq 90$\end{document } there s a transition from a region of persistent replay to a region of silence outside the green region the number of patterns exceeds the storage capacity and the retrieval fails .
when the number p of patterns exceeds the storage capacity pmax , the network responds with a self - sustained activity that is not correlated with any of the stored patterns , i.e. a spurious state . in this regime ,
5(a ) , the order parameter m(t ) is of order \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$1/\sqrt n$\end{document } for all the stored patterns ( see also raster plot in fig . 
on the contrary , in the high th regime , the network tends to become silent and unresponsive .
indeed , in the region marked with blue color , the network responds to the initial cue stimulation with a short transient and then became silent . in this case
the value of m(t ) is zero because there is no self - sustained activity at time \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$t>\bar t$\end{document } , meaning that the stored attractors become unstable when th is too high ( see raster plot in fig . 
for values of the threshold greater than \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^{\rm crit}$\end{document } , independently from p , the network activity is never persistent , as reported in fig .
5(a ) where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^{\rm crit}=90$\end{document}. at thresholds close to this critical value the network responds with a transient activity that is a short replay of the stored pattern , but not a persistent replay .
the size stot of the network response , measured as the number of spikes that follow the cue stimulation , is reported in fig .
5(b ) as a function of th , for a network with   =  3 hz , n  =  3000 and p  =  1 . in
we focused on the dependence of the frequency of collective oscillations during replay on the model parameters . figure 6(a ) shows the collective frequency of replay as a function of the frequency  of the patterns stored in the learning stage with n  =  3000 .
the red dots in the figure refer to the frequencies of oscillations observed during retrieval at the optimal spiking threshold ( where the maximum storage capacity occurs ) , while the bar indicates the available range of frequencies of replay , accessible through a change in the spiking threshold .
important to note is that , in most of the cases , the frequency of the stored pattern and the collective replay frequency do not coincide , since the pattern is replayed compressed ( or dilated ) in time , on a time scale dependent on the network parameters .
we observe that for the chosen parameters m , s of the network , and the given shape of a( ) , the replay occurs on a compressed time scale for all stored patterns of frequency lower then 25 hz , while the two time scale coincide when     25 hz .
the dependence of the frequency of the collective oscillations on the spiking threshold is shown in fig .
besides , for low stored frequencies ( 14 hz ) the frequency of the replay is very sensitive to the threshold value , changing from 6 hz at high spiking threshold to 30 hz at low threshold . 

6(a ) frequency of the collective dynamics during replay as a function of the frequency of stored patterns in the network with n  =  3000 units .
( b )  frequencies of the collective dynamics during replay as a function of the spiking threshold and for different stored frequencies ( see colors legend ) .
pattern is replayed on a time scale which becomes faster if we decrease threshold th for the most of the frequencies . the dependence is much stronger for     4 hz .
( c )  the number of spikes per cycle as a function of the spiking threshold th in networks with different frequency  of stored pattern .
( d )  frequency of the dynamics during replay as a function of ratio between spiking threshold and network size n. red dots are results for a network with n  =  10000 units , while blue squares are results for a network with n  =  1000 units .
size of the symbols refers to the stored frequency , small symbols ( on the top of the picture ) correspond to stored frequency   =  10  hz , medium size symbols correspond to stored frequency   =  3  hz , and large symbols ( bottom ) correspond to stored frequency   =  1  hz ( a ) frequency of the collective dynamics during replay as a function of the frequency of stored patterns in the network with n  =  3000 units .
( b )  frequencies of the collective dynamics during replay as a function of the spiking threshold and for different stored frequencies ( see colors legend ) .
pattern is replayed on a time scale which becomes faster if we decrease threshold th for the most of the frequencies . the dependence is much stronger for     4 hz .
( c )  the number of spikes per cycle as a function of the spiking threshold th in networks with different frequency  of stored pattern .
( d )  frequency of the dynamics during replay as a function of ratio between spiking threshold and network size n. red dots are results for a network with n  =  10000 units , while blue squares are results for a network with n  =  1000 units .
size of the symbols refers to the stored frequency , small symbols ( on the top of the picture ) correspond to stored frequency   =  10  hz , medium size symbols correspond to stored frequency   =  3  hz , and large symbols ( bottom ) correspond to stored frequency   =  1  hz we also investigate the frequency s dependence on network size n. in fig .
6(d ) red dots are results for a network with n  =  10000 units , while blue squares are results for a network with n  =  1000 units .
if what counts is only the time lag between the single units consecutive in the sequence , then one expects that result with 1 hz stored at n  =  10000 would be similar to 10 hz stored at n  =  1000 , and this is not the case .
we see that when we use a storage frequency equal to 10 hz ( which corresponds to different time lag between cells depending on n ) , then the oscillation frequency during replay is around 30 hz in both networks ( both n  =  1000 , and n  =  10000 ) , while , on the other hand , if we have a storage frequency equal to 1 hz the oscillation frequency during replay may span a large range ( 5 hz25 hz ) in both networks .
figure 6(d ) also shows that frequency of replay depends on the ratio between spiking threshold and network size n , and that the high sensitivity on spiking threshold value holds , when stored frequency is low ( 14 hz ) , also at different values of the network size .
this open the possibility to govern the oscillation frequency of the collective replay activity via neuromodulators which change the excitability and therefore the spiking threshold of the neurons . since in our model ( see eqs .  ( 1 ) and ( 2 ) ) a change of the threshold is equivalent to a change in the scale factor of all synaptic connections , a similar effect might be achieved also by simply driving the cells more due to increased synaptic input .
importantly , the sensitivity of collective oscillation frequency on spiking threshold is not a sensitivity of the single unit but of the collective behavior , since , as discussed in section 6 , if we change the spiking threshold of few units the collective rhythm is still unique for the whole population .
the replay frequency depends on the average threshold among units , but all the units have the same oscillation frequency during replay . moreover , for networks with     10 hz ,
whose replay frequency does not considerably change with spiking threshold , the replay dynamics is still affected by the spiking threshold .
indeed , in this case , the number of spikes per cycle increases with lowering of the spiking threshold .
the raster plots show the same pattern replayed in three networks having different values of the spiking threshold th : a burst of activity takes place within each cycle , with phases aligned with the pattern , with a number of spikes per cycle dependent on the value of th .
6(c ) where the number of spikes per cycle is reported as a function of spiking threshold , at different values of stored frequencies .
therefore , by lowering the spiking threshold the replay activity occurs with more than one spike per cycle , or on a faster time scale ( see fig . 
7modulation of spiking threshold changes the number of spikes per cycle , keeping preserved the phase relationship among units .
recall of the pattern   =  1 for networks of n  =  3000  units ,   =  20 hz and different values of spiking threshold th  =  80,65,40 is shown respectively in ( a)(c ) .
depending on the value of the spiking threshold th , the phase - coded pattern is replayed with a different number of spikes per cycle .
spike of the cue stimulation are shown in pink , while the response of the network in black . for clarity , the raster plot of only 50 ( randomly chosen ) units are shown , sorted according to increasing value of phase i of the stored pattern modulation of spiking threshold changes the number of spikes per cycle , keeping preserved the phase relationship among units .
recall of the pattern   =  1 for networks of n  =  3000  units ,   =  20 hz and different values of spiking threshold th  =  80,65,40 is shown respectively in ( a)(c ) .
depending on the value of the spiking threshold th , the phase - coded pattern is replayed with a different number of spikes per cycle .
spike of the cue stimulation are shown in pink , while the response of the network in black . for clarity ,
the raster plot of only 50 ( randomly chosen ) units are shown , sorted according to increasing value of phase i of the stored pattern the behavior of the output oscillation frequency suggests that a parameters region exists where the network always responds with one spike per cycle . in this region an increase of the excitability produces a growth of the frequency of oscillation up to a plateau value .
differently , for higher excitability the frequency does not increase while the number of spikes per cycle grows .
this means that the different frequencies , in addition to the information coded in the phase relationship , can code other information in relationship with the level of spiking threshold : at high frequency the threshold changes the number of spikes per cycle , while at low frequency the threshold changes the frequency of the collective oscillations during the replay . this open the possibility to have a coding scheme in which , while the phases encode pattern s information , a change in frequency or a change in rate in each cycle represents the strength and saliency of the retrieval or it may encode another variable ( lengyel and dayan 2007 ) .
the recall of the same phase - coded pattern with a different number of spikes per cycle is particularly interesting at the light of recent observations of huxter et al .
( 2003 ) in hippocampal place cells , showing the occurrence of the same phases with different rates .
the authors prove that the phase of firing and firing rate are dissociable and can represent two independent variables , e.g. the animal location within the place field and its speed of movement through the field .
notably , the recall of the same phase coded pattern with different frequencies of oscillation is also relevant and accords well with the need to have stable precise phase relationship among cells with frequency of oscillation modulated by parameters such as the speed of the animal ( blair et al .
the value of the frequency of collective activity during the replay clearly is related not only to the threshold and the stored frequency , but also to the shape of the learning window a( ) and on the two characteristic times of the model s , m .
a systematic study of the dependence of the replay time scale on the shape of stdp and the characteristic times of the neuron model has not yet done in a spiking model , however a dependence on the asymmetry of a(t ) has been analytically found in a simple model with analog neurons and a single characteristic time ( yoshioka et al .
while in the hopfield model the patterns are static , and information is coded in a binary pattern \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\bold s^\mu= s_1^\mu,\ldots , s_n^\mu$\end{document } , with \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$s_i^\mu \in \{\pm 1\}$\end{document } , here , in this study , the patterns are time dependent , and information is coded in the phase pattern \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\bold{\phi}^\mu=\phi_1^\mu, ... ,\phi_n^\mu$\end{document } with \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi_i^\mu \in [ 0,2\pi]$\end{document } , where the value \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi_i^\mu /(2\pi\nu^\mu)$\end{document } represents the time shift of the spike of unit i with respect to the collective rhythm , i.e the time delay among units .
however , as for the hopfield model , the patterns stored in the network are attractors of the dynamics , when p do not exceeds storage capacity , and the dynamics during the retrieval is robust with respect to noise .
input noise , i.e when a poissonian noise i(t ) is added to the postsynaptic potential hi(t ) given in eq . 
. the total postsynaptic potential of each neuron i is then given by 
 10\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ h_i(t)=\eta_i(t)+\sum\limits_{j}j_{ij}\sum\limits_{{\hat t}_j > { \hat t}_i } \epsilon\left(t-\hat{t_j}\right ) $ $ \end{document}where i(t ) is modelled as 
 11\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \eta_i(t)=j_{\text{noise}}\sum\limits_{\hat{t}_{\text{noise}}>\hat{t}_i } \epsilon\left(t-\hat{t}_{\text{noise}}\right ) .
$ $ \end{document}the times \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\hat{t}_{\text{noise}}$\end{document } are randomly extracted for each neuron i , and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$j_{\text{noise}}$\end{document } are random strengths , extracted independently for each neuron i and time \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\hat{t}_{\text{noise}}$\end{document}. the intervals between times \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\hat{t}_{\text{noise}}$\end{document } are extracted from a poissonian distribution \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p(\delta t ) \propto
e^{-\delta t/(n\tau_{\text{noise}})}$\end{document } , while the strength \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$j_{\text{noise}}$\end{document } is extracted from a gaussian distribution with mean \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\bar j}_{\text{noise}}$\end{document } and standard deviation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\sigma(j_{\text{noise}})$\end{document}. the network dynamics during the retrieval of a pattern in presence of noise is shown in fig . 8 with different levels of noise ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				
\usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\tau_{\text{noise}}=10$\end{document } ms , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\bar j=0 $ \end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\sigma(j_{\text{noise}})=0,10,20,30$\end{document } in a , b , c and d , respectively ) .
results show that when the noise is not able to move the dynamics out of the basin of attraction , the errors do not sum up , and the phase relationship is preserved over time ( see fig .  8(a)(c ) ) . if the input noise is very high , as in the example of fig .
raster plots show that , when the pattern retrieval is triggered , network s spikes continue to have phase alignments resembling the pattern even in presence of noise .
errors do not sum up until the system is in the basin of attraction of the phase - coded pattern , as in ( a)(c ) .
different levels of noise are used in ( a)(d ) ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\sigma(j_{\text{noise}})=0,10,20,30$\end{document } respectively ) , and pattern is triggered with m  =  n/10 as in previous cases . only in ( d )
the level of noise is too high and the system goes out of the basin of attraction .
( e ) ,  ( f )  for comparison , the dynamics , when the retrieval is not triggered ( m  =  0 ) , is shown in subplot ( ( e ) and ( f ) ) in presence of the same noise used in ( c ) and ( d ) .
figure ( e ) shows that the noise used in ( c ) usually affects strongly the dynamics of the network , however if the collective oscillation is retrieved the system is robust wrt noise .
thresholds in all figures are th  =  80 , n  =  3000 , and synaptic connections jij are build learning p  =  2 phase - patterns at   =  3 hz ( a)(d ) robustness wrt noise .
raster plots show that , when the pattern retrieval is triggered , network s spikes continue to have phase alignments resembling the pattern even in presence of noise .
errors do not sum up until the system is in the basin of attraction of the phase - coded pattern , as in ( a)(c ) .
different levels of noise are used in ( a)(d ) ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\sigma(j_{\text{noise}})=0,10,20,30$\end{document } respectively ) , and pattern is triggered with m  =  n/10 as in previous cases . only in (
d ) the level of noise is too high and the system goes out of the basin of attraction .
( e ) ,  ( f )  for comparison , the dynamics , when the retrieval is not triggered ( m  =  0 ) , is shown in subplot ( ( e ) and ( f ) ) in presence of the same noise used in ( c ) and ( d ) .
figure ( e ) shows that the noise used in ( c ) usually affects strongly the dynamics of the network , however if the collective oscillation is retrieved the system is robust wrt noise .
thresholds in all figures are th  =  80 , n  =  3000 , and synaptic connections jij are build learning p  =  2 phase - patterns at   =  3 hz in order to see the effects of input noise level used in fig 8(c ) , ( d ) , we report in fig .
8(e ) , ( f ) the network dynamics when the pattern retrieval is not initiated ( m  =  0 ) . in particular , fig .
8(c ) is strong enough to generate spontaneous random activity in absence of the initial triggering , but is not sufficient to destroy the attractive dynamics during a successful retrieval . as in the hopfield model , errors
do not sum up and the dynamics spontaneously goes back to the retrieved phase - coded pattern for all the perturbations that leave the system inside the basin of attraction .
this analysis can be carried out by using a different value \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^i$\end{document } of spiking threshold for each neuron i : 
 12\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \theta_{\rm th}^i= ( 1 + z
\zeta_i ) \theta_{\rm th } $ $ \end{document}where i is a random number extracted from a uniform distribution in [    1,1 ] , and z is the degree of heterogeneity . even with high degree of heterogeneity , the emergence of the retrieval collective dynamics forces all neurons to have exactly the same frequency of oscillation and to keep a precise phase relationship , in a very robust manner .
figure  9(a ) , ( b ) shows the dynamics with threshold heterogeneity z  =  0.2 , 0.5 respectively , while the remaining parameters are set to the values of fig .
raster plots show that , when the pattern retrieval is triggered , units participate to the network collective oscillation , showing all the same frequency and phase alignments resembling the pattern , even in presence of threshold values heterogeneity among units .
spiking thresholds of neuron i are distributed according to \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^i=\theta_{\rm th } ( 1+z \zeta)$\end{document } with average th  =  80 and z  =  0.2,0.5 in ( a ) and ( b ) respectively .
8(a ) ( n  =  3000 , m  =  n/10 , and synaptic connections jij are build learning p  =  2 phase - patterns at   =  3 hz ) robustness wrt heterogeneity of spiking threshold values .
raster plots show that , when the pattern retrieval is triggered , units participate to the network collective oscillation , showing all the same frequency and phase alignments resembling the pattern , even in presence of threshold values heterogeneity among units . spiking thresholds of neuron
i are distributed according to \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^i=\theta_{\rm th } ( 1+z \zeta)$\end{document } with average th  =  80 and z  =  0.2,0.5 in ( a ) and ( b ) respectively .
8(a ) ( n  =  3000 , m  =  n/10 , and synaptic connections jij are build learning p  =  2 phase - patterns at   =  3 hz ) the above analysis shows that a unique collective frequency emerges , which is the frequency corresponding to the mean value of the \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\theta_{\rm th}^i$\end{document}. this is also evident by comparing fig .
=  0.5 it can be seen an additional small phase shift proportional to the value of  of each neuron , but collective activity is preserved .
clearly if z is too high and threshold values are distributed out of the region of successful retrieval the network is unable to retrieve the pattern and failure happens , as already discussed in section 5 .
recently , path - integration system and the hippocampal - entorinal cortex circuit have been deeply investigated ( okeefe and recce 1993 ; okeefe and burgess 2005 ; burgess et al .
2007 ; jeewajee et al . 2008 ; blair et al . 2008 ; euston et al . 2007 ; geisler et al . 2007 ; mcnaughton et al . 2006 ) , showing that place cells and grid cells form a map in which precise phase relationship among units play a central role to generate the spatial tuning .
a number of models of the spatial firing properties of place and grid cells were offered .
generally two main categories are distinguished : models which focus on continuous attractor mechanisms and models which use interference between oscillators at dynamically modulated frequencies ; however , deeper computational principles may exist that unify the different cases of neural integration ( issa and zhang 2012 ) . in oscillatory interference models ( okeefe and burgess 2005 ;
2008 ; welday et al . 2011 ; zilli and hasselmo 2010 ) ( see also giocomo et al .
2011 , for a review ) the total synaptic input to a neuron ( such as a grid cell or a place cell ) is a weighted sum of the activities of n oscillatory inputs , whose oscillation frequency is modulated by the rat velocity and head direction .
grid and place cells , according to these models , derive their temporal and spatial properties simply by detecting synchrony among such velocity - modulated oscillatory inputs .
the oscillatory interference theory is one possible hypothesized mechanism of path integration and it has not been conclusively accepted or rejected .
it is supported by recent studies suggesting that the predicted velocity - modulated oscillators exist as theta cells ( interneurons found throughout the septo - hippocampal circuit ) whose inter - burst frequency shows a cosine modulation by running direction and a linear increase with running velocity ( welday et al .
such velocity - modulated oscillatory input was hypothesized to come from single oscillatory neurons ( okeefe and burgess 2005 ) , or networks such as subcortical  ring attractors  generating velocity - modulated theta oscillations ( blair et al .
the properties of modulation and stability of frequency , and stability of multiple phase relationships , make our circuit a possible mechanism to build the velocity - modulated oscillators of the oscillatory interference theory .
indeed our circuit has a collective oscillation frequency , which depends on the frequency stored in the connectivity matrix , and that can be modulated by changing the parameters such as th .
each neuron in the circuit has a phase determined by its position in that network , i.e. determined by the phase \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi_i^\mu$\end{document } of the stored phase - pattern .
if the parameter th is modulated by animal speed , then the collective oscillation frequency of the circuit is modulated by the animal s speed , while the neurons preserve stable phase relationship among them .
the other persistent - firing models ( blair et al . 2008 ; zilli and hasselmo 2010 ; welday et al .
2011 ) of the oscillators needed by the oscillatory interference theory suffer from problems related to robustness , as those encountered by the single - cell oscillatory models ( zilli et al .
2009 ) , due to the variability in the frequency of persisting spiking ( giocomo et al . 2011 ) .
indeed the oscillatory interference models impose strict constraints upon the dynamical properties of the velocity - modulated oscillatory inputs , which have to preserve robust velocity - modulated frequency and stable phase relationships among them on relevant time scales in a manner robust to noise ( many seconds , or dozens of theta cycle periods ) ( blair et al . 2008 ; welinder et al . 2008 ; zilli et al . 2009 ; giocomo et al .
notably , in our model , since connections jij among units in the circuit are fixed by the learning rule ( eq .  ( 7 ) ) , dynamics is a retrieval process and neurons preserve stable precise phase relationship among units and stable frequency even in noisy conditions , at least when the dynamics is in the basin of attraction of that phase relationship . even in the more recent spiking models ( zilli and hasselmo 2010 ; welday et al .
2011 ) of the oscillatory interference principle , in which many problems related to noise are solved , the heterogeneity of parameters of the cells which participate to the ring oscillator is not taken into account . here
we show that the circuit level interactions among units make the oscillation frequency and the phase - relationship of the system robust even with respect to heterogeneity of the spiking thresholds of the units ( see section 6 ) .
this robustness , due to the proposed coupling which forces all the units of the circuit to have exactly the same period of oscillation and to have precisely the same phase relationships of stored pattern , may be useful in all cases of sequence coding .
moreover , our circuit can be easily programmed to cycle in different phase orders , by storing more than one phase - pattern as attractors .
the circuit is a robust phase - shuffling ring oscillator , since the network has the capability of shuffling the order in which its neurons fire , by storing a variety of different phase - patterns within the connectivity .
if the oscillators predicted by the interference theory can generate more than one phase sequence , as in the model presented here , then this could provide a potential mechanism to explain the phenomenon of hippocampus remapping ( colgin et al .
one of the more interesting discovery of place cells behavior is indeed the remapping of the place cell representation of space in response to a changes in sensory or cognitive inputs , i.e. place cells change their firing properties ( place cells can appear disappear or move to other unpredictable locations ) .
this change may be abrupt and similar to the switch from one attractor to another ( wills et al .
if the place cell will fire at a specific place where its inputs become synchronized ( euston et al . 2007 ; blair et al .
2011 ) , by recalling a different phase - coded pattern among the ones stored in our circuit , it will change the phases of theta cells that are the inputs of the place cell , and it will change the specific  place  where the inputs become synchronized , and therefore it changes the place cell representation of the space . finally we note that even thou our model is not a continuous attractor , it shares many similarity with such a class of models .
our model is a circuit with many distinct attractors , one for each phase - relationship stored in the network , and the number of different attractor states is set by the maximal storage capacity studied here .
furthermore each attractor is a phase - coded pattern , replayed with a collective resonant frequency that can be modulated by changing for example the spiking threshold of the units . even thou during exploration the activity of place cells
may be explained by the superposition of velocity - controlled oscillators inputs , the recurrent connections inside the place cells network may have anyway a relevant role . during sleep , in absence of external input , the role of recurrent connections increases , probably due to an increase of excitability via neuromodulators or other mechanisms , and the spontaneous activity of the network show temporarily short replay of stored patterns , probably initiated simply by noise .
so the pattern activated repeatly during experience , is stored in the connectivity , and then activated during sleep when the network is near a critic point and noise is able to initiate short replay sequences .
replay of phase - coded patterns of neural activity during sleep has been observed in hippocampus and neocortex ( schwindel and mcnaughton 2011 ) .
notably in our model the time scale of reactivation is different from the time scale of storing , depends on the collective frequency which emerges from the connectivity , and may be accelerated or slowed down changing parameters such as spiking thresholds .
therefore , our model might be also relevant for replay in prefrontal cortex ( pfc ) or other cortical areas in which replay is accelerated with respect to awake activity . in the hippocampus , spikes representing adjacent place fields occur in rapid succession within a single theta cycle during behavior .
therefore , relative to this within theta cycle rate , reactivation during sleep in hippocampus is not accelerated .
however , reactivation in rat pfc is clearly compressed five to eight times relative to the waking state ( euston et al . 2007 ; schwindel and mcnaughton 2011 ) .
indeed , while in hippocampus one may think that the coding sequence is the within - theta cycle , in prefrontal cortex it is clearly seen that the cross - correlation among cells during sleep replay is time compressed compared with the cross - correlation during waking state .
the playback speed declines over time as does the strength of the replay , which is consistent for example with a simple increases of spiking threshold in our model .
we studied the temporal dynamics , including the storage and replay properties , in a network of spiking integrate and fire neurons , whose learning mechanism is based on the spike - timing dependent plasticity .
the temporal patterns we consider are periodic spike - timing sequences , whose features are encoded in the relative phase shifts between neurons .
the importance of oscillations and precise temporal patterns has been pointed out in many brain structure , such as cortex ( buzsaki and draguhn 2004 ) , cerebellum ( dangelo et al . 2009 ; dangelo and de zeeuw 2009 ) , or olfactory system ( gelperin 2006 ) . the proposed associative memory approach , with selective replay of stored sequence , can be a method for recognize an item , by activating the same memorized pattern in response to a similar input .
another possibility is to have a way to transfer a memorized item to another area of the brain , such as for memory consolidation during sleep . during sleep , indeed , few spikes with the right phase relationship may initiate the retrieval of one of the patterns stored in the network and this reactivation may be useful for memory consolidation . the stored pattern is an attractor of the network dynamics , that is the dynamics spontaneously goes back to the retrieved phase - coded pattern for all the perturbations which leave the system within the basin of attraction .
therefore phase errors do not sum up , and the phase relationships may be transferred and kept stable over long time scales . the time scale of the pattern during retrieval , i.e. the period of oscillation 1/ , depends on ( 1 ) the time constants of the single neuron m and s , ( 2 ) the spiking threshold th of the neurons , and ( 3 ) the connectivity , through the stdp learning shape a( ) and the time scale of the pattern during learning mode 1/. different areas of the brain may have different shape of stdp to subserve different oscillation frequencies and different functional role . here
we fix the shape of a( ) to the one observed in hippocampal cultures ( see fig .  1 and bi and poo 1998 ) and focus on the dependence on the spiking threshold th of the neurons .
the spiking threshold can modulate the frequency of the collective oscillation , leaving unaffected the phase relationships among the units .
this opens a possible way to govern the frequency of collective oscillation via neuromodulators , and to encode information ( such as velocity of the animal ) in the frequency of the oscillations , in addition to the information encoded in the phase relationships .
notably , in a particular range of frequencies , the spiking threshold does not affect the frequency of oscillation but changes the number of spike per cycle during the retrieval dynamics .
this means that information can be encoded via the number of spikes per cycle , independently from the information coded in the phase relationship among units , in agreement with the observations of independent rate and phase coding in hippocampus ( huxter et al .
important to note , the phase relationships and the frequency of the collective oscillation are both robust with respect to noise and to heterogeneity of the spiking threshold of the units .
a systematic study of the retrieval capacity of the network is proposed as a function of two parameters of the model : the frequency of the input pattern and the spiking threshold .
the storage capacity , evaluated as pmax / n , is always lower than the storage capacity of the hopfield model .
however , the information content of a single pattern in our dynamical model with n units is higher than the information content of a pattern in the hopfield model with n - units .
indeed , an hopfield pattern is a set of n binary values while our phase - coded pattern is a set of n real number \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi^\mu_j \in [ 0 , 2\pi]$\end{document}. the role of stdp in the formation of sequences has been recently investigated in verduzco - flores et al .
these studies have shown how it s possible to form long and complex sequences , but they did not concern themselves on how it s possible to learn and store not only the order of activation in a sequence , but the precise relative times between spikes in a closed sequence , i.e. a phase - coded pattern . in our model
not only the order of activation is preserved , but also the precise phase relationship among units . the tendency to synchronization of units
is avoided in our model , without need to introduce delays or adaptation , due to the balance between excitation and inhibition that is in the connectivity of large networks when the phase coded pattern with random phases is learned using the rule in eqs .
all the connections , both positive and negative , scale with the time of presentation of patterns , keeping always a balance . indeed , since ( 1 ) the stored phase are uniformly distributed in [ 0,2 ) and ( 2 ) the learning window has the property \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\int a(\tau ) dt = 0$\end{document } , then the connectivity matrix in eq .  ( 7 ) has the property that the summed excitation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(1/n)\sum_{i , j_{ij}>0 } j _ { ij}$\end{document } and the summed inhibition \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(1/n)\sum_{i , j_{ij}<0 } j _ { ij}$\end{document } are equal in the thermodynamic limit ( indeed they are of order unity , while their difference is of order \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$1/\sqrt{n}$\end{document } ) . under this conditions ,
we investigate how multiple phase - coded patterns can be learned and selectively retrieved in the same network , as a function of time scale of patterns and network parameters .
the task of storing and recalling phase - coded memories has been also investigated in lengyel et al .
( 5 ) in a network of if neurons , the paper lengyel et al . ( 2005b ) studies this problem from a normative theory of autoassociative memory , in which real variable xi of neuron i represents the neuron i spike timing with respect to a reference point of an ongoing field potential , and the interaction h(xi , xj ) among units is mediated by the derivative of the synaptic plasticity rule used to store memories .
the model proposed here is a mechanism which combines oscillatory and attractor dynamics , which may be useful in many models of path - integration , as pointed out in section 7 .
our learning model offers a if circuit able to keep robust phase - relationship among cells participating to a collective oscillation , with a modulated collective frequency , robust with respect to noise and heterogeneities .
notably the frequency of the collective oscillation in our circuit is not sensible to the single value of the threshold of each unit , but to the average value of the threshold of all units , since all units participate to a single collective oscillating pattern which is an attractor of the dynamics .
recently there is renewed interest in reverberatory activity ( lau and bi 2005 ) and in cortical spontaneous activity ( ringach 2009 ; luczak and maclean 2012 ) whose spatiotemporal structure seems to reflect the underlying connectivity , which in turn may be the result of the past experience stored in the connectivity .
similarity between spontaneous and evoked cortical activities has been shown to increase with age ( berkes et al .
interestingly , in our if model , in order to induce spontaneous patterns of activity reminiscent of those stored during learning stage , few spikes with the right phase relationship are sufficient .
it means that , even in absence of sensory stimulus , a noise with the right phase relationships may induce a pattern of activity reminiscent of a stored pattern .
therefore , by adapting the network connectivity to the phase - coded patterns observed during the learning mode , the network dynamics builds a representation of the environment and is able to replay the patterns of activity when stimulated by sense or by chance .
this mechanism of learning phase - coded patterns of activity is then a way to adapt the internal connectivity such that the network dynamics have attractors which represent the patterns of activity seen during experience of environment .
 we study the collective dynamics of a leaky integrate and fire network in which precise relative phase relationship of spikes among neurons are stored , as attractors of the dynamics , and selectively replayed at different time scales . using an stdp - based learning process , we store in the connectivity several phase - coded spike patterns , and we find that , depending on the excitability of the network , different working regimes are possible , with transient or persistent replay activity induced by a brief signal . 
 we introduce an order parameter to evaluate the similarity between stored and recalled phase - coded pattern , and measure the storage capacity . 
 modulation of spiking thresholds during replay changes the frequency of the collective oscillation or the number of spikes per cycle , keeping preserved the phases relationship . 
 this allows a coding scheme in which phase , rate and frequency are dissociable . 
 robustness with respect to noise and heterogeneity of neurons parameters is studied , showing that , since dynamics is a retrieval process , neurons preserve stable precise phase relationship among units , keeping a unique frequency of oscillation , even in noisy conditions and with heterogeneity of internal parameters of the units . 
Introduction
The model
Designing the connections of the network
Emerging of collective patterns in the neural dynamics of the network
Storage capacity
Effects of noise and robustness of collective oscillation frequency and phase relationships
Relationship of this model with theories of path-integration
Discussion