Search tips
Search criteria 


Logo of frontneurosciLink to Publisher's site
Front Neurosci. 2017; 11: 91.
Published online 2017 February 28. doi:  10.3389/fnins.2017.00091
PMCID: PMC5328953

Double-Barrier Memristive Devices for Unsupervised Learning and Pattern Recognition


The use of interface-based resistive switching devices for neuromorphic computing is investigated. In a combined experimental and numerical study, the important device parameters and their impact on a neuromorphic pattern recognition system are studied. The memristive cells consist of a layer sequence Al/Al2O3/NbxOy/Au and are fabricated on a 4-inch wafer. The key functional ingredients of the devices are a 1.3 nm thick Al2O3 tunnel barrier and a 2.5 mm thick NbxOy memristive layer. Voltage pulse measurements are used to study the electrical conditions for the emulation of synaptic functionality of single cells for later use in a recognition system. The results are evaluated and modeled in the framework of the plasticity model of Ziegler et al. Based on this model, which is matched to experimental data from 84 individual devices, the network performance with regard to yield, reliability, and variability is investigated numerically. As the network model, a computing scheme for pattern recognition and unsupervised learning based on the work of Querlioz et al. (2011), Sheridan et al. (2014), Zahari et al. (2015) is employed. This is a two-layer feedforward network with a crossbar array of memristive devices, leaky integrate-and-fire output neurons including a winner-takes-all strategy, and a stochastic coding scheme for the input pattern. As input pattern, the full data set of digits from the MNIST database is used. The numerical investigation indicates that the experimentally obtained yield, reliability, and variability of the memristive cells are suitable for such a network. Furthermore, evidence is presented that their strong IV non-linearity might avoid the need for selector devices in crossbar array structures.

Keywords: memristive devices, synaptic plasticity, spiking neuron, neural network, neuromorphic systems, unsupervised learning


The brains of humans, mammals, and even simple living species like invertebrates are well-adapted to permanently changing environments. Their nervous systems exhibit remarkable interactions with their surroundings—a result of millions of years of evolution and explicable by Darwinism (Shanahan, 2004). As a result, biological systems are presently unmatched in the efficient way in which they are able to perform cognitive tasks, such as pattern recognition, with extremely low power consumption. It is therefore no surprise that attempts have been made to develop bio-inspired computing systems, so-called neuromorphic systems, with the goal of reaching the performance and power efficiency of biological systems (Liu et al., 2002; Chicca et al., 2014). Machine learning dates back to the early days of serial, binary computation based on the von Neumann architecture, but today's artificial neural systems are still only partially able to mimic biological systems. The huge power dissipation, the long computational time, and the need for large datasets are among the major problems faced in attempting to realize artificial neural networks (ANNs).

In this context, analog very large-scale integration (VLSI) based on silicon complementary metal-oxide-semiconductor (CMOS) technology (Liu et al., 2002; Chicca et al., 2014) might provide advantages to software-dominated neuroinformatics (Amit, 1989; Würtz, 2008). This field has gained new momentum with the advent of the concept of memristive devices (memristors) (Strukov et al., 2008; Jo et al., 2010; Jeong et al., 2013). The memristor is a device whose resistance depends on its history of applied potentials. It was predicted in 1971 by Chua (1971). In the memristor model, the resistance (or memristance) M(x, V, t) of the device can be expressed by an internal state variable x(t), which depends on the applied voltage V(t) and on itself:


Here, f describes the dynamics of the updating process of x(t). In many memristive devices, f describes the process of ion migration due to an externally applied voltage.

Recently, Hebbian learning as an important biological concept has been realized with single memristive devices by emulating spike-timing-dependent plasticity (STDP; Jo et al., 2010; Zamarreño-Ramos et al., 2011), long-term potentiation, and long-term depression (Ohno et al., 2011). These properties are important cellular mechanisms of memory and learning in neural networks (Ziegler et al., 2015).

Although the structure of memristive devices is rather simple (consisting of a capacitor-like metal-insulator-metal sequence in the simplest case), the underlying physical mechanism is often unclear. The function of many memristive devices is based on the presence of conductive filaments, which results in poor switching reproducibility and high inter-device variability and often requires an initial and individual electrical forming step (Szot et al., 2006; Waser et al., 2009; Ha and Ramanathan, 2011; Yang et al., 2013; Dirkmann et al., 2015). Furthermore, memristive devices have mainly been considered for application in future resistive random-access memories (RRAMs; Itoh et al., 2006), whose system architecture and functionality are similar to those of digital memories. Here, switching times in the nanosecond range, long data retention times (10 years), low device variability, and good fatigue performance are essential requirements. Memristive devices for use in neuromorphic systems are subject to different requirements, and other resistive switching concepts can be of interest (Ziegler et al., 2015), such as interface-based memristive devices.

