# Spike-timing-dependent plasticity learning of coincidence detection with passively integrated memristive circuits

M. Prezioso, M. R Mahmoodi, F Merrikh Bayat, H. Nili, H. Kim, Adrien F.

Vincent, D. B Strukov

► To cite this version:

M. Prezioso, M. R Mahmoodi, F Merrikh Bayat, H. Nili, H. Kim, et al. Spike-timing-dependent plasticity learning of coincidence detection with passively integrated memristive circuits. Nature Communications, 2018, 9 (1), 10.1038/s41467-018-07757-y. http://dx.doi.org/10.1038/s41467-018-07757-y

HAL Id: hal-02487839

https://hal.science/hal-02487839

Submitted on 21 Feb 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-018-07757-y

OPEN

# Spike-timing-dependent plasticity learning of coincidence detection with passively integrated memristive circuits

M. Prezioso<fnref n="1" />, M.R. Mahmoodi<fnref n="1" />, F. Merrikh Bayat<fnref n="1" />, H. Nili<fnref n="1" />, H. Kim<fnref n="1" />, A. Vincent<fnref n="1" /> & D.B. Strukov<fnref n="1" />

Spiking neural networks, the most realistic artificial representation of biological nervous systems, are promising due to their inherent local training rules that enable low-overhead online learning, and energy-efficient information encoding. Their downside is more demanding functionality of the artificial synapses, notably including spike-timing-dependent plasticity, which makes their compact efficient hardware implementation challenging with conventional device technologies. Recent work showed that memristors are excellent candidates for artificial synapses, although reports of even simple neuromorphic systems are still very rare. In this study, we experimentally demonstrate coincidence detection using a spiking neural network, implemented with passively integrated metal-oxide memristive synapses connected to an analogue leaky-integrate-and-fire silicon neuron. By employing spike-timing-dependent plasticity learning, the network is able to robustly detect the coincidence by selectively increasing the synaptic efficacies corresponding to the synchronized inputs. Not surprisingly, our results indicate that device-to-device variation is the main challenge towards realization of more complex spiking networks.

<footnote n="1">¹Electrical and Computer Engineering Department, UC Santa Barbara, Santa Barbara, CA 93106-9560, USA. These authors contributed equally: M. Prezioso, M. R. Mahmoodi, F. Merrikh Bayat. Correspondence and requests for materials should be addressed to D.B.S. (email: strukov@ece.ucsb.edu)</footnote>
The development of spiking neural networks (SNNs) has been largely driven by the hope of replicating the extremely high energy efficiency of biological systems. The main hypothesis behind the energy efficiency is the SNN's considerably higher information capacity (e.g., in terms of Vapnik-Chervonenkis dimension), compared to the more popular firing-rate-type neural networks such as multi-layer perceptron¹⁻³. Moreover, training for the firing-rate networks typically relies on a backpropagation algorithm, the efficient implementation of which is challenging owing in part to the centralized method for computing weight updates and the requirement for large high-precision memory. By contrast, most popular SNNs' weight updates rules are local, requiring only information from pre- and post-synaptic neurons (e.g., see Eq. 2 below), which could be a significant advantage for compact and low power implementations of real-time training, and scaling towards more complex networks⁴.

In the simplest SNN models, the information is encoded in spike-time correlations¹,², while the network functionality is defined by the neural link strengths, i.e., the synaptic efficacies, which are adjusted based on the relative timing of spikes that are passed via synapses. A key network element is a leaky-integrate-and-fire (LIF) neuron, whose operation is described by the equation:

$$C \frac{dU(t)}{dt} = \sum_{i} G_{i} V_{i}^{\text{pre}}(t) - \frac{1}{R_{L}} U(t), \quad (1)$$

where $U$, $C$, and $R_L$ mimic the membrane potential, capacitance, and its leakage resistance, correspondingly. The dot product in Eq. (1) is a current flowing into the neuron, which is proportional to the sum of the individual synaptic currents, i.e., the products of pre-synaptic spike voltages $V^{\text{pre}}(t)$ and corresponding synaptic conductances (weights) $G$. The second term on the right-hand side of Eq. (1) is a membrane leakage current. Note that $R_L C$ is membrane characteristic time, which determines the integration rate of the LIF neuron. Additionally, Eq. (1) does not explicitly include the resting potential. Instead, the resting potential is used as a reference for $U(t)$.

