Abstract
In this paper, a Hewlett-Packard (HP) memristor model with a new window function and its versatile characteristics are presented. SPICE behaviors of the linear and nonlinear memristor model are studied through PSpice simulation. High flexibility is demonstrated for emulating the behaviors of the practical HP memristors. Furthermore, the characteristics of the composite SPICE behaviors are both investigated when two memristors are connected in series and in parallel. The polarity of each memristor is also taken into consideration. The relationships among flux, charge, voltage, current, and memristance of the double memristor circuits are simulated and analyzed.
1 Introduction
There are three fundamental two-terminal passive circuit elements, i.e., resistors R, capacitors C and inductors L, which can be found and introduced from any basic circuit books. Chua first theoretically postulated the existence of the fourth passive electrical circuit element in 1971 [1], named as memristor, whose memristance (i.e., an acronym for memory resistance) M can be defined by the relation of magnetic flux and electrical charge q, namely, . Together with the definitions of three commonly used basic circuit elements R, C and L, and definitions of the electric current and Faraday’s law, memristance closes the loop among four variables v, i, and q with six equations. Under the linear situation of memristor, M is seen as a constant and memristance is the same as the resistance. In fact, the memristance depends on the amount of charge that passed through the memristor, resembling a resistor with time-varying resistance.
No great progress was developed in the research field of memristor until the group of Williams at Hewlett-Packard (HP) lab successfully fabricated the first nanoscale memristive device based on titanium dioxide thin film in 2008 [2]. Since then, many researchers or scholars have focused on the memristor research [3–7] and its numerous potential applications in nonvolatile random access memories [8, 9], neuromorphic systems [10, 11], chaotic circuits [12], reconfigurable logics [13, 14] and RF/microwave devices [15, 16].
In order to emulate these memristor applications, various memristor nonlinear drift models using SPICE have been presented before fabricated realization [17–24]. Window functions need to be used to tackle the boundary conditions of the memristors [17–22]. In ref. 17, a switching window function was proposed, where the current i is involved, but it has limited scalability. Two different window functions were presented to mimic nonlinear effects in refs. 18 and 20. However, as the boundary between the doped and undoped layers approaches either end of the memristive device, it cannot be excited by external stimulus. In refs. 20 and 22, the piecewise functions were utilized to construct the window functions, which were relatively complicated when applying in the memristive devices.
In this paper, a memristor model using a new window function, called cosine window function, is proposed along with its versatile characteristics. To illustrate the significance of the proposed window function, main characteristics of state-of-the-art ones in recent years are listed in Table 1 for comparisons. Then, the SPICE behaviors of the linear and nonlinear memristor model are investigated through PSpice simulation. Furthermore, the characteristics of the composite SPICE behaviors are both studied when two memristors are connected in series and in parallel. The relationships among flux, charge, voltage, current, and memristance are simulated and analyzed using PSpice.
TABLE 1
| Function | Strukov ref. 2 | Biolek ref. 17 | Joglekar ref. 18 | Prodromakis ref. 19 | Yu ref. 20 | Zha ref. 21 | Tukey ref. 22 | Proposed |
|---|---|---|---|---|---|---|---|---|
| Symmetric | Yes | Yes | Yes | Yes | Yes | Yes | Yes | Yes |
| Resolve boundary conditions | Yes | Discontinuities | Yes | Yes | Yes | Discontinuities | Yes | Yes |
| Impose non-linear drift over entire D | Yes | Yes | Yes | Yes | Yes | Yes | Yes | Yes |
| Linkage with linear dopant drift model | No | Yes | Yes | Yes | Yes | Yes | Yes | Yes |
| Scalable 0 ≤ fmax(x) ≤ 1 | No | Limited | No | Yes | No | Yes | Yes | Yes |
| Flexibility (control parameter) | No | Yes | Yes | Yes | Yes | Yes | Yes | Yes |
Comparisons of some reported window function in memristor models.
2 Memristor Model Using Cosine Window Function
The concept of memristive devices or systems was generalized by Chua and Kang using the following equations [25],where v is the voltage, i is the current, and is the instantaneous memristance that depends on the internal state variable x and current i of the memristor. The state variable x, ranging from 0 to 1, is the normalized width of the memristor doped region , with D being the whole switching bilayer thickness of the memristor, as seen in Figure 1.
FIGURE 1
At time t, the width of the doped region w depends on the amount of charge that passed through the memristor. Thus, the time derivative of w is proportional to the current i, which can be expressed bywhere is the average dopant mobility, denotes the low resistance at high dopant concentration region. If we set , Eq. 3 can be changed as follows,The Eqs. 3, 4 is actually an expression of a linear drift model of memristor, but it does not take the boundary effects of the memristor into account. Nonlinear dopant drift phenomena exist at the boundaries of the physically realized nanoscale memristor, therefore, an appropriate window function f(x) need to be introduced in ref. 4, i.e.,where f(x) ensures zero drift at the boundaries, i.e., f (0) = f (1) = 0.
On the other hand, the memristance can be expressed bywhere denotes the high resistance at low dopant concentration region. In this paper, a new cosine window function is proposed as follows,where p is the control parameter. This new and simple cosine window function is demonstrated in Figure 2 for various values of p. As the parameter p rises, a unity window function f(x) will occur, and thus the model tends to a linear model. Here, the parameter p is set with a large range of variation, because the cosine window function is not changed significantly near the region of boundary.
FIGURE 2
Using this new cosine window function, the behaviors of the memristor can be observed in Figures 3–6 with versatile characteristics. Here, we set = 100 Ω, = 16 kΩ, w0 = 5 nm, D = 10 nm, and = 10–14 m2/Vs, thus the initial state variable x0 is equal to 0.5 at t = 0. When the voltage biased at two sides of memristor is set as , i.e., the frequency of the AC voltage f = 1 Hz, the calculated dynamical properties of the normalized doped state variable with different values of p are illustrated in Figure 3A. For instance, when p = 16, the state variable periodically saturates between the normalized boundary “1” from 0.295 + N to 0.705 + N s, where N is a natural number starting from 0. On the other hand, according to (6), the curves of the corresponding memristances can be obtained as seen in Figure 3B. When the memristor biased voltage is set as , the calculated dynamical properties of the normalized doped state variable with different values of p is shown in Figure 3C, and the corresponding memristance curves are plotted in Figure 3D.
FIGURE 3
Actually, the minimum state variable corresponding to maximum memristance is dependent on the frequency of the biased AC voltage. From Figure 3C, we can see that the state variable does not arrive the normalized boundary “0”. Through decreasing the frequency of the biased voltage, two edge states, i.e., saturation (w = D) and depletion (w = 0), can be reached when the memristor device is biased with positive voltage sinusoidal stimulus and reversed polarity stimulus, respectively. Figure 4 shows the influence of the frequency variation of the biased voltage on the state variable x and memristance M. As the frequency of the AC voltage f decreases from 1 Hz to 1/3 Hz, the minimum value of state variable x will approach 0 and the corresponding maximum memristance will be 16 kΩ, under the condition of p = 8.
FIGURE 4
Figure 5 illustrates a nonlinear calculated current-voltage (i-v) curve with p set as 8, where a periodic pinched hysteresis loop with hard switching emerges. This memristive model has zero-crossing property in a form of i-v Lissajous figure, which represents no current through the system when the voltage value is zero. Figure 6 shows the versatility of our model, where the i-v characteristic curves are plotted for eight different values of p. This model demonstrates high flexibility for emulating the behaviors of the practical HP memristors, which can be beneficial in characterizing memristive devices before fabrication.
FIGURE 5
FIGURE 6
3 SPICE Behaviors of Double Memristor Circuits
The HP memristor model using cosine window function and its characteristics have been addressed in the above section. To further investigate the relationships among charge, flux, and memristance of composite memristors connected in different topologies, simulated SPICE behaviors of single and double memristor circuits will be presented in details in this section.
3.1 SPICE Behaviors of Single Memristor
The single memristor circuit using the above-mentioned cosine window function has been proposed, and the corresponding SPICE code was written for use in PSpice [26], as seen in Figure 7. The single memristor’s SPICE behaviors can be seen in Figures 8, 9, 10, where the initial state of the linear and nonlinear memristor models are both M(0) = 15.9 kΩ with a sinusoidal input signal of . The behaviors in Figures 9, 10 are plotted together with those of two serial memristor circuits for comparisons, which will be introduced in subsection B.
FIGURE 7
FIGURE 8
FIGURE 9
FIGURE 10
As shown in Figure 8A, when the value of p increases, the curves of the nonlinear models will be toward that of linear model, which means that the nonlinearity will be weakened. At the middle area of the whole charging process, the slopes of all the curves with various p are identical including that of the linear model. In Figure 8B, as the value of p increases for a fixed electrical charge, the magnetic flux will decrease, i.e., the nonlinearity will be weakened.
3.2 SPICE Behaviors of Two Serial Memristor Circuits
Then, we construct a circuit with two nonlinear memristors (p = 8), in series with same polarities in PSpice as shown in Figure 11A. The initial states of these two memristors are both M(0) = 11.25 kΩ. Due to the movement of the doping ions caused by the bias voltage, the memristance associated with the doping boundary in the HP model definition will be changed accordingly. Therefore, when the magnitude Vm of the sinusoidal input signal is positive, the memristance will keep decreasing, and vice versa. Simultaneously, the value of Vm will influence the migration rate of doping ions.
FIGURE 11
Figure 9A shows the relationships between voltage and current of the single memristor (M1 or M2) and two serial memristors with same polarities in Figure 11A, where the magnitude of the sinusoidal input signal Vm = 2.3V. Due to the two memristors connected in series, equal current flows through both memristors. Therefore, the behaviors of each single memristor are identical. Moreover, the applied voltage is distributed across each memristor identically at any instant of time due to same polarities of these two memristors M1 and M2. Consequently, the slope of pinched hysteresis loop of the composite memristor circuit is half of that of the single memristor.
Figures 9B,C illustrate the curves of memristance versus charge (M-q) and flux versus charge (-q), respectively. In Figure 9B, the stored charge q in the memristor will be gradually increased as the memristance M decreases. The slope of the M-q curve of the single memristor is half of that of the composite memristor circuit. Meanwhile, the memristance M of the single memristor is always half of that of the composite memristor circuit at any time. When the memristance M decreases to the value of , namely, approaching the doping boundary of the memristor at high dopant concentration region, it will keep “ON” of this memristor even if the charge increases. Similarly, the slope of the -q curve of the single memristor is also half of that of the composite memristor circuit, as shown in Figure 9C.
For a two-serial-memristor circuit with opposite polarities in PSpice shown in Figure 11B, the current versus voltage curves of the individual memristors and two serial memristors are plotted in Figure 10A. The initial states of the two memristances are both set as M(0) = 11.25 kΩ in PSpice software. In contrast to the case of the same polarities in Figures 9A, 10A has three separate hysteresis loops, denoting the individual memristor M1, M2 and composite M of the two memristors. Due to the two memristors connected with opposite polarities, they have the opposite variations. The specific expression is that when the memristance M1 decreases, the other memristance M2 will increase, and vice versa.
Figure 10B shows the relationship between memristance and charge. The initial memristances of M1 and M2 are both M(0) = 11.25 kΩ. Therefore, the memristance of composite M is 22.5 kΩ at first. In the linear region, the amplitude changes of the two memristors are identical, but the trends are opposite. Therefore, the memristance of composite M remains unchanged. When the shift is close to the boundary, memristors are in the nonlinear region. Thus, the memristance of composite M will decrease accordingly. This phenomenon is caused from the boundary effects in physical memristor model. At this time, the two memristors reach the minimum and maximum memristances, respectively, resulting in the composite memristance of approximately 16 kΩ. For the relationship between flux and charge in Figure 10C, it is also different from that of Figure 9C. The magnetic fluxes of M1 and M2 are not the same any more due to the two memristors connected with opposite polarities. In Figure 10C, the slope of composite flux in black is always equal to the arithmetic sum of the two slopes of individual memristor fluxes. Moreover, the composite flux is varied more linearly as a function of charge, compared with those of the individual memristor fluxes.
3.3 SPICE Behaviors of Two Parallel Memristor Circuits
Similarly, the parallel circuits of two nonlinear memristors (p = 8) are constructed in PSpice, where the initial states of the memristiors are M1(0) = M2 (0) = 11.25 kΩ. Since the two memristors are connected in parallel with same and opposite polarities as shown in Figures 12A,B, respectively, equal voltage passes through the two memristors, where the magnitude Vm of the sinusoidal input signal is set as 1.15V.
FIGURE 12
The relationships between voltage and current of the single memristor (M1 or M2) and two memristors in parallel with same polarities are shown in Figure 13A. Since the initial states of the two memristors are identical, the applied voltage is distributed across each memristor equally at any instant of time. Hence, the behaviors of each single memristor are identical. In addition, the slope of pinched hysteresis loop of the composite memristor circuit is twice as big as that of the single memristor.
FIGURE 13
Figure 13B exhibits the variations of memristance versus flux (M-φ). In contrast to composite memristor circuit in series, the memristance of composite memristor circuit in parallel is half of that of the individual memristor at any instant of flux. As the memristance decreases, the magnetic flux will rise accordingly and its variations (i.e., slopes) for single memristor and composite memristor circuit are both increasingly large until the memristance reaches the minimum value. Figure 13C plots the charge versus flux (q-φ) curves of the single memristor and two memristors in parallel with same polarities. Seen from Figure 13C, the slope of the q-φ curve of the single memristor is also half of that of the composite memristor circuit.
For a double memristor circuit connected in parallel with opposite polarities shown in Figure 12B, the curves of current versus voltage are plotted in Figure 14A. It has three different hysteresis loops in contrast to the case of the same polarities in Figure 13A. The inset of Figure 14A shows that the curve of M2 has a relatively linear hysteresis loop compared with those of the other two. However, the composite current is always equal to the sum of the individual branch currents of the two memristors at any time whatever the applied voltage is.
FIGURE 14
The curves of memristance versus flux (M-φ) and charge versus flux (q-φ) are plotted in Figures 14B,C, respectively. As seen in Figure 14B, the variation trends of the two individual memristors M1 and M2 are opposite, due to opposite polarities of these two memristors. Additionally, the composite memristance is always lower than the memristance of each single memristor at any time. Seen from Figure 14C, it can be found that the q-φ curves of the two individual memristors with opposite polarities are symmetrical about the initial origin point. Moreover, the slope of composite charge is always equal to the arithmetic sum of the two slopes of individual memristor charges.
4 Conclusion
A new window function has been created and applied into the HP memristor model in this paper. To validate its advantages, versatile characteristics and SPICE behaviors of the memristor have been demonstrated. Moreover, we also investigate the composite SPICE behaviors when two memristors are connected in series and in parallel with same or opposite polarities, and make comparisons with those of the individual memristors. Their relationships among flux, charge, voltage, current, and memristance are analyzed using PSpice simulation.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This work was supported in part by the FY2019 JSPS Postdoctoral Fellowship for Research in Japan under Grant P19350, Grant-in-Aid for JSPS Research Fellow under Grant JP19F19350, National Natural Science Foundation of China under Grant 61601390, Opening project of Key Laboratory of Cognitive Radio and Information Processing under Grant CRKL180202, Natural Science Foundation of Guangxi under Grant 2019GXNSFFA245002, and Dean Project of Guangxi Key Laboratory of Wireless Wideband Communication and Signal Processing under Grant GXKL06190118.
Conflict of interest
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.
References
1.
ChuaL. Memristor-The missing circuit element. IEEE Trans. Circuit Theory (1971) 18(5):507–19. 10.1109/tct.1971.1083337
2.
StrukovDBSniderGSStewartDRWilliamsRS. The missing memristor found. Nature (2008) 453(7191):80–3. 10.1038/nature06932
3.
ChanthboualaAGarciaVCherifiROBouzehouaneKFusilSMoyaXet alA ferroelectric memristor. Nat Mater (2012) 11(10):860–4. 10.1038/nmat3415
4.
YakopcicCTahaTMSubramanyamGPinoRERogersS. A memristor device model. IEEE Electron Device Lett (2011) 32(10):1436–8. 10.1109/led.2011.2163292
5.
AdhikariSPSahMPKimHChuaLO. Three fingerprints of memristor. IEEE Trans. Circuits Syst (2013) 60(11):3008–21. 10.1109/tcsi.2013.2256171
6.
RzigaFOMbarekKGhediraSBesbesK. The basic I–V characteristics of memristor model: simulation and analysis. Appl Phys A (2017) 123:288. 10.1007/s00339-017-0902-9
7.
SinghJRajB. Comparative analysis of memristor models and memories design. J Semicond (2018) 39(7):074006. 10.1088/1674-4926/39/7/074006
8.
MeijerGI. Materials science. Who wins the nonvolatile memory race?Science (2008) 319:1625–6. 10.1126/science.1153909
9.
PalSBoseSKiW-HIslamA. Design of power- and variability-aware nonvolatile RRAM cell using memristor as a memory element. IEEE J Electron Devices Soc (2019) 7:701–9. 10.1109/jeds.2019.2928830
10.
JoSHChangTEbongIBhadviyaBBMazumderPLuW. Nanoscale memristor device as synapse in neuromorphic systems. Nano Lett (2010) 10(4):1297–301. 10.1021/nl904092h
11.
LiuNYangGHeYMaGChenAChenQet alRealization of synapse behaviors based on memristor and simulation study with KMC method. IEEE J Electron Devices Soc (2020) 8:981–5. 10.1109/jeds.2020.3023015
12.
ChenJ-JYanD-WDuanS-KWangL-D. Memristor-based hyper-chaotic circuit for image encryption. Chin Phys B (2020) 29(11):110504. 10.1088/1674-1056/abbbfe
13.
XiaQRobinettWCumbieMWBanerjeeNCardinaliTJYangJJet alMemristor-CMOS hybrid integrated circuits for reconfigurable logic. Nano Lett (2009) 9(10):3640–5. 10.1021/nl901874j
14.
WeiFCuiXCuiX. An improved iMemComp OR gate and its applications in logic circuits. IEEE J Electron Devices Soc (2020) 8:57–61. 10.1109/jeds.2019.2962822
15.
XuKDZhangYHWangLYuanMQFanYJoinesWTet alTwo memristor SPICE models and their applications in microwave devices. IEEE Trans Nanotechnology (May. 2014) 13(3):607–16. 10.1109/tnano.2014.2314126
16.
XuKDLiXSGuoYJLiuQH. Simple memristive SPICE macro-models and reconfigurability in filter and antenna. Radioengineering (2016) 25(4):700–6. 10.13164/re.2016.0700
17.
BiolekZBiolekDBiolkovaV. SPICE model of memristor with nonlinear dopant drift. Radioengineering (2009) 18(2):210–4.
18.
JoglekarYNWolfSJ. The elusive memristor: properties of basic electrical circuits. Eur J Phys (Jul. 2009) 30(4):661–75. 10.1088/0143-0807/30/4/001
19.
ProdromakisTPehPBPapavassiliouCToumazouC. A versatile memristor model with nonlinear dopant kinetics. IEEE Trans. Electron Devices (Sep. 2011) 58(9):3009–105. 10.1109/ted.2011.2158004
20.
YuJTMuXMXiXMWangSN. A memristor model with piecewise window function. Radioengineering (2013) 22(4):969–74.
21.
ZhaJHuangHLiuY. A novel window function for memristor model with application in programming analog circuits. IEEE Trans. Circuits Syst (2016) 63(5):423–7. 10.1109/tcsii.2015.2505959
22.
TakahashiYSekineTYokoyamaM. SPICE model of memristive device using Tukey window function. IEICE Electron Express (2015) 12(5):20150149. 10.1587/elex.12.20150149
23.
ThomasSPrakashSAn accurate analytical memristor model for SPICE simulators. IEICE Electronics Express (2018) 15(18). p. 20180724. 10.1587/elex.15.20180724
24.
JangJTMinJKimDParkJChoiS-JKimDMet alA highly reliable physics-based SPICE compact model of IGZO memristor considering the dependence on electrode metals and deposition sequence. Solid-State Electronics (2020) 166:107764. 10.1016/j.sse.2020.107764
25.
ChuaLOSung Mo KangS. Memristive devices and systems. Proc IEEE (Feb. 1976) 64(2):209–23. 10.1109/proc.1976.10092
26.
PSpice. Cadence © [online]. Available: https://www.pspice.com/.
Summary
Keywords
memristive devices, memristor circuits, simulation program with integrated circuit emphasis model, window function, memristor
Citation
Xu K-D, Li D, Jiang Y and Chen Q (2021) SPICE Behaviors of Double Memristor Circuits Using Cosine Window Function. Front. Phys. 9:648737. doi: 10.3389/fphy.2021.648737
Received
01 January 2021
Accepted
21 January 2021
Published
03 March 2021
Volume
9 - 2021
Edited by
Gang Zhang, Nanjing Normal University, China
Reviewed by
Wen Pin, Nanchang University, China
Di Wu, Shenzhen University, China
Updates
Copyright
© 2021 Xu, Li, Jiang and Chen.
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(s) 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.
*Correspondence: Yannan Jiang, ynjiang@guet.edu.cn
This article was submitted to Radiation Detectors and Imaging, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.