In interface-based devices, uniform interface effects lead to homogeneous rather than abrupt changes in resistance. Furthermore, such devices do not suffer from the randomness generated by electroforming or filament growth (Kohlstedt et al., 1993; Baikalov et al., 2003; Meyer et al., 2008; Park et al., 2008; Sawa, 2008; Baik and Lim, 2010; Hu et al., 2011; Jeong et al., 2011; Aoki et al., 2014; Mikheev et al., 2014). In some designs for these devices, the resistive switching results from changes at a Schottky-like contact (Baikalov et al., 2003; Mikheev et al., 2014). In another approach, the electron tunneling probability is varied when a memristive layer is in contact with a tunnel barrier (Kohlstedt et al., 1993; Meyer et al., 2008; Baik and Lim, 2010; Jeong et al., 2011). Recently, we have been able to combine both concepts into a single device by sandwiching an ultra-thin memristive layer between a tunnel barrier and a Schottky-like contact (Hansen et al., 2015). This offers several benefits: the tunnel barrier defines the lower resistance boundary and limits the current through the device. Further, both the tunnel and the Schottky barrier define chemical barriers to ion migration, leading to improved data retention compared with single-barrier concepts. However, the possible benefits for bio-inspired neuromorphic circuits are not that obvious and hence have to be explored.

The purpose of this work is to provide a thorough analysis of interface-based memristive devices and their potential use in bio-inspired neuromorphic systems. With this aim in mind, double-barrier memristive devices with a layer sequence Al/Al2O3/NbxOy/Au have been fabricated. The thickness of the Al2O3 tunnel barrier is ~1.3 nm and that of the memristive NbxOy layer is 2.5 nm. The use of an ultra-thin memristive layer allows the required electrical field strength for oxygen ion migration to be achieved, reducing the resistance variation due to interfacial processes, and also allows interference between electron tunneling and the Schottky barrier (Hansen et al., 2015). To determine the important device characteristics and necessary electrical conditions for use in neuromorphic systems, automated measurements are undertaken, providing valid statistical data from 84 single devices. The resistive switching properties are investigated using voltage pulses with different lengths and amplitudes. The yield, reliability, and device variability are determined. To obtain a realistic mathematical description of the device behavior from the obtained experimental data for network-level applications, the memristive plasticity model of Ziegler et al. (2015) is used. This model allows us to explore functional consequences of the use of individual memristive devices for the emulation of Hebbian plasticity. As guidelines for network simulation, previously reported neural-network-based paradigms are adopted (Querlioz et al., 2011, 2013; Sheridan et al., 2014; Zahari et al., 2015). In particular, a neural network for pattern recognition is simulated consisting of inhibitorily linked output neurons within a winner-takes-it-all architecture and a homeostasis-like rule for the spiking-neurons thresholds. The network is trained on the MNIST database. On the basis of this pattern recognition network, essential requirements for the development of neuromorphic circuits with interface-based memristive devices are discussed.

The paper is organized as follows: Section Materials and Methods describes the experimental techniques used and the fabrication and electrical characterization of the memristive devices. It also presents the phenomenological synaptic learning model that we use to model the experimental data. The section closes with a brief description of the network architecture used for the pattern recognition simulation. Section Results contains the results of our investigation, namely, plasticity measurements, and network performance. The behavior of the fabricated devices is discussed with respect to their reliability and variability in Section Discussion.

Materials and methods

Electrical measurements

All measurements were performed using an Agilent E5260 source measurement unit. Current–voltage measurements (IV curves) were obtained by sweeping the applied voltage and measuring the current simultaneously. For synaptic plasticity measurements on single devices, rectangular voltage pulses with different amplitudes, polarities, and pulse durations were applied to the devices. Positive voltage is defined as positive voltage on the top electrode and ground on the bottom electrode.

Device fabrication and characterization

The memristive devices were fabricated on 4-inch Si wafers with 400 nm SiO2 (thermally oxidized) using a standard optical lithography process (a schematic of the material stack is shown in Figure Figure1A).1A). The junctions were arranged in 1 × 1 mm cells across the wafer, each containing six different contact sizes ranging from 70 to 2300 μm2. The devices were fabricated using the following procedure: First, the multilayer (including top and bottom electrodes) was deposited using DC magnetron sputtering, without breaking the vacuum. The Al2O3 tunnel barrier was fabricated by depositing Al, which was afterwards partially oxidized in situ. The NbxOy layer was deposited by reactive sputtering in an O2/Ar atmosphere. Following the subsequent lift-off, the junction area was defined by wet chemical etching of the Au top electrode using a potassium iodide solution. The etched parts were then covered with thermally evaporated SiO to insulate the bottom electrode from the subsequently deposited Nb wiring to contact the top electrode.

Figure 1
Double-barrier memristive device: (A) Schematic cross-section of the Al/Al2O3/NbxOy/Au double-barrier memristive device. (B) Current density J as function of the applied bias voltage. In the upper graph in a linear scale for J was used, while in the lower ...

Figure Figure1B1B shows a typical recorded I–V curve of a double-barrier memristive device, exhibiting the typical pinched bipolar hysteresis of a memristive device (to reveal the I–V non-linearity more clearly, in the lower panel, the curve is plotted with a logarithmic scale on the vertical axis). The applied voltage was varied between −2 and 2.8 V, while the current was measured simultaneously. The voltage was ramped from 0 to 2.8 V to switch the device from its initial high-resistance state (HRS) to the low-resistance state (LRS). To reset the device, the voltage was decreased to −2 V. The asymmetry between positive and negative current can be attributed to the Schottky-like NbxOy/Au contact, while the gradual resistance change indicates a non-filamentary resistance switching mechanism. This results from homogeneously changed interface properties (Hansen et al., 2015; Dirkmann et al., 2016). Neither an initial forming procedure nor current compliance was used. This is important for the integration of such devices into crossbar architectures, as we will discuss below.

