Noisy Softplus: an activation function that enables SNNs to be trained as ANNs

We extended the work of proposed activation function, Noisy Softplus [12], to fit into training of layered up spiking neural networks (SNNs). Thus, any ANN employing Noisy Softplus neurons, even of deep architecture, can be trained simply by the traditional algorithm, for example Back Propagation (BP), and the trained weights can be directly used in the spiking version of the same network without any conversion. Furthermore, the training method can be generalised to other activation units, for instance Rectified Linear Units (ReLU), to train deep SNNs off-line. This research is crucial to provide an effective approach for SNN training, and to increase the classification accuracy of SNNs with biological characteristics and to close the gap between the performance of SNNs and ANNs.

Keywords: activation function, spiking neural network, leaky-integrate-and-fire, ReLU, Noisy Softplus

DNNs are the most promising research field in computer vision, even exceeding human-level performance on image classification tasks [8]. To investigate whether brains might work similarly on vision tasks, these powerful DNN models have been converted to SNNs. In addition, spiking DNNs offer the prospect of neuromorphic systems that combine remarkable performance with energy-efficient training and operation.

Theoretical studies have shown that biologically-plausible learning, e.g. Spike-Timing-Dependent Plasticity (STDP), could approximate a stochastic version of powerful machine learning algorithms such as Contrastive Divergence [15], Markov Chain Monte Carlo [1] and Gradient Descent [16]. It is the stochasticity, in contrast to the continuously differentiable functions used by ANNs, that is intrinsic to the event-based spiking process, making network training difficult. In practice, ANNs use neuron and synapse models very different from biological neurons, and it remains an unsolved problem to develop SNNs with equivalent performance.

Conversely, the offline training of an ANN, which is then mapped to an SNN, has shown near loss-less conversion and state-of-the-art classification accuracy. This research aims to prove that SNNs are equally capable as their non-spiking rivals of pattern recognition, and at the same time are more biologically realistic and energy-efficient. Jug et al. [10] first proposed the use of the Siegert function to replace the sigmoid activation function in training Restricted Boltzmann Machine (RBM). The Siegert units map incoming currents driven by Poisson spike trains to the response firing rate of a Leaky Integrate-and-Fire (LIF) neuron. The ratio of the spiking rate to its maximum is equivalent to the output of a sigmoid neuron. A spiking Deep Belief Network (DBN) [21] of four layers of RBMs was implemented on neuromorphic hardware, SpiNNaker [6], to recognise hand written digits in real time.

However, cortical neurons seldom saturate their firing rate as sigmoid neurons. Thus ReLU were proposed to replace sigmoid neurons and surpassed the performance of other popular activation units thanks to their advantage of sparsity [7] and robustness towards the vanishing gradient problem. Recent developments on ANN-trained SNN models has focused on using ReLU units and converting trained weights to fit in SNNs. Better performance [2,4] than Siegerttrained RBM has been demonstrated in Spiking ConvNets, but this employed simple integrate and fire (IF) neurons without leakage. The training used only ReLUs and zero bias to avoid negative outputs, and applied a deep learning technique, dropout [19], to increase the classification accuracy. Normalising the trained weights for use on an SNN employing IF neurons only was relatively straightforward and maintained classification accuracy. This work was extended to a Recursive Neural Network (RNN) [5] and run on the TrueNorth [14] neuromorphic hardware platform. Except for the popular, simplified version of ReLU,  max(0, � wx), the otherimplementation of log(1 +  ex), “Softplus”, is more biologically realistic. Recent work [9] proposed the Soft LIF response function for training SNNs, which is equivalent to Softplus activation of ANNs.

We have discussed the difficulty of SNN training, and considered existing state-of-the-art solutions. We then propose an activation function, Noisy Softplus, and demonstrate how it fits the network dynamics composed of spiking neurons in Section 2. Section 3 will illustrate a complete SNN training method which employs Noisy Softplus activations and can be generalised to other activation functions. To validate the classification accuracy, a convolutional network (ConvNet) is trained on the MNIST database following this training mechanism and tested directly on an SNN in Section 4.

