Comparing memristor models.

## Abstract

The memristor has been theoretically investigated as one of the fundamental electrical elements by Pr. Leon Chua in 1971. Meanwhile, its electrical characteristics are not yet fully understood. The nonlinear characteristics and the ability to examine large-scale amounts of storing data of this device reveal an interesting development in emerging electronic systems. Research on memristor modeling based on SPICE tools has grown rapidly. This leads us to study the behavior of such devices. Our aim is to simulate different types of memristor behavior. The adjustment of the model is based on the implementation of several parameters, which enables the switching of this device. In this chapter, we prove the flexibility and the correlation of memristor model with different memristive characterization data, by applying different voltage bias, sinusoidal and with a repetitive sweeping. Moreover, we demonstrate the memristor behavior as four types of switching. This includes bipolar switching, unipolar switching, bipolar switching with forgetting effect, and a reversible process between bipolar and unipolar switching. In order to validate this study, we compare our simulation results with experimental data and we prove a good agreement. The SPICE model used in our simulations shows a special advantage for its flexibility and simplicity.

### Keywords

- memristor
- I-V characteristics
- SPICE model
- switching behavior
- hysteresis loop
- bipolar behavior
- unipolar behavior

## 1. Introduction

Significant interest has been focused on the development of memristor-based systems. It has been first developed on symmetry consideration by Prof. Leon Chua in 1971 [1, 2]. In addition, it has been admitted physically by the HP Labs Team in 2008 [3]. This device does have a great potential to be the future memory cell, due to the small feature size and ability to retain the content (nonvolatile). The identity of such device is obvious on the *I-V* characteristics, i.e., its “pinched hysteresis loop.” Thus, the choice of the model and the structure are necessary to achieve better endurance and performance. Hence, the correlation of one model to other memristive devices is an interesting development to further research.

In the literature, memristor models studied in [3, 4, 5, 6, 7, 8, 9, 10] have been published for basic mathematical functioning properties of the memristor, which have been proposed by HP Labs in [3]. Other models [11, 12, 13, 14] focus on extracting the *I-V* characteristic of the model with other mathematical method using boundary conditions. They differ in complexity, materials, and accuracy. Thus, since our interest is in the behavior of the memristor, we choose to explore and investigate a simple SPICE model, which has been proposed by the present authors in [11]. The main differences are displayed in the implementation of parameters such as the state variable of the device. However, so far, no SPICE model could be correlated to several characterizations data of memristive devices. Our goal is to use a memristor model to analyze its functioning for different voltage bias. We study the dynamical behavior of memristor and we demonstrate that this model accounts for four different types of a memristor cited as the following: the bipolar behavior of memristor, the unipolar behavior, the bipolar with forgetting effect, and the reversible process between the bipolar and the unipolar behavior of memristor. Those types of memristor change under distinct stimulus such as sinusoidal, triangular, and repetitive DC sweeping voltage.

## 2. Theoretical principles