The underlying physical mechanism of the resistive switching process has recently been studied in a combined experimental and theoretical investigation (Dirkmann et al., 2016). This investigation identified the transport of oxygen ions within the NbxOy as the key mechanism for the resistive switching process. During the set process, oxygen ions (within the NbxOy) move under positive voltage toward the Au interface and affect essential interfacial parameters (e.g., the density of states, the local barrier height, and the barrier thickness) at the Al2O3/NbxOy and NbxOy/Au (Schottky) interfaces simultaneously. By applying a negative bias voltage, the original ion distribution is restored. As a consequence, the electron transport is altered in accordance with the local ion distribution. This leads to the observed memristive IV curve. For further information about the device performance and technology, the reader is referred to Hansen et al. (2015).

Simulation model

To investigate device performance at the network level, we used a computing scheme for pattern recognition similar to Querlioz et al. (2011, 2013), Sheridan et al. (2014), Zahari et al. (2015). As a numerical model of the resistive switching under voltage pulsing, we used the plasticity model of Ziegler et al. (2015). This model is compatible with advanced biophysical plasticity models that account for experimental data on STDP. Furthermore, the model is suitable for describing plasticity emulation using memristive devices. We should mention here that the plasticity model provides a behavioral description of the memristive device, rather than an explanation of the underlying physical mechanisms. In the following, the plasticity model and network structure for pattern recognition are presented.

Plasticity model

The resistive switching of double-barrier memristive devices in neural circuits are described in the framework of the phenomenological learning model of Ziegler et al. (2015). This model links the change in conductance in a memristive device and the applied voltage pulse. The limited weight growth and the weight-dependent (memristive) learning rate make this model appropriate for memristive device and synapse emulations. In the plasticity model, the updating process for the synaptic weights ω is given by


where β is the weight-dependent learning rate and ωmax is the maximum synaptic weight. The resistive switching dynamic of the memristive device is manifested in β, which depends on the current device conductance and electrical stimuli. In particular, β depends on the switching mechanism of the memristive device and can lead to very different learning behaviors (Querlioz et al., 2015; Ziegler et al., 2015). In determining β for the double-barrier memristive device, two particularly important points should be noted: First, it is necessary to use distinct learning rates for potentiation (βp) and depression (βd) owing to the asymmetry between positive and negative voltages in the memristive device (cf. Figure Figure1B).1B). Second, the synaptic weight ω is associated with the conductance G of a memristive device, which depends on height ΔV, length Δt, and number n of applied electrical pulses. Hence, we assume that βp and βd depend on ΔV, Δt, and n, so that


Here kp, kd, and γ are positive constants, while α and λ depend on the height ΔV and length Δt of the applied voltage stimulus. The functions α and λ account for the nonlinearity of the switching dynamics of the memristive device and have to be determined for the particular memristive device under consideration (see Section Results).

Network architecture for pattern recognition

The network structure is a two-layer feedforward network and is shown schematically in Figure Figure2.2. The memristive devices are arranged in a crossbar array to which input (blue circles) and output (red circles) neurons are connected. The individual pixels of the input pattern are coded by voltage pulses within the input layer. Leaky integrate-and-fire neurons are used in the output layer. These neurons are laterally coupled within an inhibitory winner-takes-it-all network including adaptive thresholds for the spiking, as proposed in Querlioz et al. (2011). The network can be briefly described as follows: The output neuron's spiking depends on the applied input pattern and the particular resistance of the memristive devices. The conductances of the memristive devices are changed by associative learning, and the voltage across each device depends on the input pattern and the activity of the particular output neuron. This enables unsupervised learning, because every output neuron creates its own specific receptive field during learning. Afterwards, in the recognition phase, each of the output neurons will spike in accordance with the previously learned pattern for a varying input.

Figure 2
Schematic of the simulated neural network: Positive and negative voltage pulses are applied to the input of the network, representing the intensity of the individual values of the 28 × 28 pixel MNIST image. Stochastic coding of the input data ...

The input patterns were taken from the MNIST database, which consists of 60,000 handwritten digits from 250 different writers. Each digit is stored in a 256-level grayscale image with 28 × 28 pixels (LeCun et al., 1998). In the network implemented here, the images are rearranged into a 784-row input vector applied to the input neurons (one neuron per image pixel). The value of each pixel is applied as a voltage pulse, which, in combination with the output neuron's pulse, causes a change in the conductance of a particular memristive cell in the cross-bar array. The device conductance is most strongly affected when the input and output pulses match, as sketched in Figure Figure2.2. In particular, the overlap of a negative input pulse with an output neuron's pulse leads to a decrease in the particular device conductance, while the combination of a positive input pulse with an output pulse increases the conductance of the memristive devices connected to the specific neurons. Memristive devices in input rows and output columns where there is no overlap of voltage pulses remain unaffected.