Another important SNN feature is spike-timing-dependent plasticity (STDP), which is a timing-dependent specialization of Hebbian learning¹,²,⁵,⁶. A typical goal of STDP learning is to strengthen the synaptic efficiency when two events happen in the expected causal temporal order, and to weaken it otherwise. STDP rule formally describes the change in synaptic weight as a specific function $f_{\text{STDP}}$ of a difference in firing times between pre-synaptic ($t^{\text{pre}}$) and post-synaptic ($t^{\text{post}}$) spikes, i.e.,

$$\Delta G = f_{\text{STDP}}(t^{\text{pre}} - t^{\text{post}}). \quad (2)$$

For example, a very common STDP rule found in biological synapses of neocortex layer 5 is shown in Fig. 1c⁶, according to which weight is increased (decreased) when pre-synaptic spike reaches the synapse before (after) the post-synaptic one.

Compact, low power hardware implementation of STDP learning is of particular importance, because practically valuable neural networks require a massive number of synapses. Owing to the extremely high density and analogue functionality, resistive switching memory devices ('memristors') are considered as one of the most attractive candidates for artificial synapses⁷. Indeed, crossbar circuit arrays with passively integrated memristors⁸,⁹ enable the most efficient hardware for dot-product computation¹⁰. Moreover, STDP learning in memristors can be implemented by applying specific bipolar pulses (see Fig. 1)¹¹.

There have been many experimental demonstrations of various STDP rules¹²⁻¹⁸ and more advanced synaptic functionalities based on memristors¹⁹⁻²¹ (also see comprehensive reviews in refs. ²²,²³). However, memristor implementations of even simple SNNs, which rely on STDP learning, are so far very scarce²⁴⁻²⁸. For example, Milo et al. used nine discrete memristors, connected in 1 transistor + 1 resistor (1T1R) configuration, to demonstrate learning of 3 × 3 patterns²⁴. In ref. ²⁵, Kim et al. reported SNN implementation based on CMOS-integrated 2T1R phase-change memory devices. Their demonstration did not include LIF neurons, so that the post-spikes were triggered externally after post-processing the neuron input. The focus of ref. ²⁶ was on emulating homeostasis using organic-based devices immersed in an electrolyte. However, the demonstrated coincidence detection, i.e. detection of the occurrence of temporally close but spatially distributed input signals, was performed without any learning. Ambrogio et al. reported the experimental STDP dynamics for hafnium-oxide memristors and used these results to simulate an SNN based on 1T1R discrete memristors²⁷. A more advanced demonstration of STDP-based spatial pattern recognition was demonstrated more recently by the same group, though still based on the same discrete 1T1R memristors²⁸.

Here, we experimentally demonstrate operation and STDP learning in SNN implemented with the most prospective, passively integrated (0T1R) memristive synapses connected to a silicon LIF neuron. Similar to ref. ²⁶, the focus of our work is on coincidence detection, i.e., the task of identifying correlated spiking activity. Coincidence detection functionality is found in various parts of nervous systems, such as auditory²⁹,³⁰ and visual³¹ cortices and is generally assumed to play a very important role in the brain³²⁻³⁴. In fact, entire neural systems have been hypothesized based purely on synchrony³⁵.

## Results

**Experimental setup.** The implemented SNN features 20 input neurons connected via 20 memristive crossbar-integrated synapses to a single LIF neuron (Fig. 2). Figure 3 shows particular voltage spike shapes employed in the coincidence detection experiment. (More details on the experimental setup and spike parameters are reported in Methods section.) The specific amplitudes and durations of voltage pulses were chosen to ensure balanced STDP windows, with approximately 50–200% maximum changes in conductance, for a majority of the devices. The post-synaptic spike has a larger amplitude negative pulse, which was needed owing to larger absolute values of reset switching thresholds compared to those of set transition.