The existing work of modelling the response firing activity of LIF neurons using Siegert [10] function has suggested an approach to modelling activation functions for spiking neurons whose input is a synaptic current generated by spike arrivals and whose output is the firing rate of a sequence of spikes. This section will start from the neural science background to demonstrate those physical quantities used in the Siegert formula. Then we will compare and show the differ-ence between analytical estimation and practical simulations of spiking neurons. Consequently, we propose the new activation function, Noisy Softplus, to replace Siegert function for modelling the practical spiking activities of LIF neurons.

2.1 Neural Science Background

The LIF neuron model follows the following membrane potential dynamics:


The membrane potential V changes in response to the input current I, starting at the resting membrane potential  Vrest, where the membrane time constant is  τm =RmCm, Rmis the membrane resistance and  Cmis the membrane capacitance. The central idea in converting spiking neurons to activation units lies in the response function of a neuron model. Given a constant current injection I, the response function, i.e. firing rate, of the LIF neuron is:


otherwise the membrane potential cannot reach the threshold  Vthand the output firing rate is zero. The absolute refractory period  trefis included, during which period synaptic inputs are ignored.

However, in practice, a noisy current generated by the random arrival of spikes, rather than a constant current, flows into the neurons. The noisy current is typically treated as a sum of a deterministic constant term,  Iconst, and a white noise term,  Inoise. Thus the value of the current is Gaussian distributed with  mImean and  sI2variance. The white noise is a stochastic process  ξ(t) with mean 0 and variance 1, which is delta-correlated, i.e., the process is uncorrelated in time so that a value  ξ(t) at time t is totally independent on the value at any other time  t′. Therefore, the noisy current can be seen as:


and accordingly, Equation (1) becomes:


We then multiply the both sides of Equation (4) by a short time step dt, the stochastic differential equation of the membrane potential satisfies an Ornstein-


Fig. 1: Response function of the LIF neuron with noisy input currents with different standard deviations.

Uhlenbeck process:


The last term dWtis a Wiener process, that  Wt+dt − Wtobeys Gaussian distribution with mean 0 and variance dt. The instantaneous mean  µand variance  σ2of the change in membrane potential characterise the statistics of V in a short time range, and they can be derived from the statistics of the noisy current:


The response function [17,11] of the LIF neuron to a noisy current, also known as Siegert formula [18], is driven by the  µand  σ:


Figure 1 shows the response curves (Equation (7)) of a LIF neuron driven by noisy currents with Gaussian noise of  mImean and  sIstandard deviation. The parameters of the LIF neuron are all biologically valid (see the listed values in Table 1), and the same parameters are used throughout this paper. The solid (zero noise) line in Figure 1 illustrates the response function of such an LIF

Table 1: Parameter setting for the current-based LIF neurons using PyNN.


neuron injected with constant current, which inspired the proposal of ReLUs. As noise increases, the level of firing rates also rises, and the Softplus activation function approximates the response firing activity driven by current with Gaussian white noise added. Softplus units only represent a single level of noise that, for example, the dotted line in Figure 1 is drawn when  sI = 1.

Although the use of Siegert function opened the door for modelling response activities of LIF neurons as activation function used in ANNs, there are several drawbacks of this method:

most importantly, the numerical analysis on an LIF response function is far from accurate in practice. “Practice” in the paper means SNN simulations using LIF neurons. Thus the inaccurate model generates error between the estimation and the practical response firing rate. This issue will be addressed in Section 2.2.

the high complexity of the Siegert function causes much longer training time and more energy, let alone the high-cost computation on its derivative.

the training uses Siegert function to estimate the output firing rate in the forward path, but uses the derivative of sigmoid function during back propagation for lower complexity. However, due to the difference between Siegert and sigmoid functions, error is introduced.

neurons have to fire at high frequency (higher than half of the maximum firing rate) to represent the activation of a sigmoid unit; as a result the network activity results in high power dissipation.

better learning performance has been reported using ReLU, so modelling ReLU-like activation function for spiking neurons is needed.

Therefore, we propose the Noisy Softplus function which provides solutions to the drawbacks of Siegert unit.

2.2 Mismatch of Siegert Function to Practice

To verify the response function in practice, simulation tests were carried out using PyNN [3] to compare with the analytical results. The noisy current was


Fig. 2: Recorded response firing rate of a LIF neuron driven by NoisyCurrentSource sampled at every (a) 1 ms and (b) 10 ms. Averaged firing rates of 10 simulation trails are shown in bold lines, and the grey colour fills the range between the minimum to maximum of the firing rates. The analytical LIF response function, Siegert formula (Equation (7)), is drawn in thin lines (shown in Figure 1) to compare with the practical simulation.