It is important for the operation of the network that the input data is stochastically coded (Querlioz et al., 2011, 2013; Sheridan et al., 2014; Zahari et al., 2015), and this has been achieved by the following steps: The image pixels are normalized to the interval −1 to 1. Following Sheridan et al. (2014), the grayscale value pi of each pixel i is best normalized by


where pmean is the mean grayscale value of all pixels for an individual image. The absolute value of pinorm denotes the probability of an input spike generation. Therefore, a random number r for each pixel is generated at each iteration step. If the condition r |pinorm| is satisfied, an input pulse is generated. The sign of pinorm corresponds to the polarity of the input voltage (see Figure Figure22).

Leaky integrate-and-fire neurons are used as output neurons and are arranged in a laterally coupled inhibitory network (see Figure Figure2).2). To guarantee that individual input patterns are learned by different output neurons, the winner-takes-it-all approach is used, in which the first spiking neuron resets the integration of all other neurons. This allows unsupervised learning with the network structure that is used. Crucial for unsupervised learning is an adjustable neuron firing threshold, which guarantees that all output neurons participate equivalently in the learning phase. This can be motivated by considering the process of homeostasis in biological systems (Querlioz et al., 2011). Therefore, the threshold of a neuron to fire is increased whenever the spike number (activity) of a neuron is above the desired activity, and vice versa. Following Querlioz et al. (2011), this can be achieved by using


for the threshold voltage adaptation. Here, γth, Aavg, and Atar are respectively the gain factor, the mean activity of an individual neuron, and the target activity (γth = 7.5 × 10−3 and Atar = 10).


Synaptic plasticity measurements

The basis for the emulation of synaptic functionalities with memristive devices is the fact that the synaptic weight ω between individual neurons can be related to the state variable x of a memristive device according to Equation (1). Zamarreño-Ramos et al. showed that the conductance of a memristive device is proportional to ω in the ideal voltage-driven memristor model (see Equation 1) (Zamarreño-Ramos et al., 2011). Hence, biological plasticity mechanisms can be emulated by changes in the device conductance under suitable voltage pulses. To study the capability of synaptic plasticity emulations using the Al/Al2O3/NbxOy/Au double-barrier memristive device of Figure Figure1,1, voltage pulses with different pulse widths Δt and amplitudes ΔV were applied to 86 individual devices.

Figure Figure3A3A shows a typical plasticity measurement (conductance vs. pulse number). We found that the double-barrier memristive device exhibits a gradually changing conductance under voltage pulsing, as is desired for plasticity emulation. To study the pulse dependence, 1000 equivalent positive voltage pulses (potentiation pulses) followed by 1000 equivalent negative voltage pulses (depression pulses) were used, as shown schematically in the inset of Figure Figure3A.3A. For the emulation of potentiation, voltage pulses of ΔV = 3.9 V and Δt = 1 ms in length were chosen, while for depression the voltage pulse height was reduced to ΔV = −2.5 V. To measure the device conductance, voltage sweeps below the threshold voltage of the device up to 0.48 V were applied and the current was measured after every 100 potentiation/depression pulses. For better illustration, the obtained conductances G(n) were normalized by the average maximum conductance Gmax (100 nS).

Figure 3
Synaptic plasticity emulation: (A) Typical plasticity measurement for 1000 potentiation pulses of +3.9 V and 1000 depression pulses of −2.5 V with pulse lengths of 1 ms. Each point represents the conductance after 100 pulses. In (B–D) ...

In order to study in greater detail the variations of the device conductance as the amplitude and length of the applied voltage pulses are changed, potentiation pulses with amplitudes ranging from 2.4 to 3.7 V and pulse lengths ranging from 1 to 30 ms were applied. The results are shown in Figures 3B,C. For the variation with amplitude, the pulse length was fixed at 1 ms, while for the variation with pulse length, the amplitude was fixed at 3.9 V. The data points reflect the device conductance after 1000 voltage pulses. For an amplitude of 2.4 V, the conductance only increases by 0.1–0.5% from the initial value G0. However, if amplitude is set to 3.7 V, there is a 10-fold increase (see Figure Figure3B).3B). Of further interest is that the increase in conductance saturates for pulses above 10 ms (see Figure Figure3C).3C). This is in accordance with the diffusion times of oxygen ions within the memristive layer, which are at the heart of the resistance switching process (Dirkmann et al., 2016).

The recorded data for synaptic depression are shown in Figure Figure3D.3D. A voltage train of 1000 single voltage pulses of 30 ms length with amplitudes ranging from −1.4 to −2.6 V were applied to memristive cells that had previously been set to the high-conductance (low-resistance) state. The unchanged conductance for voltages above −1.4 V shows the threshold behavior for the reset process. In contrast, pulses with amplitudes of −2.6 V nearly reset the device conductance completely (see Figure Figure3D3D).

To adapt the experimental findings for the plasticity model described above, the functions α(ΔV) and λ(Δt) from Equation (3) have to be determined. The functions α(ΔV) and λ(Δt) account for the nonlinear dependencies of the learning rates βp and βd on the pulse amplitude and length. Good agreement with the experimental data (the black points in Figure Figure3)3) could be achieved using the following expressions (the red lines in Figure Figure33):