In the coincidence detection task, STDP learning mechanism is employed for training a single neuron to fire an output spike when it receives simultaneous (i.e., correlated) spikes from several input neurons. Figure 2a shows an example of a spike pattern applied to the rows of the crossbars in the experiment. The correlated input spikes, i.e., rows 11 through 15 in Fig. 2a, were always synchronized and repeated every 840 ms, which represents one 'frame' of spiking input. Furthermore, in some experiments, the position of each synchronous spike was randomly shifted by $t_{\text{jitter}}$ to emulate temporal jitter. The uncorrelated ('noise') inputs were approximated by randomly generating spikes, with a probability $p_{\text{spike}}$ for each row in each frame. The position of these spikes within the frame was also generated probabilistically, by adding a random offset $\Delta_{\text{offset}}$ to the specific reference time, corresponding to the shifted spike phase—see Supplementary Figure 1 for more details.

**Coincidence detection results.** At the beginning of coincidence detection experiments, conductances of all memristors were set to approximately 16 µS (Figs. 4e–h, 5b) using an automatic tuning algorithm³⁶, which roughly corresponds to the middle of the device's dynamic range⁸,⁹. Figure 4 summarizes evolution of the
<figure><img src="image_3.png" /><figcaption>**Fig. 1** STDP implementation. Relative timing of the pre- and post-synaptic spikes for performing **a** potentiation and **b** depression of the synaptic weight for the memristor-based implementation¹¹, resulting in **c** the specific STDP, which is commonly observed in biology</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Spiking neural network setup. **a** Example of pattern of spikes applied to the inputs of the crossbar circuit. **b** Scanning electron top-view image of the utilized crossbar circuit. Scale bar is 2 µm. **c** LIF neuron implementation. See Methods section for more details on the crossbar circuit and neuron implementations</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 3** Signal waveforms. An example of measured oscilloscope signals in the lower-noise pattern coincidence experiment, showing **a** a pre-synaptic synchronous spike, that along with other synchronous inputs, led to the generation of the postsynaptic spike, and **b** an out-of-phase pre-synaptic synchronous spike with correspondingly increasing and then gradually leaking membrane potential</figcaption></figure>

system for the first set of experiments. The left panels show specific frames of spiking input, while the right panels show the trajectory of each synaptic weight (during 20 epochs or a total of 1.68 s of spiking activity for Fig. 4e–g and 60 epochs for Fig. 4h), with thick/thin lines corresponding to the synchronized/