produced by NoisyCurrentSource in PyNN which is a constant current of  mIadded to a Gaussian white noise of zero mean and  s2Ivariance. The noise was drawn from the Gaussian distribution in a time resolution of dt. We chose dt = 1 ms which is the finest resolution for common SNN simulators, and dt = 10 ms for comparison. For a given pair of  mIand  s2I, a noisy current was injected into a current-based LIF neuron for 10 s, and the output firing rate was the average over 10 trials.

The curves in Figures 2 illustrate output firing rate driven by different level of noise The differences relative to the analytical results (thin lines) is due to the time resolution, dt, of the NoisyCurrentSource. The sampled signals are shown in Figures 3 (a) and (b). The discrete sampling of the noisy current introduces time step correlation to the white noise, shown in Figures 3 (e) and (f), where the value remains the same within a time step dt. Although both current signals followed the same Gauss distribution, see Figures 3 (g) and (h), the current is a white noise when dt = 1 ms, but a coloured noise, e.g. increased Power Spectral Density (PSD) at lower frequency, when dt = 10 ms, see Figures 3 (c) and (d). Thus the Siegert formula, Equation (7), can only approximate the LIF response of noisy current with white noise.

A more realistic simulation of a noisy current can be generated by 100 Poisson spike trains, where the mean and variance of the current are given by:


where  τsynis the synaptic time constant, and each Poisson spike train connects to the neuron with a strength of  wiand a firing rate of  λi. Two populations of Poisson spike sources, for excitatory and inhibitory synapses respectively, were


Fig. 3: NoisyCurrentSource samples noisy currents from Gauss distribution at every 1 ms (left) and 10 ms (right). The signals are shown in time domain in (a) and (b), and spectrum domain in (c) and (d). The autocorrelation of both current signals are shown in (e) and (f). The distribution of the discrete samples are plotted in bar chart to compare with PDF of the original Gauss distribution, shown in (g) and (h).


Fig. 4: Recorded response firing rate of a LIF neuron driven by synaptic current with two synaptic time constants: (a)  τsyn=10 ms and (b) τsyn=10 ms. Averagedfiring rates of 10 simulation trails are shown in bold lines, and the grey colour fills the range between the minimum to maximum of the firing rates. The firing rates recorded driven by NoisyCurrentSource, are drawn in thin lines (shown in Figure 3) to compare with.

connected to a single LIF neuron to mimic the noisy currents. The firing rates of the Poisson spike generators were determined by the given  mIand  sI. Figures 4 illustrate the recorded firing rates responding to the Poissoin spike trains compared to the mean firing rate driven by NoisyCurrentSource in Figure 3. Note that the estimation of LIF response activity requires noisy current with white noise, however in practice the release of neurotransmitter takes time (τsyn >>0)

and the synaptic current the decays exponentially  Isyn = I0e

τsyn. Figures 5 (a) and (b) show two examples of synaptic current of 0 nA mean and 0.2 standard deviation driven by 100 neurons firing at the same rate and with the same synaptic strength (half excitatory, half inhibitory), but of different synaptic time constant. Therefore, the current at any time t during decaying period is dependant to the value at previous time step, which makes the synaptic current a coloured noise, see Figures 5 (c) and (d).

We observe in Figure 4 (a) that the response firing rate to synaptic current is higher than the NoisyCurrentSource for all the 10 trials. This is caused by the coarse resolution (1 ms) of the spikes, thus the standard deviation of the current is larger than 0.2, shown in Figure 5 (g); and the  τsyn, even short as 1 ms, adds coloured noise instead of white noise to the current. However, Figure 5 (h) shows a similar firing rate of both the synaptic driven current and the NoisyCurrentSource, since both of the current signals have similar distribution and time correlation. Nevertheless, the analytical response function, Siegert formula, cannot approximate either of the practical simulations.


Fig. 5: Noisy currents generated by 100 Poinsson spike trains to a LIF neuron with synaptic time constant  τsyn=1 ms (left) and τsyn=10 ms (right). The cur-rents are shown in time domain in (a) and (b), and spectrum domain in (c) and (d). The autocorrelation of both current signals are shown in (e) and (f). The distribution of the generated samples are plotted in bar chart to compare to the expected Gauss distribution, shown in (g) and (h).