α(Δ V) ={α1 ΔV for ΔV>0, exp[α2(ΔVα3)]for ΔV<0,


Here α1, 2, 3 and λ1.2.3 are positive constants, which are listed in Table Table11 together with the other parameters of the fitting procedure. Owing to the differences between potentiation and depression, two different definitions for α(ΔV) are necessary: for potentiation, a linear dependence on ΔV is in good agreement with experiment, while an exponential function of ΔV is necessary to reproduce the experimentally recorded data for the emulation of synaptic depression. A linear function is best at reproducing the pulse length dependence (see Equation 7).

Table 1
Parameters for the plasticity model.

Pattern recognition

The network performance of the double-barrier memristive device was investigated by training the pattern recognition network with the complete MNIST dataset of 60,000 handwritten digits. To improve the learning performance of the network, the complete training set was applied three times to the network. After training, the network was tested by using the MNIST test dataset, which contains 10,000 digits that are not included in the training set. The performance of the network can be evaluated by counting the number of correctly identified digits. The recognition rate thereby obtained increases with the number of output neurons used. For 10, 20, 50, and 100 output neurons, recognition rates of 65, 70, 77, and 82%, respectively, were obtained. These rates are in good agreement with previously published investigations (Querlioz et al., 2011, 2013; Sheridan et al., 2014; Zahari et al., 2015). For visualization, each 784-value vector, above the particular output neuron, is rearranged into a 28 × 28 pixel image representing the receptive fields. A typical set of receptive fields for a simulation with 50 output neurons is shown in Figure Figure4.4. We found that the implemented network structure is able to learn all of the 10 input digits and the different details of each digit. We should mention here that the obtained recognition rates are much lower than those from other spiking network architectures (for an overview of different network architectures, the reader is referred to Diehl and Cook, 2015). However, the aim of our investigation was to study network requirements for the use of memristive devices in neural networks rather than to improve pattern recognition computing schemes.

Figure 4
Receptive fields: Obtained receptive fields after unsupervised learning in the case of 50 output neurons. The white pixels correspond to the maximum conductance (strong synaptic weight), while the black pixels represent minimum conductance values (weak ...

Appropriately chosen voltage pulses with lengths Δt and amplitudes ΔV are important for the performance of the neural network, to precisely adjust the conductance of the double barrier memristive devices. In the case of matching input (Vpre) and output (Vpost) pulses, the device conductance must be always affected. The voltage pulses used here are depicted in Figure Figure5.5. Based on the experimental data of Figure Figure3,3, suitable amplitudes and lengths of the voltage pulses were chosen for the simulations. For the input neuron, the pulse height Vpre was set to either +0.6 or −0.6 V (see the top row of Figure Figure5).5). The voltage pulses Vpost generated by the leaky integrate-and-fire output neurons contain a positive part of 2.9 V followed by a negative part of −2.3 V (see the middle row of Figure Figure5).5). The combination of Vpost with a negative input pulse (Vsum = −2.9 V) decreases the conductance, while the overlap of Vpost with a positive input pulse (Vsum = +3.5 V) increases the conductance of the device connected to the specific neurons (see the bottom row in Figure Figure5).5). An important aspect of the double-barrier memristive device is its strong IV nonlinearity (see Figure Figure1B).1B). According to Figure Figure3,3, an increase in pulse height from 3 to 3.8 V leads to an increase in device conductance by a factor of 10, while below 1 V the device conductance remains unaffected. This provides a threshold voltage for the change in conductance, which is a very important requirement for the learning behavior of this crossbar-based computing scheme. Without this threshold voltage, the voltage drop across neighboring devices during the set or reset could be sufficient to change also the conductance of these neighboring devices. This problem is commonplace with devices that have a symmetrical IV curve, requiring selector devices that access individual devices. Owing to the strong asymmetry and nonlinearity in the current of double-barrier memristive devices, there is no need for selector devices (Kügeler et al., 2009).

Figure 5
Pulse forms used for the simulations: For the coding of the grayscale input images, positive and negative voltage pulses (Vpre) with ±0.6 V in amplitude and 6 ms in width have been used. The voltage pulses of the output neurons (Vpost) were set ...


Three different scenarios have been investigated to determine the reliability, variability, and yield of the fabricated double-barrier memristive devices in the framework of the implemented network model: In the first scenario, a device-to-device variability was modeled. In this approach, the learning rates of all devices were initially changed once to represent a Gaussian distribution, but were kept constant during the simulation. The second scenario reflects the variability of individual memristive devices. In this case, the local learning rate of all memristive devices was varied in each iteration step. We also investigated a combination of the first and second scenarios, since this is the most likely case in reality. In all cases, the local learning rates βp and βd of the memristive devices (see Equation 3) were varied according to a Gaussian distribution.


where x=βp(d)x-βp(d)0 (here βx is the variable and β0 is the undisturbed learning rate defined in Equation 3) and σ denotes the standard deviation (see the inset of Figure Figure6A).6A). The standard deviation of the device-to-device variability (first case) is denoted by σdevice, while σiteration denotes the variability for each iteration step (second case). The simulation results are shown in Figure Figure6.6. The full set of 60,000 MNIST training images was applied three times to a network with 10 output neurons for the simulation. The presented data were averaged over three total simulation runs.