uncorrelated inputs. (One epoch is defined as an application of ten frames. In each frame the spikes have the same reference positions with different random jitter and offsets.) Specifically, for all these experiments, $p_{spike} = 1$ so that the average spiking activity for all inputs was 1.19 Hz. To investigate the impact of
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Coincidence detection of higher-noise pattern. **a–d** Example of specific spiking input frame with **a** ($t_{\text{jitter}}_{\text{max}} = 0$, $\Delta_{\text{offset}}_{\text{max}} = 0$, **b** ($t_{\text{jitter}}_{\text{max}} = 27$ ms, $\Delta_{\text{offset}}_{\text{max}} = 27$ ms, **c** ($t_{\text{jitter}}_{\text{max}} = 27$ ms, $\Delta_{\text{offset}}_{\text{max}} = 216$ ms, and **d** ($t_{\text{jitter}}_{\text{max}} = 27$ ms, $\Delta_{\text{offset}}_{\text{max}} = 432$ ms and **e–h** corresponding results for time evolution of all devices’ conductances involved in the experiment. A single data point on right panels corresponds to the conductance measured at the end of a particular epoch. For all experiments, $p_{\text{spike}} = 1$. All inputs and corresponding conductances are shown using unique colours. Thicker lines show synchronous inputs/weights</figcaption></figure>

noise inputs, ($\Delta_{\text{offset}}_{\text{max}}$) was first set to 0 (Fig. 4a), which represents the simplest and more reproducible case with non-overlapping spikes. The ($\Delta_{\text{offset}}_{\text{max}}$) was then increased up to 432 ms (Fig. 4d), i.e., approximately 750% of the duration of pre-synaptic voltage spikes, in which case all 15 uncorrelated input spikes can fully overlap in time by chance and provide a current three times higher than the maximum current from the all-synchronous spikes.

Figure 5 shows coincidence detection results for another study, which demonstrate that the same neuron can learn a new pattern, while forgetting the previously learnt one. In this case, to simplify the experiment, $p_{\text{spike}}$ was set to 0.2, so that the average spiking activity was 0.238 Hz for each of the 15 noise inputs (Fig. 5a). ($\Delta_{\text{offset}}_{\text{max}}$) was set to 500 ms, so that, again, all of the noise spikes can overlap in time by chance. Note that the relative spiking activity of the noise inputs with respect to the synchronous ones was still much higher, and hence the considered task is more challenging, compared to those studied in ref. 28.

Figure 3 provides more details on the dynamics of the system by showing two representative voltage waveforms recorded by the oscilloscope. The recorded data in Fig. 3a shows that the pre-synaptic synchronous pulse and post-synaptic pulse overlapped in time. The former pulse was sent earlier in time, causing membrane potential to reach a firing threshold and, subsequently, generate the post-synaptic spike, which would be required to strengthen the synaptic efficacy. Obviously, for the considered task, this potentiation condition has happened more frequently and consistently for the synchronized inputs, so that their synaptic efficacies were reinforced. Figure 3b also shows an example of the synchronous pre-synaptic spike that was out of phase with other synchronous inputs. In this case, the membrane potential was below the firing threshold and leaked over time.

Finally, the additional results show that the coincidence detection experiment can be successfully performed multiple times using the same column (Fig. 6) and using different columns in the crossbar (Fig. 7).
<figure><img src="image_7.png" /><figcaption>**Fig. 5** Coincidence detection of two lower-noise patterns. **a** A spiking input with $p_{spike} = 0.2$, $(t_{jitter})_{max} = 5$ ms, $(\Delta_{offset})_{max} = 500$ ms. The first pattern of synchronous spikes was applied to the rows #11 through #15 in the first 30 epochs, while the second pattern was applied to rows #1 through #5 in the last 30 epochs. **b** Time evolution of all devices' conductances involved in the experiment. **c** Post-processed results from panel **b**, showing time evolution of the total and the averaged conductances for all devices and the averages of subsets corresponding to the synchronized (pattern #1 and pattern #2) and noise inputs. The colour coding in panels **a** and **b** are similar to those of Fig. 4</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 6** Trial statistics. Coincidence detection experiment repeated 15 times when connecting LIF neuron to **a** column #2 with $p_{spike} = 1$, with no jitter or offset added; and **b** column #19 with $p_{spike} = 1$, $(t_{jitter})_{max} = 25$ ms, and $(\Delta_{offset})_{max} = 200$ ms. The noise spikes were generated randomly and were unique for each frame in all trials</figcaption></figure>

## Discussion

In both sets of experiments, the synaptic weights did not always evolve in the proper direction e.g., synchronized input #15 in Fig. 4. Additionally, the conductance change is not monotonic, with clearly noticeable jumps for some of synapses. This is partially due to random post-synaptic spikes generated by the

output neuron in response to noise, but also due to non-idealities in memristors, most importantly their device-to-device variations in dynamic current-voltage characteristics. Nevertheless, the results in all experiments clearly show the progressive cumulative potentiation of the synapses that correspond to the synchronized inputs.
For example, Fig. 5c shows that when the second pattern was first applied (beginning of the epoch 31) the average synaptic strength corresponding to the synchronized inputs was initially

<figure><img src="image_9.png" /><figcaption>**Fig. 7** Column statistics. The relative change in device conductances after 20 epochs of applying synchronous input to rows #11 through #15, when connecting LIF neuron to a specified column in the crossbar. The experiments were performed separately for each column. $p_{spike} = 1$, $(t_{jitter})_{max} = 5$ ms, $(\Delta_{offset})_{max} = 27$ ms. The noise spikes were always generated randomly for each row in all frames</figcaption></figure>

well below the average across all of the devices and more than two times lower than that of the uncorrelated inputs. However, as the training continued, the network was able to learn to discriminate those synchronized inputs by increasing their weights, eventually reaching and surpassing the average strength of uncorrelated inputs. In fact, by the end of the training, even the cumulative strength for the uncorrelated inputs was substantially lower compared to the synchronized ones for both patterns (Fig. 5c). In addition, as expected, the convergence took longer for the truly random noise input (Fig. 4h), because of the higher probability of firing a post-synaptic spike based on the spiking activity on the uncorrelated inputs (cf. with Fig. 4e).

Our experimental work has confirmed one of the main challenges for SNN implementation with memristors—their device-to-device variations³⁷,³⁸. Indeed, given similar pre- and post-synaptic spikes, the conductance updates should ideally follow the same STDP window for all synapses. In reality, the updates vary significantly across memristive synapses (Fig. 8), primarily due to variations in switching thresholds of the devices. Indeed, SNNs typically rely on fixed-magnitude spikes. In SNN's crossbar circuit implementation, the fixed-magnitude spikes are intended to update in parallel weights of multiple devices. In our experiment, the specific spike amplitudes were selected according to the average switching thresholds across the devices. Consequently, the change in the conductances for the devices with

<figure><img src="image_10.png" /><figcaption>**Fig. 8** STDP variations. Measured STDP device-to-device variations for 20 memristors. The vertical axis shows the relative change in the conductance for a particular time difference between pre and post-synaptic spikes. (Note the different vertical axis scale for the panels.) The top label shows the device index, counted from the top of the crossbar circuit. Similar to ref. ¹², the device conductance was always tuned to approximately 16 µS before measuring each data point</figcaption></figure>
larger switching thresholds is naturally smaller, e.g., device #18 in Fig. 8, and alternatively, larger for those with smaller threshold (device #19). For the same reason, some STDP windows appeared to be noisy, e.g., that of device #15 in Fig. 8. In this respect, memristor implementation of the simpler ex situ trained firing-rate neural networks is much less challenging, as the write amplitude voltages in these networks can be adjusted uniquely for each device based on the feedback information during conductance tuning⁹,³⁸.

The STDP window for many memristors is also not balanced with either dominant potentiation (such as device #6 in Fig. 8) or depression (device #16). This is in part due to device-to-device variations in the on/off dynamic range—see, e.g., impact of the initial state on STDP function shape in ref. ¹². However, these variations only result in slight variation in the synaptic equilibrium currents, i.e., the memristor states reached after application of random spiking input¹², and hence, are not very critical for the network operation.

Recently, several groups reported promising devices with very uniform switching characteristics³⁹,⁴⁰, that would naturally alleviate the device variation issue. An even better solution would be to holistically optimize training algorithms, circuits, and architectures with respect to device variations. For example, learning rules based on stochastic binary synapses⁴¹,⁴² would be more robust to the non-ideal analogue switching. The device variations could also be addressed by considering more complex implementation of synapses²²,⁴³ and/or having redundant devices in the network⁴⁴, though at the expense of larger area and power consumption.

In fact, our choice of the rectangular-shaped pulses utilized for pre- and post-synaptic voltage spikes is one example of the holistic optimization. Such shapes were selected in part for their simplicity of implementation but, more importantly, to ensure more reproducible STDP windows. In principle, the use of triangular pulses (e.g., those shown in Fig. 1a, b and implemented in ref. ¹²) would cause smoother STDP curves, i.e., with more gradual changes in conductance. However, we found that such schemes resulted in less reproducible STDP behaviour when there are significant device variations. This is likely due to the specific switching dynamics of the considered devices (which is representative for many metal-oxide memristors⁷,³⁶). Specifically, when the spikes are implemented with triangular pulses, the dropped voltage over the device also has a triangular waveform (see, e.g., Fig. 2a from ref. ¹²), and hence maximum voltage is applied only over a short period of time. By contrast, with rectangular pulse implementation, the maximum voltage is applied over a longer time, which is comparable to pulse duration when pre- and post-synaptic spikes are fully overlapped in time. Because of quasi DC operation, i.e., relatively long pulse durations compared to the intrinsic switching times, and exponential switching dynamics with respect to the applied voltage, which saturates in time, the latter scheme seems to be more tolerant to memristor variations.

In conclusion, we have developed a complete SNN in which artificial synapses were implemented with passive integrated metal-oxide memristors, while LIF neuron was realized with discrete conventional semiconductor circuits. A specific but representative task for SNN, i.e., coincidence detection, which involves training a neuron to discriminate between synchronized inputs from uncorrelated ones, was successfully demonstrated by employing spike-timing-dependent plasticity learning. Our work showed that device-to-device variation in memristors' switching thresholds is the major challenge. Finally, we discussed several approaches for overcoming this challenge to build more practical, efficient memristor-based SNNs.

## Methods

**Spiking neural network implementation.** Artificial synapses are implemented with Pt/Al₂O₃/TiO₂⁻/Pt memristors, which were passively integrated in a 20 × 20 crossbar array (Fig. 2b). (A detailed discussion of the fabrication methods and electrical characterization results can be found in refs. ⁸,⁹.) The crossbar array was packaged (Supplementary Figure 2c) and inserted in a custom-printed circuit board (Supplementary Figure 2a) that allows connecting the crossbar lines either to the input/output neurons during network operation or to a switch matrix, which in turn is connected to a Keysight B1500A parameter analyser, for device forming, testing, and conductance tuning.

The input neurons, one for each row of the crossbar array, were implemented using off-the-shelf digital-to-analogue converter circuits. The output LIF neuron was connected to the third column of the crossbar array for the experiments reported in Figs. 3–5. (Other columns of the crossbar were grounded.) The neuron was implemented with a combination of another dedicated custom-printed circuit board (Supplementary Figure 2b), whose circuitry is shown in Fig. 2c, and an arbitrary waveform generator. The design of the neuron allowed fine-tuning its characteristic time, membrane potential threshold and the scale of synaptic current by adjusting the neuron's variable resistors. Operational amplifiers, TL074CN, were used in all three stages, while typical neuron circuit parameters are as following: $R_1 = 5.6 \text{ K}\Omega$, $R_2 = 2.5 \text{ K}\Omega$, $R_L = 900 \text{ K}\Omega$, $V_{b1} = 10 \text{ mV}$, $V_{b2} = 0.8 \text{ V}$ and $C = 0.1 \text{ }\mu\text{F}$.

Upon reaching a threshold, an LIF neuron triggers an arbitrary waveform generator. The latter was used to produce the post-synaptic spikes and, by using its 'sync' output, which is a transistor-transistor logic pulse synchronized with the signal output, to temporarily connect its signal output to the crossbar column and reset the neuron membrane potential. The entire experimental setup was controlled from a personal computer using a custom C program. The program controlled the generation of the pre-synaptic spikes by digital-to-analogue converters, an additional oscilloscope which was used to measure the inputs, post-synaptic spikes, and neuron membrane potential, and a switch matrix and parameter analyser.

**Spike pulse details.** For the experimental results shown in Figs. 3–7, the utilized pre-synaptic voltage spike consists of 5-ms ramp to +0.45 V, 7-ms 0.63 V plateau, 1-ms ramp down to -0.4 V, 20-ms -0.4 V plateau and then 5-ms ramp up to 0 V (Fig. 2). The post-synaptic voltage spike was 5-ms ramp from -0.15 V to +0.65 V, 23-ms +0.6 V plateau, 1-ms ramp down to -1.15 V, 22-ms -1.15 V plateau, and 5-ms ramp up to -0.15 V (Fig. 2a). The constant -0.15 V in the post-synaptic spike was for convenience of neuron implementation and not essential to the network operation. The total durations of pre-synaptic and post-synaptic pulses were 38 and 57 ms, respectively.

The results in Fig. 8 were obtained using a different 20 × 20 crossbar circuit, which was implemented using the same fabrication process. For this particular experiment, the pre-synaptic voltage spike shape was 5 ms ramp to +0.63 V, 10 ms to 0.63 V plateau, 2 ms ramp down to -0.54 V, 20 ms to -0.54 V plateau and then 5 ms ramp up to 0 V. The post-synaptic voltage spike shape was 5 ms ramp to +0.60 V, 23 ms to +0.6 V plateau, 2 ms ramp down to -0.5 V, 22 ms to -0.5 V plateau, and 5 ms ramp up to 0 V. The total durations of pre-synaptic and post-synaptic pulses were 42 and 57 ms, respectively.

## Data availability

The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.

Received: 12 June 2018 Accepted: 25 November 2018

Published online: 14 December 2018

## References

1. Gerstner, W. & Kistler, W. *Spiking Neural Models* (Cambridge University Press, Cambridge, MA, 2002).

2. Maass, M. Networks of spiking neurons: the third generation of neural network models. *Neural Netw.* **10**, 1659–1671 (1997).

3. Maass, M. & Schmitt, M. On the complexity of learning for spiking neurons with temporal coding. *Inf. Comput.* **153**, 26–46 (1999).

4. Merolla, P. A. et al. A million spiking-neuron integrated circuit with a scalable communication network and interface. *Science* **345**, 668–673 (2014).

5. Caporale, N. & Dan, Y. Spike timing-dependent plasticity: a Hebbian learning rule. *Annu. Rev. Neurosci.* **31**, 25–46 (2008).

6. Bi, G. & Poo, M. Synaptic modifications in cultured hippocampal neurons: dependence on spike timing, synaptic strength, and postsynaptic cell type. *J. Neurosci.* **18**, 10464–10472 (1998).

7. Yang, J. J., Strukov, D. B. & Stewart, D. R. Memristive devices for computing. *Nat. Nanotechnol.* **8**, 13–24 (2013).

8. Prezioso, M. et al. Training and operation of an integrated neuromorphic network based on metal-oxide memristors. *Nature* **521**, 61–64 (2015).
9. Merrikh Bayat, F. et al. Implementation of multilayer perceptron network with highly uniform passive memristive crossbar circuits. *Nat. Commun.* **9**, 2331 (2018).

10. Ceze, L. et al. Nanoelectronic neurocomputing: status and prospects. In *Proc. Device Research Conference 1-2* (IEEE and MRS, Newark, DE, 2016).

11. Linares-Barranco, B. & Serrano-Gotarredona, T. Memristance can explain spike-timing-dependent plasticity in neural synapses. *Nat. Preced.* http://precedings.nature.com/documents/3010/version/1 (2009).

12. Prezioso, M., Merrikh Bayat, F., Hoskins, B. D., Likharev, K. K. & Strukov, D. Self-adaptive spike-time-dependent plasticity of metal-oxide memristors. *Sci. Rep.* **6**, 21331 (2016).

13. Jo, S. H. et al. Nanoscale memristor device as synapse in neuromorphic systems. *Nano Lett.* **10**, 1297-1301 (2010).

14. Shimeng, Y., Yi, W., Jeyasingh, R., Kuzum, D. & Wong, H. S. P. An electronic synapse device based on metal oxide resistive switching memory for neuromorphic computation. *IEEE Trans. Electron. Devices* **58**, 2729-2737 (2011).

15. Kuzum, D., Jeyasingh, R. G., Lee, B. & Wong, H.-S. P. Nanoelectronic programmable synapses based on phase change materials for brain-inspired computing. *Nano Lett.* **12**, 2179-2186 (2011).

16. Alibart, F. et al. An organic nanoparticle transistor behaving as a biological spiking synapse. *Adv. Funct. Mater.* **20**, 330-337 (2010).

17. Li, S. et al. Synaptic plasticity and learning behaviours mimicked through Ag interface movement in an Ag/conducting polymer/Ta memristive system. *J. Mater. Chem. C* **1**, 5292-5298 (2013).

18. Subramaniam, A. et al. Spike-timing-dependent plasticity using biologically realistic action potentials and low-temperature materials. *IEEE Trans. Nanotechnol.* **12**, 450-459 (2013) (2015).

19. Ohno, T. et al. Short-term plasticity and long-term potentiation mimicked in single inorganic synapses. *Nat. Mater.* **10**, 591-595 (2011).

20. Kim, S. et al. Experimental demonstration of a second-order memristor and its ability to bioreistically implement synaptic plasticity. *Nano Lett.* **15**, 2203-2211 (2015).

21. Wang, Z. et al. Memristors with diffusive dynamics as synaptic emulators for neuromorphic computing. *Nat. Mater.* **16**, 101-108 (2017).

22. Serrano-Gotarredona, T. et al. STDP and STDP variations with memristors for spiking neuromorphic learning systems. *Front. Neurosci.* **7**, 2 (2013).

23. Saighi, S. et al. Plasticity in memristive devices for spiking neural network. *Front. Neurosci.* **9**, 51 (2015).

24. Milo, V. et al. Demonstration of hybrid CMOS/RRAM neural networks with spike time/rate-dependent plasticity. In *Proc. International Electron Devices Meeting* 16.8.1-16.8.4 (IEEE, San Francisco, CA, 2016).

25. Kim, S. et al. NVM neuromorphic core with 64k-cell (256-by-256) phase change memory synaptic array with on-chip neuron circuits for continuous in-situ learning. In *Proc. International Electron Devices Meeting* 17.1.1-17.1.4 (IEEE, Washington, DC, 2015).

26. Gkoupidenis, P., Koutsouras, D. A. & Malliaras, G. G. Neuromorphic device architectures with global connectivity through electrolyte gating. *Nat. Commun.* **8**, 15448 (2017).

27. Ambrogio, S. et al. Neuromorphic learning and recognition with one-transistor-one-resistor synapses and bistable metal oxide RRAM. *IEEE Trans. Electron. Devices* **63**, 1508-1515 (2016).

28. Pedretti, G. et al. Memristive neural network for on-line learning and tracking with brain-inspired spike timing dependent plasticity. *Sci. Rep.* **7**, 5288 (2017).

29. Agmon-Snir, H., Carr, C. E. & Rinzel, J. The role of dendrites in auditory coincidence detection. *Nature* **393**, 268-272 (1998).

30. Oertel, D., Bal, R., Gardner, S. M., Smith, P. H. & Joris, P. X. Detection of synchrony in the activity of auditory nerve fibers by octopus cells of the mammalian cochlear nucleus. *Proc. Natl. Acad. Sci. USA* **97**, 11773-11779 (2000).

31. König, P., Engel, A. K. & Singer, W. Integrator or coincidence detector? The role of the cortical neuron revisited. *Trends Neurosci.* **19**, 130-137 (1996).

32. Softky, W. R. & Koch, C. The highly irregular firing of cortical cells is inconsistent with temporal integration of random EPSPs. *J. Neurosci.* **13**, 334-350 (1993).

33. Mainen, Z. F. & Sejnowski, T. J. Reliability of spike timing in neocortical neurons. *Science* **268**, 1503-1506 (1995).

34. Stevens, C. F. & Zador, A. M. Input synchrony and the irregular firing of cortical neurons. *Nat. Neurosci.* **1**, 210-217 (1998).

35. Shastri, L. & Ajjanagadde, V. From simple associations to systematic reasoning: a connectionist representation of rules, variables and dynamic bindings using temporal synchrony. *Behav. Brain Sci.* **16**, 417-451 (1993).

36. Alibart, F., Gao, L., Hoskins, B. & Strukov, D. B. High-precision tuning of state for memristive devices by adaptable variation-tolerant algorithm. *Nanotechnology* **23**, 075201 (2012).

37. Dongale, T. D. et al. Investigation of process parameter variation in the memristor based resistive random access memory (RRAM): Effect of device size variations. *Mater. Sci. Semicond. Process.* **35**, 174-180 (2015).

38. Strukov, D. B. Tightening grip. *Nat. Mater.* **17**, 293-295 (2018).

39. Choi, S. et al. SiGe epitaxial memory for neuromorphic computing with reproducible high performance based on engineered dislocations. *Nat. Mater.* **17**, 335-340 (2018).

40. Goswami, S. et al. Robust resistive memory devices using solution-processable metal-coordinated azo aromatics. *Nat. Mater.* **16**, 1216-1224 (2017).

41. Ma, X. & Likharev, K. K. Global reinforcement learning in neural networks. *IEEE Trans. Neural Netw.* **18**, 573-577 (2006).

42. Suri, M. et al. Bio-inspired stochastic computing using binary CBRAM synapses. *IEEE Trans. Electron Devices* **60**, 2402-2409 (2013).

43. Indiveri, G., Corradi, F. & Qiao, N. Neuromorphic architectures for spiking deep neural networks. In *Proc. IEEE International Electron Devices Meeting* 4.2.1-4.2.4 (IEEE, Washington, DC, 2015).

44. Lee, J. H. & Likharev, K. K. Defect-tolerant nanoelectronic pattern classifiers. *Int. J. Circ. Theor. Appl.* **35**, 239-264 (2007).

## Acknowledgements

This work was supported by ARO under contract #W91NF-16-1-0302 and DARPA under contract HR00113-13-C-0051 UPSIDE via BAE Systems, Inc.

## Author contributions

M.P., F.M.B., and D.S. conceived the original concept and initiated the work. M.P., H.N., and H.K. fabricated memristors. F.M.B., M.P., and A.V. developed the experimental setup. M.M., M.P., and H.K. performed measurements. D.S. and M.P. wrote the manuscript. All authors discussed the results.

## Additional information

Supplementary Information accompanies this paper at https://doi.org/10.1038/s41467-018-07757-y.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Journal peer review information:** *Nature Communications* thanks the anonymous reviewers for their contribution to the peer review of this work.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018