2.3 Noisy Softplus

Due to the limited time resolution of common SNN simulators and the time taken for neurotransmitter release,  τsyn, mismatches exist between the analytical response function, the Siegert formula, and practical neural activities. Consequently, a unified activation function is required to model the practical responses of a LIF neuron. Inspired by the set of response functions triggered by different levels of noise, we propose the Noisy Softplus activation function:


where x refers to the mean current, y indicates the strength of the output firing rate,  σplays an important role to define the noise level, and k, which is determined by the neuron parameters, controls the shape of the curves. Note that the novel activation function we propose contains two parameters, the current and its noise; both are naturally obtained in spiking neurons. Figure 6 shows the activation function in curve sets corresponding to different discrete noise levels, and in a 3D plot.

The derivative is the logistic function scaled by  kσ:


which could be easily applied to back propagation in any network training. However, such a derivative function of low complexity does not present in Siegert function.

The activation function can be scaled up by a factor, S, to represent the firing rate  λof a LIF neuron driven by a noisy current x.


Noisy Softplus fits well to the practical response firing rate of the LIF neuron with suitable calibration of k and S, see Figure 7. The parameter pair of (k, S) is curve-fitted with the triple data points of (λ, x, σ). The fitted parameter was set to (k, S) = (0.19, 208.76) for the practical response firing rate driven by synaptic noisy current with  τsyn = 1 ms and was set to (k, S) = (0.35,201.06) when  τsyn = 10 ms. The calibration currently is conducted by linear least squaresregression; numerical analysis is considered however for future work to express the factors with biological parameters of a LIF neuron.

We have discussed modelling response firing activity of a LIF neuron with a uni-fied activation function, Noisy Softplus. However the demonstration of mapping the physical activity to numerical ANN calculations is still required to train the layered-up deep network.


Fig. 6: Noisy Softplus fits to the response function of the LIF neuron. Noisy Softplus in (a) curve sets and (b) 3D.

3.1 Equivalent Input and Output

Neurons in ANNs take inputs from their previous layer, and feed the weighted sum of their input,  netj = �i wijxi, to the activation function. The transformed signal then forms the output of an artificial neuron,  yj = f(netj), see Figure 8. However, Noisy Softplus takes physical quantities of current, and firing rate as input/output, thus an extra step is still needed to map the firing rate to numerical values in ANNs. According to Equation (8), the mean of the current


Fig. 7: Noisy Softplus fits to the response firing rates of LIF neurons.


Fig. 8: Artificial neuron model in ANNs.

feeding into a spiking neuron is equivalent to net of artificial neurons, where


The noise level of Noisy Softplus,  σ2, is the variance of the current, which also can be seen as a weighted sum of the same input x but with different weights:


Noisy Softplus transforms the noisy current with parameters of (netj, σj) to the equivalent ANN output  yj, where it can be scaled up by the factor S to the firing rate of SNNs. Note that the calculation of noise level is not necessary for activation functions other than Noisy Softplus, for example, it can be set to a


Fig. 9: Artificial spiking neuron takes scaled firing rates as input, then transforms weighted sum in some activation unit to its output which can be scaled-up to the firing rate of an output spike train.

constant for Softplus or 0 for ReLU. We name the neuron model ‘artificial spiking neurons’ which takes firing rates of spike trains as input and output. The entire artificial spiking neuron model is then generalised to any ReLU/Softplus-like activation functions, See Figure 9.

3.2 Layered-up Network

Referred to Figure 9, if we move the left end process of  ×τsynto the right end after  λj, Figure 9 forms the same model structure as artificial neurons shown in Figure 8: neurons take x as input and outputs y, and this conversion is illustrated in Figure 10. The process within such an artificial neuron is divided into weighted summation and activation, which also applies to SNN modelling by combining the scaling factor S and the synaptic time constant  τsynto activation functions. Thus the combined activation function for modelling SNNs should be:


and its derivative function which is used when back propagates is:


Thus, using this method of ANN-trained SNNs, the activation functions are of lower complexity than the Siegert formula, and their corresponding derivative functions can be directly used for back propagation. Furthermore, the method enables ReLU-like activation functions for SNN training, thus improving the recognition accuracy while keeping a relative lower firing rate compared to sigmoid neurons. Most significantly, the ANN-trained weights are ready for use in SNNs without any transformation, and the output firing rate of a spiking neuron can be estimated in the ANN simulation.


Fig. 10: Transforming artificial spiking neurons to artificial neurons for SNN modelling. The combined activation links the firing activity of a spiking neuron to the numerical value of ANNs.

3.3 Fine Tuning

There are two aspects to the fine tuning which makes the ANN closer to SNNs: Firstly, using Noisy Softplus activation functions in a whole trained network operates every single neuron running in a similar noise level as in SNNs, thus the weights trained by other activation functions will be tuned to fit closer to SNNs. Secondly, the output firing rate of any LIF neuron is greater than zero as long as noise exists in their synaptic input. Thus adding up a small offset on the labels directs the model to approximate to practical SNNs.

The labels of data are always converted to binary values for ANN training. This enlarges the disparities between the correct recognition label and the rest to train the network for better classification capability. Consequently, we train the network as stated in Section 3.2 with any activation function and then fine-tune it with Noisy Softplus to take account of both accuracy and practical network activities of SNNs. However, we add a small number, for example 0.01, to all the binary values of the data labels. Doing so helps the training to loosen the strict objective function to predict exact labels with binary values. Instead, it allows a small offset to the objective. An alternative method is to use Softmax function at the top layer, which aims to map real vectors to the range of (0, 1) that add up to 1. However, without a limit on the input of Softmax, it will be easy to reach or even exceed the highest firing rate of a spiking neuron. The result of fine tuning on a Convnet will be demonstrated in Section 4.2.

A convolutional network model was trained on MNIST, a popular database in neuromorphic vision, using the ANN-trained SNN method stated above. The architecture contains 28  ×28 input units, followed by two convolutional layers 6c5-2s-12c5-2s, and 10 output neurons fully connected to the last pooling layer to represent the classified digit.

The training only employed Noisy Softplus units that all the convolution, average sampling, and the fully-connected neurons use Noisy Softplus function with no bias. The parameters of the activation function were calibrated as, (k = 0.30, S = 201), for LIF neurons (see Tablel 1) of τsyn = 5 ms. The in-put images were scaled by 100 Hz to present the firing rates of input spikes. The weights were updated using a decaying learning rate, 50 images per batch and 20 epochs. The ANN-trained weights were then directly applied in the corresponding convolutional SNN without any conversion for recognition tasks.

4.1 Neural Activity

To validate how well the Noisy Softplus activation fits to the response firing rate of LIF neurons in a real application, we simulated the model on NEST using the Poisson MNIST dataset [13] and the neurons of a convolutional map were observed.

Figure 11 shows a small test of ten MNIST digits presented in Poisson spike trains for 1 s each. A trained 5×5 kernel (Figure 11(c)) was convolved with these input digits, and the fourth digit ‘0’ (Figure 11(b)) and its convolved output of the feature map was shown in Figure 11(d) as firing rate. The output firing rate was recorded during a real-time SNN simulation on NEST, and compared to the modelled activations of Equation (14) in ANNs.

The input x of the network was calculated as Equation (12):  xi = λiτsyn, and so as the weighted sum of the synaptic current (see Equation (12)),  netjand its variance (see Equation (13)),  σ2j. With three combined activation functions as Equation (14):


we compare the output to the recorded SNN simulations. ReLU assumes a nonnoise current, and Softplus takes a static noise level thus  σjis not used for either of them, meanwhile Noisy Softplus adapts to noise automatically with σj. The experiment took the sequence of 10 digits shown in Figure 11(a) to the same kernel and the estimated spike counts using Noisy Softplus fit to the real recorded firing rate much more accurately than ReLU and Softplus, see 12. The Euclidean distance,��j(yj/τsyn − λj), between the spike counts and the predicted firing rates by Noisy Softplus, ReLU and Softplus was 184.57, 361.64 and 1102.76 respectively. We manually selected a static noise level of 0.45 for Softplus, whose estimated firing rates located roughly on the top slope of the real response activity. This resulted in longer Euclidean distance than using ReLU, since most of the input noisy currents were of relatively low noise level in this experiment. Hence, the firing rate driven by lower noise level is closer to ReLU curve than Softplus.


