Polarity Reversal Effect of a Memristor From the Circuit Point of View and Insights Into the Memristor Fuse

As the memristor device is asymmetrical in nature, it is not a bilateral element like the resistor in terms of circuit functionality. Thus, it causes hindrance in some memristor-based applications such as in cellular nonlinear network neighborhood connections and in some application areas where its orientation is essentially expected to act as a bilateral circuit element reliable for bidirectional communication, for example, in signal and image processing or in electrical synapse devices. We introduce a memristor-based network for each purpose where we replace the conventional series resistances by memristors. The memristor asymmetry is described from the circuit point of view allowing us to observe its interaction within the network. Moreover, a memristor fuse is proposed in order to achieve the memristive effect with symmetry, which is formed basically by connecting two memristors antiserially. We, therefore, analyze the memristor fuse from its basic principle along with the theoretical analysis and then observe the response from the circuit point of view.


INTRODUCTION
The memristor was predicted in 1971 (Chua) by observing the symmetrical nature of the three known basic circuit elements, resistor R, capacitor C, and inductor L, with respect to the four circuit variables, namely, electric voltage v, electric current i, electric charge q, and magnetic flux ϕ, see Figure 1. As stated in the work of Chua (1971) for the sake of completeness, there should be a fourth passive circuit element describing the relationship between magnetic flux ϕ and electric charge q, hence named memristor. A broader class of this device known as the memristive system is given by Chua and Kang (1976). Memristor (M) is the short form of memory resistor, this name being due to the fact that the device remembers its previous history (resistance), hence the memory effect, and is analogous to a resistor with memory. Depending on the type of excitation, a memristor can be described as charge controlled ϕ h(q) or flux controlled q h(ϕ) (Chua (2015)) which are preferably described as memristance and memductance, having unit as Ohms (Ω) and Siemens (S), respectively.
For more than 3 decades, memristor remained a mystery until in 2008 (Strukov et al. (2008)) a group of researchers from the HP laboratory announced the successful realization of the first solidstate memristor in a device form (Stanley Williams (2013)). This recent discovery of the HP lab allured many scientists, engineers, and researchers to explore the feasible applications of memristor in discrete and crossbar array configurations Mazumder et al. (2012) and more possible device technologies.
Since the invention by Strukov et al. (2008), many memristor technologies emerged which are basically adhered to the principle of bipolar resistance switching between two extreme values, namely, R on and R off that, respectively, correspond to the lowest and highest resistance state of the device. Note that another commonly used memristor technology is the selfdirected channel device (Campbell (2017)) whose conductivity is based on the formation and dissolution of ionic bridges that result in low and high resistance states, respectively, but we will mainly refer to the TiO 2 memristor in the remaining part of this article. Figure 2 shows the formation of the titanium-oxide memristor with the TiO 2 doped with some positive oxygen vacancies. Hence, the TiO 2 memristor is an example of MIM devices, i. e., metal-insulator-metal, in which a thin bilayer of TiO 2 film is sandwiched between platinum electrodes (labeled as Pt 1 and Pt 2 ). The small positive spots in the doped region refer to the positive charges due to the oxygen vacancies (Stanley Williams (2013); Paris and Taioli (2016)).
The mathematical description of the titanium-oxide memristor when a positive voltage V(t) is applied to the 2port device shown in Figure 2A corresponds to a current i(t) flowing into the memristor, while the voltage is shared in two parts: a voltage V(t) R on w/R on w + R off (D − w) across the doped region (the left part of Figure 2A) and the complementary part V(t)R off (D − w)/R on w + R off (D − w) present across the undoped region (the right part of Figure 2A). Let us consider the oxygen vacancies in the doped region. These charge carriers with mass m and charge q are accelerated by the electric field E V(t)R on /R on w + R off (D − w) and are broken as they collide together. They finally reach a limit speed V l q.E/mΓ, where Γ is the mean time between two consecutive collisions, leading to their averaged velocity < v oxygen vacancy > μ v E, where μ v is the mobility of the oxygen vacancies. They expand the doped region toward the right (w↑), whose boundary increases with a positive current i(t) such that dw/dt μ v E μ v R on /Di(t). Finally, with the normalized form x w/D, the behavior of the memristor is given by where M(x) is the memristance and δR R off − R on is the difference between R off and R on . When integrating (1c) for x from 0 to 1, where q d w D w 0 i(t)dt D 2 /u v R on is the charge required to move completely the doped/undoped boundary from w 0 to w D. Then, Eq. 1c can be rewritten as A window function f (x) is often introduced as a factor in the right-hand side of Eq. 3 for nonlinear dopant drift modeling, i.e., to avoid x from taking values outside of the interval [0, 1] (Joglekar and Wolf (2009) Eq. 1 characterizes a bipolar memristor where the resistance switching depends on the voltage polarity (Strachan et al. (2011), Krzysteczko et al. (2009), Teixeira et al. (2009). However, there are other reported memristors exhibiting symmetry in polarity, such as unipolar, nonpolar, and complementary resistive switching memristors (Yoshida et al. (2008), Huang et al. (2010), Wang et al. (2017), Linn et al. (2010). Here, the resistance switching between R off and R on (and vice versa) can be completed in the same voltage polarity. As such, unipolar memristors are important elements in memory arrays and logic circuit implementation (Yin et al. (2020)).
The interesting features of memristors, such as connection flexibility, nanoscaleability, memory capability, and conductance modulation, are essential properties affirming the reliability of the memristor in neuromorphic networks, especially as synaptic function. The memristor is studied in the coupling mode between two neuron cells where the synchronization phenomena is investigated numerically and theoretically (Ascoli et al. (2015), Zhang and Liao (2017), Xu et al. (2018), Bao et al. (2020)). Unidirectional coupling and bidirectional or mutual coupling are the commonly used coupling modes for nonlinear chaotic systems (Volos et al. (2015)). The synchronization and chaos between two neuron cells is also investigated by using unidirectional and bidirectional coupling (Zhang and Liao (2017)).
The main application of our work is to use the memristor as a synaptic link between neurons in electronic models, as for example, in hybrid technologies with neuronal electronic prosthesis between real neurons. The network is initially composed of a linear capacitor and a nonlinear resistance in each cell and a linear resistor in series (Comte et al. (2001)), see Figure 3A. The equivalent memristorbased network is shown in Figure 3B where the series resistance is replaced by a memristor. The memory circuit element being intrinsically asymmetric (Di Ventra et al. (2009)), the TiO 2 memristor crossbar is used to visualize the nature of current flowing through the device with respect to the polarities of the applied voltage. Using a memristor for the image processing technique was also reported by Prodromakis and Toumazou (2010), where a memristive grid is employed to perform edge detection. In a first step to implement 2D-memristor-based cellular nonlinear networks for signal and image processing purposes or for modeling a neural network with memristors as synapses, we rather focus here on the interaction of the memristor between pixel cells by considering a system of two cells in order to assess the behavior of the memristor quantitatively and qualitatively.
The underlining task is accompanied by observing the pinched hysteresis loop (PHL), memristance transition, voltage evolution of the cells, and the explanation of the memristor's lack of bilaterality. We derive analytically a second-order nonlinear differential equation characterizing the interaction of the memristor between the two cells bidirectionally. The system is studied in the phase plane allowing visualizing the memristor asymmetry. The memristance variation of a bipolar memristor with respect to the direction of flowing current affects its reliability in some potential applications where sensitivity in direction is important, for example, the memristive grid for neuromorphic application and image processing, hence the need for the so-called memristor fuse (Jiang and Shi (2009)). The formation is achieved by connecting two identical memristors antiserially. We describe the memristor fuse and obtain some results in accordance to our application.

DESCRIPTION OF THE POLARITY REVERSAL
Figures 4A,B show two identical memristors M 1 and M 2 connected in parallel across a voltage source V with terminals for a direct polarization for M 1 but reversed in the case of M 2 , both memristors having the same initial condition. The schematic is shown in a way to illustrate the trending of the mobile charge carriers under the influence of external bias. The currents through M 1 and M 2 are measured as I 1 and I 2 , respectively, and V(t) I 1 (t)M 1 I 2 (t)M 2 . Although the memristors are identical, we found that |I 1 | ≠ |I 2 |; hence, the conductivity differs if the polarity is reversed, even though with the same initial condition and voltage excitation. Hence, the device offers low resistance path with the orientation of M 1 and high resistance path for that of M 2 .
The schematic is shown in Figure 2A, where initially, the width of TiO 2 altogether is D and the width of the doped (TiO 2-e ) region is w and the undoped one is (D − w). Figure 4A shows that the positive charges in the doped region are repelled by the positive terminal of the power supply, thereby making the width of the doped region to expand such that w → D, as illustrated by the width trending w p . If the terminals of the applied voltage are reversed ( Figure 4B), the negative terminal from the power supply attracts the positive charges in the doped region, thereby causing the contraction of the doped region such that w → 0 with the width trending illustrated as w n . Figure 4C shows the comparison of the current flowing through M 1 and M 2 with respect to the applied voltage source. The result is obtained using a sine input voltage source and a window function by Joglekar and Wolf (2009). The current I 1 (t) and I 2 (t) have different absolute values. The current-voltage graph of I 2 (t) falls in the second and fourth quadrants due to the reversed terminals of V(t).
It follows that the conductivity of a memristor can be compared to that of a diode in terms of terminal polarity. However, unlike the diode, the memristor conducts electricity in both directions but the conductivity increases if its higher polarity terminal is connected to the positive terminal of the applied voltage source and decreases if its lower polarity terminal is connected to the positive terminal of the applied input voltage source.
In a nanoscale device, even small voltages can generate large electric field required to cause current to flow through the device. The smaller the device, the higher the electric field developed, and hence, more current flows through the device (Strukov et al. (2008)). The electroforming process forms the oxygen vacancies which cause a high-conducting channel (TiO 2-e ) shunting the bulk of the insulation film TiO 2 (Yang et al. (2008), Pickett et al. (2009)). Depending on the nature of the bipolar input source, the conducting channel is affected by the variation of the tunneling barrier width (D − w). Therefore, the device conductivity can be described by the Simmons tunneling barrier model (Simmons   (2009)). Respectively, the overall process affects the device conducting channel widths w p and w n , and it gives the trend on how w approaches 0 or D, (C) the currents I 1 and I 2 through M 1 and M 2 , respectively, and the corresponding I-V curves for sine input voltage with 1V amplitude. Reversing the polarity of the memristor also affects the absolute value of the flowing current. The current I 2 (t) falls in the second and fourth quadrants due to the reversed polarity of the input voltage V(t).
Frontiers in Communications and Networks | www.frontiersin.org July 2021 | Volume 2 | Article 647528 (1963a), Simmons (1963b), Simmons (1971)), devoted to the conduction of a material in a given medium. Depending on the magnitude of the input source, the boundary moves back and forth proportionally to the concentration of the dopant, thus setting the resistance of the memristor. Note that, in Figure 2, in reality, the boundary can never be outside of the interval [0, D] because there is always doped and undoped material present; in other words, the doped region can only expand or contract but never ceases to exist.

MEMRISTOR ASYMMETRY FROM THE CIRCUIT POINT OF VIEW
To vividly visualize the effect of memristor asymmetry, we consider two identical RC cells shown in Figure 5, which is the simplified setup of the system in the work of Isah et al. (2020b). The cells are labeled as cell-1 and cell-2 having potentials V 1 and V 2 , respectively, coupled together by a memristor M with its orientation as shown. Two tests are carried out which allow to observe the interaction of the memristor bidirectionally: In the former, the direction of I(t) is as shown in Figure 5; meanwhile, in the latter, the direction is reversed. The voltage across the memristor is for Cond-2. By taking into account the history of the memristor, we have where q 0 is the amount of charge already flowed through the device from its last usage, and thus, it becomes the initial charge at time t 0. Therefore, we consider the same memristor M with the same previous history, characterized by the initial charge q 0 .
In Cond-1, it is placed in one way as shown in Figure 5 and in Cond-2, on the opposite way. Figure 5 is simulated in SPICE using the memristor model by Biolek et al. (2009), which can easily be implemented experimentally. The setup is activated by closing the switches S 1 and S 2 simultaneously. For each Cond-1 and Cond-2, we considered the initial charge q 0 38μC; meanwhile, the low and high resistance limits of the memristor are 100Ω and 16KΩ, respectively (Strukov et al. (2008)). Given the initial conditions of the cells, that is, for Cond-1, V 1 0 1V, V 2 0 0V and for Cond-2, V 1 0 0V, V 2 0 1V, the result is shown in Figure 6. It is, however, important to note that the initial voltages can have any numerical values. Figure 6A shows the memristance transition for Cond-1 and Cond-2, illustrated, respectively, by the solid and dash curves. Figure 6B shows the corresponding time evolution of V 1 (t) and V 2 (t). The results show that Cond-1 and Cond-2 lead to quite different scenarios.
In Cond-1, the memristance decreases corresponding to the expansion of the doped region; meanwhile, for Cond-2, the memristance increases corresponding to the contraction of the doped region. In both Cond-1 and Cond-2, the memristance transition eventually flattens as the cells stabilize, that is, at a time when V 1 (t) V 2 (t). The memristance transition depends on the initial conditions of the cells. For example, Figure 7 shows the case where the initial conditions of the cells are changed with two different initial charges as q 0 1 46μC and q 02 31μC. Figure 7A and 7B, respectively, show the memristance transition and the corresponding voltage evolution of the cells. Furthermore, Figure 6B and 7B show that the orientation of the memristor according to Cond-1 or Cond-2 affects the time taken for the system to stabilize. Let V m (t) represent the voltage across memristor; hence, V m (t) V 1 (t) − V 2 (t) for Cond-1 and V m (t) V 2 (t) − V 1 (t) Cond-2. No current flows through the memristor when V 1 (t) V 2 (t) because the voltage across the memrisor is also zero even though V 1 (t) V 2 (t) ≠ 0, see Figure 8. The combined evolution of V 1 (t) and V 2 (t) eventually stabilizes to zero due to the resistive nature of the cells.
Furthermore, from Figure 5, one can deduce the following equations: where M(q) is defined in the work of Isah et al. (2020a) as follows: with M(q) R off if q(t) ≤ 0 and M(q) R on if q(t) ≥ q d . The dynamics of the memristor between the two cells can be expressed analytically. Eqs. 6-8 are simplified to give FIGURE 5 | Memristor asymmetry from the circuit point of view: two charged RC cells coupled together by a memristor. The circuit is invoked by the switches S 1 and S 2 . The cells are at different potentials so that the current i(t) will flow through the memristor. The test is performed for V 1 > V 2 and then V 1 < V 2 , namely, Cond-1 and Cond-2, respectively. For example V 10 1V, V 20 0V and then V 10 0V, V 20 1V. For Cond-2, the voltage across the memristor becomes [V 2 (t) − V 1 (t)].
Frontiers in Communications and Networks | www.frontiersin.org July 2021 | Volume 2 | Article 647528 where τ c RC is the time constant of the cell. Substituting Eq. 9 into 11 and taking the time derivative on the left-hand side give Eq. 12 can be expressed in a normalized form as where τ t/τ c is the normalized time, X q/q d is the normalized charge (Joglekar and Wolf (2009), Biolek et al. (2012)), Y dX/dτ is the first derivative of X with respect to τ, _ Y dY/dτ dX 2 /dτ 2 is the second derivative of X with respect to τ, and M M/δR is the normalized form of (10), thus rewritten as follows: with R R off δR . Substituting (14) into (13) and posing c 1 R δR + R 2 and c 2 R 2 , Eq. 13 is better studied when replaced by the set of equations allowing to study the system in the phase plane (X,Y) which facilitates the observation of cell interaction in relation to the memristive effect under different initial conditions. Furthermore, Eq. 13 requires a continuous first derivative of M with respect to X at X 0 and X 1, and it is achieved perfectly by using Eq. 14.
The history of a memristor marks its initial value in ohms, determined by the last amount of electric charge passed through it. However, it is reported in the work of Isah et al. (2020b) and Chua (2015)that the initial memristance is unknown.
FIGURE 7 | System evolution using two different initial charges q 01 46μC and q 02 31μC with V 10 1V, V 20 0V for Cond-1 (solid curves) and V 10 0V, V 20 1V for Cond-2 (dash curves). (A) Memristance transition showing the variation effect of the initial charge; (B) the corresponding evolution of V 1 (t) and V 2 (t).
Frontiers in Communications and Networks | www.frontiersin.org July 2021 | Volume 2 | Article 647528 It is important to note that the state X 0 is not fixed (in other words, it is unknown) and it strongly depends on the history of the device. Moreover, from Eq. 15, we get where H is a conservative expression for the equation set Eq. 15, only depending on the initial conditions X 0 , V 10 , V 20 , and Y 0 . Hence, X 0 q0 q d , being simply the normalized form of q 0 , and we obtained Y 0 from Eq. 9 as follows: Recall that Y dX dτ ; then, Eq. 16 can be expressed as and the analytical relation between the normalized time τ and the normalized charge X becomes Here, X 1 and X 2 are the real roots in the denominators of Eq. 18; meanwhile, the coefficientsα 1 ,α 2 ,α 3 ,α 4 ,β 1 ,β 2 , and τ 1 areα No current flows through the memristor when V 1 (t) V 2 (t), and the voltage across the memristor is also zero. V 1 (t) and V 2 (t) eventually decay to zero due to the resistive nature of the cells.
Frontiers in Communications and Networks | www.frontiersin.org July 2021 | Volume 2 | Article 647528 Furthermore, the equilibrium point of equation system Eq. 15 is met when dY dτ 0, that is, Y 0; in other words, V 1 (t) V 2 (t); then, from (15), X 3 − 3 2 X 2 + c 1 0, having at least one real root corresponding to the value of X(τ) at the equilibrium point Y 0. The singularity point of the system is where the derivative dY dτ does not exist, that is, dY dτ ∞; hence, we get from (15) X 3 − 3 2 X 2 + c 2 0, having at least one real root corresponding to the singular line, for any given c 2 . Depending on the initial conditions, Y evolves positively according to Cond-1 and negatively according to Cond-2. Recall that the value of X 0 is not known; however, we considered all the possible occurrences as shown in Figure 9. Each trajectory begins with the corresponding value of Y 0 . Therefore, depending on X 0 , we observed different evolution patterns in the phase portraits.
The phase portraits show the families of curves for different initial conditions. The results are obtained for R on 100Ω and R off 16KΩ. Note that the memristance is unchanged for X ≤ 0 and X ≥ 1 and is, respectively, given by R off and R on as depicted by the parallel evolution of the curves outside the interval [0, 1]. Therefore, different possibilities are considered that take into account the case where X 0 0 or one and beyond. The lack of symmetry is highly observable as the curves evolve from left to right for Y > 0 and then from right to left for Y < 0 according to Cond-1 and Cond-2, respectively. Furthermore, the time evolution of cells 1 and 2, respectively, can be obtained from Eqs. 6-9 as follows:

MEMRISTOR FUSE
The lack of bilaterality manifested in a bipolar memristor device is challenging in terms of its usage for certain applications, such as communication link in bidirectional applications (Comte et al. (2001)). As shown in Figure 6 and Figure 7, using a memristor to link two possible sources of information communicating together bidirectionally is not advisable owing to its resistance dependency on the amount and direction of flowing current. To convert it, a memristor fuse is proposed and then demonstrated in the work of Shi (2009), Gelencser et al. (2012), and Serb et al. (2016). It is basically formed by connecting two memristors antiserially in order to avoid the lack of bilaterality (Yildirim et al. (2018)). The memristor fuse is reported to be useful in a memristive grid network for CNN neighborhood connection and image processing (Pershin and Di Ventra (2011), Gelencser et al.
). In general (Fouda et al. (2013)), Figure 10 shows the four possible ways to form series connections of two memristors with respect to their polarities. As shown in Figure 10, cases 1 and 2 refer to a serial connection of two memristors and the memristive effect is retained for these branches. Meanwhile, cases 3 and 4 are identical in structure and refer to antiserial connection of two memristors, thus forming a memristor fuse. The memristive effect for branches in cases 3 and 4 could be suppressed (Joglekar and Wolf (2009)). However, case 3 is the commonly adopted formation of a memristor fuse (Gelencser et al. (2012)). Note that cases 1 and 2 resemble FIGURE 9 | Phase portraits showing the charge evolution from left to right for Y 0 > 0 and from right to left for Y 0 < 0 under different initial conditions. The lack of symmetry is noticeable within the bulk of the device.
FIGURE 10 | Four possible series connections of two memristors with respect to the input source. The memristive effect is retained for cases 1 and 2 whereas it is balanced for cases 3 and 4 (Joglekar and Wolf (2009)). Although cases 3 and 4 are identical and both form a memristor fuse, only case 3 is commonly considered as memristor fuse formation (Gelencser et al. (2012)).
Frontiers in Communications and Networks | www.frontiersin.org July 2021 | Volume 2 | Article 647528 Figures 4A,B, respectively, with the exception that two memristors are involved. It is important to note that the resistance of the memristor fuse is the sum of the resistance of each of the individual memristor because the equivalent memristance is additive in a serially connected memristor. In addition, this could be a disadvantage to the desired amount of current and it also affects the dynamic features of the memristor to the extent that the current-voltage graph is merely linear; hence, the formation resembles normal resistor. Therefore, the resistance limits of the memristor fuse must be the same as those of memristor if acting alone.
The pinched hysteresis loop is one of the most distinguished fingerprints of a memristor (Adhikari et al. (2013), Chua (2014)) and is a reflection of its memory effect. As pointed out in the work of Chua (1971) and Chua (2015), without the memory, a memristor is nothing different from a resistor. A verification test is performed to compare the memristor fuse with a standalone memristor, as demonstrated in Figure 11. M p and M n are set with their polarity reversed in parallel with a memristor fuse M f , all connected across the same voltage source V(t). Note that the orientation of M p and M n is, respectively, similar to M 1 and M 2 illustrated in Figures 4A,B. Therefore, for the same resistance limits, the memristance of a FIGURE 12 | Schematic of a memristor fuse formed by two TiO 2 memristors M 1 and M 2 . The positive spot signifies the positive oxygen vacancies, which serve as the higher conducting part of each memristors. w 1 and w 2 represent the instantaneous width of the doped region, meanwhile the black and red arrows describe the trending expansion and contraction of the doped region, respectively. Hence, for any input voltage, say V(t), an increase in the width w 1 corresponds to a decrease in w 2 and vice versa.
Frontiers in Communications and Networks | www.frontiersin.org July 2021 | Volume 2 | Article 647528 memristor fuse is higher than the memristance of a standalone memristor, as can be observed from the respective slopes of Figure 11. Figure 12 shows the schematic of a memristor fuse formed by antiseries connection of two memristors M 1 and M 2 with the instantaneous dopant width denoted by w 1 and w 2 , respectively. During positive bias, w 1 expands while w 2 contracts, and during negative bias, w 1 contracts while w 2 expands. As w 1 tends to D, w 2 tends to 0 and the converse gives the opposite. Moreover, w 1 and w 2 could be represented in normalized forms as x 1 and x 2 , where x 1 w 1 D and x 2 w 2 D . Then, the drift speeds of the respective dopant are expressed as The same current flows through a series connection of two memristors; hence, i(t) is the same for both M 1 and M 2 . Here, we consider the expression of the memristance given by Eq. 1b rather than the one given in Eq. 10 because it is simpler and already investigated in the work of Fouda et al. (2013). From Eqs. 1b, 23, the rates of change of the instantaneous memristance of M 1 and M 2 are, respectively, obtained to be where 9 1 − δR q d 1 and 9 2 δR q d 2 . Although the current i(t) flowing through them is the same, however, the rate of change of memristance for one differs from the other. Therefore, from 24a and 24b, the rate of change of M 1 with respect to that of M 2 is given by where 9 − q d 2 q d 1 is called the mismatch factor (Fouda et al. (2013)) describing the increase of M 1 with respect to the decrease of M 2 and vice versa. Note that q d1 ≠ q d2 , each possibly able to depend on the dimension (D), the mobility of charge carriers, and the value of the lowest resistance (R on ) for M 1 and M 2 , respectively. Thus, FIGURE 13 | Comparing the system evolution using the circuit of Figure 5 by considering memristor (M) and then memristor fuse (M f ). Mc1 and Mc2 correspond, respectively, to a single memristor used according to Cond-1 and Cond-2, while Mf1 and Mf2 correspond, respectively, to a memristor fuse used according to Cond-1 and Cond-2. The results show the response for each case until the system stabilized: (A) the voltage evolution as V 1 (t) and V 2 (t), (B) the flowing currents through the memristor (I e1 and I e2 for Cond-1 and Cond-2, respectively) and memristor fuse (I f1 and I f2 for Cond-1 and Cond-2, respectively), and (C) memristance transition. Given the same initial condition, the results show that the memristor fuse conducts equally in both directions. For Mc1, Mc2, Mf1, and Mf2, c1, c2, f1, and f2 are subscripts denoting Cond-1 and Cond-2, accordingly.
Frontiers in Communications and Networks | www.frontiersin.org July 2021 | Volume 2 | Article 647528 the mismatch factor is determined by the mobilities and velocities of charge carriers in both memristors. Given the initial memristance of M 1 and M 2 as M 1 0 and M 2 0 , respectively, integrating Eq. 25 gives where δM 0 M 10 − 9M 20 . The net memristance is additive in a series connection of memristors. The instantaneous memristance of the memristor fuse M f (t) is given by M f (t) M 1 (t) + M 2 (t); thus, The instantaneous memristance of the memristor fuse depends then on the mismatch factor (9). Figure 13 shows the comparison of the circuit response for the memristor and memristor fuse using the setup shown in Figure 5. Then, for Cond-1, V 1 0 1V and V 2 0 0V, while for Cond-2, V 10 0V and V 20 1V with q 0 45.58μC, C 1μF, and R 100KΩ in each case. Furthermore, M e1 and M e2 represent memristor according to Cond-1 and Cond-2, respectively. Similarly, M f1 and M f2 represent the memristor fuse according to Cond-1 and Cond-2, respectively. For each case, the system evolves and eventually stabilizes when V 1 (t) V 2 (t). Figure 13A shows the evolution of V 1 (t) and V 2 (t) for M e1 , M e2 , M f1 , and M f2 . The results of M e1 and M e2 show a shift difference during the transient state while there is no such shift between the curves of M f1 and M f2 , showing that the memristor fuse behaves equally in both Cond-1 and Cond-2. Figure 13B shows the currents through the memristor as i e1 and i e2 according to Cond-1 and Cond-2, respectively, and then through the memristor fuse as i f1 and i f2 according to Cond-1 and Cond-2, respectively. Furthermore, the results show that no current is flowing through the memristor when V 1 (t) V 2 (t) as similarly observed in the analytical result of Figure 9, that is, Y 0 when V 1 (t) V 2 (t). Figure 13C shows the flowing current through the memristor and the memristor fuse. The results show the differences in the memristor responses according to Cond-1 and Cond-2, but the memristor fuse behaves indifferently in both conditions.

CONCLUSION
We introduced the application of memristor in a nonlinear network, focusing specifically on the behavior of the memristor with respect to the polarity reversal effect of the input signal. Our target is the implementation of memristorbased 2D nonlinear networks for versatile applications, such as signal processing and electronic prosthesis for a synaptic link between real neurons. Here, we investigate the interaction of a bipolar memristor between two RC cells communicating together bidirectionally. In this way, the interaction of the memristor within the network is studied qualitatively and quantitatively. We have shown from the circuit point of view and the analytical solution that the conductivity of the memristor depends on the polarity of the applied input signal, thus affecting the mobility of its charge carriers, this property being due to the intrinsic nature of the device. It is an inevitable nature of a bipolar memristor, irrespective of its device technology. Hence, the memristive effect changes according to the connection mode and the amount of current flowing through it, showing that the memristor is not a bilateral circuit element as verified by our study.
To achieve the memristive effect with symmetry, a memristor fuse is proposed. We present the detail analytical interpretation of the memristor fuse. We also authenticate the memristor fuse prior to apply it in a circuit, and the results show that the memristor fuse behaves like a standalone memristor under high input frequency. Although connecting two memristors antiserially to form a memristor fuse lets the dynamic of the two state variables system become more intricate, as well as the dynamics of the resistance switching (Serb et al. (2016)), terminals' asymmetry is resolved as confirmed by the results shown in Figure 13. The symmetry displayed by the memristor fuse suggests it to be a promising element useful as a memristive grid in neighborhood connections and it could become an important concept for the ongoing study of our memristorbased network.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.