Versatile CMOS Analog LIF Neuron for Memristor-Integrated Neuromorphic Circuits
thanks: We acknowledge financial support from the EU: ERC-2017-COG project IONOS (# GA 773228) and the CHIST-ERA UNICO project. FA thanks the MEIE for support through the chair in neuromorphic engineering.

Nikhil Garg 3IT-LN2, Université de Sherbrooke
Sherbrooke, Canada
IEMN-CNRS, Université de Lille
Lille, France
   Davide Florini Patrick Dufour 3IT-LN2, Université de Sherbrooke
Sherbrooke, Canada
   Eloir Muhr Mathieu C. Faye Marc Bocquet IM2NP-CNRS, Aix-Marseille Université
Marseille, France
   Damien Querlioz C2N-CNRS, Université Paris-Saclay
Paris, France
   Yann Beilliard Dominique Drouin 3IT-LN2, Université de Sherbrooke
Sherbrooke, Canada
   Fabien Alibart 3IT-LN2, Université de Sherbrooke
Sherbrooke, Canada
IEMN-CNRS, Université de Lille
Lille, France
   Jean-Michel Portal IM2NP-CNRS, Aix-Marseille Université

Marseille, France
Abstract

Heterogeneous systems with analog CMOS circuits integrated with nanoscale memristive devices enable efficient deployment of neural networks on neuromorphic hardware. CMOS Neuron with low footprint can emulate slow temporal dynamics by operating with extremely low current levels. Nevertheless, the current read from the memristive synapses can be higher by several orders of magnitude, and performing impedance matching between neurons and synapses is mandatory. In this paper, we implement an analog leaky integrate and fire (LIF) neuron with a voltage regulator and current attenuator for interfacing CMOS neurons with memristive synapses. In addition, the neuron design proposes a dual leakage that could enable the implementation of local learning rules such as voltage-dependent synaptic plasticity. We also propose a connection scheme to implement adaptive LIF neurons based on two-neuron interaction. The proposed circuits can be used to interface with a variety of synaptic devices and process signals of diverse temporal dynamics.

Index Terms:
Neuromorphic computing, Analog circuits, In-memory computing, memristors, ASIC

I Introduction

Analog circuits and devices can efficiently emulate neural dynamics in real-time by utilizing physical mechanisms such as Kirchoff’s law through memristive synapses for signal transmission and capacitive charging on CMOS devices for temporal integration of current [1]. Such strategies are inherited from the seminal work of C. Mead [2] to implement silicon neurons with rich temporal dynamics at low power. More recently, nanoscale non-volatile memristive devices have been considered for synaptic function implementation and offer several interesting features, such as non-volatile synaptic weight storage, analog programming, and low-power operation. Such devices are also CMOS compatible and could allow 3D integration [4] for high-density synaptic arrays to implement in-memory computing architectures [5]. The development of spiking neural networks based on these technologies still needs to address several challenges, such as synapses/neuron impedance matching and operation stability validation.

Refer to caption
Fig. 1: (A) Optical micrograph of the realized ”UNICO” ASIC with an arrow denoting the signal chain from IO pad to neurons of the input layer (LIFin𝐿𝐼subscript𝐹𝑖𝑛LIF_{in}italic_L italic_I italic_F start_POSTSUBSCRIPT italic_i italic_n end_POSTSUBSCRIPT) through a resistive synaptic column and then to output layer (LIFout𝐿𝐼subscript𝐹𝑜𝑢𝑡LIF_{out}italic_L italic_I italic_F start_POSTSUBSCRIPT italic_o italic_u italic_t end_POSTSUBSCRIPT) through a memristive synaptic array. (B) The signal chain for spike transmission consists of a synaptic resistance, a low dropout regulator (LDO) for regulating the voltage at other synaptic terminals, and a current attenuator for down-scaling the current by a factor of ’K’. The leaky integrate and fire (LIF) neuron integrates the read current. The neuron membrane potential tracks the neuron state and is reset when crossing a threshold voltage level. Consequently, a spike is transmitted to the next layer of neurons through the synaptic array. (C) Test structures with connection pads probed for electrical characterization of implemented circuit blocks.

This study presents an architecture and circuits that enable stable reading of the memristive synapse into an analog spiking neuron circuit. These circuit blocks, including a low-dropout regulator (LDO) and current attenuator, were implemented alongside a LIF neuron in a 130nm CMOS integrated circuit (IC). The measurement results demonstrate the sensitivity of the neuron’s activity for a range of synaptic resistance, excitation voltage levels, and pulse widths representative of the memristive synapse operation. The modulation of threshold, leak rate, and refractory period showcase the configurability of neuron dynamics that could be adapted to match different synaptic array sizes, device conductance range, and the time scale of the deployed application. Finally, an architecture to configure generic memristive LIF neurons to an adaptive variant is proposed.

Neurons with slow leak rates and long-time-scale dynamics can be critical for emulating bio-realistic dynamics in real time. Specialized low-pass filtering circuits ([10], [11]) are often used to implement such slow dynamics with small capacitance. We implemented such an integrator circuit with modifications to realize controllable bi-directional leakage for local learning with voltage-dependent synaptic plasticity (VDSP) [15, 16]. In contrast to previous studies [8, 9], the interfacing between the integrator operating at low-current levels and the memristive device was enabled by a current attenuator implemented in the signal chain. VDSP-based synaptic plasticity can further exploit such dynamics to detect spike patterns in longer time windows and enable efficient local learning.

Refer to caption
Fig. 2: Key circuit blocks. (A) This schematic shows the synaptic resistance, Low Dropout Regulator (LDO), and current attenuator arranged to read the synaptic current to the neuron. (B) The LIF neuron is depicted as comprising a leaky integrator, fire and reset block, and buffer. (C) The schematic displays the synaptic array architecture on the left, and on the right, a micrograph of the realized chip is zoomed in around the signal chain from the synaptic array to the neuron.

The manuscript is organized as follows. The Methods section details the realized ASIC architecture with key circuit blocks to illustrate the signal chain. The results section presents the electrical characterization of the neuron subject to various input synaptic currents. The membrane threshold level, and refractory period impact on spiking activity is further evaluated. Next, the characterization results in response to spike-based (pulsed) stimuli are presented with an analysis of leak rate modulation. In the end, a connection topology for the dynamic configuration of generic LIF neurons to adaptive LIF neurons is presented, and hardware overheads for different circuit functionalities are discussed.

II Materials and Methods

II-A Implemented Design

The application-specific integrated circuit (ASIC) was realized for hardware implementation of SNN on a monolithic chip through CMOS analog neurons and back end of the line (BEOL) integrated memristive [17] synaptic devices (Fig. 1(A). The memristive devices are integrated in this design in a 1T1R configuration. In such a configuration, the output terminal of the pre-neuron sends voltage spikes (pulses) to the gate of the 1T1R downstream synaptic array that will enable input current to the post-neurons. The signal chain of the synapse and neuron is illustrated in Fig. 1(B). For neuron characterization, the presynaptic input is a voltage applied to the resistor’s first terminal (IN), emulating the memristive synapse. LDO clamps the second terminal at Vrefsubscript𝑉𝑟𝑒𝑓V_{ref}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_f end_POSTSUBSCRIPT, and the Isynsubscript𝐼𝑠𝑦𝑛I_{syn}italic_I start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT current weighted by the resistance Rsynsubscript𝑅𝑠𝑦𝑛R_{syn}italic_R start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT is scaled by a factor of K by the current attenuator. The LIF neuron membrane voltage performs temporal integration of the current and resets to Vresetsubscript𝑉𝑟𝑒𝑠𝑒𝑡V_{reset}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_e italic_t end_POSTSUBSCRIPT when Vmemsubscript𝑉𝑚𝑒𝑚V_{mem}italic_V start_POSTSUBSCRIPT italic_m italic_e italic_m end_POSTSUBSCRIPT is higher than the Vth threshold voltage. Test structures with connection pads implemented on the ASIC were probed for testing (Fig. 1C). The analog bias voltages and input signals were generated through an FPGA-controlled PCB, and signals from ASIC were measured through an oscilloscope.

The full SNN integrates two layers of neurons in which the output spike is transmitted to the gates of all the downstream synapses of the memristive array. The output of the synaptic array is connected to the second layer of LIF neurons, which present again the successive blocks of LDO and current attenuator.

TABLE I: Bias voltages used in circuits
Name Typical (V) Purpose
Vddsubscript𝑉𝑑𝑑V_{dd}italic_V start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT 3.3 Power supply
Vrefsubscript𝑉𝑟𝑒𝑓V_{ref}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_f end_POSTSUBSCRIPT 2.4 LDO reference
Vopasubscript𝑉𝑜𝑝𝑎V_{opa}italic_V start_POSTSUBSCRIPT italic_o italic_p italic_a end_POSTSUBSCRIPT 2.4 OPAMP bias (LDO)
Vgainsubscript𝑉𝑔𝑎𝑖𝑛V_{gain}italic_V start_POSTSUBSCRIPT italic_g italic_a italic_i italic_n end_POSTSUBSCRIPT 2.1 Gain modulation
Vtaunsubscript𝑉𝑡𝑎𝑢𝑛V_{taun}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_n end_POSTSUBSCRIPT 1.2 Leak rate (Down)
Vtaupsubscript𝑉𝑡𝑎𝑢𝑝V_{taup}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_p end_POSTSUBSCRIPT 1.2 Leak rate (UP)
Vrestsubscript𝑉𝑟𝑒𝑠𝑡V_{rest}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_t end_POSTSUBSCRIPT 0.6 Resting potential
Vthrsubscript𝑉𝑡𝑟V_{thr}italic_V start_POSTSUBSCRIPT italic_t italic_h italic_r end_POSTSUBSCRIPT 1.2 Neuron threshold level
Vbcompsubscript𝑉𝑏𝑐𝑜𝑚𝑝V_{bcomp}italic_V start_POSTSUBSCRIPT italic_b italic_c italic_o italic_m italic_p end_POSTSUBSCRIPT 2.4 Bias for COMP (LIF)
Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT 1 Pulse width modulation
Refer to caption
Fig. 3: Characterization of neuron’s sensitivity to synaptic resistance. (A) Measured output spike rate of neuron with respect to Rsynsubscript𝑅𝑠𝑦𝑛R_{syn}italic_R start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT for different Vreadsubscript𝑉𝑟𝑒𝑎𝑑V_{read}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_a italic_d end_POSTSUBSCRIPT levels between 100mV and 400mV. The respective Vreadsubscript𝑉𝑟𝑒𝑎𝑑V_{read}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_a italic_d end_POSTSUBSCRIPT was applied for 1s for each experiment, and the neuron’s response was recorded. Comparison of neuron’s membrane voltage and output response is shown for Rsynsubscript𝑅𝑠𝑦𝑛R_{syn}italic_R start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT of 10kΩΩ\Omegaroman_Ω (B) and 50kΩΩ\Omegaroman_Ω (C) with read voltage (Vreadsubscript𝑉𝑟𝑒𝑎𝑑V_{read}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_a italic_d end_POSTSUBSCRIPT) of 250mV.

II-B Circuits

The circuit schematics of LDO and current mirror are shown in Fig. 2(A). External Input voltage (Vinsubscript𝑉𝑖𝑛V_{in}italic_V start_POSTSUBSCRIPT italic_i italic_n end_POSTSUBSCRIPT) is applied to the first terminal of the synaptic resistance (Rsynsubscript𝑅𝑠𝑦𝑛R_{syn}italic_R start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT). The other terminal of the synapse is connected to the (+) input terminal of the operational amplifier (OPAMP) in LDO. The other terminal (-) of OPAMP is connected to Vrefsubscript𝑉𝑟𝑒𝑓V_{ref}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_f end_POSTSUBSCRIPT. Whenever a voltage greater than Vrefsubscript𝑉𝑟𝑒𝑓V_{ref}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_f end_POSTSUBSCRIPT is applied on IN, the output of OPAMP rises, and transistor M1 is turned on. The current read from Rsynsubscript𝑅𝑠𝑦𝑛R_{syn}italic_R start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT is transmitted via the feedback branch, and the voltage at the positive input of OPAMP is pulled to Vrefsubscript𝑉𝑟𝑒𝑓V_{ref}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_f end_POSTSUBSCRIPT. The output current of LDO is the accumulation of weighted current from all upstream synapses, which is fed to the downstream current attenuator.

The current attenuator is implemented with a cascoded current mirror. This topology is beneficial for a wide swing [18] across the resistance range of upstream memristive synapse. The attenuation factor was tuned by adjusting the W/L ratio of each transistor pair (e.g., M1 and M2, M3 and M4, etc.). Additionally, the sizing was performed to ensure a constant attenuation factor of 500 throughout the read current range. Such topology of LDO and current mirror is also referred to as active current mirror in previous studies [19].

The LIF neuron (Fig. 2B) performs temporal integration of signals from the current attenuator and consists of a low pass filtering based on difference pair integrator (DPI) circuit [11] for implementing the integration of synaptic current on the capacitor (C1𝐶1C1italic_C 1). The resting state potential bias voltage controls the level to which the membrane potential leaks (Vrestsubscript𝑉𝑟𝑒𝑠𝑡V_{rest}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_t end_POSTSUBSCRIPT) and was set higher than the reset potential (0V). A transistor (M3) was added to the DPI circuit to realize bi-directional leakage of membrane potential, and the leak rate is controlled by Vtaupsubscript𝑉𝑡𝑎𝑢𝑝V_{taup}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_p end_POSTSUBSCRIPT and Vtaunsubscript𝑉𝑡𝑎𝑢𝑛V_{taun}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_n end_POSTSUBSCRIPT bias voltages.

The fire and reset block is composed of a comparator for membrane voltage threshold crossing detection with an externally supplied threshold level (Vthrsubscript𝑉𝑡𝑟V_{thr}italic_V start_POSTSUBSCRIPT italic_t italic_h italic_r end_POSTSUBSCRIPT). This block is biased by Vbiascompsubscript𝑉𝑏𝑖𝑎𝑠𝑐𝑜𝑚𝑝V_{biascomp}italic_V start_POSTSUBSCRIPT italic_b italic_i italic_a italic_s italic_c italic_o italic_m italic_p end_POSTSUBSCRIPT. The neuron’s pulse width and refractory period are configured by modulating the discharging rate of capacitance C2𝐶2C2italic_C 2 through bias voltage Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT. A reset transistor (M8) discharges the membrane capacitor to the ground on a spike event, and the generated spike is transmitted to the downstream synapse through a buffer. As the reset transistor is activated throughout the spike generation period, Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT also controls the neuron’s refractory period.

The architecture of the memristive synaptic array is shown in Fig. 2(C)(left). The neurons in the input layer (LIFI1𝐿𝐼subscript𝐹𝐼1LIF_{I1}italic_L italic_I italic_F start_POSTSUBSCRIPT italic_I 1 end_POSTSUBSCRIPT to LIFI16𝐿𝐼subscript𝐹𝐼16LIF_{I16}italic_L italic_I italic_F start_POSTSUBSCRIPT italic_I 16 end_POSTSUBSCRIPT) transmit a generated spike to the next layer by applying the spike output to the transistor’s gate of the 1T1R synaptic cell. The spike amplitude is set to Vddsubscript𝑉𝑑𝑑V_{dd}italic_V start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT to minimize the voltage drop access transistor during spike transmission. All the synaptic cells in a column are connected to the corresponding output neuron through LDO and the current attenuator. The circuits were implemented in 130nm CMOS technology, and a picture from a microscope is shown to visualize the layout of the signal chain in Fig. 2(C)(right). The typical values and short descriptions of bias voltages are summarized in Table I.

Refer to caption
Fig. 4: Characterization of neuron’s transfer function with DC excitation. (A) Measured spike rate of neuron with respect to input current (Isynsubscript𝐼𝑠𝑦𝑛I_{syn}italic_I start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT) and spiking threshold bias level (Vthrsubscript𝑉𝑡𝑟V_{thr}italic_V start_POSTSUBSCRIPT italic_t italic_h italic_r end_POSTSUBSCRIPT). Comparison of the measured temporal response of membrane voltage and output for different threshold levels obtained with Isynsubscript𝐼𝑠𝑦𝑛I_{syn}italic_I start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT of 10μ𝜇\muitalic_μA (B) and 40μ𝜇\muitalic_μA (C). (D) Neuron’s spike rate with respect to input current for different Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT bias levels. (E) Pulse width (Tpwsubscript𝑇𝑝𝑤T_{pw}italic_T start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT) of the generated output spike plotted with respect to Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT bias voltage.

III Results and Discussion

III-A Synaptic resistance reading

For characterizing the synaptic resistance reading, an external resistor (Rextsubscript𝑅𝑒𝑥𝑡R_{ext}italic_R start_POSTSUBSCRIPT italic_e italic_x italic_t end_POSTSUBSCRIPT) was connected in series with the on-chip resistor of 10kΩΩ\Omegaroman_Ω (RICsubscript𝑅𝐼𝐶R_{IC}italic_R start_POSTSUBSCRIPT italic_I italic_C end_POSTSUBSCRIPT) to emulate a synaptic resistance, resulting in read current Isyn=(VinVref)/(Rext+RIC)subscript𝐼𝑠𝑦𝑛subscript𝑉𝑖𝑛subscript𝑉𝑟𝑒𝑓subscript𝑅𝑒𝑥𝑡subscript𝑅𝐼𝐶I_{syn}=(V_{in}-V_{ref})/(R_{ext}+R_{IC})italic_I start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT = ( italic_V start_POSTSUBSCRIPT italic_i italic_n end_POSTSUBSCRIPT - italic_V start_POSTSUBSCRIPT italic_r italic_e italic_f end_POSTSUBSCRIPT ) / ( italic_R start_POSTSUBSCRIPT italic_e italic_x italic_t end_POSTSUBSCRIPT + italic_R start_POSTSUBSCRIPT italic_I italic_C end_POSTSUBSCRIPT ), where Vrefsubscript𝑉𝑟𝑒𝑓V_{ref}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_f end_POSTSUBSCRIPT, the reference voltage supplied to LDO, was set to 1V, and Rext+RICsubscript𝑅𝑒𝑥𝑡subscript𝑅𝐼𝐶R_{ext}+R_{IC}italic_R start_POSTSUBSCRIPT italic_e italic_x italic_t end_POSTSUBSCRIPT + italic_R start_POSTSUBSCRIPT italic_I italic_C end_POSTSUBSCRIPT represents the net synaptic resistance (Rsynsubscript𝑅𝑠𝑦𝑛R_{syn}italic_R start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT). The relationship between the measured output spike rate and input resistance (Rsynsubscript𝑅𝑠𝑦𝑛R_{syn}italic_R start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT) is shown in Fig. 3(A). The input resistance was varied from 10kΩΩ\Omegaroman_Ω to 1MΩΩ\Omegaroman_Ω, and experiments were repeated for five different read voltages ranging between 100mV and 400mV. Such resistance range matches the one observed in our memristive devices. The applied voltage is also compatible with the read voltage range from the memristive synapses (i.e., read without weight disturbance).

To illustrate the temporal response of the neuron, the measured membrane voltage and output of the neuron are compared in Fig. 3(B-C) for synaptic resistance of 10kΩΩ\Omegaroman_Ω and 50kΩΩ\Omegaroman_Ω and read voltage (Vreadsubscript𝑉𝑟𝑒𝑎𝑑V_{read}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_a italic_d end_POSTSUBSCRIPT) of 250mV.

Sensitivity to the large range of synaptic resistance makes the neuron suitable for various memristive technologies. In the above experiments, the output spike rate was observed to vary between 8Hz and 25kHz. The ability to read and differentiate between low resistances (with high firing rates) can enable reading from several presynaptic resistances in parallel.

III-B Neuron transfer characteristics

Refer to caption
Fig. 5: Illustrative plot to show the behavior of a neuron with bi-directional leakage. The downward and upward leakage occurs when membrane voltage is greater or less than resting state potential (Vrestsubscript𝑉𝑟𝑒𝑠𝑡V_{rest}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_t end_POSTSUBSCRIPT), respectively. A leak in the upward direction occurs when the neuron is reset after the spike event (t1subscript𝑡1t_{1}italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT) and enables estimation of time elapsed (in blue). Similarly, the occurrence of a spike event in the near future (t2subscript𝑡2t_{2}italic_t start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT) can be predicted through a high value of membrane voltage (in red).
Refer to caption
Fig. 6: Characterization of temporal dynamics. (A) Input pulses of width 10μ𝜇\muitalic_μs were applied at 100Hz, with the neuron’s membrane potential and output shown across time. (B) The input pulse width was varied between experiments (7μ𝜇\muitalic_μs and 15μ𝜇\muitalic_μs), comparing the charging events of membrane voltage. (C) The experiment varied Vtaunsubscript𝑉𝑡𝑎𝑢𝑛V_{taun}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_n end_POSTSUBSCRIPT across experiments to compare the neuron’s membrane voltage evolution across time. In (A-C), the Isynsubscript𝐼𝑠𝑦𝑛I_{syn}italic_I start_POSTSUBSCRIPT italic_s italic_y italic_n end_POSTSUBSCRIPT for each pulse was of magnitude 40μ𝜇\muitalic_μA. (D) A single input pulse was applied to charge the membrane voltage close to the neuron’s threshold, and the downward leak rate is controlled by Vtaunsubscript𝑉𝑡𝑎𝑢𝑛V_{taun}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_n end_POSTSUBSCRIPT bias level and varied across experiments. (E) The magnitude of the applied spike was increased to induce neuron firing. The leakage in the upward direction (from Vresetsubscript𝑉𝑟𝑒𝑠𝑒𝑡V_{reset}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_e italic_t end_POSTSUBSCRIPT to Vrestsubscript𝑉𝑟𝑒𝑠𝑡V_{rest}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_t end_POSTSUBSCRIPT) is compared by varying Vtaupsubscript𝑉𝑡𝑎𝑢𝑝V_{taup}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_p end_POSTSUBSCRIPT bias level between experiments.

To characterize the excitation of neurons for a range of synaptic currents and bias voltages, the excitation current was varied from 10μ𝜇\muitalic_μA to 200μ𝜇\muitalic_μA. The read voltage was applied through an on-chip resistor of 10kΩΩ\Omegaroman_Ω for a long duration (100ms). The rate of output spikes generated by the neuron is plotted with respect to the injected current, and the neuron’s threshold voltage level (Vthrsubscript𝑉𝑡𝑟V_{thr}italic_V start_POSTSUBSCRIPT italic_t italic_h italic_r end_POSTSUBSCRIPT) is shown in Fig. 4(A). For the threshold level of 1.8V, the spike rate was 419Hz and 59kHz for input current of 10μ𝜇\muitalic_μA and 200μ𝜇\muitalic_μA, respectively. Similarly, for a threshold level of 0.8V, the spike rate varied between 800Hz and 68kHz. The measured membrane voltage and output spike are compared for different threshold voltages and input currents of 10μ𝜇\muitalic_μA (B) and 40μ𝜇\muitalic_μA (C).

The threshold impacts the range of firing rate (activity) that could be observed with respect to input current. A lower threshold could increase the sensitivity to low current (high synaptic resistance), but the spike rate saturate early. Conversely, a high threshold led to lower firing rates for small excitation currents, and a variation in spike rate could be observed even for high input currents. The maximum spike rate of the neuron is limited by the output pulse width or refractory period controlled by Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT bias level. The neuron’s response(rate) to different (Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT) is shown in Fig. 4(D). The maximum firing rate was 20kHz and 92kHz for (Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT) of 0.8V and 1.1V, respectively. The Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT bias voltage modulates the output spike’s pulse width and the neuron’s refractory period by controlling the discharging rate of the C2𝐶2C2italic_C 2 capacitor in Fig. 2(C). The pulse width was measured as the time difference between the crossing of OUT and Vdd/2subscript𝑉𝑑𝑑2V_{dd}/2italic_V start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT / 2 and is plotted for Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT in Fig. 4(E). The resultant pulse width of the output spike varied between 20ms and 10μ𝜇\muitalic_μs for Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT of 0.45V and 1V, respectively.

III-C Temporal dynamics

The neuron was characterized by applying short pulses at frequent intervals to access charging with respect to input pulse width and dynamics in the absence of excitation. Spikes of magnitude 40μ𝜇\muitalic_μA and width 10μ𝜇\muitalic_μs were fed to the neuron by applying voltage pulses (Vreadsubscript𝑉𝑟𝑒𝑎𝑑V_{read}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_a italic_d end_POSTSUBSCRIPT=400mV) across the 10kΩΩ\Omegaroman_Ω on-chip resistor at the rate of 100Hz. The measured response of the neuron membrane voltage and output spikes is plotted in Fig. 6(A) and results in a neuron spike rate of 8Hz. The above experiment was repeated for different pulse widths of input spikes between 5μ𝜇\muitalic_μs to 15μ𝜇\muitalic_μs. The measured response of the neuron’s membrane voltage across time is plotted in Fig. 6(B) to compare the magnitude of the increment in membrane voltage for every charging event. A pulse of 7μ𝜇\muitalic_μs led to a small increment in membrane voltage at every spike event, resulting in a firing rate of 2Hz. Whereas a 15μ𝜇\muitalic_μs input pulse could fully charge the neuron with 5 input spikes, resulting in an output spike rate of 18Hz. Since the pulse width influences the number of spikes, the postsynaptic neuron integrates before firing the pulse width of the presynaptic neuron spikes, which can be adapted for the crossbar size. For example, suppose a neuron in the postsynaptic layer is expected to integrate signals from four times input neurons in parallel. In that case, the pulse width of input neurons should be accordingly reduced by increasing the Vpwsubscript𝑉𝑝𝑤V_{pw}italic_V start_POSTSUBSCRIPT italic_p italic_w end_POSTSUBSCRIPT bias level to maintain similar temporal dynamics in the output layer. A shorter pulse also leads to a granular increment in membrane voltage, allowing down-sizing of the membrane capacitance. In the absence of input current, the membrane voltage leaks to a resting state potential (Vrestsubscript𝑉𝑟𝑒𝑠𝑡V_{rest}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_t end_POSTSUBSCRIPT). Since the reset level was set to 0V, the neuron exhibits leakage in upward direction after the refractory phase to Vrestsubscript𝑉𝑟𝑒𝑠𝑡V_{rest}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_t end_POSTSUBSCRIPT (set to 600mV). This bi-directional mechanism was implemented to differentiate between idle neurons and those in the refractory period, as shown in Fig. 5, and is beneficial for implementing local synaptic learning.

An important aspect of memristive synapses / CMOS neuron co-integration for SNN is to allow local learning based on the activity of the adjacent neurons  [12]. In local learning rules such as Spike Timing Dependent Plasticity (STDP), the spike time difference between the pre and post-neurons is converted into a programming voltage by overlapping of slow decaying voltage pulses [13]; however, this approach can consume a significant fraction of area-energy budget [14] and implementing pulses with complex shapes can become challenging. More recently, Voltage-dependent synaptic plasticity (VDSP) [15, 16] was proposed to implement learning efficiently in hardware without requiring pulse shaping circuits. In this approach, the recent pre neuron’s activity can be estimated through the neuron membrane potential (i.e., low membrane potential being associated with a recent firing event and high membrane potential associated with an imminent spiking event). Mapping of this concept to memristive devices programming requires the implementation of more complex membrane dynamics. In this study, we present a neuron circuit with bi-directional leakage to enable learning of memristive weights through VDSP.

For characterizing the tunability of leak rate in the downward direction, spikes of magnitude 40μ𝜇\muitalic_μA and width 10μ𝜇\muitalic_μs were fed to the neuron at 100Hz, and Vtaunsubscript𝑉𝑡𝑎𝑢𝑛V_{taun}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_n end_POSTSUBSCRIPT was varied between experiments. The measured membrane voltage response is compared in Fig. 6(C). For Vtaunsubscript𝑉𝑡𝑎𝑢𝑛V_{taun}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_n end_POSTSUBSCRIPT of 1.2V (typical), the output spike rate was 8Hz, which could be increased to 12Hz by lowering the leak rate. Conversely, increasing the leak rate could decrease the output spike rate to 2Hz.

For evaluation of Vtaunsubscript𝑉𝑡𝑎𝑢𝑛V_{taun}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_n end_POSTSUBSCRIPT, the neuron was excited by a single input pulse of higher magnitude to charge the membrane voltage to just below the spiking threshold, and the response was recorded for 2s. The membrane voltage is compared for values of Vtaunsubscript𝑉𝑡𝑎𝑢𝑛V_{taun}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_n end_POSTSUBSCRIPT between 1V and 1.35V in Fig. 6(D) for time intervals of 2s and 200ms. The neuron could fully leak from the threshold level to resting state potential in 8 ms for a high leak rate. Lowering the leak rate could modulate this duration to more than 2s.

The leak rate impacts the neuron’s memory window and can be adjusted to match the dynamics of the input signals. For instance, in scenarios where the neuron receives sparse spikes at a low frequency, reducing the leak rate can help preserve the neuron’s memory over a large temporal window. Conversely, increasing the leak rate when processing high-frequency input signals can regulate the neuron’s output firing rate.

Next, we characterized the upward modulation of the leak rate. The neuron was stimulated with a single high-magnitude spike to trigger an output spike event, and its membrane potential subsequently leaked from Vresetsubscript𝑉𝑟𝑒𝑠𝑒𝑡V_{reset}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_e italic_t end_POSTSUBSCRIPT to Vrestsubscript𝑉𝑟𝑒𝑠𝑡V_{rest}italic_V start_POSTSUBSCRIPT italic_r italic_e italic_s italic_t end_POSTSUBSCRIPT, as illustrated in Fig. 6(E). The biasing voltage (Vtaupsubscript𝑉𝑡𝑎𝑢𝑝V_{taup}italic_V start_POSTSUBSCRIPT italic_t italic_a italic_u italic_p end_POSTSUBSCRIPT) was adjusted between 1.1V and 1.6V. We monitored the membrane voltage over 2-second periods to assess the upward leak rate. The ability to tune the upward leak rate and the resting state voltage is beneficial for adjusting the learning dynamics by modifying the probabilities of potentiation and depression.

Refer to caption
Fig. 7: (A) Connection scheme to realize an adaptive LIF neuron with a pair of LIF neurons. The regulator neuron integrates the output signals of the primary neuron, and its membrane voltage (MEMreg𝑀𝐸subscript𝑀𝑟𝑒𝑔MEM_{reg}italic_M italic_E italic_M start_POSTSUBSCRIPT italic_r italic_e italic_g end_POSTSUBSCRIPT) is used as the threshold of the primary neuron. (B) The primary neuron was stimulated at a 1kHz rate and connected to a regulator neuron. The membrane potentials of both neurons and the output of the primary neuron are plotted using results from a SPICE circuit simulation. (C) Single regulator neuron (Nregsubscript𝑁𝑟𝑒𝑔N_{reg}italic_N start_POSTSUBSCRIPT italic_r italic_e italic_g end_POSTSUBSCRIPT) shared by a sub-population of generic LIF neurons (N1,,Nmsubscript𝑁1subscript𝑁𝑚N_{1},...,N_{m}italic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_N start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT).

III-D Configurability to adaptive neuron

The spiking neurons’ information processing and memory capacity can be further enhanced through threshold adaptation [20, 21, 22]. The designed neuron block can be used as a regulator neuron to monitor the spiking activity of a primary neuron. As depicted in Fig. 7(A), the output of the primary neuron (OUT) is connected to the input terminal of the regulator neuron (INreg𝐼subscript𝑁𝑟𝑒𝑔IN_{reg}italic_I italic_N start_POSTSUBSCRIPT italic_r italic_e italic_g end_POSTSUBSCRIPT). This connection tracks the neuron activity through the membrane potential of the regulator neuron (MEMreg𝑀𝐸subscript𝑀𝑟𝑒𝑔MEM_{reg}italic_M italic_E italic_M start_POSTSUBSCRIPT italic_r italic_e italic_g end_POSTSUBSCRIPT), which in turn sets the threshold (THR𝑇𝐻𝑅THRitalic_T italic_H italic_R) for the primary neuron.

The designed circuit blocks, connected as described in Fig. 7, were simulated in SPICE. The primary neuron, excited by 1 kHz pulses, and the corresponding changes in membrane potential for both the primary and regulator neurons are depicted in Fig. 7(B). During this simulation, the threshold voltage of the primary neuron (MEMreg𝑀𝐸subscript𝑀𝑟𝑒𝑔MEM_{reg}italic_M italic_E italic_M start_POSTSUBSCRIPT italic_r italic_e italic_g end_POSTSUBSCRIPT) increased from 600mV to 2V, resulting in the inter-spike interval of the output spikes increasing from 2ms to 5ms. A chain consisting of a 10kΩΩ\Omegaroman_Ω synaptic resistance, followed by an LDO and a current attenuator, was used to transmit output voltage pulses from the primary to the regulator neuron. Alternatively, this synaptic resistance can be replaced with a memristor, which can be programmed to introduce an additional level of plasticity at the neuron level.

Additionally, a sub population of neurons (N1subscript𝑁1N_{1}italic_N start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, N2subscript𝑁2N_{2}italic_N start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, …, Nmsubscript𝑁𝑚N_{m}italic_N start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT) could share a common regulator neuron (Nregsubscript𝑁𝑟𝑒𝑔N_{reg}italic_N start_POSTSUBSCRIPT italic_r italic_e italic_g end_POSTSUBSCRIPT), as illustrated in Fig. 7(C). Having a pool of LIF with the same regulator neuron can improve the system’s scalability by sharing resources. This architecture can configure a fully analog LIF neuron chip to an arbitrary population of LIF and ALIF neurons in runtime. Such a heterogeneous population has been shown to enhance the learning capabilities of SNNs  [23, 24, 25]. Furthermore, monitoring neuron activity over time through traces can also enhance local learning [26], providing a dynamic feedback mechanism to improve performance and adaptability.

III-E Overhead of different functionalities

The power consumption of the implemented LDO and current attenuator was estimated through SPICE simulations of implemented circuits. The static power consumption, measured when no signal was present at the input, was 10µW for the LDO and 10pW for the current attenuator. During the operation of reading from the resistive synapse, the dynamic power consumption was estimated to be 18µW for the LDO and 20 µW for the current attenuator. The static power dissipation of LDO can be accounted for leakage due to the high biasing current used to obtain fast response time and high output current levels (up to 200µA). This quick activation is necessary when using short-reading pulses (neuron output spikes). As discussed in the earlier section, short pulses lead to gradual charging of (post-synaptic) neurons, thus helpful for scaling up crossbar size or the number of pre-synaptic neurons. Moreover, power dissipation through the memristive device during reading can be lower with short-read pulses. As these circuits are shared by synapses in a row or column, energy utilization can be balanced.

The LIF neuron’s static energy and energy per spike were estimated to be 5µW and 200 pJ/spike. The high static power can be attributed to leakage due to biasing the comparator. With a lower bias current, the static power can be reduced to 17nW, but the energy per spike increases to 7nJ. This trade-off occurs because a slow comparator can lead to high leakage during the firing phase. Feedback mechanisms [27] can potentially reduce the energy per spike of the neuron with a low comparator bias current. Overheads associated with reading memristive states are also present in hardware Artificial Neural Network (ANN) implementations, where trans-impedance amplifiers and analog-to-digital converters (ADCs) sense the current from the synaptic column for subsequent computational processing and account for a major fraction of energy consumption. However, analog neurons are more suitable for interacting with analog memristive synapses. SNNs benefit from events’ sparse activity, and gating the power supply [28] of synapse reading blocks around spike events can reduce static power dissipation through leakage.

IV Conclusion

We propose a versatile CMOS circuit to integrate memristive synapses into the signal chain of analog neuromorphic circuits. This integration pathway includes a Low-Dropout Regulator (LDO), a current attenuator, and an analog Leaky Integrate-and-Fire (LIF) neuron. The circuit blocks were implemented on a 130nm CMOS ASIC. The chip features multiple instances of these blocks alongside a synaptic array composed of 1T1R cells, facilitating the integration of memristive synaptic devices. The circuit building blocks of the signal chain in the implemented ASIC were characterized. We demonstrated the neuron’s sensitivity to the synaptic state by testing resistances ranging from 10kΩΩ\Omegaroman_Ω to 1MΩΩ\Omegaroman_Ω. As a result, the neuron’s firing rate was observed to vary between 8 Hz and 25 kHz. Additionally, the neuron’s activity was characterized across a range of applied read voltages. By increasing the spiking threshold and pulse width of the neuron, the neuron can be fine-tuned to lower and higher firing rates, respectively. Excitation with a sparse train of pulses was used to measure the neuron’s temporal response. The neuron’s charging was controllable by changing the input pulse width. Consequently, the presynaptic neuron spike’s width can be adjusted per the postsynaptic neuron’s fan-in characteristics. In the absence of an input signal, the neuron retains its memory through leakage, with an adjustable leak time constant. The bi-directional leakage enables the estimation of neurons’ recent activity and should benefit online learning through local synaptic learning.

We experimentally validated the modulation of the leak rate across short and long-term intervals, which can be helpful in processing signals with various temporal dynamics. The topology of pair neurons for configuration with adaptive neurons could further enhance the capabilities of the implemented neuron. Future work will detail the architecture and characterization of the integrated CMOS-RRAM ASIC with in-situ learning.

Author Contributions

N.G. designed the circuits with contributions from E.M. and M.C.F., under the direction of J.M.P.; D.F. performed post-processing of the chip to expose the pads for probing. N.G. performed the on-chip experimental measurements, with contributions from P.D.; N.G. wrote the initial version of the manuscript. J.M.P., F.A., D.D., Y.B., D.Q., and M.B. directed the project and edited the manuscript. All authors discussed the results and reviewed the manuscript.

Acknowledgment

We acknowledge the support from CMC Microsystems for fabrication services using the 130-nanometre CMOS technology from TSMC. The authors thank Bastien Imbert for assisting in the characterization setup.

References

  • [1] Joubert, A., Belhadj, B., Temam, O. & Héliot, R. Hardware spiking neurons design: Analog or digital?. The 2012 International Joint Conference On Neural Networks (IJCNN). pp. 1-5 (2012)
  • [2] Mead, C. Neuromorphic electronic systems. Proceedings Of The IEEE. 78, 1629-1636 (1990)
  • [3] Sarpeshkar, R. Analog versus digital: extrapolating from electronics to neurobiology. Neural Computation. 10, 1601-1638 (1998)
  • [4] Shulaker, M., Hills, G., Park, R., Howe, R., Saraswat, K., Wong, H. & Mitra, S. Three-dimensional integration of nanotechnologies for computing and data storage on a single chip. Nature. 547, 74-78 (2017)
  • [5] Sebastian, A., Le Gallo, M., Khaddam-Aljameh, R. & Eleftheriou, E. Memory devices and applications for in-memory computing. Nature Nanotechnology. 15, 529-544 (2020)
  • [6] Amirsoleimani, A., Alibart, F., Yon, V., Xu, J., Pazhouhandeh, M., Ecoffey, S., Beilliard, Y., Genov, R. & Drouin, D. In-Memory Vector-Matrix Multiplication in Monolithic Complementary Metal–Oxide–Semiconductor-Memristor Integrated Circuits: Design Choices, Challenges, and Perspectives. Advanced Intelligent Systems. 2, 2000115 (2020)
  • [7] Rao, M., Tang, H., Wu, J., Song, W., Zhang, M., Yin, W., Zhuo, Y., Kiani, F., Chen, B., Jiang, X. & Others Thousands of conductance levels in memristors integrated on CMOS. Nature. 615, 823-829 (2023)
  • [8] Wu, X., Saxena, V., Zhu, K. & Balagopal, S. A CMOS Spiking Neuron for Brain-Inspired Neural Networks With Resistive Synapses and In Situ Learning. IEEE Transactions On Circuits And Systems II: Express Briefs. 62, 1088-1092 (2015)
  • [9] Lecerf, G., Tomas, J., Boyn, S., Girod, S., Mangalore, A., Grollier, J. & Saïghi, S. Silicon neuron dedicated to memristive spiking neural networks. 2014 IEEE International Symposium On Circuits And Systems (ISCAS). pp. 1568-1571 (2014)
  • [10] Chicca, E., Stefanini, F., Bartolozzi, C. & Indiveri, G. Neuromorphic Electronic Circuits for Building Autonomous Cognitive Systems. Proceedings Of The IEEE. 102 pp. 1367-1388 (2014)
  • [11] Indiveri, G. A low-power adaptive integrate-and-fire neuron circuit. Proceedings Of The 2003 International Symposium On Circuits And Systems, 2003. ISCAS ’03.. 4 pp. IV-IV (2003)
  • [12] Querlioz, D., Bichler, O., Dollfus, P. & Gamrat, C. Immunity to device variations in a spiking neural network with memristive nanodevices. IEEE Transactions On Nanotechnology. 12, 288-295 (2013)
  • [13] Saïghi, S., Mayr, C., Serrano-Gotarredona, T., Schmidt, H., Lecerf, G., Tomas, J., Grollier, J., Boyn, S., Vincent, A., Querlioz, D., La Barbera, S., Alibart, F., Vuillaume, D., Bichler, O., Gamrat, C. & Linares-Barranco, B. Plasticity in memristive devices for spiking neural networks. Frontiers In Neuroscience. 9 pp. 51 (2015), https://www.frontiersin.org/article/10.3389/fnins.2015.00051
  • [14] Ambrogio, S., Ciocchini, N., Laudato, M., Milo, V., Pirovano, A., Fantini, P. & Ielmini, D. Unsupervised learning by spike timing dependent plasticity in phase change memory (PCM) synapses. Frontiers In Neuroscience. 10 pp. 56 (2016)
  • [15] Garg, N., Balafrej, I., Stewart, T., Portal, J., Bocquet, M., Querlioz, D., Drouin, D., Rouat, J., Beilliard, Y. & Alibart, F. Voltage-dependent synaptic plasticity: Unsupervised probabilistic Hebbian plasticity rule based on neurons membrane potential. Frontiers In Neuroscience. 16 pp. 983950 (2022)
  • [16] Goupy, G., Juneau-Fecteau, A., Garg, N., Balafrej, I., Alibart, F., Frechette, L., Drouin, D. & Beilliard, Y. Unsupervised and efficient learning in sparsely activated convolutional spiking neural networks enabled by voltage-dependent synaptic plasticity. Neuromorphic Computing And Engineering. 3, 014001 (2023)
  • [17] El Mesoudy, A., Lamri, G., Dawant, R., Arias-Zapata, J., Gliech, P., Beilliard, Y., Ecoffey, S., Ruediger, A., Alibart, F. & Drouin, D. Fully CMOS-compatible passive TiO2-based memristor crossbars for in-memory computing. Microelectronic Engineering. 255 pp. 111706 (2022)
  • [18] Bruun, E. & Shah, P. Dynamic range of low-voltage cascode current mirrors. Proceedings Of ISCAS’95-International Symposium On Circuits And Systems. 2 pp. 1328-1331 (1995)
  • [19] Miguez, M., Gak, J., Oliva, A. & Arnaud, A. Active current mirrors for low-voltage analog circuit design. Circuits, Systems, And Signal Processing. 36, 4869-4885 (2017)
  • [20] Laughlin, S. A simple coding procedure enhances a neuron’s information capacity. Zeitschrift Für Naturforschung C. 36, 910-912 (1981)
  • [21] Naud, R., Marcille, N., Clopath, C. & Gerstner, W. Firing patterns in the adaptive exponential integrate-and-fire model. Biological Cybernetics. 99 pp. 335-347 (2008)
  • [22] Benda, J. & Herz, A. A Universal Model for Spike-Frequency Adaptation. Neural Computation. 15 pp. 2523-2564 (2003), https://api.semanticscholar.org/CorpusID:380557
  • [23] Gutierrez, G. & Denève, S. Population adaptation in efficient balanced networks. ELife. 8 pp. e46926 (2019)
  • [24] Salaj, D., Subramoney, A., Kraisnikovic, C., Bellec, G., Legenstein, R. & Maass, W. Spike frequency adaptation supports network computations on temporally dispersed information. Elife. 10 pp. e65459 (2021)
  • [25] Patterson, C., Wissig, S. & Kohn, A. Distinct effects of brief and prolonged adaptation on orientation tuning in primary visual cortex. Journal Of Neuroscience. 33, 532-543 (2013)
  • [26] Mitra, S., Fusi, S. & Indiveri, G. Real-time classification of complex patterns using spike-based learning in neuromorphic VLSI. IEEE Transactions On Biomedical Circuits And Systems. 3, 32-42 (2008)
  • [27] Nair, M. & Indiveri, G. An ultra-low power sigma-delta neuron circuit. 2019 IEEE International Symposium On Circuits And Systems (ISCAS). pp. 1-5 (2019)
  • [28] Chundi, P., Wang, D., Kim, S., Yang, M., Cerqueira, J., Kang, J., Jung, S., Kim, S. & Seok, M. Always-on sub-microwatt spiking neural network based on spike-driven clock-and power-gating for an ultra-low-power intelligent device. Frontiers In Neuroscience. 15 pp. 684113 (2021)