Fig. 11: Images presented in spike trains convolve with a weight kernel. (a) The 28  ×28 Poisson spike trains in raster plot, representing 10 digits in MNIST. (b) The firing rate of all the 784 neurons of the fourth image, digit ‘0’, is plotted as a 2D image. (c) One out of six of the trained kernels (5  ×5 size) in the first convolutional layer. (d) The spike trains plotted as firing rate of the neurons in the convolved 2D map.

Figure 13 demonstrates the output firing rates of the 10 recognition neurons when tested with the digit sequence. The SNN successfully classified the digits where the correct label neuron fired the most. We trained the network with binary labels on the output layer, thus the expected firing rate of correct classification was 1/τsyn = 200 Hz according to Equation (10). The firing ratesof the recognition test fell to the valid range around 0 to 200 Hz. This shows another advantage of the proposed ANN-trained method that we can constrain the expected firing rate of the top layer, thus preventing SNN from exceeding its maximum firing rate, for example 1000 Hz when time resolution of SNN simulation set to 1 ms.


Fig. 12: Noisy Softplus fits to the neural response firing rate in an SNN simulation. The recorded firing rate of the same kernel convolved with 10 images shown in Figure 11 in SNN simulation, comparing to the prediction of activations of Softplus, ReLU, and Noisy Softplus.

4.2 Recognition Performance

Here we focus on the recognition performance of the proposed ANN-trained SNN method. Before looking into the recognition results, it is significant to see the


Fig. 13: Output firing rates for recognising 10 hand written digits.

learning capability of the proposed activation function, Noisy Softplus. We compared the training using ReLU, Softplus, and Noisy Softplus by their loss during training averaged over 3 trials, see Figure 14. ReLU learned fastest with the lowest loss, thanks to its steepest derivative. In comparison, Softplus accumulated spontaneous firing rates layer by layer and its derivative may experience vanishing gradients during back propagation, which result in a more difficult training. Noisy Softplus performance lay between these two in terms of loss and learning speed. However, the loss stabilised fastest, which means a possible shorter training time.

The recognition test took the whole testing dataset of MNIST which contains 10, 000 images. At first, all trained models were tested on the same artificial neurons as used for training in ANNs, and these experiments were called ‘DNN’ test since the network had a deep structure (5 layers). Subsequently, the trained weights were directly applied to SNN without any transformation, and these ‘SNN’ experiments tested their recognition performance on the NEST simulator. The LIF neurons had the same parameters as in training. The input images were converted to Poisson spike trains and presented for 1 s each. The output neuron which fired the most indicated the classification of an input image. Moreover, a ‘Fine tuning’ test took the trained model for fine tuning, and the tuned weights were tested on the same SNN environment. The tuning only ran for one epoch, 5% cost of the ANN training (20 epochs), using Noisy Softplus neurons with labels shifted for +0.01.

The classification errors for the tests are investigated in Table 2 and the averaged classification accuracy is shown in Figure 15. From DNN to SNN, the classification accuracy declines by 0.80%, 0.79% and 3.12% on average for Noisy softplus, ReLU and Softplus The accuracy loss was caused by the mismatch between the activations and the practical response firing rates, see example in


Fig. 14: Comparisons of Loss during training using Noisy Softplus, ReLU and Softplus activation functions. Bold lines show the average of three training trials, and the grey colour illustrates the range between the minimum and the maximum values of the trials.

Figure 12, and the strict binary labels for Noisy Softplus and Softplus activations. Fortunately, the problem is alleviated by fine tuning which increased the classification accuracy by 0.38%, 0.19% and 2.06%, and resulted in the total loss of 0.43%, 0.61%, and 1.06% respectively. The improvement of ReLU is not as great as the others, because there is no problem of strict labels during training. Softplus benefits the most from fine tuning, since not only the huge mismatch of response firing rate is greatly corrected, but also the offset on the labels helps the network to fit SNNs.

Table 2: Comparisons of classification accuracy (in %) of ANN-trained convo- lutional neural models on original DNN, NEST simulated SNN, and SNN with fine-tuned (FT) model.



Fig. 15: Classification accuracy compared among trained weights of Noisy Softplus, ReLU, Softplus on DNN, SNN and fine-tuned SNN.