Figure 6
Device variation: To test the reliability of the network model under device variation, the local learning rates of the memristive devices were changed to represent a Gaussian distribution [var phi](x). Each curve shows the obtained recognition rates for ...

The recognition rates for the first scenario (increased device-to-device variability) are shown in Figure Figure6A.6A. We found that, up to 50%, the recognition rate is not affected by the device-to-device variability and decreases only slightly with a further increase in device variability. The outcome of the second scenario is presented in Figure Figure6B.6B. We found that the network is more robust under variations within individual devices (see the black line in Figure Figure6B).6B). Even if we combine both, a device-to-device and individual device variability (see the red and blue curves in Figure Figure6B),6B), a relatively robust network performance is obtained. Even if we combine both, a device-to-device and individual device variability (see the red and blue curves in Figure Figure6B),6B), a relatively robust network performance is obtained. Therefore, a device-to-device variability up to 100% might be acceptable for a suitable network performance. Hence, our numerical investigation provides evidence that the most striking problems are introduced by a constant device-to-device variability. This can be explained by the fact that a variation in each iteration step is averaged over the total number of iteration steps. However, a constant variation of the local learning rate between the memristive devices, i.e., a device-to-device variability, leads to an ever-increasing variation with every iteration.

To estimate the performance of double-barrier devices within the network model considered here, it is necessary to analyze the experimentally obtained device variations in greater detail. Considering the error bars in Figure Figure3,3, the variation seems to be too high to fulfill the obtained network requirement of device-to-device variation of <100%. However, these error bars were determined from the final conductance after 1000 potentiation/depression pulses. A better representation of the device performance could be obtained by measuring the conductance variations after every 100th potentiation pulse. The obtained results are shown in Figure Figure6C.6C. For a suitable comparison between simulated and measured data, the measured data were normalized by the highest recorded conductance Gmax total (see Figure Figure6C).6C). Furthermore, the solid lines in Figure Figure6C6C indicate the variation range of learning rates with, respectively, 40, 80, and 100% device-to-device variability. In particular, the experimentally recorded device-to-device variability is in the required variation interval. This suggests that the device variability can be assumed to be sufficient for the network learning process. Nevertheless, the device-to-device variability is one of the most challenging parameters for the realization of the pattern recognition computing scheme depicted in Figure Figure22.

Another possible negative impact on network performance can be introduced by defective devices. Two scenarios can be distinguished here: devices with a low resistance (e.g., devices that are shorted) and devices with a high resistance (e.g., those that are not properly connected during the final metallization step). In either case, during the simulation, the resistance was assumed to be constant. The most problematic situation arises from devices that are initially in the low-resistance state and are not able to perform resistive switching under voltage pulsing. Owing to the increased conductivity of the device and the subsequent greater weight during the training phase, the receptive field for a given column of devices will be distorted. To study the influence of such devices on overall network performance, the dependence of the recognition rates on the number of defective devices was calculated. The results for a network with 10 output neurons are shown in Figure Figure7A,7A, where the blue curve represents the situation where only one of the 10 receptive fields contains defective devices, while the red curve was obtained from a simulation where all 10 receptive fields contain defective devices. If only one receptive field contains defective devices, the overall recognition rate of the circuit is nearly constant up to 784 defective devices. However, if all 10 receptive fields above the output neurons contain defective devices, the recognition rate is drastically decreased. For defective devices with an unchangeable high resistance, the recognition rate shows significantly different behavior (Figure (Figure7B).7B). In comparison with defective low-resistance devices, defects in high-resistance devices do not drastically change the formation of the receptive fields. However, the observed influence of defective devices on network performance can be assumed to play a minor role owing to the high quality of the fabrication process. Figure Figure88 shows a map of a typical recorded resistance distribution across a wafer. We should mention here that only the center of the 4-inch wafer is depicted, because the 4-inch sputter targets limit the usable area on these wafers. Each square represents the resistance of one memristive device (although the actual density of devices is larger because only one in six devices was measured). In this wafer map, a total of 966 devices were measured. Besides the working devices, 14 devices (black squares) are defective in terms of an initial low resistance (i.e., are shorted), while five devices (yellow squares) have an unexpected high resistance (e.g., as a result of problems with the final metallization step). Yellow squares at borders (e.g., the lower right corner) are devices with a significantly larger resistance, which can be attributed to side effects of the sputtering process. In the interesting center region, the yield is thus about 98%, which is typical for our wafers. Comparing the number of defective devices with Figure Figure7,7, it is obvious that the recognition rates will remain nearly unaffected.

Figure 7
Impact on the recognition rate of defective devices: Recognition rate as function of the number of defective devices with low and high resistance per receptive field. The blue curve represents the situation where only one of the 10 receptive fields contains ...
Figure 8
Resistance distribution and yield: Wafermap showing the resistance of the memristive devices in the center of the wafer. Each square represents one device with the color-coded resistance. Black squares are shorted devices; individual yellow squares are ...

