Review of Population Activity
Information from spatiotemporal pattern:

Temporal average: mean firing rate of the network:
where $\Delta t$ is the time binning. We can take the limit $\Delta t\to 0$. For spiking neuron models,
Assume the $j$th neuron fires at $t_j^{(f)}$ (where $f$ means neuron $j$ can fire at $t_j^{(1)}$, $t_j^{(2)}$, etc), we know which neuron would fire at any time $t$ by calculating
If the argument is zero for neuron $j$, we know this neuron fires at time $t$. Thus the total number of active neurons at time $t$ is
The summation of $f$ counts all the contributions from neuron $j$ since it can fire multiple times.

Temporal correlations between pulses of different neurons;
Define Homogeneous Network
Homogeneous:
 Each and every neurons are the same;
 Neurons take the same external input;
 Weight between arbitrary neuron $i$ and arbitrary neuron $j$ are the same $w_{ij}=J_0/N$.
Modeling Using Leaky Integrateandfire Model
where the input on neuron $i$ at time $t$ ($I_i(t)$) includes both external and internal signal
where $\alpha(tt_j^{(f)})$ is the time course after a $f$th spike in neuron $j$.
QUESTION HERE:
Assuming that we start from a homogeneous initial condition that all neurons have the same input will the network evolve to an inhomogeneous state under tiny perturbation? This should be called the stability of the network.
The method I can think of to check the instability is to add a field of perturbation, i.e., a perturbation to each input and check the evolution of these perturbations. If all the perturbations increase with time, we obtain a hint that the system could be unstable. Then we can investigate using numerical methods or linear stability analysis to find out exactly how.
COMMENT: This perturbation method is similar to the one discussed in Section 6.2.1 as stochastic input.
Regardless of the question above, we assume that no perturbation is included so that the system evolve homogeneously. Input current on arbitrary neuron is the same as each other
Due to homogeneity, we can plug in the definition of population activity
for every neurons. At time $s$, we know the activity is calculated as , the input current at time $t$ which is later than $s$ depends on all the history of population activity, which means we include all the history to find the current input current, i.e.,
Modeling Using SRM
SRM
Using an argument similar to integrateandfire model, we can plug in the population activity to equation (\ref{srm}).
Recall that spikes occur as the potential $u_i$ reaches the threshold $\theta$.
Density Equations
Intergrateandfire and stochastic spike arrival
Each neuron is described by its membrane potential. We can find the percentage of neurons within a range of membrane potential $[u,u+du]$, which is
where $p(u,t)$ is the membrane potential (distribution) density.

Due to the normalization $1/N$, we expect integral over all the possible potentials gives us 1, i.e.,

Define the flux of neutrons going through the threshold by inspecting the population increase at potential $u_{\mathrm{reset}}$.
About the stochastic input, we assume that a total number of M types of synapses are introduced. For each type $k$ we have a input rate $\nu_k$. The equation for the $p(u,t)$ is written down in Section 6.2.1 as equation (6.14). In general, we can define the flux going across a potential by checking the time derivative of $p(u,t)$,
where
where $dJ_{\mathrm{drift}}$ rate of neurons that has a potential continuously evolved into $u$ during time $[t,t+dt]$, while $dJ_{\mathrm{jump}}$ stands for the rate of neurons jumping from all different potential $uw_k$’s into $u$ due to stochastic input. READ Fig. 6.2
For a better understanding of $dJ$, please refer to Eq. (\ref{continuityequation}).
Flux means the percentage of neurons that goes across a specific potential $u$ at time $t$ per unit potential. In other words, it’s the changing rate of the membrane potential density at potential $u$ during time interval $[t,t+dt]$.
Equation 6.18 shows the population activity.
Continuity Equation
Continuity equation describes the rate of change. In this case, continuity is the rate of change in percentage of neurons at time $t$ within a range of potentials $[u_0,u_1]$, so that
An example of continuity equation is the conservation of charge. The change of charge equals to the charge coming in and the charge going out.
The we know the change in the percentage of neurons within potential range $[u_0,u_1]$ is the percentage of neurons moving within the potential range $J(u_0,t)$ and the percentage of neurons going out $J(u_1,t)$. In the language of density, we have
where $u\neq u_{\mathrm{reset}}$ and $u\neq \theta$.
Diffusion Approximation
Expand equation 6.14 in terms of $w_k$ and keep only $\mathscr{O}(w_k^2)$,
Plug into equation 6.14
The first two terms are perfect derivatives
The equation we need to deal with is
Fokker–Planck equation
$A(t)$ is the flux at threshold $\theta$. To solve it we simply set $u=\theta$,
To expand $A(t)$ around $u=\theta$, we need
so that where $p(\theta,t)=0$ since this is the threshold and no neurons can stay at this potential. Eq.(\ref{eqputexpandtheta}) and Eq.(\ref{eqintputexpandtheta}) are plugged back into the express for $A(t)$. The simplified threshold flux is
I obtained a different sign. In the text book $A$ is negative.
Stationary Solution
Why stationary solution? Experiments.
Define
Stationary state means $p(u,t)$ is time translation invariant, i.e., $p$ only depends on $u$. The FokkerPlanck equation becomes
where $h=R I^{\mathrm{ext}}+\tau_m \sum_k v_k w_k$.
$A$ is the change of percentage of neurons at threshold (some kind of flux), dimension analysis shows the terms left is of dimension $[\partial_u J]$.
Dimension of $\delta(uu_r)$ is $[1/u]$.
Indeed we define a flux
Educated guess
comments powered by Disqus