The wide variety in memristor structure and composition has led to the development of many different memristor modeling techniques. Some of them have been designed to represent a specific device for a specific type of application, such as AHaH [12], ANN [15, 16], Slime mold [17], and neuromorphic applications [5]. Implementation of the memristor could be generated on several tools of simulation, such as SPICE [18, 19, 20, 21, 22, 23, 24, 25], Matlab [26, 27, 28, 29, 30, 31, 32], Verilog-A [33], and VHDL-AMS [34, 35, 36, 37]. Resistive switching behavior is one of the fundamental properties showed in memristors; the well-known HP lab model of a memristor [3] shown in Figure 1(a) consists of a thin TiO_{2} double-layer of width D between a pair of platinum electrodes. One of the TiO_{2} layers of width w is doped with oxygen vacancies. The second, undoped layer of width *w*-*D* has insulating properties. As a result of complex processes in the device, the width w of the doped layer varies by applying a voltage or current to the electrodes of the memristor, and there will be dramatic changes in resistance. Therefore, the boundary, defined as the state variable *x* = *w*/*D*, between the two layers moves simultaneously. The well-known characteristic of the memristor is shown in Figure 1(b), the pinched hysteresis loop, which indicates the switching behavior of memristive devices. An application of a positive bias voltage to the electrodes of the device leads to the switching between Off and On states, and this switching is labeled SET. A RESET switching corresponds to the exchange between On and Off states. As current flows through the device, the cross section between the regions moves. As a result, the doped and the undoped regions have resistance R_{on} and R_{off} when each of them reaches the (D-w) and the full-length D, respectively. Also, the width w of the doped region of the memristor increases by applying a positive voltage bias, which causes the total resistance of the device to decrease. The same process is carried out by applying a negative voltage to the opposite side of the device. Moreover, there are two methods of the behavior of resistive switching for memristors: static and dynamic switching.

### 2.1 Static characteristics

The characteristic of a static switching behavior is obtained with a slow sweep of the voltage applied to the terminals of the device between the minimum and maximum values eligible (typically a triangular signal).

### 2.2 Dynamic characteristics

In dynamic switching, voltage pulses are applied to the device and the current rises under the constant voltage bias during the pulsing interval.

A comprehensive mathematical illustration of a SPICE memristor model has been reported in [11], which will be used later on for our simulation results. This model can illustrate the static and dynamic switching behavior, which will be studied in the next section. Thus, this model is based on the assumption that the switching behavior of the memristor is small or fast, below or above a threshold voltage V_{SET} or V_{RESET}, respectively, which is considered as the minimum voltage required to impose a change on the physical structure and thus the memristance of the device. This assumption is encapsulated in the use of the multiple implemented parameters, which are included in the set of equations below:

The relationship between the memristor voltage and the memristor current is given by Eq. (1), and it comprises three main parameters: *a*_{1}, *a*_{2}, and *b*. These parameters are responsible for the modeling of the nonlinear phenomenon of the pinched hysteresis loop. *a*_{1} and *a*_{2} are the magnitude parameters that vary according to the polarity of the input voltage; it is also related to the thickness of the dielectric layer of the memristor. Meanwhile, *b* is defined as the control parameter, which refers to the amount of oxygen deficiencies presented in the device, and it controls the conductivity of the device. The main voltage equation is defined by the relation *g*(*t*) defined below:

Equation (2) incorporates the threshold voltage with *V*_{p} and *V*_{n} which refers to the positive and negative polarizations, respectively, which makes a change in the switching behavior for value below the external voltage of the memristor. *A*_{p} and *A*_{n} are fitting parameters that affect the conductivity of the device. Accurately, it controls the speed of the oxygen deficiencies motion. The demonstration of the linearity of the model is described by parameters included in the following equations:

The physical parameters *x*_{p} and *x*_{n} have been defined in Eqs. (3) and (4); it represents the value of the state variable, which is responsible for the linearity of the device. Fitting parameters *α*_{p} and *α*_{n} are also included in these equations; are responsible for the linearity of the device; and they determine the degree of motion including the amortization of the state variable. The parameters w_{n} and w_{p} are defined by Eqs. (5) and (6), respectively. Those functions are used to shape the intensity of the state variable dynamics, i.e., the rate of memristance change.

Equation (7) represents the modeling function of the state variable. The fitting parameter *η* represents the direction of the movement of the state variable depending on the polarity of the input voltage. When *η* = 1, a positive voltage greater than the threshold voltage will increase the value of the state variable; and when *η* = −1, a positive voltage will decrease the value of the state variable.

Each pair of the parameters indicates the variation in the positive and negative region of the polarization. These multiple parameters make it possible for this device to be adaptable to a variety of characterization data of memristive devices, which we will discuss in the next section.

## 3. Analysis of the *I-V* characteristics

In our previous work [38], we illustrated a methodology for a simple memristor model to automatically adjust other behaviors of memristive devices. It effectively demonstrates the basic *I-V* characteristics of a memristive device. In addition, it acts differently in the positive and negative regions of the applied voltage, and the implemented parameters of the device take account on this, which makes the analysis of the pinched hysteresis loop simple and coherent for the positive and the negative regions independently. Therefore, we analyze the fundamental fingerprint of the macromodel and its memristance switching behavior. In the simulation results of Figure 2(a), we used a sinusoidal voltage 0.46 V with a frequency of 100 Hz, Figure 2(b) shows the resultant pinched hysteresis loop, which correlates the characterization data of the proposed model [11]. Thereby, our results agree well with the experimental results already published in [11], and we prove the linearity property of the device for a higher value of frequency. The next simulation results reveal the richness of memristor’s switching behavior confirming the usefulness of the specific design approach. The effect of memristive switching is inspected by varying the implemented parameters of the model such as the magnitude of the voltage bias, the initial charge, and the state variable.

This changes the operating regime so the memristance value may not remain constant and the memristor operates in different segments or takes different memristance values. This sudden jump of memristance is called “memristive switching” or “resistive switching.”

In this case, memristive switching depends on the bias of the applied voltage across the device, which is represented in Figure 3(a), the curve of the state variable motion at memristor boundaries. We consider the memristor in an Off state, as an initial state of the device, switching the device to On state, requires a positive bias across the device. While switching it to Off state requires negative bias. Figure 3(b) represents the curve of the memristance or the resistance of the memristor on the On and Off states. From this results, R_{on} and R_{off}’s value estimated by 1.3 *k*Ω and 1.1 *k*Ω, respectively.

## 4. Correlation of different memristive devices

In this section, we present our simulation results, in which we implement the different values of parameters on PSPICE, to fit a set of memristive devices studied for different types of applications. Those results describe the static and dynamical characteristics of the model. Thus, we prove that the SPICE model fits well with the characterization data of memristors defined in [8, 15, 17, 25, 39, 40, 41, 42, 43, 44]. The polarization voltages studied are either sinusoidal pulses or repetitive DC sweeping voltage to represent the different switching resistive levels of the memristor. The simulation results of the proposed model [11] shown in Figure 4, which indeed shows the characterization data of several memristive devices. In these simulation results, we adjust the different implement parameters on the SPICE model to fit the experimental results presented in [11]. Thus, it describes the *I*-*V* characteristics for devices defined in [8, 9, 40, 41, 42, 43], which has been correlated by the SPICE model.

Figure 4(a) describes the simulation results of the device published by the State University of Boise in [40].

Figure 4(b) describes the simulation results of the device published by the Tel Aviv University in [9, 41].

Figure 4(c) describes the simulation results of the device published by the University of Michigan in [8].

Figure 4(d) describes the simulation results of the device published by the state University of Iowa in [42].

Figure 4(e) describes the simulation results of the device published by the University of Michigan in [43].

### 4.1 Memristive device of the laboratory of slime mold

We adjust the implemented parameters to find the appropriate shape of *I*-*V* characteristics of the Slime mold device [17] shown in Figure 5, our results fits well the experimental results described in [17]. The application of Slime mold is a group of bacteria that lives mainly in the soil, which has the ability to change its shape by sliding every 50 s (by extension and retraction). This outcome contributes to the development of bioelectronics circuits of self-growth. The *I*-*V* characteristics curve for a DC voltage and a repetitive sweeping and the curve of the resistance of the device shown in Figure 5(a)–(c), respectively. These results present the functionality of this model by applying a repetitive DC sweeping voltage to present the various resistance switching states.

### 4.2 Memristive device of Strachan of the HP laboratory

Another memristive device based on TaO*x* was proposed by the team of HP Labs in [44]; we adjust the fitting parameters with the characterization data of this device. The results are shown in Figure 6; it agrees with the experimental results represented in [44]. The *I-V* characteristic curve for a DC voltage and a repetitive sweeping and the curve of the resistance of the device are shown in Figure 6(a)–(c), respectively. The simulation results present the functionality of this model by applying a repetitive DC sweeping voltage to present the several resistance switching states.

### 4.3 Memristor of Nugent for AHaH applications

Our model also fits well a learning AHaH application accomplished by Nugent in [12]. After the application of a sinusoidal signal to the memristor with amplitude 0.25 V for a period 10 *m*s, we found the resultant *I*-*V* characteristics shown in Figure 7(a), which seems compatible with the *I*-*V* characterization profile and it fits well the experimental results revealed in [12]. This model works well for a repetitive DC sweeping voltage which is represented by the *I*-*V* characteristics curves and the curve of the resistance of the device shown in Figure 7(b), (c). The simulation results present the functionality of this model by applying a repetitive DC sweeping voltage to present the several resistance switching states.

### 4.4 Memristive device of the University of Pittsburgh

The device represented by Zhang in [39] is based on TaO*x* material. The simulation results of this device are represented in Figure 8(a) by applying a triangular voltage with a sweeping of the magnitude, 0.74 V for the positive region, and −1.25 V for the negative region. After adjustment of the parameters, we have a nonlinear *I*-*V* curve which seems to be compatible with the *I*-*V* characterization profile recorded in the experiments of HP labs in [39]. In addition, we proved the functioning of this device with a repetitive DC sweeping voltage to present several resistance switching states shown in Figure 8(b) and (c).

### 4.5 Memristive device for ANN learning application

Meanwhile, we also simulated the model with a square wave excitation shown in Figure 9. This excitation method is used as a learning method, which presents the behavior of this model as an artificial neural network ANN [15]. We follow the learning experience carried out in [16]. As shown in **Figure 9(c)**, the memristance of the device increases along with the applied voltage. However, this behavior response is different from the other previous excitation, and this depends on the type of excitation and the followed current. The current curve decreases with each pulse of the excitation voltage positive and negative.

In conclusion, we notice that this SPICE model is a general model that can be applied in multiple domains. Furthermore, according to the simulation results of both devices Slime mold and HP Labs, we notice that the hysteresis loop of the memristor maintains its nonlinear shape even for a remarkable range of the values of the parameters. The speed of movement for the memristive devices of Slime mold and Nugent are faster compared to the other memristive devices since A_{p} and A_{n} have higher values. Thus, the memristive devices of Slime mold and Zhang application have the lowest values of b, which decrease its conductivity. The ANN learning application presents another type of excitation, which largely affects the dynamic of the memristor’s behavior. Thus, we will deal later with the use of different memristor switching behavior, we demonstrate not only bipolar, but also the unipolar switching behavior of memristors, which differs from bipolar memristors in the fact that only the magnitude of the voltage across the device determines the change in the resistance.

## 5. Behavior of the SPICE model for different types of memristors

Memristor models, in literature, have different responses, which are generated for four different types of a memristor, i.e., bipolar, bipolar with a forgetting effect, unipolar and reversible behavior between the bipolar, and the unipolar memristor. For our simulation results, we used a SPICE model that can not only describe the basic memory ability of memristor, but also be able to capture all of the four types of memristor switching behavior.

Models with bipolar switching [45, 46] distinguishable by the memristance which increases and decreases by different polarity voltages. Models with unipolar switching behavior [45, 46] are distinguishable by its memristance, which can increase and decrease by the same polarity voltage. The bipolar with forgetting effect [47, 48] is distinguishable by its memristance which increases and decreases by a different polarity voltage, but memristance can spontaneously decrease at the same time, even with no voltage applied. The reversible bipolar and unipolar switching behavior [49], here the memristor will behave like a bipolar memristor at first, but after a few iterations, it will turn to a unipolar memristor.

In the same context, we use different polarization voltages, either sinusoidal or repetitive DC sweeping voltage, exploited in order to present the different states of resistance of the memristor, and thus it shows the behavior of the four different types of a memristor. Those results reveal the richness of memristor’s dynamical behavior confirming the usefulness of the specific model approach.

To verify the memristive characteristics and the coexistence of different switching behavior of our proposed model, we employed different excitations presented in the following figures of the rest of the paper. In fact, to characterize different types of memristors, we need to verify the behavior of the model for the well-known fundamental switching behavior in both bipolar and unipolar switching behavior. In this case, we have adapted our model according to the experimental results demonstrated in [50]. These observations are consistent and in very good qualitative agreement with the experimental results of the memristor switching behavior already published in [50].

In addition, Figure 10 shows the dynamical characteristics of a bipolar memristor behavior. The sweeping voltage bias approaches a set value v = 1.8 V with current value I = 7 *m*A. Reversing the voltage polarity, the device switches to a reset value at v = −1.8 V with current value I = −4.5 *m*A. The attained pinched hysteresis *I-V* curves are shown in Figure 10(b), (c), which are the typical fingerprint of bipolar resistive switching. The corresponding resistance response is illustrated in Figure 10(d), which was measured according to voltage sweeping with a maximum value of 1.2 *k*Ω and a lower value of 0.9 *k*Ω in the negative and positive voltage application respectively. In addition, Figure 10(e) confirms that hysteresis takes place in both *I-V* and *C-V* relationships of the device, and it shows a closed switching cycle in the bipolar switching behavior for a maximum value of 25 *m*S.

Memristor behavior for a bipolar switching with forgetting effect is shown in Figure 11. An obvious overlap of the *I-V* curve is shown in Figure 11(b), (c), due to the repetitive sweeping of the applied voltage. The sweeping voltage bias approaches a set value v = 1.2 V with a maximum current value I = 6 *m*A. Reversing the voltage polarity, the device switches to a reset value at v = −1.2 V with a maximum current value I = −5 *m*A. Also, it can be seen from these curves, an accumulation of the current on each pulse. The corresponding resistance response is illustrated in Figure 11(d), which was measured according to voltage sweeping with a maximum value of 1.8 *k*Ω decreasing to a lower value of 0.3 *k*Ω in the positive voltage application, on the opposite side of the negative voltage application the resistance response shows an increase from 0.3 to 0.9 *k*Ω. The curve in Figure 11(e) shows five switching cycles for a maximum value of 35 *m*S. This *C-V* curve shows that the memristance not only increases and decreases by a different polarity voltage, but it also can spontaneously decrease at the same time, even with no voltage applied, and this is a unique switching behavior of memristor. In fact, these curves show the operation of the model as a bipolar memristor with forgetting effect.

Furthermore, the simulation results for the unipolar behavior of memristor are shown in Figure 12, which show that another switching behavior is characterized by the memory devices and also that the memristance of the device can increase and decrease by the same polarity of the voltage. For this type of memristor, we use a positive voltage excitation for a value of 2 V and maximum current value I = 17 *m*A, which is shown in the curve (Figure 12(a)), we notice a slight accumulation of the current on each pulse. The characteristics shown in Figure 12(d), (e) describe the resistance and the conductance curve of the memristor. The resistance was measured with a maximum value of 2.8 *k*Ω decreasing to a lower value of 0.25 *k*Ω. However, the conductance curve (Figure 12(e)) shows three switching cycles related to voltage sweeping for a maximum value of 80 *m*S. In fact, this *C-V* curve shows that the conductance change in response to three positive pulses, it initially increases (during each pulse stimulus) and subsequently decays toward its original value (between stimuli).

In the end, we represent the results of memristor under a reversible state between the bipolar and the unipolar behavior in Figure 13. The sweeping voltage bias, shown in Figure 13(a), approaches a set value v = 5 V with a maximum current value I = 60 *m*A. Reversing the voltage polarity, the device switches to a reset value at v = −5 V with a maximum current value I = −10 *m*A.

An obvious overlap of the *I-V* curve in Figure 13(b),(c) occurs due to the repetitive sweeping of the applied voltage. The corresponding resistance response is illustrated in Figure 13(d), which occurs in a different switching behavior; the rise and fall of the resistance exist but with a large gap between high and low values. The conductance curve in Figure 13(e) shows four switching cycles for a maximum value of 100 *m*S. This *C-V* curve shows that the first cycle of the switching behavior differs to the other cycles of the switching behavior of the memristor. In fact, the first cycle that can be seen from these curves shows bipolar operation, but after the second pulse, it automatically turned to a unipolar memristor behavior. These characteristic curves are shown, respectively, in Figures 10, 11, 12, 13. The operation of the memristor model as a bipolar memristor behavior is shown in Figure 10. The response to a bipolar with forgetting effect is shown in Figure 11. The response to a unipolar memristor behavior is shown in Figure 12. And, the memristor model response to a reversible bipolar and unipolar behavior is shown in Figure 13. We can conclude that our simulation results are consistent and in very good qualitative agreement with the results already published in [15]. A detailed comparison between our work model and other popular memristor models (the Chua [1], the Strukov (HP) [3], Vourkas [45], and the Chen [15]) is shown in Table 1. We can notice that the SPICE model gets a special advantage on describing various memristors models with an average number of parameters, and for its flexibility, and low complexity.

## 6. Conclusion

The consideration of the SPICE memristor model as a simple and flexible model was proved to explain the memristor switching, not only processing the general memristor properties, but also catching the different types of memristor behavior: the bipolar, unipolar, the bipolar with forgetting effect, and the reversible process between the bipolar and the unipolar behavior. Our simulation results demonstrate that for the bipolar memristor, a regular hysteresis curve can be obtained. For the bipolar memristor with forgetting effect, an obvious overlap between the neighbor loops of the *I-V* curve, and for the unipolar memristor, a positive voltage is applied, but the conductance will increase only when the voltage is over 1 V. Also, for the reversible process between bipolar and unipolar behavior, the memristor firstly behaves as a bipolar switching, and its conductance increases and decreases according to the polarity of the voltage. However, after applying a second pulse, it will turn to behave as a unipolar switching. This chapter provides a practical memristor model that can be simulated with different types of stimulus, and further studies are aimed at integrating the memristor model into a computing design with complementary metal-oxide-semiconductor (CMOS) circuits that can perform the necessary functions on a chip.