The most efficient training in terms of both classification accuracy and algorithm complexity, takes ReLU for ANN training and Noisy Softplus for fine tuning. Softplus does not exhibit better classification capability and more importantly the manual selected static noise level hugely influences the mismatch between the predicted firing rates and the real data. Although Noisy Softplus shows the least classification drop from ANNs to SNNs, the training performance is still worse than ReLU.

The best classification accuracy achieved by SNN was 98.85%, a 0.20% drop from ANN test (99.05%), which was trained with ReLU and fine-tuned by Noisy Softplus. The network structure was the same with the state-of-the-art model which reported the best classification accuracy of 99.1% [4] in ANN-trained SNNs: 12c5-2s-64c5-2s-10fc. Their nearly loss-less conversion from ANNs to SNNs was achieved by using IF neurons, while our network performs the best among SNNs consisted of LIF neurons to our knowledge.

As it is a major concern in neuromorphic vision, the recognition performance over short response times is also estimated in Figure 16. After fine tuning, Softplus significantly reduced the mismatch since the randomness among the three trials shrinks to a range similar to other experiments. More obviously, fine tuning improved its classification accuracy and the response latency. Notice that all of the networks trained by three different activation functions showed a very similar stabilisation curve against time, which means they all reached an accuracy close to their best by only taking 300 ms of test.


Fig. 16: The classification accuracy of 3 trials (averaged in bold lines, grey shading shows the range between minimum to maximum) over short response times, with (a) trained weights before fine tuning, and (b) after fine tuning.

4.3 Power Consumption

Noisy Softplus can easily be used for energy cost estimation for SNNs. For a single neuron, the energy consumption of the synaptic events it triggers is:


where  λjis the output firing rate,  Njis the number of post-synaptic neurons it connects to, T is the testing time, and  Esynis the energy cost for a synaptic event of some specific neuromorphic hardware, for example, about 8 nJ on SpiNNaker [20]. Thus to estimate the whole network, we can sum up all the synaptic events of all the neurons:


Thus, it may cost SpiNNaker 0.064 W, 192 J running for 3, 000 s with synaptic events of 8  ×106/sto classify 10, 000 images (300 ms each) with an accuracy of 98.02%. The best performance reported using the larger network may cost SpiNNaker 0.43 W operating synaptic event rate at 5.34 ×107/s, consume 4271.6 J to classify all the images for 1 s each.

Most significantly, we proposed the Noisy Softplus activation function which accurately models response firing rate of LIF neurons and overcomes the drawbacks of Siegert units.

Moreover, we proposed complete SNN modelling method by using artificial neurons of combined activation; this method can be generalised to activation units other than Noisy Softplus. The training of an SNN model is exactly the same as ANN training, and the trained weights can be directly used in SNN without any transformation. This method is simpler and even more straightforward than the other ANN offline training methods which requires an extra step of converting ANN-trained weights to SNN’s.

In terms of classification/recognition accuracy, the performance of ANNtrained SNNs is nearly equivalent as ANNs, and the performance loss can be partially solved by fine tuning. The best classification accuracy of 98.85% using LIF neurons in PyNN simulation outperforms state-of-the-art SNN models of LIF neurons and is very close to the result using IF neurons [4].

The research leading to these results has received funding from the European Research Council (FP/2007-2013) / ERC Grant Agreement n. 320689 and from the EU Flagship Human Brain Project (FP7-604102). Yunhua Chen received funding from the Natural Science Foundation of Guangdong Province, China (No: 2016A030313713) and also from the Natural Science Foundation of Guangdong Province, China (No: 2014A030310169). Qian Liu personally thanks to the funding from the National Natural Science Foundation of China (61662013,U1501252), the Guangxi Natural Science Foundation (2014GXNSFDA118036), and The High Level of Innovation Team of Colleges and Universities in Guangxi and Outstanding Scholars Program Funding.

1. Buesing, L., Bill, J., Nessler, B., Maass, W.: Neural dynamics as sampling: a model for stochastic computation in recurrent networks of spiking neurons. PLoS Comput Biol 7(11), e1002211 (2011)

2. Cao, Y., Chen, Y., Khosla, D.: Spiking deep convolutional neural networks for energy-efficient object recognition. International Journal of Computer Vision 113(1), 54–66 (2015)

