Fundamentals of Low-Power Neuromorphic Circuit Design With Spiking Neural Networks (SNNS)
Fundamentals of Low-Power Neuromorphic Circuit Design With Spiking Neural Networks (SNNS)
1 Background
400
Pre-synaptic 300 AMPA
0
–100 GABAB
–200 GABAA
Post-synaptic –300
0 100 200 300 400 500
Time (ms)
a b
Fig. 1 (a) The dynamics of neurotransmitters. (b) Multiple synaptic currents
Figure 2 depicts a fundamental neuron circuit that integrates, fires and resets once the
output reaches the voltage potential ϴ. Despite its straightforward nature, this model
serves as the keystone for various neuron circuits, playing an essential role in
understanding the circuits of spiking neural networks (SNNs) and their operational
attributes.
The equivalent circuit behaviour can be described as:
d
τ u = -ðu - urest Þ þ RI ðt Þ ð1Þ
dt
In Fig. 3, the output fire and reset mechanisms are depicted for both the linear
region, where the potential has not yet reached the threshold value ϴ, and the point at
which the potential reaches ϴ.
Input spikes
The overall synaptic current is determined by the disparity between the instanta-
neous voltage u and the synaptic reversal potential Esyn. The conductance gsyn(t)
varies over time and can be conceptualized as an exponential pulse for ease of
interpretation, aligning with the synaptic current output AMPA (Fig. 1):
t -t k
gsyn ðt Þ = gsyn e τ Θðt - t k Þ ð3Þ
For t > tk, where the spike initiates at tk, the exponential curve described by Eq. 3
commences at time tk and exponentially decays with a time constant τ. When a
sequence of spikes is received, their outputs exhibit additive behaviour. Conse-
quently, the overall synaptic conductance is obtained by summing these individual
contributions:
t -t k
gsyn ðt Þ = gsyn e τ Θðt - t k Þ ð4Þ
k
2 Fundamentals of SNN Circuits and Mathematical Models 83
The synaptic conductance, shown in Eq. 1 and multiplied by the drive, yields
the synaptic output. To enhance the model’s realism, it is beneficial to incorporate
the rise time, as biological neurons typically exhibit two distinct time constants, with
the rising phase typically faster than the decaying phase. A slight adjustment to Eq. 4
could accommodate the inclusion of the rise time, as depicted below:
t -t k t -t k
-τ
gsyn ðt Þ = gsyn e τ 1-e rise Θðt - t k Þ ð5Þ
k
du
C = -gNa m3 hðuENa Þ - gK n4 ðu - EK Þ - gl ðu - E l Þ þ I stim ðt Þ ð6Þ
dt
functioning of the synapse circuit within this subthreshold region, as shown in Eq. 7,
where VT is the thermal voltage.
.
kðV 1 -V 2 Þ
I0 = Iie 2V T
ð7Þ
I1
I0 = Ii ð8Þ
I2
The scaling property observed in the circuits depicted in Fig. 6a, b can be applied
to synapse circuits requiring adaptation. An adaptive circuit utilizing the scaling
circuit is exemplified in Fig. 7, as documented by [3]. In this circuit, the input Ii
yields an output I0, which represents a scaled version of Ii. A supplementary input to
3 Methodology for Neuromorphic Circuit Simulation 85
Fig. 6 (a) Adjusting the current output by utilizing the two voltage controls, V1 and V2
(b) Adjustment of current output through the utilization of I1 and I2
the circuit, denoted by a pulse input to the transistor gate, regulates the scaling of the
output current to exceed that of the input current. When the supplementary input is
inactive, no scaling occurs. Additionally, the voltage Vl establishes an upper limit for
the current scaling. The supplementary input serves to facilitate the implementation
of learning algorithms, wherein it can be replaced with a logical combination of
various inputs.
In this chapter, a synapse circuit was devised and executed for software simulation
employing LTspice electronic simulator software. LTspice is freely available for
download and used by students and researchers, offering a range of analysis options,
86 Fundamentals of Low-Power Neuromorphic Circuit Design with Spiking. . .
including DC and transient circuit analysis, which were utilized in crafting the circuit
discussed herein. The LTspice library encompasses numerous transistor options,
permitting the adjustment of parameters to suit specific needs. The parameters of the
transistors employed in this study are detailed in Tables 1, 2 and 3. Notably, a pulse-
type input voltage source serves as the input to the circuit, with the LTspice
facilitating the customization of various parameters such as the period, duty cycle,
rise time, fall time, delay and number of cycles in alignment with the circuit’s
requirements. Tables 1 and 2 delineate the parameter values for both the NMOS
and PMOS transistors, while Table 3 provides insight into the dimensions of the
transistors utilized.
Numerous synapse circuits have been documented in the literature, with many
relying on the current mirror synapse due to its simplicity, albeit plagued by
concerns about gain, tau and other parameters. An alternative model introduced by
[4] employs exponential decay via log domain filtering [5], endowing the circuit
with characteristics akin to an LTI (linear time invariant) system. Additionally, the
CMOS synapse proposed by [4] shows independent controllability over the leakage
current, synaptic gain and decay time constant. This chapter aims to provide a
comprehensive exposition of the operation of the circuit, complemented by detailed
analysis and simulation outcomes.
In this chapter, the synaptic circuit employed offers distinct control over the leakage
current, synaptic weight and decay time constant, thereby enabling precise manip-
ulation of synaptic behaviour. Utilizing current feedback, the circuit achieves a first-
order response characteristic. The architecture of five PMOS transistors, three
NMOS transistors, and a capacitor is delineated, as depicted in Fig. 8. The arrival
of an action potential is simulated by applying an inverted logic pulse with a narrow
width to the gate of transistor M1. Concurrently, the voltage Vτ governs the magni-
tude of current passing through transistor M7, denoted as Iτ, thereby dictating the
decay time constant of the output. The voltage Vw, on the other hand, regulates
synaptic strength. Transistor M6 serves to convert the voltage across the capacitor
into a corresponding current. The resulting current is channelled through a current
mirror formed by transistors M4 and M5, subsequently passing through a source
follower composed of transistors M3 and M4. The inhibitory synaptic current output
is derived from the drain current of transistor M8. Notably, all NMOS transistor bulk
terminals are grounded, while the bulk terminals of the PMOS transistors are
connected to Vdd, except for M3.
3
Value 3 4.8 m 3.6 m 0.8 0.7 0.7n 0.25n 1n 0.36n 0.1u 6m 20 12 m 18n
87
88
4 Circuit Simulations
The transistors within the simulated circuit function within the subthreshold regime,
also known as the weak inversion region. In this operational mode, the relationship
between the current and voltage exhibits an exponential characteristic, succinctly
encapsulated by the drain current equation depicted in Eq. 9. This mode of operation
is pivotal in achieving the desired synaptic behaviour, facilitating nuanced control
over synaptic currents and ensuring efficient utilization of power resources:
where V T = kTq is the thermal voltage and S = WLnn . To determine the technology
current Ion, the gate and source terminals are connected to the ground, with the bulk
terminal being shorted to the source terminal. This configuration ensures a specific
operating condition where the transistor’s behaviour can be accurately characterized,
and its technology current can be reliably estimated. This step is fundamental in the
comprehensive analysis of transistor performance within the circuit context, provid-
ing essential insights into its operational dynamics and facilitating informed design
decisions where vbs = vgs = 0, and the dimensions are set to W = L = 10μ so that
S = 1. This results in the drain current equation being simplified to Eq. 10:
90 Fundamentals of Low-Power Neuromorphic Circuit Design with Spiking. . .
Fig. 9 Technology current for the NMOS. Approximately 1 nA when the threshold voltage VT
exceeds 0.1 V
V ds
ids = I 0n 1 - e - V T ð10Þ
The decay time constant of the synaptic output is governed by the leakage current
passing through M7. The LTspice simulation of the circuit is depicted in Fig. 11,
which provides a visual representation of its functionality and characteristics. Fig-
ure 12 illustrates the synaptic output obtained from the drain of M6 in conjunction
with the input potential.
4 Circuit Simulations 91
Fig. 10 Technology current for the PMOS at Vds = 1 V. Ion≈18 nA when VT > 0.1 V
Fig. 12 The simulation results depicting the synapse current output alongside the pulse input for
Vτ = 0.2 and a pulse width of 10 ms
k n v c ðt Þ
isyn ðt Þ = S8 I on e VT
ð12Þ
0:7ð0:63555862Þ
isyn ðt Þ = 10 -9 e 0:02586492 ffi 0:295μA
This result closely matches the value of 0.299 μA observed in Fig. 12.
Figures 13 and 14 illustrate how the decay time constant changes proportionally
with Vτ while maintaining a constant VW at 3.75 V, with the pulse width set at 10 ms.
Figures 12, 13 and 14 illustrate the direct impact of Vτ on the time constant, with a
slight increase in magnitude observed concurrently. Meanwhile, the voltage VW
governs the synapse gain. As depicted in Fig. 15 below, with Vτ held constant at
0.25 V, the VW is varied. Notably, while the magnitude (gain) decreases with
increasing weight, the time constant remains consistent throughout.
This visual representation offers insight into how the synapse circuit handles the
accumulation of multiple spike inputs over time, highlighting any limitations or
saturation effects that may arise under specific conditions.
94 Fundamentals of Low-Power Neuromorphic Circuit Design with Spiking. . .
Fig. 15 Vw = 3.6 V
Fig. 16 This figure shows the input (blue) spikes and corresponding synaptic outputs (green)
The analysis of the proposed design was simplified under two key assumptions.
First, it was assumed that all transistors operated solely in the subthreshold region.
However, this assumption does not hold for transistor M7 (as depicted in Fig. 11),
which occasionally operates in the triode region when Vc(t) < 4VT, which is
equivalent to 0.103 V. Consequently, the synapse current response may be faster
than initially anticipated. Second, the analysis disregarded all parasitic capacitances,
whose impact on the circuit’s output at low currents can be significant due to a
sluggish feedback response. Specifically, when Vc(t) is equivalent to 0 and an input
potential is introduced, the transconductances of the feedback loop transistors
typically range from 10-12 to 10-15, causing a delayed response from the loop.
Consequently, the voltage across the capacitor and hence the synaptic output
deviates from the behaviour of a true first-order system. To address this issue, a
small bias voltage was introduced at the source of M7, ensuring that the minimum
5 Addressing Nonidealities and Optimizing Design 95
value of Vc(t) corresponds to the bias voltage. This approach results in a small
leakage current for the synaptic output when no input potential is present. Alterna-
tively, this characteristic can be rectified by directly connecting the source of M6 to
the bias voltage. The modified circuit design is illustrated in the schematic provided
in Fig. 17.
The schematic in Fig. 17 illustrates the modified circuit design after the incorpo-
ration of a small bias voltage. This adjustment addresses the deviation from the ideal
behaviour observed in the original circuit. By introducing a bias voltage at the source
of transistor M7, the minimum value of Vc(t) is now fixed to the bias voltage,
ensuring consistent behaviour even in the absence of an input potential. Conse-
quently, the synaptic output current exhibits a small leakage current when no input
potential is applied, thus rectifying the previously observed nonideal characteristics.
This modification enhances the circuit’s performance, aligning it more closely with
the desired behaviour of a true first-order system.
96 Fundamentals of Low-Power Neuromorphic Circuit Design with Spiking. . .
Fig. 18 The integration of the bias current. Upon complete discharge of the capacitor, it diminishes
to 100 mV, aligning with the bias voltage value. The depicted period represents the entire discharge
cycle
However, integrating the bias voltage effectively regulated the quiescent current,
as shown in Fig. 18.
As depicted in Fig. 19, setting the bias voltage Vs to 100 mV results in the
synaptic current exhibiting both an exponential increase and a decay. This config-
uration ensures minimal delay in the feedback loop, thereby providing a more
precise representation of the desired first-order dynamic characteristic.
5 Addressing Nonidealities and Optimizing Design 97
Fig. 20 Schematic view of circuit modifications for independent control over rising and falling
time constants
As depicted in Fig. 20, the time constant for the synapse output is solely determined
by the current passing through M7. Nonetheless, biological synapses exhibit dual
time constants, with the rising phase typically considerably swifter than its decaying
counterpart. This poses a significant challenge in SNNs, particularly given the
narrow pulse width of the spike concerning the spike period. The subsequent circuit
design overcomes this obstacle by enabling independent regulation of both the rising
and falling time constants, thus enhancing the model’s fidelity to biological systems.
A supplementary current branch composed of transistors M8 and M9 runs parallel
to M7, providing a combined current of ITr + Iτ. In the absence of any input signal,
M8 and M9 act as open circuits, operating similarly to the configuration depicted in
Fig. 17. When the current passing through M7 increases, the steady-state value isyn
(1) diminishes because the additional current branch reduces the flow of current
responsible for charging the capacitor. Figure 21 shows the autonomous regulation
of both the rising and falling time constants. Moreover, the expression for the decay
time constant is derived from Eqs. 14 and 15:
98 Fundamentals of Low-Power Neuromorphic Circuit Design with Spiking. . .
Fig. 21 The impact of altering Vt2 on the decay time constant, with varying values for Vt2: 0.1, 0.25
and 0.4, arranged from top to bottom. The blue trace represents the input to M1, while the red trace
represents the input to M9
CV T
τd = ð14Þ
kn I τ
CV T
τr = ð15Þ
kn ðI τr þ I τ Þ
6 Discussion
This book chapter presents an insightful exploration into the realm of energy-
efficient circuit design tailored for spiking neural networks (SNNs). While this
book chapter highlights several promising aspects of the proposed circuit design,
critical discussion is warranted to assess its strengths, limitations and potential
implications. One notable strength of the chapter lies in its emphasis on leveraging
analogue complementary metal-oxide-semiconductor (CMOS) technology for VLSI
implementation, aligning with the growing demand for low-power, high-perfor-
mance neuromorphic computing solutions. By capitalizing on the nonidealities of
real transistors, particularly their minimal leakage currents at subthreshold voltages,
the circuit demonstrates commendable power efficiency, a crucial attribute for
7 Future Tasks for Students/Researchers 99
The following activities are highly recommended and provide a comprehensive hands-on
experience for students to explore and experiment with circuit modelling of
neuromorphic circuits as low-power building blocks. It is highly recommended to
implement and simulate the following circuit for further extension. Some of the work
presented in [6–10] is highly recommended, and readers interested in accessing the
software code or seeking hands-on implementation are encouraged to contact the author.
(continued)
100 Fundamentals of Low-Power Neuromorphic Circuit Design with Spiking. . .
V2
V11 V5
4.4
4.4 4.4
PULSE(0 4.7 0.001 1p 1p 0.9 1)
PULSE(0 4.7 0.001 1p 1p 0.99 1)
V4
M24 M8
AO6407
AO6407 AO6407
M7
M30 M14
V1
V9 V3
AO6407
AO6407 AO6407
M2 3.7
M25 3.75 M9 3.75
M5
M28 M12
AO6408
AO6408 AO6408
.tran 3
.tran 5 C1 .tran 5
C4 M6 C2
M29 M13
Vt6 M31 AO6408 0.1µ
AO6408 0.35µ AO6408 0.35µ
AO6408
0.25
Vt5
0.25
Vt2
0.3
Vs2
0.05
Acknowledgements The authors would like to thank M.Sc. thesis student(s) for contributing to
the experimental work conducted in the IoT laboratory within the School of Computing, Electronics
and Maths, Coventry University, UK. The work reported in this chapter was solely supervised by
Dr. Arfan Ghani at Coventry University, UK.
References