<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Neurosci.</journal-id>
<journal-title>Frontiers in Neuroscience</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Neurosci.</abbrev-journal-title>
<issn pub-type="epub">1662-453X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fnins.2018.00057</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Spike-Timing Dependent Plasticity in Unipolar Silicon Oxide RRAM Devices</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Zarudnyi</surname> <given-names>Konstantin</given-names></name>
</contrib>
<contrib contrib-type="author">
<name><surname>Mehonic</surname> <given-names>Adnan</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/246152/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Montesi</surname> <given-names>Luca</given-names></name>
</contrib>
<contrib contrib-type="author">
<name><surname>Buckwell</surname> <given-names>Mark</given-names></name>
</contrib>
<contrib contrib-type="author">
<name><surname>Hudziak</surname> <given-names>Stephen</given-names></name>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Kenyon</surname> <given-names>Anthony J.</given-names></name>
<xref ref-type="author-notes" rid="fn001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/141040/overview"/>
</contrib>
</contrib-group>
<aff><institution>Department of Electronic and Electrical Engineering, University College London</institution>, <addr-line>London</addr-line>, <country>United Kingdom</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Gert Cauwenberghs, Jacobs School of Engineering and Institute for Neural Computation, University of California, San Diego, United States</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Sadique Sheik, BioCircuits Institute, University of California, San Diego, United States; Duygu Kuzum, Department of Electrical and Computer Engineering, University of California, San Diego, United States</p></fn>
<fn fn-type="corresp" id="fn001"><p>&#x0002A;Correspondence: Anthony J. Kenyon <email>a.kenyon&#x00040;ucl.ac.uk</email></p></fn>
<fn fn-type="other" id="fn002"><p>This article was submitted to Neuromorphic Engineering, a section of the journal Frontiers in Neuroscience</p></fn></author-notes>
<pub-date pub-type="epub">
<day>08</day>
<month>02</month>
<year>2018</year>
</pub-date>
<pub-date pub-type="collection">
<year>2018</year>
</pub-date>
<volume>12</volume>
<elocation-id>57</elocation-id>
<history>
<date date-type="received">
<day>11</day>
<month>07</month>
<year>2017</year>
</date>
<date date-type="accepted">
<day>23</day>
<month>01</month>
<year>2018</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2018 Zarudnyi, Mehonic, Montesi, Buckwell, Hudziak and Kenyon.</copyright-statement>
<copyright-year>2018</copyright-year>
<copyright-holder>Zarudnyi, Mehonic, Montesi, Buckwell, Hudziak and Kenyon</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>Resistance switching, or Resistive RAM (RRAM) devices show considerable potential for application in hardware spiking neural networks (neuro-inspired computing) by mimicking some of the behavior of biological synapses, and hence enabling non-von Neumann computer architectures. Spike-timing dependent plasticity (STDP) is one such behavior, and one example of several classes of plasticity that are being examined with the aim of finding suitable algorithms for application in many computing tasks such as coincidence detection, classification and image recognition. In previous work we have demonstrated that the neuromorphic capabilities of silicon-rich silicon oxide (SiO<sub>x</sub>) resistance switching devices extend beyond plasticity to include thresholding, spiking, and integration. We previously demonstrated such behaviors in devices operated in the unipolar mode, opening up the question of whether we could add plasticity to the list of features exhibited by our devices. Here we demonstrate clear STDP in unipolar devices. Significantly, we show that the response of our devices is broadly similar to that of biological synapses. This work further reinforces the potential of simple two-terminal RRAM devices to mimic neuronal functionality in hardware spiking neural networks.</p></abstract>
<kwd-group>
<kwd>resistive switching</kwd>
<kwd>resistance switching</kwd>
<kwd>STDP</kwd>
<kwd>RRAM</kwd>
<kwd>machine learning</kwd>
<kwd>neuromorphic systems</kwd>
</kwd-group>
<contract-num rid="cn001">EP/K01739X/1</contract-num>
<contract-sponsor id="cn001">Engineering and Physical Sciences Research Council<named-content content-type="fundref-id">10.13039/501100000266</named-content></contract-sponsor>
<counts>
<fig-count count="5"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="35"/>
<page-count count="8"/>
<word-count count="4858"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Non-von Neumann computing architecture, inspired by the neuronal architecture of the brain, is receiving considerable interest as a more efficient way to tackle a range of machine learning tasks such as classification, sorting, and image recognition (Indiveri et al., <xref ref-type="bibr" rid="B9">2006</xref>; Smith, <xref ref-type="bibr" rid="B28">2006</xref>; Izhikevich and Edelman, <xref ref-type="bibr" rid="B12">2008</xref>; Ananthanarayanan et al., <xref ref-type="bibr" rid="B3">2009</xref>; Jo et al., <xref ref-type="bibr" rid="B13">2010</xref>). Such work is not limited to purely software implementations: hardware &#x0201C;neuromorphic&#x0201D; systems, which mimic specific functionalities of biological neurons, most often synaptic plasticity, have been increasingly studied in recent years (Mead, <xref ref-type="bibr" rid="B15">1990</xref>; Indiveri, <xref ref-type="bibr" rid="B8">2000</xref>; Le Masson et al., <xref ref-type="bibr" rid="B14">2002</xref>; Vogelstein et al., <xref ref-type="bibr" rid="B31">2008</xref>; Mitra et al., <xref ref-type="bibr" rid="B23">2009</xref>; Indiveri et al., <xref ref-type="bibr" rid="B10">2011</xref>).</p>
<p>The biological synapse can be naively modeled as a two-terminal structure that permits communication between two neurons via electrical or chemical signals (Chua, <xref ref-type="bibr" rid="B6">1971</xref>; Strukov et al., <xref ref-type="bibr" rid="B29">2008</xref>; Schacter et al., <xref ref-type="bibr" rid="B26">2011</xref>) and whose connection strength depends on electrical history, encoded in the dynamics of ion diffusion. Likewise, Resistive Random Access Memory devices (Waser and Aono, <xref ref-type="bibr" rid="B33">2007</xref>), often described as memristors (Chua, <xref ref-type="bibr" rid="B6">1971</xref>), which are a promising novel non-volatile memory technology (Torrezan et al., <xref ref-type="bibr" rid="B30">2011</xref>; Chen et al., <xref ref-type="bibr" rid="B5">2012</xref>; Mehonic et al., <xref ref-type="bibr" rid="B18">2012a</xref>), are also two-terminal devices whose conductance can be modified dynamically by changes in both ion and electron current flow.</p>
<p>Spike timing dependent plasticity is a form of Hebbian learning rule that results in synaptic connection strength modification based on the relative timing of voltage spikes between neurons (Hebb, <xref ref-type="bibr" rid="B7">1949</xref>; Bi and Poo, <xref ref-type="bibr" rid="B4">1998</xref>; Abbott et al., <xref ref-type="bibr" rid="B1">2000</xref>), more details of which are given below. In recent years, there have been multiple reports of RRAM devices showing STDP behavior in a range of bipolar devices and materials (Jo et al., <xref ref-type="bibr" rid="B13">2010</xref>; Yu et al., <xref ref-type="bibr" rid="B35">2011</xref>, <xref ref-type="bibr" rid="B34">2013</xref>; Ambrogio et al., <xref ref-type="bibr" rid="B2">2013</xref>; Indiveri et al., <xref ref-type="bibr" rid="B11">2013</xref>; Sa&#x000EF;ghi et al., <xref ref-type="bibr" rid="B25">2015</xref>; Serb et al., <xref ref-type="bibr" rid="B27">2016</xref>). In most studies, memristive devices are operated purely in bipolar mode, in which potentiation and depression of synaptic strength are governed by positive and negative voltage pulses, respectively. However, in order to implement other neuronal functions, including thresholding, spiking and integration, more complex CMOS circuitry is required in such systems to interface with the RRAM &#x0201C;synapse.&#x0201D; Nevertheless, in our previous work we have demonstrated how this broader class of functions can be realized in unipolar SiO<sub>x</sub> RRAM devices (Mehonic et al., <xref ref-type="bibr" rid="B16">2015</xref>). In the majority of other reported examples of the use of RRAM devices in neuromorphic systems, the role of the RRAM element is to model the behavior of the biological synapse&#x02014;that is, to provide a programmable synaptic weight in the form of an adjustable conductance. As is the case with the biological system, this has most conveniently been achieved by operating in the bipolar mode. However, oxide-based RRAM devices exhibit interesting and useful behavior when operated in a unipolar mode. In the case of SiO<sub>x</sub>, devices biased with a constant current can enter a metastable state in which field-driven filament formation competes with current-driven filament (partial) dissolution. Adjusting the bias correctly enables us to generate voltage spikes with a current-controlled frequency, removing the need for much of the CMOS electronics used in conventional systems to model neuronal function. While much work remains to be done in this field, the reduction of CMOS circuitry footprint in neuromorphic systems could yield savings in valuable silicon chip real-estate. Other benefits such as power reduction are more uncertain, and remain so until further work has been performed on algorithm development and a fuller understanding of how devices can be designed to operate at very low current and/or duty cycle. Nevertheless, a reduction in circuit complexity by exploiting simple two-terminal devices is an attractive goal. Building on our previous work, here we report, for the first time, spike timing dependent plasticity in such unipolar devices&#x02014;a result that, when combined with neuronal functionality, promises a greatly reduced component count for RRAM hardware spiking neural networks. Further, we present a statistical analysis of STDP conductance modification mimicking the response of a biological synapse. RRAM devices in our experiments are based on SiO<sub>x</sub> switching layers; this ensures full compatibility with CMOS technology (Mehonic et al., <xref ref-type="bibr" rid="B18">2012a</xref>). This study serves as a proof of principle of device functionality; energy efficiency will be studied in our future work.</p></sec>
<sec sec-type="materials and methods" id="s2">
<title>Materials and methods</title>
<p>All experiments were carried out on metal/SiO<sub>x</sub>/metal devices whose structure is shown schematically in <bold>Figure 2A</bold>. One hundred nanometers of thick titanium nitride (TiN) electrodes sandwich a 37 nm thick SiO<sub>x</sub> resistance switching layer. All layers were sputter deposited (the SiO<sub>x</sub> layer was co-sputtered from Si and SiO<sub>2</sub> targets). The SiO<sub>x</sub> active switching layer was amorphous and sub-stoichiometric (x &#x02248; 1.3). The top TiN layer was patterned using conventional lithography techniques into square electrode structures with edges ranging from 100 to 400 &#x003BC;m (for more details of fabrication see Mehonic et al., <xref ref-type="bibr" rid="B16">2015</xref>; Montesi et al., <xref ref-type="bibr" rid="B24">2016</xref>).</p>
<p>Devices were operated in ambient conditions, exhibiting intrinsic switching through the formation of oxygen vacancy filaments, and may therefore be classed as valence change memory (VCM) systems (Mehonic and Kenyon, <xref ref-type="bibr" rid="B20">2015</xref>). Pristine devices require an electroforming step to generate a nanometer-scale conductive oxygen vacancy filament, after which devices may be cycled many times between multiple resistance states (Waser, <xref ref-type="bibr" rid="B32">2012</xref>). During electroforming, current was limited to avoid destructive breakdown of the oxide. Following electroforming, the device enters a low resistance state (LRS). To make a transition from the LRS to a high resistance state (HRS) a lower bias voltage, with no current compliance, is applied to re-oxidize a section of the conductive filament by increasing oxygen ion mobility via Joule heating. To reinstate the filament connection between the anode and cathode a localized process similar to electroforming can be driven by the application of an appropriate field, limited to avoid currents sufficiently high to damage the oxide through catastrophic dielectric breakdown. Further details of the resistance switching mechanism are found in Mehonic et al. (<xref ref-type="bibr" rid="B18">2012a</xref>), and a study of the role of oxygen movement in electroforming and switching is in Mehonic et al. (<xref ref-type="bibr" rid="B17">2016</xref>).</p>
<p>While characterization of RRAM devices is often carried out using voltage sweeps, in normal operation devices are programmed using short voltage pulses, which may be of the order of nanoseconds in duration. In this case, the major control factors for pulsed electroforming are duration and voltage amplitude of the pulse. We used a Keithley 4200-SCS Semiconductor Characterization System with a Pulse Measure Unit (PMU) to create custom pulse shapes. The PMU does not allow manual control of current compliance, which we implemented instead by carefully controlling pulse voltage and duration. We used custom Matlab scripts to control the system either directly or through custom C drivers. In all cases we contacted devices with tungsten needle probes on a Signatone probe station. Devices were tested in ambient conditions in an isolated low-light environment.</p></sec>
<sec id="s3">
<title>Results and discussion</title>
<sec>
<title>Unipolar resistance switching</title>
<p>Details of resistance switching in intrinsic oxide RRAM devices may be found elsewhere (Mehonic et al., <xref ref-type="bibr" rid="B18">2012a</xref>,<xref ref-type="bibr" rid="B19">b</xref>, <xref ref-type="bibr" rid="B16">2015</xref>, <xref ref-type="bibr" rid="B17">2016</xref>, <xref ref-type="bibr" rid="B22">2017</xref>) but to briefly summarize, application of an initial electroforming field, limited to prevent catastrophic dielectric breakdown, across a pristine, highly insulating, oxide generates a filament of conductive oxygen vacancies bridging the oxide film. Subsequent application of a field, either of the opposite polarity (bipolar operation) or similar polarity (unipolar) without a current limit resets the device to a high resistance state. Applying a field of the original polarity (bipolar operation) or the same polarity with a current limit (unipolar) sets the device back to the low resistance state. Devices may typically be cycled between the high and low resistance states many times.</p>
<p>Figure <xref ref-type="fig" rid="F1">1A</xref> shows a typical unipolar resistance switching of the devices we use in this study. The electroforming process typically occurs around 6 V, while the set and the reset voltages are typically 4 and 2.5 V, respectively. For the electroforming and the set process, a current limit of 3 mA is used to prevent hard breakdown. Similar to other RRAM devices, SiO<sub>x</sub> unipolar devices show a variability of the switching voltages. This is shown in the cumulative plots in Figure <xref ref-type="fig" rid="F1">1B</xref>, in which data is obtained from more than 1,000 I/V sweeps (top contacts size 200 &#x000D7; 200 &#x003BC;m).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p><bold>(A)</bold> Typical I/V curve for devices used in this study showing resistance switching in the unipolar mode. The current compliance limit is labeled as CC. <bold>(B)</bold> Cumulative probability plots of set and reset voltages obtained from more than 1,000 switching cycles (during the I/V sweeps).</p></caption>
<graphic xlink:href="fnins-12-00057-g0001.tif"/>
</fig></sec>
<sec>
<title>Setting and resetting dependence on pulse shape</title>
<p>Devices may be set or reset depending on the shape of voltage pulses. Devices may be set by an ascending sequence of 300 ns voltage pulses, which, over a sufficiently long time, look like the leading edge of a voltage pulse, and then they may be reset in response to a descending sequence of pulses that look like the trailing edge of a triangular or sawtooth pulse. The gradual resistance change can be tracked by applying short read voltage pulses (0.6 V, 300 ns) that do not induce any switching (shown in in <bold>Figure 4</bold>). These pulses were applied in between each of the ascending or descending pulses that make up the sawtooth voltage ramp (Figures <xref ref-type="fig" rid="F2">2B,C</xref>&#x02014;shown as a continuous blue line for clarity) while the current supplied by the Keithley system was measured, and hence the device resistance obtained. The transition from LRS to HRS begins to occur at around 2.7 V, a voltage that agrees well with the statistical measurements we carried out shown in Figure <xref ref-type="fig" rid="F2">2B</xref>. The ascending pulse train can then set the device to a LRS (Figure <xref ref-type="fig" rid="F2">2C</xref>). Resistance gradually decreases as the voltage bias increases, which demonstrates the system&#x00027;s capacity to exhibit plasticity. An important feature is that the resistance value for any SiO<sub>x</sub> device can be incrementally adjusted by controlling the initial resistance state, voltage starting position and cut-off voltage of the pulse train.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p><bold>(A)</bold> Schematic illustration of a SiO<sub>x</sub> RRAM device (not to scale). TiN layers (electrodes) are 100 nm thick (yellow). The active SiO<sub>x</sub> layer is &#x0007E;37 nm thick (blue). The gray needles represent tungsten probes. Controlling device resistance by SET and RESET with leading and trailing edges of voltage pulses. <bold>(B)</bold> Descending train of 300 ns voltage pulses (6.4&#x02013;0 V: blue line for clarity) increases device resistance from &#x0007E;100 &#x003A9; to &#x0007E;20 k&#x003A9; (orange data points). <bold>(C)</bold> Ascending train of 300 ns voltage pulses (0&#x02013;5 V: blue line for clarity) decreases device resistance from &#x0007E;4 k to &#x0007E;200 &#x003A9; (orange data points).</p></caption>
<graphic xlink:href="fnins-12-00057-g0002.tif"/>
</fig>
</sec>
<sec>
<title>Spike-timing dependent plasticity (STDP)</title>
<p>That the sign of voltage pulse gradient control both setting and resetting processes in our devices suggests that combining suitable tailored pulses and applying these to individual devices may allow us to emulate STDP in biological synapses. In this study, external circuitry is used to provide programming voltage pulses, but it should be borne in mind that they may in principle be generated by optimized unipolar devices, as we have previously demonstrated (Mehonic and Kenyon, <xref ref-type="bibr" rid="B21">2016</xref>).</p>
<p>STDP is a form of Hebbian learning rule that results in synaptic modification based on the relative timing of pre- and post-synaptic voltage spikes (Hebb, <xref ref-type="bibr" rid="B7">1949</xref>; Bi and Poo, <xref ref-type="bibr" rid="B4">1998</xref>; Abbott et al., <xref ref-type="bibr" rid="B1">2000</xref>; Jo et al., <xref ref-type="bibr" rid="B13">2010</xref>). That is to say, the relative arrival time of two spikes (&#x0201C;pre-synaptic&#x0201D; and &#x0201C;post-synaptic&#x0201D;) at the synaptic gap determines whether the synaptic strength will be increased, decreased, or remain unaltered. The change of state in a biological synapse is represented the change in excitatory postsynaptic current amplitude; in a RRAM device this change is represented by the change in device conductance (Chua, <xref ref-type="bibr" rid="B6">1971</xref>).</p>
<p>In order to control both setting and resetting processes we need to design suitable voltage pulses and to vary the delay between &#x0201C;pre-synaptic&#x0201D; and &#x0201C;post-synaptic&#x0201D; pulses. We chose two different combinations of pulse shape to test the STDP properties of our devices, with programming pulses being the sum of distinct pulses with variable inter-pulse delays. Figure <xref ref-type="fig" rid="F3">3</xref> illustrates the concept of changing device conductance by tailoring pulse shape. In the first case, we used square and triangle pulses as the pre- and post-synaptic spikes, respectively (both 1.5 V, Figure <xref ref-type="fig" rid="F3">3A</xref>). To ensure that the transition point between setting and resetting is at zero the triangle pulse had different rise and fall times (Figure <xref ref-type="fig" rid="F3">3A</xref>). The leading edge of the triangle pulse was set to be equal to the duration of the square pulse. When the applied pulse, or combination of pulses, is below a threshold voltage, V<sub>th</sub>, there is no change to the device conductance. However, for pulses above V<sub>th</sub>, the device conductance will change dependent on the sequence of voltages applied to the device. In the case of a slowly rising pulse, the device will SET into a low resistance state once the threshold is reached. Rapidly reducing the voltage has little effect, and the device remains in the LRS. However, if the pulse has a sharp leading edge and a slow trailing edge, the integral of current flowing through the conductive filament can become high enough to generate sufficient Joule heating to rupture the filament. The device therefore enters the HRS. In Figure <xref ref-type="fig" rid="F4">4</xref> we present an example of 12-cycle set of combined pulses for a non-identical pulse configuration. Each subsequent pulse has the delay between pre- and post-synaptic pulses reduced by 90 &#x003BC;s, and hence the resultant pulse shape changes from square followed by triangle to triangle followed by square. At intermediate delays we can see the device under test initially setting to a LRS and then resetting to a HRS. Device resistance is thus dependent on the time difference between the two pulses. As was predicted, a rising edge above V<sub>th</sub> lowers the resistance of the device under test, while a gradually descending voltage edge above V<sub>th</sub> increases it.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Illustration of a non-identical STDP-mimicking pulse set up. If a square pulse and a triangular pulse are below threshold there is no change in device resistance if a single pulse is applied; however, if the sum of these two pulses is above the threshold it is possible to adjust device resistance. In these examples the square pulse is a pre-synaptic spike and triangle pulse is a post-synaptic spike. <bold>(A)</bold> Pre-synaptic spike arrives earlier than the post-synaptic spike; the resulting sum is a slow leading edge above the threshold. This leads to decrease in resistance (increase in conductance&#x02014;SET process). <bold>(B)</bold> The post-synaptic spike arrives earlier than the pre-synaptic spike; the resulting sum is a slow trailing edge above the threshold. This leads to an increase in resistance (decrease in conductance&#x02014;RESET process).</p></caption>
<graphic xlink:href="fnins-12-00057-g0003.tif"/>
</fig>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p><bold>(A)</bold> Experimental device resistance measurements for a sequence of 12 programming pulses with varying pre-post (square-triangle) delays. Square 0.6 V 300 ns reading pulses (indistinguishable on this time scale, so present as a blue band) between shaped programming pulses yield 450 resistance measurements. We see clear evidence of potentiation and depression depending on programming pulse shape. <bold>(B)</bold> Experimental resistance measurement for a sequence of 11 pulses with varying pre-post delays in which the post-synaptic spike is altered with an emulated capacitor. We see again evidence of potentiation and depression.</p></caption>
<graphic xlink:href="fnins-12-00057-g0004.tif"/>
</fig>
<p>In order to generate an STDP-like response with nominally identical pre- and post-synaptic pulses we used two square voltage pulses and modified the post-synaptic pulse with a capacitor; this is a scheme that may easily be implemented in a real system. Our previous work has demonstrated how voltage spikes can be generated by individual RRAM elements, and it would be straightforward to combine these with appropriate capacitive elements. The result is a modified square pulse that looks more like a triangular pulse. This allows us to exploit the previously-shown device response to different leading and trailing pulse edge slopes. To achieve the required gradient we emulated the response of a 3 &#x003BC;F capacitor, but it should be noted that much smaller capacitance values will be sufficient for devices operating in the nanosecond regime rather than the microsecond regime studied here for demonstration. The square pulse peak voltage remained the same at 1.5 V, and the pulse width was 450 &#x003BC;s. The shift that was previously required for a non-identical configuration was no longer needed due to the modification of the post-synaptic pulse. Figure <xref ref-type="fig" rid="F4">4B</xref> shows an example of a sequence of 11 pulses. The DUT shows similar behavior to the previous setup, and we are able to demonstrate that this pulse setup mimics biological STDP (Figure <xref ref-type="fig" rid="F5">5</xref>).</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p><bold>(A)</bold> Plotted percentage occurrence of expected operation (decrease or increase in conductance for negative and positive time delays, respectively) for 306 independent pulse sets. The occurrence rate is above 50% for most time difference configurations, indicating that expected operations are statistically favorable. <bold>(B)</bold> Median percent change in synaptic weight (conductance) vs. time delay between pre- and post-synaptic spikes for 306 independent identical pulse sets with 11 cycles in each set. <bold>(C)</bold> Mean percent change in synaptic weight (conductance) vs. time delay between pre- and post-synaptic spikes. Data were limited to events that were classified as successes in the statistical analysis shown in <bold>(A)</bold>. Multiple memristive devices were used in this experiment. There is a distinct shape resemblance to biological STDP. Dotted lines are guides for the eye.</p></caption>
<graphic xlink:href="fnins-12-00057-g0005.tif"/>
</fig></sec>
<sec>
<title>Statistical analysis</title>
<p>We conducted a statistical analysis of the second pulse configuration (square pulse and modified square pulse). The results in Figure <xref ref-type="fig" rid="F5">5</xref> are based on 306 pulse cycles identical to those represented in Figure <xref ref-type="fig" rid="F4">4</xref>. Before applying the pulse trains the devices were electroformed to a preconfigured resistance of &#x0007E;2 k&#x003A9;. The pre-configuration routine is necessary to force a conductive filament into a metastable state whose conductance can be either increased or decreased by the appropriate voltage pulse combination. If the filament were in a very high resistance state, for example, further increases in resistance would not be possible. To investigate the statistical significance of timing-dependent changes in conductance, we determined the success rate of expected occurrences (Figure <xref ref-type="fig" rid="F5">5A</xref>). Here we measured for the same 306 pulse cycles the probability of a negative difference in pre- to post-synaptic pulse timing resulting in a positive change in synaptic weight (conductance), and vice-versa for positive differences in spike timing. We find that nearly all instances had a success rate &#x0003E;50%, indicating the statistical favorability of STDP-like behavior. The median of synaptic weight change (i.e., change in device conductance) is plotted against the relative time delay between pre- and post-synaptic pulses shows a striking similarity to biological STDP (Figure <xref ref-type="fig" rid="F5">5B</xref>). There is a transition at small positive delays that mimics biological behavior. Figure <xref ref-type="fig" rid="F5">5C</xref> shows the mean percent change in synaptic weight (conductance) versus time delay between pre- and post-synaptic spikes. In this case, we selected those instances of successful change in synaptic weight and constructed an STDP response curve. We have thus a proof-of-principle demonstration that not only it is possible to achieve STDP behavior, but also it is a statistically favorable process. We note at this point that the unsuccessful STDP events in our non-optimized devices are primarily related to the variability in set and reset voltages. This stochasticity may be greatly reduced through material optimization and device design, as we have recently demonstrated in RRAM devices for non-volatile memory applications. We are therefore confident that further work will improve substantially the success rates demonstrated here.</p></sec></sec>
<sec sec-type="conclusions" id="s4">
<title>Conclusion</title>
<p>We have shown that it is possible to achieve an STDP response from unipolar SiO<sub>x</sub> RRAM devices. Furthermore, we have demonstrated that STDP behavior is a statistically favorable occurrence, even in non-optimized devices. This extends and adds to our previous work showing spiking, integration and thresholding in unipolar SiO<sub>x</sub> RRAM devices (Yu et al., <xref ref-type="bibr" rid="B35">2011</xref>), and allows us to speculate that such devices can be used in hardware spiking neural networks to mimic a broad range of rich neuronal behavior, greatly simplifying the design of hardware neuromorphic systems. Further work is certainly required to establish design rules around neuromorphic systems based on RRAM devices, and the work we report is very much proof-of-concept at this stage. In particular, work is required to determine the flexibility of RRAM devices in terms of control they offer over metrics such as spiking frequency, power dissipation, conductance, and so on. We know from previous work that the electrical response of RRAM devices can be tailored through appropriate programming&#x02014;particularly by the initial electroforming step&#x02014;but work is needed to explore the limits of this approach. Further work is also needed to optimize both material and device structure to increase reliability. Nevertheless, combined with excellent CMOS compatibility, silicon oxide unipolar memristive devices are an interesting prospect for future neuromorphic systems.</p></sec>
<sec id="s5">
<title>Author contributions</title>
<p>KZ: Performed the STDP experiments, designed the study, interpreted the data, and wrote the first draft of the manuscript; AM: Assisted in designing the study, interpreted the data, and assisted with drafting the manuscript; LM and MB: Assisted with experimental design and measurement; SH: Assisted with experiments; AK: Oversaw the work, co-designed the study, and assisted with interpreting the results; All authors contributed to the manuscript.</p>
<sec>
<title>Conflict of interest statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p></sec>
</sec>
</body>
<back>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Abbott</surname> <given-names>L. F.</given-names></name> <name><surname>Song</surname> <given-names>S.</given-names></name> <name><surname>Miller</surname> <given-names>K. D.</given-names></name></person-group> (<year>2000</year>). <article-title>Competitive Hebbian learning through spike-timing-dependent synaptic plasticity</article-title>. <source>Nat. Neurosci</source>. <volume>3</volume>, <fpage>919</fpage>&#x02013;<lpage>926</lpage>. <pub-id pub-id-type="doi">10.1038/81453</pub-id><pub-id pub-id-type="pmid">10966623</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ambrogio</surname> <given-names>S.</given-names></name> <name><surname>Balatti</surname> <given-names>S.</given-names></name> <name><surname>Nardi</surname> <given-names>F.</given-names></name> <name><surname>Facchinetti</surname> <given-names>S.</given-names></name> <name><surname>Ielmini</surname> <given-names>D.</given-names></name></person-group> (<year>2013</year>). <article-title>Spike-timing dependent plasticity in a transistor-selected resistive switching memory</article-title>. <source>Nanotechnology</source> <volume>24</volume>:<fpage>384012</fpage>. <pub-id pub-id-type="doi">10.1088/0957-4484/24/38/384012</pub-id><pub-id pub-id-type="pmid">23999495</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Ananthanarayanan</surname> <given-names>R.</given-names></name> <name><surname>Esser</surname> <given-names>S.</given-names></name> <name><surname>Simon</surname> <given-names>H.</given-names></name> <name><surname>Modha</surname> <given-names>D.</given-names></name></person-group> (<year>2009</year>). <article-title>The cat is out of the bag: cortical simulations with 109 neurons, 1013 synapses</article-title>&#x0201D; in <source>Proceedings of 2009 IEEE/ACM Conference High Performance Networking Computing</source> (<publisher-loc>Portland</publisher-loc>).</citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bi</surname> <given-names>G. Q.</given-names></name> <name><surname>Poo</surname> <given-names>M. M.</given-names></name></person-group> (<year>1998</year>). <article-title>Synaptic modifications in cultured hippocampal neurons: dependence on spike timing, synaptic strength, and postsynaptic cell type</article-title>. <source>J. Neurosci</source>. <volume>18</volume>, <fpage>10464</fpage>&#x02013;<lpage>10472</lpage>. <pub-id pub-id-type="pmid">9852584</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>H. Y.</given-names></name> <name><surname>Yu</surname> <given-names>S.</given-names></name> <name><surname>Gao</surname> <given-names>B.</given-names></name> <name><surname>Huang</surname> <given-names>P.</given-names></name> <name><surname>Kang</surname> <given-names>J.</given-names></name> <name><surname>Wong</surname> <given-names>H.-S. P.</given-names></name></person-group> (<year>2012</year>). <article-title>HfOx based vertical resistive random access memory for cost-effective 3D cross-point architecture without cell selector</article-title>, in <source>2012 International Electron Devices Meeting</source>, <fpage>20.7.1</fpage>&#x02013;<lpage>20.7.4</lpage>. <publisher-loc>San Francisco</publisher-loc>.</citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chua</surname> <given-names>L.</given-names></name></person-group> (<year>1971</year>). <article-title>Memristor-The missing circuit element</article-title>. <source>IEEE Trans. Circ. Theory</source> <volume>18</volume>, <fpage>507</fpage>&#x02013;<lpage>519</lpage>. <pub-id pub-id-type="doi">10.1109/TCT.1971.1083337</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Hebb</surname> <given-names>D. O.</given-names></name></person-group> (<year>1949</year>). <source>The Organization of Behavior: A Neuropsychological Theory</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Wiley</publisher-name>.</citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Indiveri</surname> <given-names>G.</given-names></name></person-group> (<year>2000</year>). <article-title>Modeling selective attention using a neuromorphic analog VLSI device</article-title>. <source>Neural Comput</source>. <volume>12</volume>, <fpage>2857</fpage>&#x02013;<lpage>2880</lpage>. <pub-id pub-id-type="doi">10.1162/089976600300014755</pub-id><pub-id pub-id-type="pmid">11112258</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Indiveri</surname> <given-names>G.</given-names></name> <name><surname>Chicca</surname> <given-names>E.</given-names></name> <name><surname>Douglas</surname> <given-names>R.</given-names></name></person-group> (<year>2006</year>). <article-title>A VLSI array of low-power spiking neurons and bistable synapses with spike-timing dependent plasticity</article-title>. <source>IEEE Trans. Neural Netw.</source> <volume>17</volume>, <fpage>211</fpage>&#x02013;<lpage>221</lpage>. <pub-id pub-id-type="doi">10.1109/TNN.2005.860850</pub-id><pub-id pub-id-type="pmid">16526488</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Indiveri</surname> <given-names>G.</given-names></name> <name><surname>Linares-Barranco</surname> <given-names>B.</given-names></name> <name><surname>Hamilton</surname> <given-names>T. J.</given-names></name> <name><surname>van Schaik</surname> <given-names>A.</given-names></name> <name><surname>Etienne-Cummings</surname> <given-names>R.</given-names></name> <name><surname>Delbruck</surname> <given-names>T. S.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Neuromorphic silicon neuron circuits</article-title>. <source>Front. Neurosci</source>. <volume>5</volume>:<fpage>73</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2011.00073</pub-id><pub-id pub-id-type="pmid">21747754</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Indiveri</surname> <given-names>G.</given-names></name> <name><surname>Linares-Barranco</surname> <given-names>B.</given-names></name> <name><surname>Legenstein</surname> <given-names>R.</given-names></name> <name><surname>Deligeorgis</surname> <given-names>G.</given-names></name> <name><surname>Prodromakis</surname> <given-names>T.</given-names></name></person-group> (<year>2013</year>). <article-title>Integration of nanoscale memristor synapses in neuromorphic computing architectures</article-title>. <source>Nanotechnology</source> <volume>24</volume>:<fpage>384010</fpage>. <pub-id pub-id-type="doi">10.1088/0957-4484/24/38/384010</pub-id><pub-id pub-id-type="pmid">23999381</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Izhikevich</surname> <given-names>E. M.</given-names></name> <name><surname>Edelman</surname> <given-names>G. M.</given-names></name></person-group> (<year>2008</year>). <article-title>Large-scale model of mammalian thalamocortical systems</article-title>. <source>Proc. Natl. Acad. Sci. U.S.A</source>. <volume>105</volume>, <fpage>3593</fpage>&#x02013;<lpage>3598</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.0712231105</pub-id><pub-id pub-id-type="pmid">18292226</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jo</surname> <given-names>S. H.</given-names></name> <name><surname>Chang</surname> <given-names>T.</given-names></name> <name><surname>Ebong</surname> <given-names>I.</given-names></name> <name><surname>Bhadviya</surname> <given-names>B. B.</given-names></name> <name><surname>Mazumder</surname> <given-names>P.</given-names></name> <name><surname>Lu</surname> <given-names>W.</given-names></name></person-group> (<year>2010</year>). <article-title>Nanoscale memristor device as synapse in neuromorphic systems</article-title>. <source>Nano Lett</source>. <volume>10</volume>, <fpage>1297</fpage>&#x02013;<lpage>1301</lpage>. <pub-id pub-id-type="doi">10.1021/nl904092h</pub-id><pub-id pub-id-type="pmid">20192230</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Le Masson</surname> <given-names>G.</given-names></name> <name><surname>Renaud-Le Masson</surname> <given-names>S.</given-names></name> <name><surname>Debay</surname> <given-names>D.</given-names></name> <name><surname>Bal</surname> <given-names>T.</given-names></name></person-group> (<year>2002</year>). <article-title>Feedback inhibition controls spike transfer in hybrid thalamic circuits</article-title>. <source>Nature</source> <volume>417</volume>, <fpage>854</fpage>&#x02013;<lpage>858</lpage>. <pub-id pub-id-type="doi">10.1038/nature00825</pub-id><pub-id pub-id-type="pmid">12075353</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mead</surname> <given-names>C.</given-names></name></person-group> (<year>1990</year>). <article-title>Neuromorphic electronic systems</article-title>. <source>Proc. IEEE</source> <volume>78</volume>, <fpage>1629</fpage>&#x02013;<lpage>1636</lpage>. <pub-id pub-id-type="doi">10.1109/5.58356</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mehonic</surname> <given-names>A.</given-names></name> <name><surname>Buckwell</surname> <given-names>M.</given-names></name> <name><surname>Montesi</surname> <given-names>L.</given-names></name> <name><surname>Garnett</surname> <given-names>L.</given-names></name> <name><surname>Hudziak</surname> <given-names>S.</given-names></name> <name><surname>Fearn</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Structural changes and conductance thresholds in metal-free intrinsic SiOx resistive random access memory</article-title>. <source>J. Appl. Phys</source>. <volume>117</volume>:<fpage>124505</fpage>. <pub-id pub-id-type="doi">10.1063/1.4916259</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mehonic</surname> <given-names>A.</given-names></name> <name><surname>Buckwell</surname> <given-names>M.</given-names></name> <name><surname>Montesi</surname> <given-names>L.</given-names></name> <name><surname>Munde</surname> <given-names>M. S.</given-names></name> <name><surname>Gao</surname> <given-names>D.</given-names></name> <name><surname>Hudziak</surname> <given-names>S. J.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Nanoscale transformations in metastable, amorphous, silicon-rich silica</article-title>. <source>Adv. Mater</source>. <volume>28</volume>, <fpage>7486</fpage>&#x02013;<lpage>7493</lpage>. <pub-id pub-id-type="doi">10.1002/adma.201601208</pub-id><pub-id pub-id-type="pmid">27334656</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mehonic</surname> <given-names>A.</given-names></name> <name><surname>Cueff</surname> <given-names>S.</given-names></name> <name><surname>Wojdak</surname> <given-names>M.</given-names></name> <name><surname>Hudziak</surname> <given-names>S.</given-names></name> <name><surname>Jambois</surname> <given-names>C.</given-names></name> <name><surname>Labb&#x000E9;</surname> <given-names>O.</given-names></name> <etal/></person-group>. (<year>2012a</year>). <article-title>Resistive switching in silicon suboxide films</article-title>. <source>J. Appl. Phys</source>. <volume>111</volume>:<fpage>074507</fpage>. <pub-id pub-id-type="doi">10.1063/1.3701581</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mehonic</surname> <given-names>A.</given-names></name> <name><surname>Cueff</surname> <given-names>S.</given-names></name> <name><surname>Wojdak</surname> <given-names>M.</given-names></name> <name><surname>Hudziak</surname> <given-names>S.</given-names></name> <name><surname>Labbe</surname> <given-names>C.</given-names></name> <name><surname>Rizk</surname> <given-names>R.</given-names></name> <etal/></person-group>. (<year>2012b</year>). <article-title>Electrically tailored resistance switching in silicon oxide</article-title>. <source>Nanotechnology</source> <volume>23</volume>:<fpage>455201</fpage>. <pub-id pub-id-type="doi">10.1088/0957-4484/23/45/455201</pub-id><pub-id pub-id-type="pmid">23064085</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Mehonic</surname> <given-names>A.</given-names></name> <name><surname>Kenyon</surname> <given-names>A. J.</given-names></name></person-group> (<year>2015</year>). <article-title>Resistive switching in oxides</article-title>, in <source>Defects at Oxide Surfaces</source>, eds <person-group person-group-type="editor"><name><surname>Jupille</surname> <given-names>J.</given-names></name> <name><surname>Thornton</surname> <given-names>G.</given-names></name></person-group> (<publisher-loc>Cham</publisher-loc>: <publisher-name>Springer International Publishing</publisher-name>), <fpage>401</fpage>&#x02013;<lpage>428</lpage>.</citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mehonic</surname> <given-names>A.</given-names></name> <name><surname>Kenyon</surname> <given-names>A. J.</given-names></name></person-group> (<year>2016</year>). <article-title>Emulating the electrical activity of the neuron using a silicon oxide RRAM cell</article-title>. <source>Front. Neurosci</source>. <volume>10</volume>:<fpage>57</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2016.00057</pub-id><pub-id pub-id-type="pmid">26941598</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mehonic</surname> <given-names>A.</given-names></name> <name><surname>Munde</surname> <given-names>M. S.</given-names></name> <name><surname>Ng</surname> <given-names>W. H.</given-names></name> <name><surname>Buckwell</surname> <given-names>M.</given-names></name> <name><surname>Montesi</surname> <given-names>L.</given-names></name> <name><surname>Bosman</surname> <given-names>M. A. L. A. J.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Intrinsic resistance switching in amorphous silicon oxide for high performance SiOx ReRAM devices</article-title>. <source>Microelectr. Eng</source>. <volume>178</volume>, <fpage>98</fpage>&#x02013;<lpage>103</lpage>. <pub-id pub-id-type="doi">10.1016/j.mee.2017.04.033</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mitra</surname> <given-names>S.</given-names></name> <name><surname>Fusi</surname> <given-names>S.</given-names></name> <name><surname>Indiveri</surname> <given-names>G.</given-names></name></person-group> (<year>2009</year>). <article-title>Real-time classification of complex patterns using spike-based learning in neuromorphic VLSI</article-title>. <source>IEEE Trans. Biomed. Circuits Syst</source>. <volume>3</volume>, <fpage>32</fpage>&#x02013;<lpage>42</lpage>. <pub-id pub-id-type="doi">10.1109/TBCAS.2008.2005781</pub-id><pub-id pub-id-type="pmid">23853161</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Montesi</surname> <given-names>L.</given-names></name> <name><surname>Buckwell</surname> <given-names>M.</given-names></name> <name><surname>Zarudnyi</surname> <given-names>K.</given-names></name> <name><surname>Garnett</surname> <given-names>L.</given-names></name> <name><surname>Hudziak</surname> <given-names>S.</given-names></name> <name><surname>Mehonic</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Nanosecond analog programming of substoichiometric silicon oxide resistive RAM</article-title>. <source>IEEE Trans. Nanotechnol</source>. <volume>15</volume>, <fpage>428</fpage>&#x02013;<lpage>434</lpage>. <pub-id pub-id-type="doi">10.1109/TNANO.2016.2539925</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sa&#x000EF;ghi</surname> <given-names>S.</given-names></name> <name><surname>Mayr</surname> <given-names>C. G.</given-names></name> <name><surname>Serrano-Gotarredona</surname> <given-names>T.</given-names></name> <name><surname>Schmidt</surname> <given-names>H.</given-names></name> <name><surname>Lecerf</surname> <given-names>G.</given-names></name> <name><surname>Tomas</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Plasticity in memristive devices for spiking neural networks</article-title>. <source>Front. Neurosci</source>. <volume>9</volume>:<fpage>51</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2015.00051</pub-id><pub-id pub-id-type="pmid">25784849</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Schacter</surname> <given-names>W. D.</given-names></name> <name><surname>Gilbert</surname> <given-names>D. T.</given-names></name> <name><surname>Wegner</surname> <given-names>D. M.</given-names></name></person-group> (<year>2011</year>). <source>Psychology, 2nd Edn</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Worth</publisher-name>.</citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Serb</surname> <given-names>A.</given-names></name> <name><surname>Bill</surname> <given-names>J.</given-names></name> <name><surname>Khiat</surname> <given-names>A.</given-names></name> <name><surname>Berdan</surname> <given-names>R.</given-names></name> <name><surname>Legenstein</surname> <given-names>R.</given-names></name> <name><surname>Prodromakis</surname> <given-names>T.</given-names></name></person-group> (<year>2016</year>).Unsupervised learning in probabilistic neural networks with multi-state metal-oxide memristive synapses. <source>Nat. Commun.</source> <volume>7</volume>:<fpage>12611</fpage>. <pub-id pub-id-type="doi">10.1038/ncomms12611</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>L. S.</given-names></name></person-group> (<year>2006</year>). <source>Handbook of Nature-Inspired and Innovative Computing: Integrating Classical Models with Emerging Technologies</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Springer</publisher-name>.</citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Strukov</surname> <given-names>D. B.</given-names></name> <name><surname>Snider</surname> <given-names>G. S.</given-names></name> <name><surname>Stewart</surname> <given-names>D. R.</given-names></name> <name><surname>Williams</surname> <given-names>R. S.</given-names></name></person-group> (<year>2008</year>). <article-title>The missing memristor found</article-title>. <source>Nature</source> <volume>453</volume>, <fpage>80</fpage>&#x02013;<lpage>83</lpage>. <pub-id pub-id-type="doi">10.1038/nature06932</pub-id><pub-id pub-id-type="pmid">18451858</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Torrezan</surname> <given-names>A. C.</given-names></name> <name><surname>Strachan</surname> <given-names>J. P.</given-names></name> <name><surname>Medeiros-Ribeiro</surname> <given-names>G.</given-names></name> <name><surname>Williams</surname> <given-names>R. S.</given-names></name></person-group> (<year>2011</year>). <article-title>Sub-nanosecond switching of a tantalum oxide memristor</article-title>. <source>Nanotechnology</source> <volume>22</volume>:<fpage>485203</fpage>. <pub-id pub-id-type="doi">10.1088/0957-4484/22/48/485203</pub-id><pub-id pub-id-type="pmid">22071289</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vogelstein</surname> <given-names>R. J.</given-names></name> <name><surname>Tenore</surname> <given-names>F. V. G.</given-names></name> <name><surname>Guevremont</surname> <given-names>L.</given-names></name> <name><surname>Etienne-Cummings</surname> <given-names>R.</given-names></name> <name><surname>Mushahwar</surname> <given-names>V. K.</given-names></name></person-group> (<year>2008</year>). <article-title>A silicon central pattern generator controls locomotion <italic>in vivo</italic></article-title>. <source>IEEE Trans. Biomed. Circuits Syst</source>. <volume>2</volume>, <fpage>212</fpage>&#x02013;<lpage>222</lpage>. <pub-id pub-id-type="doi">10.1109/TBCAS.2008.2001867</pub-id><pub-id pub-id-type="pmid">23852970</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Waser</surname> <given-names>R.</given-names></name></person-group> (<year>2012</year>). <source>Nanoelectronics and Information Technology</source>. <publisher-loc>Weinheim</publisher-loc>: <publisher-name>Wiley</publisher-name>.</citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Waser</surname> <given-names>R.</given-names></name> <name><surname>Aono</surname> <given-names>M.</given-names></name></person-group> (<year>2007</year>). <article-title>Nanoionics-based resistive switching memories</article-title>. <source>Nat. Mater</source>. <volume>6</volume>, <fpage>833</fpage>&#x02013;<lpage>840</lpage>. <pub-id pub-id-type="doi">10.1038/nmat2023</pub-id><pub-id pub-id-type="pmid">17972938</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>S.</given-names></name> <name><surname>Gao</surname> <given-names>B.</given-names></name> <name><surname>Fang</surname> <given-names>Z.</given-names></name> <name><surname>Yu</surname> <given-names>H.</given-names></name> <name><surname>Kang</surname> <given-names>J.</given-names></name> <name><surname>Wong</surname> <given-names>H.-S. P.</given-names></name></person-group> (<year>2013</year>). <article-title>A low energy oxide-based electronic synaptic device for neuromorphic visual systems with tolerance to device variation</article-title>. <source>Adv. Mater</source>. <volume>25</volume>, <fpage>1774</fpage>&#x02013;<lpage>1779</lpage>. <pub-id pub-id-type="doi">10.1002/adma.201203680</pub-id><pub-id pub-id-type="pmid">23355110</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>S.</given-names></name> <name><surname>Wu</surname> <given-names>Y.</given-names></name> <name><surname>Jeyasinghe</surname> <given-names>R.</given-names></name> <name><surname>Kuzum</surname> <given-names>D.</given-names></name> <name><surname>Wong</surname> <given-names>S. P. H.</given-names></name></person-group> (<year>2011</year>). <article-title>An electronic synapse device based on metal oxide resistive switching memory for neuromorphic computation</article-title>. <source>IEEE Trans. Electron Dev.</source> <volume>58</volume>, <fpage>2729</fpage>&#x02013;<lpage>2737</lpage>. <pub-id pub-id-type="doi">10.1109/TED.2011.2147791</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> We gratefully acknowledge financial support provided by the Engineering and Physical Sciences Research Council (EPSRC). AM gratefully acknowledges the Royal Academy of Engineering and the funding provided through the Research Fellowship &#x0201C;Next Generation Adaptive Electronics for Neuromorphic Engineering&#x0201D;.</p>
</fn>
</fn-group>
</back>
</article>