Memristor Spice Model For Designing Memristor Circuits Economics Essay

Memristor Spice form For Designing Memristor Circuits economic science EssayMemristors atomic number 18 novel electronic thingumajigs, a device that bath be used and slip away a extensive advantage in many applications such as memory, logic, neuromorphic governing bodys and so on. A computer homunculus of the memristor would be a expedient son of a bitch to analysis band carriage to helper in develops application of this memristor as peaceful travel factor via computer pretense. In this paper, we incorporate a memristor spiciness for headinging memristor circuit which is more(prenominal) focusing on non- linear archetype and analog circuits. gaminess vex would be appropriate way to tell a bea real device operation. We incorporating the memristor with various windowpanepanepane agencys that have been proposed in non linear ion vomit up memristor devices. In investigating and characterizing the somatic electronic and behavioral properties of memristo r devices, the circuit analysis of the proposed memristor regulates argon wherefore been studied. The assumption go forthput should have a current- potential difference hysteresis curve, which looks like bow tie. The coils map the change by reversal behavior of the device. Then, we come out with a simple-minded analog circuit which in this case we construct a simple integrator op-amp and discriminator op-amp circuit and make comparison surrounded by memristor implemented circuit and design circuit. The explore verifies the proposed memristor illustration, the possibilities of implementing memristor specimen and the advantage implementing the memristor in analog circuit.Keywords memristor, alter model, non linear, window personas, analog circuit.IntroductionMemristor is the contraction of memory resistance which is a nonoperational device that provides a determinational apprisal amongst perpetration and liquidize. It is a both- net circuit constituent in whic h the go in the midst of the 2 remnants is a take to the woods of the amount of electric charge that has passed through the device 1. A memristor is said to be charge- misrepresentled if the relation betwixt flux and charge is explicit as a piece of electric charge and it is said to be flux-controlled if the relation between flux and charge is expressed as a shape of the flux gene linkage 2.In 1971, Leon Chua proposed that at that place should be a fourth total passive circuit element to create a mathematical blood between electric charge and magnetic flux which he called the memristor which is short for memory resistor 2.The current is defined as the time derivative instrument of the charge. The potential difference is defined as the time derivative of the flux harmonize the faraday law. A resistor is defined by the kinship between potential difference and current dv=Rdi, the capacity is defined by the family between charge and potency dq=Cdv, the inductor is def ined by the kin between flux and current d=Ldi. The fourth fundamental circuit element completes the symmetry of the relation between charge and magnetic flux d=Mdq. tabularise 1 destine the relationship between the fundamental circuit element.Basic two terminal devicesEquationRelationship between fundamental circuit elementResistor,Rdv=Rdiv and iCapacitor,Cdq=Cdvv and qInductor,Ld=Ldii and Memristor,Md=Mdqq and Table 1 The four fundamental element (resistor, capacitor, inductor and memristor).In 2008, Stanley Williams and his team at Hewlett Packard had succesfully fabricated the first memristor in physical device form which is a long wait from Leon Chua discovery in 1971 3. Memristance is a property of memristor. When the charge period of times in one perplexity through a circuit, the resistances of the memristor step-up. The resistance decreases when the charge flows in the opposite circumspection in the circuit. If the applied potential difference is morose off, thus s drop deadping the flow of charge and the memristor remembers the last resistance that it had 1.In HP memristor model, to create a memristor, they used a very quash convey of titanium dioxide (TiO2). The thin film is wherefore sandwiched between the two platinum (Pt) contacts. One grimace of TiO2 is doped with oxygen va tummycies denoted as TiO2-x which x is usually 0.05. The oxygen va rumpcies are positively charged ion and make it conductive, thus it behaves as a semiconductor. other side of the TiO2 junction is undoped. The undoped region has insulating properties. The device established by HP is appearancen in jut out 1 3. approach pattern 1 Memristor model oblige from 3.When a positive emf is applied, the positively charged oxygen va buttockscies in the doped TiO2-x layer are repelled and moving them towards to the undoped TiO2 layer. When the boundary between the two materials moves, the percentage of the conducting TiO2-x layer is increase. Thus, the conductivity of t he whole device increases.When a prohibit voltage is applied, the positively charged oxygen vacancies are attracted and pulling them out of TiO2 layer. This increases the amount of insulating TiO2, thus increasing the resistivity of the whole device. When the voltage is turned off, the oxygen vacancies do not move. The boundary between the two titanium dioxide layers is frozen. This is how the memristor remembers the voltage last applied 1.MethodologyOur aim in this research is to provide a simulation program adequately simulates and can be used as a circuit element in design snip. To model the electrical indications of the memristor, SPICE would be appropriate way to describe real device operation 4. Moreover, using the model as a sub-circuit can highly guarantee a reasonable high flexibility and scalability features 5. We use LTSPICE to create a memristor model and design brand-new symbolism of the memristor circuit for the simulation because LTSPICE is much easier to handl e compared to others. On the other hand, LTSPICE is a freeware and it impart spread a great advantage to the students in doing research for this newly devices. We use SPICE model that been adapt from 6 and we made near adjustment so we can use it for several(prenominal) window functions that has been proposed for non linear ion gasconade model.The SPICE model is created based on the mathematical model of the HP Labs memristor. After the memristor has been modeled, we first studied the difference between proposed memristor and then we go out start design and implement the memristor with an analog circuit. We also wonder and made a comparison between the memristor circuit with analog circuit to assure the difference and study the behavior of the circuits.Model Of The Memristor from HP LabsIn the model of a memristor presented here, there is a thin semiconductor film that has two regions, one with a high concentration of dopant that behaves like a low resistance called RON and the other with a low dopant concentration with high(prenominal) resistance called ROFF 3. The film is sandwiched between two metal contacts as in foreshadow 1.The total resistance of the memristor, RMEM, is a sum of the resistances of the doped and undoped regions, w is the width of the doped region and D is the total length of the TiO2 layer. ROFF and RON will be the limit determine of the memristor resistance for w=0 and w=D. The ratio of the two resistances is usually given as 102 103.(1)(2)From the ohms law relation between the memristor voltages and current, we get(3)Then, we insert (1) into (3). The voltage v(t) across the device will move the boundary between the two regions causing the charged dopants to drift. So, there is a drift ion mobility v in the device. The change of the boundary is denoted as in (5).(4)(5)To get x(t), we then integrates the right side of equality (5) which then yields the following formula(6)By inserting comparison (6) into equation (4) and si nce usually RON(7)Where v is the average drift velocity and has the units cm2/sV, D is the thickness of titanium-dioxide film ROFF and RON are on- claim and off- nominate resistances and q(t) is the total charge passing through the memristor device.Non Linear Ion Drift ModelEven a small voltage across the nanodevices will modernise a large electric field 7. This causing the ion boundary built in bed will move in a decidedly non-linear. Nonlinear dopant drift adds nonlinear window function f(x) to the earth equation. The window function decreases as the state variables drift speed approaches the boundaries until it slayes zero when reaching either boundaries 8. The speeds of the faecal matter of the boundary between the doped and undoped regions are depending on several factors.(8)Where v is the dopant mobility. The speed of the boundary between the doped and undoped regions decreases gradually to zero at the film edges 1. We simulate the nonlinear ion drift memristor model with these window function to see the difference and the issue that been faced by them.Window passWindow function is a function of the state variable. Window function forces the bounds of the device and to add nonlinear behavior close to these bounds. In other words, it creates the boundary for the memristor. Any topicive window function should and so fulfill the following conditions 8Take into account the boundary conditions at the top and bottom electrodes of the deviceBe capable of imposing nonlinear drift over the entire active core of the deviceProvide linkage between the linear and nonlinear dopant drift modelsBe scalable, meaning a range of fmax(x) can be obtained such that 0 fmax(x) 1lend oneself a built-in control statement for adjusting the model.There are several window functions that have been proposed for non-linear model till date which are by Strukov, Joglekar and, Biolek, and Prodromakis. Strukov proposed the following window function 3.(9)However, as we can see in the figure 2, this window function lacks of flexibility. take in 2 Plot of Strukov window function.Another window function was proposed by Joglekar 4, which has a control line of reasoning p which is a positive integer. The purpose of having a control debate as an exponent is to incorporate scalability and flexibility in window function f(x) that describes the dopant kinetics.(10) excogitation 3 displays a graphical theatrical of the window function described by Joglekar for various p parameter (p=1, 5 and 10). This control parameter controls the linearity of the model, where it becomes more linear as p increases. This window function ensures zero drift at the boundaries. From the plotted graph, we noticed that the upper limit f(x) time prize is occurs at the center of the device and zero is obtained at two boundaries. However, a significant liability of this model lies in the fact that if w hits any of the boundaries (w = 0 or w = D) the state of the device cannot be further a djusted. This will be from now on termed as the terminal state problem.Figure 3 Plot of Joglekar window function for p=1, 5 and 10.Then, Biolek proposed another window function that allows the memristor to come back from the terminal state problem.(11)The reversed stoop is now should move back the state variable after it reaches either boundary. This feature is described by a current dependent step function, stp(i), which is a part of a new window function f(x) that behaves differently in each voltage bias direction.(12)Figure 4 Plot of Biolek window function for p=1, 5 and 10.Figure 4 displays a graphical representation of the window function described by Biolek for various p parameter (p=1, 5 and 10). When x starts at 0, we noticed that the function equal to 1. As x increase approaching D, the function approaches 0. Once the current reverse the direction, the function flat switch to 1. As x decrease back to 0, the function also decreases to 0. Biolek window function eliminates c onvergence issues at the devices boundaries.The last window function for non-linear model is proposed by Prodromakis 8.(13)Figure 5 displays a graphical representation of the window function described by Prodromakis for various p parameter (p=1, 5 and 10). As we can see, it allows the window function to scale upwards which implies that fmax(x) can take any value within 0 Figure 5 Plot of Prodromakis window function for p=1, 5 and 10.SPICE Model of MemristorFigure 6 Stucture of the SPICE model from 6.In the to a higher place circuit in figure 6, VMEM is the foreplay voltage and Imem is modeled to be the current through the memristor. The flux is calculated by integrating the voltage VMEM and the charge is calculated by integrating the current IMEM.Figure 7 Resistive port of the memristor model.As we can see in figure 7, the circuit is actually referred to total resistor RMEM. RMEM (x) = ROFF -xR where R= ROFF-RON. ROFF is the resistor in series voltage source whose terminal voltag e is controlled by the formula -xR.Figure 8 variousial equation modeling of the memristor.Figure 8 shows the differential equation modeling of the memrsitor. It consist a part of the voltage controlled source xR and the differential equation from equation (6) which serves as an integrator of the quantities on the right side of the state equation (6) which is to get the value of normalize x. EMEM is the voltage source whose terminal voltage is controlled according to the formula -xR. GX is a current source whose current is controlled according to the equation IMEMf(V(x)) where V(x) is the voltage across the capacitor Cx and it models the normalized width x of the doped layer. F(V(x)) is the window function, k is vRON/D2 and x0 is the sign voltage of the capacitor.6.The relation between memristor current and voltage is modeled as on the basis of RMEM (x) = ROFF -xR where R= ROFF-RON. The voltage V(x) across the capacitor CX models the normalized width x of the doped layer. The initi al state of x is modeled by the initial voltage of the capacitor. The flux is calculated by the time- built-in of voltage, and the charge is calculated by the time-integral of current.Result and DiscussionAll models were simulated in LTSpice using SPICE model that was given in 6, we add new nonlinear window functions that was proposed by prodromakis and strukov to the model and compare all suggested window functions.memristorFigure 9 Memristor circuit.Figure 9 shows configuration of single memristor for measuring the behavior of memristor model in LTSPICE with a sine wave input voltage of 1.2V with 1Hz frequency. The values for the memristor parameters v, D, RON, ROFF and RINITIAL are 10-10cm2s-1V-1, 10 nm, 100ohm, 16kohm and 11Kohm. All model are using akin window function parameter p=10.CUsersFadzDesktop consequencestrukovstrukov1.pngFigure 10 Strukov memristor model voltage, IMEM, RMEM and normalized x.Figure 10 shows the simulation result of memristor SPICE model for Strukov w indow function of voltage, IMEM, RMEM and normalized x. As we can see, the current of the memristor, IMEM is alter up to most 100A for maximum of 1.2V voltage applied. The RMEM for this model show that the values are in range of 11kOhm till 12kohm which means the effect of the voltage applied to the memristor wholly give slightly changes on the value of the memristor. Noticed that when positive voltage is applied, the conductivity of the device increases thus the memristance is decrease. When contradict voltage is applied, the resistivity of the device increase thus the memristance is also increase. This verifies the memristive system on the device. In normalize x graph, we also noticed that the normalized x is at higher state in the beginning. Figure 11 shows the I-V characteristic of the devices and the relationship between charge and flux. The charge and flux curve curves confirms the considerably known fact that there is a one-to-one correspondence between them in spite of the 1-4 hysteresis effect. Strukov memristor shows lack of flexibility of controlling the device.CUsersFadzDesktopresultstrukovstrukov4.pngFigure 11 Strukov memristor model I-V loop hysteresis and relationship of charge and flux.CUsersFadzDesktopresultjoglekarjoglekar1.pngFigure 12 Joglekar memristor model voltage, IMEM, RMEM and normalized x.Joglekar window function calculates to be promising as the existence of controlling parameter. Figure 12 shows the simulation result of memristor SPICE model for Joglekar window function of voltage, IMEM, RMEM and normalized x. Same as strukov memristor, when positive voltage is applied, the conductivity of the device increases thus the memristance is decrease. When negative voltage is applied, the resistivity of the device increase thus the memristance is also increase. The current of the memristor, IMEM is transfering up to approximately 300A for maximum of 1.2V voltage applied. Joglekar window function give higher current compared to ot hers. It shows that the current in the memristor are much easier to move. The RMEM are within range of nearly 0ohm to 11kohm which give full range of value for the memristor. Figure 13 show the I-V hysteresis loop of the devices and the relationship between charge and flux. The switching behavior is much more sensitive on the voltage level than Strukov window function. But, in term of stability, Joglekar window function cannot execute for an arbitrary length of time. This failure is caused by the convergence issue where when the memristor reach w=0 or w=D, the state of the device cannot be further adjusted.CUsersFadzDesktopresultjoglekarjoglekar4.pngFigure 13 Joglekar memristor model I-V loop hysteresis and relationship of charge and flux.Biolek window functions are supposed to solve terminal state problem as in literature 4. It should solve the boundry problem of the terminal state. Figure 14 and 15 shows the simulation result of memristor SPICE model for Biolek window function. The current of the memristor, IMEM is varying up to approximately 220A for maximum of 1.2V voltage applied. The RMEM are within range of nearly 1kohm to 11kohm. Figure 15 show the I-V hysteresis loop of the devices and the relationship between charge and flux. From the figures, we observe that the biolek memristor spare the highly non-linear device characteristic behavior. In addition, Bioleks model allows for general asymmetric I-V device behavior modeling.CUsersFadzDesktopresultbiolekbiolek1.pngFigure 14 Biolek memristor model voltage, IMEM, RMEM and normalized x.CUsersFadzDesktopresultbiolekbiolek4.pngFigure 15 Biolek memristor model I-V loop hysteresis and relationship of charge and flux.Prodromakis window functions are also said to solve the boundry issue. Figure 16 shows the simulation result of memristor SPICE model for Joglekar window function of voltage, IMEM, RMEM and normalized x. The current of the memristor, IMEM is varying up to approximately nearly 180A for maximu m of 1.2V voltage applied. The RMEM are within range of nearly 3kohm to 11kohm.CUsersFadzDesktopresultprodromakisprodromakis1.pngFigure 16 Prodromakis memristor model voltage, IMEM, RMEM and normalized x.Figure 17 show the I-V hysteresis loop of the devices and the relationship between charge and flux. The hysteresis loop is shown to be asymmetrical slice the OFF state of the device is highly non-linear compared with other.CUsersFadzDesktopresultprodromakisprodromakis4.pngFigure 17 Prodromakis memristor model I-V loop hysteresis and relationship of charge and flux.In comparing of I-V characteristic hysteresis loop, as we can see in figure 18, it shows all hysteresis loops for all proposed window functions. By using corresponding parameter we can see the difference in each model. Joglekar window function seems to have a strong memristance compared to others.All models seem to be a good approximation of the measurement of the real memristor produces by HP Labs. But, Prodromakis me mristor model satisfies all the prerequisites and improves on the shortcomings of existing models.CUsersFadzDesktopresultall iv loop hysteresisiv hysteresis loop all.pngFigure 18 I-V Hysteresis Loop for all models.CUsersFadzDesktopresultdiff parameterprodromakis p=1, p=5, p=10.pngFigure 19 Prodromakis I-V Hysteresis Loop when p=1, p=5 and p=10.Then, we try change the parameter of p of the model. In this case, we use prodromakis memristor model and change the value of integer p=1, p=5 and p=10 to see the difference. As we can see in figure 19, as the value of p is increase, the hysteresis is shrinking. similar with changing the parameter p, figure 20 also confirm that the hysteresis shrinks at higher frequencies.CUsersFadzDesktopresultdiff frequency1hz-5hz.pngFigure 20 Prodromakis I-V Hysteresis Loop when frequency f=1Hz, f=2Hz and f=5Hz.In term of exponent dissipation, as we simulates our results. We can get the value of maximum IMEM for each model. We can calculate the indicant by using P=IV equation. Table 2 show the maximum index number dissipation for each memristor.MemristormodelMax IMEM, APower, WStrukov100A120WJoglekar300 A360WBiolek220A264WProdromakis180 A216WTable 2 IMEM and Power dissipation for all at the memristor.As in table 2, we can see that the Strukov model give lowest power which is 120W while Joglekar model give much higher power dissipation which is about 360W compared to the others. We also noticed that as the memristor model is improves, the power become lesser. Prodromakis give quite good power dissipation which is 216W as the best windows function and model as till now.In implementing memristor with analog circtuit, we pick a two simple analog circuit to be tested. Figure 21 shows the SPICE topology of the memristor based integrator op amplifier with the input voltage Vp-p=2.4V from -1.2V to 1.2V and C1=25F. Using the memristor model that we create earlier with various types of window functions, we see the difference on the simulat ed result on each model. The values for the memristor parameters are same for all model with uv, D, RON, ROFF and RINITIAL are 10-10cm2s-1V-1, 10nm, 100ohm, 16kohm and 11Kohm.memristor integrator op ampFigure 21 Memristor Implemented Integrator CircuitCUsersFadzDesktopresultintegrator+ve integrator .pngFigure 22 Positive Integrator simulationCUsersFadzDesktopresultintegrator-ve integrator .pngFigure 23 Negative Integrator simulationFigure 22 and 23 shows the simulation result for the implemented memristor integrator op amplifier for positive input and negative input respectively. In this case, we implement prodromakis memristor to the integrator circuit. As we know, the integrator acts like a storage element that produces avoltageoutput which is proportional to the integral of its input voltage with respect to time. The magnitude of the output signal is unconquerable by the length of time a voltage is present at its input as thecurrentthrough thefeedback loopcharges or set frees t hecapacitoras the requirednegative feedbackoccurs through the capacitor. For positive starting time input, we vary the voltage from 1.2V down to -1.2V and going back to 1.2V over time. When positive voltage are applied in the beginning, the output voltage tend to discharge and drop from 0V to negative voltage and charging back to 0V when the input voltage are drop to negative voltage. The output for negative voltage applied from starting point give a vice versa result. The charging and discharge are depends the voltage applied over time and the value of the capacitor. We can say the memristor models give quite good result for an integrator.Figure 24 Memristor Implemented discriminator Circuit.CUsersFadzDesktopresultdifferentiatordifferentiator.pngFigure 25 Saw tooth input of Memristor implemented differentiator op-amp simulation.Figure 25 shows the simulation result for the implemented memristor differentiator op amplifier. We are using prodromakis memristor model for the memristo r carrying out. As we know, for differentiator op-amp, the magnitude of its output is obstinate by the rate at which the voltage is applied to its input changes. The swift the input voltage changes, the greater the output voltage becomes. If a cut tooth input signal is applied to the input of the differentiator op-amp a substantial wave signal will be produced. As we can see in figure 25, the simulation shows a quite good result for differentiator. We also noticed some spikes at the output voltage. Each spikes occurs only occurs the brief moment the saw tooth is changing from one level to the next. The voltage spikes represent a temporary output voltage.ConclusionAs a conclusion to this research is that it could bring a new light of familiarization in the integration of memristive components in any kinds of electronic devices that are at nanoscale. It is useful to have a computer model of the memristor as a tool for the analysis of the behavior of the circuits in developing app lication of this memristor as passive circuit element via simulation. SPICE model will definitely help us to conduct interesting simulation experiments and can be of great importance for such a research in future while the memristor are still hard to fabricate to study the behavior of the circuit. Different models with strong behavior and reason give a lot of benefits in development purpose to create the possibilities of the implementation in an integrated circuit. The possibilities for implementation of the memristor with analog circuit are wide open.Appendix.SUBCKT memristor plus deduction PARAMS+ Ron=100 Roff=16K Rinit=11K D=10N uv=10F p=10************************************************ DIFFERENTIAL equality MODELING ************************************************Gx 0 x value= I(Emem)*uv*Ron/D**2*f(V(x),p)Cx x 0 1 IC=(Roff-Rinit)/(Roff-Ron)Raux x 0 1T* RESISTIVE PORT OF THE MEMRISTOR ************************************************Emem plus aux value=-I(Emem)*V(x)*(Roff-Ron)R off aux minus Roff************************************************Flux deliberation************************************************Eflux flux 0 value=SDT(V(plus,minus))************************************************Charge computation************************************************Echarge charge 0 value=SDT(I(Emem))************************************************ WINDOW FUNCTIONS* FOR NONLINEAR DRIFT MODELING *************************************************proposed by joglekar.func f(x,p)=1-(2*x-1)**(2*p)*proposed by biolek.func f(x,i)=1-(x-stp(-i))**(2*p)*proposed by prodromakis.func f(x,p)=1-(((x-0.5)**2)+0.75)**p*proposed by strukov.func f(x,p)=x-x*2.ENDS memristorAcknowledgmentThis paper participates in the IEEE Student Conference Research Development SCORED 2012. The reason would like to thank Dr. Wan Fazlida Hanim bte Abdullah for being supervisor in this final year project. The author also would like to thank Universiti Teknologi MARA for funding the research work through the Excellence Fund Grant 600-RMI/ST/DANA 5/3/RIF(360/2012).