3. Davison, A.P., Br¨uderle, D., Eppler, J., Kremkow, J., Muller, E., Pecevski, D., Perrinet, L., Yger, P.: PyNN: a common interface for neuronal network simulators. Frontiers in neuroinformatics 2 (2008)

4. Diehl, P.U., Neil, D., Binas, J., Cook, M., Liu, S.C., Pfeiffer, M.: Fast-classifying, high-accuracy spiking deep networks through weight and threshold balancing. In: Neural Networks (IJCNN), 2015 International Joint Conference on. pp. 1–8. IEEE (2015)

5. Diehl, P.U., Zarrella, G., Cassidy, A., Pedroni, B.U., Neftci, E.: Conversion of artificial recurrent neural networks to spiking neural networks for low-power neuromorphic hardware. arXiv preprint (2016)

6. Furber, S.B., Galluppi, F., Temple, S., Plana, L., et al.: The SpiNNaker project. Proceedings of the IEEE 102(5), 652–665 (2014)

7. Glorot, X., Bordes, A., Bengio, Y.: Deep sparse rectifier neural networks. In: In- ternational Conference on Artificial Intelligence and Statistics. pp. 315–323 (2011)

8. He, K., Zhang, X., Ren, S., Sun, J.: Delving deep into rectifiers: Surpassing human- level performance on imagenet classification. In: Proceedings of the IEEE International Conference on Computer Vision. pp. 1026–1034 (2015)

9. Hunsberger, E., Eliasmith, C.: ’spiking deep networks with lif neurons’. arXiv preprint (2015)

10. Jug, F., Lengler, J., Krautz, C., Steger, A.: Spiking networks and their rate-based equivalents: does it make sense to use Siegert neurons? In: Swiss Soc. for Neuroscience (2012)

11. La Camera, G., Giugliano, M., Senn, W., Fusi, S.: The response of cortical neurons to in vivo-like input current: theory and experiment. Biological cybernetics 99(4-5), 279–301 (2008)

12. Liu, Q., Furber, S.: Noisy Softplus: a biology inspired activation function. In: In- ternational Conference on Neural Information Processing. pp. 405–412. Springer (2016)

13. Liu, Q., Pineda-Garc´ıa Garibaldia nd Stromatias, E., Serrano-Gotarredona, T., Furber, S.: Benchmarking spike-based visual recognition: a dataset and evaluation. Frontiers in Neuroscience (under review) (2016)

14. Merolla, P.A., Arthur, J.V., Alvarez-Icaza, R., Cassidy, A.S., Sawada, J., Akopyan, F., Jackson, B.L., Imam, N., Guo, C., Nakamura, Y., et al.: A million spikingneuron integrated circuit with a scalable communication network and interface. Science 345(6197), 668–673 (2014)

15. Neftci, E., Das, S., Pedroni, B., Kreutz-Delgado, K., Cauwenberghs, G.: Event- driven contrastive divergence for spiking neuromorphic systems. Frontiers in neuroscience 7 (2013)

16. O’Connor, P., Welling, M.: Deep spiking networks. arXiv preprint (2016)

17. Rauch, A., La Camera, G., L¨uscher, H.R., Senn, W., Fusi, S.: ’neocortical pyramidal cells respond as integrate-and-fire neurons to in vivo–like input currents’. Journal of neurophysiology 90(3), 1598–1612 (2003)

18. Siegert, A.J.: On the first passage time probability problem. Physical Review 81(4), 617 (1951)

19. Srivastava, N., Hinton, G.E., Krizhevsky, A., Sutskever, I., Salakhutdinov, R.: Dropout: a simple way to prevent neural networks from overfitting. Journal of Machine Learning Research 15(1), 1929–1958 (2014)

20. Stromatias, E., Galluppi, F., Patterson, C., Furber, S.: Power analysis of large- scale, real-time neural networks on SpiNNaker. In: Neural Networks (IJCNN), The 2013 International Joint Conference on. pp. 1–8 (2013)

21. Stromatias, E., Neil, D., Galluppi, F., Pfeiffer, M., Liu, S.C., Furber, S.: Scal- able energy-efficient, low-latency implementations of trained spiking deep belief networks on SpiNNaker. In: Neural Networks (IJCNN), 2015 International Joint Conference on. pp. 1–8. IEEE (2015)

designed for accessibility and to further open science