In conclusion, we have provided evidence that double-barrier devices are interesting candidates for use as artificial synapses in neuromorphic circuits. In particular, the gradual change in their resistance under voltage pulsing and the resulting reliability, variability, and yield of such devices might fulfill the requirements of neural networks. An experimental realization of the use of these components seems possible and may pave the way to a real-time implementation of a pattern recognition system.

Author contributions

MH prepared the samples, performed the measurements, analyzed the experimental results and co-wrote the manuscript. MZ supported the measurements and data interpretation. FZ and MZ developed the simulation model. The simulation results were discussed and interpreted between MZ, FZ, MH, and HK. HK and MZ conceived the idea, initiated, and supervised the experimental research. MH, MZ, and HK discussed the experimental results and contributed to the refinement of the manuscript.

Conflict of interest statement

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.


This work was supported by the Deutsche Forschungsgemeinschaft (DFG) through FOR2093.


  • Amit D. J. (1989). Modeling Brain Function, The World of Attractor Neural Networks. Cambridge: Cambridge University Press.
  • Aoki Y., Wiemann C., Feyer V., Kim H.-S., Schneider C. M., Ill-Yoo H., et al. . (2014). Bulk mixed ion electron conduction in amorphous gallium oxide causes memristive behaviour. Nat. Commun. 5, 3473. 10.1038/ncomms4473 [PubMed] [Cross Ref]
  • Baik S. J., Lim K. S. (2010). Bipolar resistance switching driven by tunnel barrier modulation in TiOx/AlOx bilayered structure. Appl. Phys. Lett. 97, 072109 10.1063/1.3481076 [Cross Ref]
  • Baikalov A., Wang Y. Q., Shen B., Lorenz B., Tsui S., Sun Y. Y., et al. (2003). Field-driven hysteretic and reversible resistive switch at the Ag–Pr0.7Ca0.3MnO3 interface. Appl. Phys. Lett. 83, 957 10.1063/1.1590741 [Cross Ref]
  • Chicca E., Stefanini F., Bartolozzi C., Indiveri G. (2014). Neuromorphic elecroncis circuits for building autonomous cognitive systems. Proc. IEEE 102, 1397–1388. 10.1109/JPROC.2014.2313954 [Cross Ref]
  • Chua O. L. (1971). Memristor-the missing circuit element. IEEE Trans. Circ. Theor. 18, 507–519. 10.1109/TCT.1971.1083337 [Cross Ref]
  • Diehl P. U., Cook M. (2015). Unsupervised learning of digit recognition using spike-timing-dependent plasticity. Front. Comput. Neurosci. 9:99. 10.3389/fncom.2015.00099 [PMC free article] [PubMed] [Cross Ref]
  • Dirkmann S., Hansen M., Ziegler M., Kohlstedt H., Mussenbrock T. (2016). The role of ion transport phenomena in memristive double barrier devices. Sci. Rep. 6:35686. 10.1038/srep35686 [PMC free article] [PubMed] [Cross Ref]
  • Dirkmann S., Ziegler M., Hansen M., Kohlstedt H., Trieschmann J., Mussenbrock T. (2015). Kinetic simulation of filament growth dynamics in memristive electrochemical metallization devices. J. Appl. Phys. 118:214501 10.1063/1.4936107 [Cross Ref]
  • Ha S. D., Ramanathan S. (2011). Adaptive oxide electronics: a review. J. Appl. Phys. 110, 071101 10.1063/1.3640806 [Cross Ref]
  • Hansen M., Ziegler M., Kolberg L., Soni R., Dirkmann S., Mussenbrock T., et al. . (2015). A double barrier memristive device. Sci. Rep. 5:13753. 10.1038/srep13753 [PMC free article] [PubMed] [Cross Ref]
  • Hu J., Nainani A., Yun S., Krishna S., Philip Wong H.-S. (2011). Impact of fixed charge on metal-insulator-semiconductor barrier height reduction. Appl. Phys. Lett. 99, 252104 10.1063/1.3669414 [Cross Ref]
  • Itoh K., Horiguchi M., Tanaka H. (2006). Ultra-Low Voltage Nano-Scale Memories. New York, NY: Springer Science and Business.
  • Jeong D. S., Cheng B., Kohlstedt H. (2011). Pt/Ti/Al2O3/Al tunnel junctions exhibiting electro forming-free bipolar resistive switching behaviour. Solid State Electron. 63, 1 10.1016/j.sse.2011.05.028 [Cross Ref]
  • Jeong D. S., Kim I., Ziegler M., Kohlstedt H. (2013). Towards artificial neurons and synapses: materials point of view. RSC Adv. 3, 3169–3183. 10.1039/c2ra22507g [Cross Ref]
  • Jo S. H., Chang T., Ebong I., Bhadviya B. B., Mazumder P., Lu W. (2010). Nanoscale memristor device as synapse in neuromorphic systems. Nano Lett. 10, 1297. 10.1021/nl904092h [PubMed] [Cross Ref]
  • Kohlstedt H., Kuriki S., Gundlach K.-H. (1993). Electro forming and telegraph noise in Pb/Bi/InOx/Pb/Bi tunnel junctions. J. Appl. Phys. 73, 2564 10.1063/1.353066 [Cross Ref]
  • Kügeler C., Meier M., Rosezin R., Gilles S., Waser R. (2009). High density 3D memory architecture based on the resistive switching effect. Solid State Electron. 53, 1287–1292. 10.1016/j.sse.2009.09.034 [Cross Ref]
  • LeCun Y., Bottou L., Bengio Y., Haffner P. (1998). Gradient-based learning applied to document recognition. Proc. IEEE 86, 2278–2324. 10.1109/5.726791 [Cross Ref]
  • Liu S. C., Kramer J., Indiveri G., Delbrück T., Douglas R. (2002). Analog VLSI: Circuits and Principles. Cambridge: MIT Press.
  • Meyer R., Schloss J., Brewer J., Lambertson R., Kinney W., Sanchez J., et al. (2008). Oxide dual-layer memory element for scalable non-volatile cross-point memory technology. NVMTS 9, 1–5. 10.1109/nvmt.2008.4731194 [Cross Ref]
  • Mikheev E., Hoskins B. D., Strukov D. B., Stemmer S. (2014). Resistive switching and its suppression in Pt/Nb:SrTiO3 junctions. Nat. Commun. 5, 3990. 10.1038/ncomms4990 [PMC free article] [PubMed] [Cross Ref]
  • Ohno T., Hasegawa T., Tsuruoka T., Terabe K., Gimzewski J. K., Aono M. (2011). Short-term plasticity and long-term potentiation mimicked in single inorganic synapses. Nat. Mat. 10, 591. 10.1038/nmat3054 [PubMed] [Cross Ref]
  • Park C., Seo Y., Jung J., Kim D.-W. (2008). Electrode-dependent electrical properties of metal/Nb-doped SrTiO3 junctions. J. Appl. Phys. 103, 054106 10.1063/1.2872707 [Cross Ref]
  • Querlioz D., Bichler O., Dollfus P. (2013). Immunity to device variations in a spiking neural network with memristive nanodevices. IEEE Trans. Nanotechnol. 12, 288 10.1109/TNANO.2013.2250995 [Cross Ref]
  • Querlioz D., Bichler O., Gamrat C. (2011). Simulation of a memristor-based spiking neural network immune to device variations, in Neural Networks (IJCNN), The 2011 International Joint Conference on (San Jose, CA: IEEE; ), 1775–1781.
  • Querlioz D., Bichler O., Vincent A. F., Gamrat C. (2015). Bioinspired programming of memory devices for implementing an inference engine. Proc. IEEE 103, 1398–1416. 10.1109/JPROC.2015.2437616 [Cross Ref]
  • Sawa A. (2008). Resistive switching in transition metal oxides. Mater. Today 11, 28 10.1016/S1369-7021(08)70119-6 [Cross Ref]
  • Shanahan T. (2004). The Evolution of Darwinism: Selection, Adaptation and Progress in Evolutionary Biology. Cambridge: Cambridge University Press.
  • Sheridan P., Ma W., Lu W. (2014). Pattern recognition with memristor networks, in Circuits and Systems (ISCAS), 2014 IEEE International Symposium on (IEEE), 1078–1081.
  • Strukov D. B., Snider G. S., Stewart D. R., Williams R. S. (2008). The missing memristor found. Nature 453, 80–83. 10.1038/nature06932 [PubMed] [Cross Ref]
  • Szot K., Speier W., Bihlmayer G., Waser R. (2006). Switching the electrical resistance of individual dislocations in single-crystalline SrTiO3. Nat. Mater. 5, 312. 10.1038/nmat1614 [PubMed] [Cross Ref]
  • Waser R., Dittmann R., Staikov G., Szot K. (2009). Redox-based resistive switching memories - nanoionic mechanisms, prospects, and challenges. Adv. Mater. Weinheim. 21, 2632 10.1002/adma.200900375 [Cross Ref]
  • Würtz R. P. (2008). Organic Computing, Understanding Complex Systems. Berlin; Heidelberg: Springer.
  • Yang J. J., Strukov D. B., Stewart D. R. (2013). Memristive devices for computing. Nat. Nanotechnol. 8, 13. 10.1038/nnano.2012.240 [PubMed] [Cross Ref]
  • Zahari F., Hansen M., Mussenbrock T., Ziegler M., Kohlstedt H. (2015). Pattern recognition with TiOx-based memristive devices. AIMS Mater. Sci. 2, 203–216. 10.3934/matersci.2015.3.203 [Cross Ref]
  • Zamarreño-Ramos C., Camu-as-Mesa L. A., Pérez-Carrasco J. A., Masquelier T., Serrano-Gotarredona T., Linares-Barranco B. (2011). On spike-timing-dependent-plasticity, memristive devices, and building a self-learning visual cortex. Front. Neurosci. 5:26. 10.3389/fnins.2011.00026 [PMC free article] [PubMed] [Cross Ref]
  • Ziegler M., Riggert C., Hansen M., Bartsch T., Kohlstedt H. (2015). Memristive hebbian plasticity model: device requirements for the emulation of hebbian plasticity based on memristive devices. IEEE Trans. Biomed. Circ. Syst. 9, 197–206. 10.1109/TBCAS.2015.2410811 [PubMed] [Cross Ref]

Articles from Frontiers in Neuroscience are provided here courtesy of Frontiers Media SA