CONTENTS
21.1 Introduction
21.2 Describing Function
21.2.1 Memristor
21.2.2 Nanodevices Based on DNA or CNTs
21.3 The Liénard Equation
21.3.1 Conditions for a Limit Cycle
21.3.2 Quadratic Memristor
21.3.3 Nonlinear Nanodevice
21.4 Simulations
21.4.1 Limit Cycle
21.4.2 Stable Fixed Point
21.5 Conclusions
References
Nanotechnology is a promising research area with applications in different disciplines. In electrical and electronic engineering, nanodevices are new building blocks for circuits that increase the designer possibilities beyond traditional elements. In this chapter, two methods to study the behavior of nonlinear circuits and systems are illustrated for memristors, and nanodevices based on deoxyribonucleic acid (DNA) or carbon nanotubes (CNTs).
The two points of view are: (i) frequency analysis with the describing function (DF) [1,2], and (ii) time analysis with the Liénard equation [3]. The purpose of these tools is to predict the behavior of circuits and systems with nonlinear nanodevices or even to guide the design of nanodevices to synthesize any desired dynamics. Also, the tunable nanodevices could increase the complexity and applications of new circuits.
In the literature, it is common to find experimental v−i characteristics for nanodevices; these plots show nonlinear behavior that can be described by piecewise linear equations. Some typical nanodevices are based on oxide films [4, 5 and 6], DNA [7], or CNTs [8]. In other nanodevices, the nonlinear v−i relationship can be expressed as formulas, v = φ(i), ϕ = f(q); the current is related to the charge, i = dq/dt, and the voltage is related to the flux, v = dϕ/dt.
This chapter is divided as follows. Section 21.2 presents the definition of DF and its calculation, step by step, for a general nonlinearity. Then, the general formula is applied to memristors, and DNA or CNTs-based nanodevices. Section 21.3 introduces the Liénard equation for series RLC circuits with nanodevices; two types of nonlinearities are considered v = φ(i) and ϕ = f(q); also, five conditions for the existence of limit cycles are formulated. Section 21.4 presents the simulations to support the theoretical predictions for self-sustained oscillations. Finally, Section 21.5 includes some conclusions and future work.
In nonlinear control theory, the DF is used to predict the behavior of a nonlinear device in closed loop with a linear system, see Figure 21.1. The characteristic polynomial (denominator) of the closed-loop transfer function produces two equations; the solution yields the frequency and amplitude (ω,a) of the limit cycle or self-sustained oscillation [1]. This approach is called analysis in the frequency domain; the input voltage is periodic with known amplitude.
The closed-loop transfer function obeys the expression,
(21.1) |
where Ψ(a) is the DF of the nonlinear nanodevice and G(jω) is the linear system transfer function in the frequency domain. The closed loop has a limit cycle when,
(21.2) |
The function G(jω) is a complex quantity with real and imaginary parts; if the DF is real, there are two equations,
(21.3) |
(21.4) |
Solving Equation 21.4 produces the limit cycle frequency and Equation 21.3 yields the oscillation amplitude [2]. When a closed-loop system like Figure 21.1 is given, in a practical situation, the first step is to find the DF Ψ(a) for the nonlinearity.
The definition of a real DF Ψ(a) for odd, time invariant, memoryless nonlinearities is the following [1]:
(21.5) |
This analysis is valid for the input voltage,
(21.6) |
Figure 21.2 shows a voltage–current characteristic for a general nanodevice where (k0,k1,k2) are slopes; this v−i function can be studied using the DF as follows.
To evaluate integral (21.5), it is necessary to find the current i(θ) for different intervals following the input voltage (21.6),
(21.7a) |
where
(21.7b) |
Writing again Equation 21.5 and using Equation 21.7 yields,
(21.8) |
Integrating each term [9], and after some algebra, the DF for the nonlinearity of Figure 21.2 obeys,
(21.9) |
Formula (21.9) can be used to find the DF of known nanodevices by replacing the value δ0, δ1, k0, k1, k2 from the corresponding v−i characteristics.
Figure 21.3 shows the v−i response for a memristor [4, 5 and 6].
To reduce Figure 21.2 to Figure 21.3,
Replacing these relationships in Equation 21.9,
(21.10) |
Equation 21.10 is the memristor DF when the frequency of the input voltage is constant.
21.2.2 NANODEVICES BASED ON DNA OR CNTS
Other types of nanodevices use the molecule of life (DNA) or CNTs between a pair of electrodes [7,8]. For these electronic devices, the general DF (21.9) can be applied to the v−i characteristic, see Figure 21.4.
To reduce Figure 21.2 to Figure 21.4,
Replacing these parameters in Equation 21.9,
(21.11) |
Equation 21.11 is the DF for the typical nonlinearity of DNA or CNT-based nanodevices.
The second approach to study circuits and systems with nonlinear nanodevices is known as time domain analysis and involves nonlinear differential equations; here, numerical solutions are required when the differential equation is not standard.
The forced Liénard equation [3] has the structure,
(21.12) |
where a1(x); a0(x) are nonlinear functions; u is the forcing input.
21.3.1 CONDITIONS FOR A LIMIT CYCLE
The solution of Equation 21.12 presents a limit cycle or self-sustained oscillation if the following five conditions are satisfied,
1. The coefficients a1(x); a0(x) are continuously differentiable for all x.
2. Coefficient a0(x) has odd symmetry,
(21.13) |
3. Coefficient a0(x) is positive for x positive,
(21.14) |
4. Coefficient a1(x) has even symmetry,
(21.15) |
5. The function F(x) has one positive zero at x = a,
(21.16) |
F(x) < 0, for 0 < x < a;
F(x) > 0 and nondecreasing for x > a;
F(x) → ∞ as x → ∞;
The Liénard equation 21.12 appears naturally in series RLC circuits with memristors or nanodevices with nonlinear v–i characteristics.
Consider the circuit shown in Figure 21.5 with a memristor ϕ = f(q).
The corresponding differential equation,
(21.17) |
Here, the memristance M = df(q)/dq is quadratic in relation with the charge,
(21.18) |
where (k,α) are constants; replacing Equation 21.18 by Equation 21.17 with i = dq/dt yields,
(21.19) |
This Liénard equation, for the memristive circuit shown in Figure 21.5, has coefficients a0(x) and a1(x) that satisfy the limit cycle conditions 1–4. Equation 21.19 can be written as Equation 21.12 with coefficients (21.20).
(21.20) |
Condition five for a limit cycle (21.16) requires the integration of a1(x),
(21.21) |
Finding the positive zero,
(21.22) |
Formula (21.22) is dependent on the memristor parameters (k,α) and the resistive part R of the circuit, that is, the elements that damp oscillations or dissipative elements.
Figure 21.6 is a series RLC circuit with a general nonlinear nanodevice v = φ(i).
The differential equation in this case includes the function φ(i),
(21.23) |
If the input u is constant, then the derivative with respect to time of Equation 21.23 follows the expression,
(21.24) |
In some cases, from the nanodevice experimental v–i characteristic, it could be possible to interpolate a continuous function v = φ(i).
In this section, two series RLC circuits with quadratic memristors are simulated [10] to verify the presence of limit cycles if the five conditions of the Liénard equation are satisfied [3]. The parameters for the simulation are: R = 4 Ω, L = 2 H, and C = 0.5 F.
In the first case, the constants for M(q) are (k, α) = (10, 0.5); the positive zero is located at x = 0.55. Integral (21.16) is the cubic equation (21.25) shown in Figure 21.7; notice that F(x) satisfies condition five for a limit cycle.
(21.25) |
Writing the unforced Liénard equation for this memristive circuit (k,α) = (10, 0.5) with initial conditions (21.26),
(21.26) |
Figure 21.8 is the solution of Equation 21.26; there is a limit cycle because a1(x); a0(x) satisfy all the required conditions.
In the second case, the constants for M(q) are (k,α) = (2, 0.5); there is no positive zero; instead, x = ± j 2.12. Integral (21.16) is the cubic equation (21.27) shown in Figure 21.9; notice that F(x) does not satisfy condition five for a limit cycle.
(21.27) |
Writing the unforced Liénard equation for this memristive circuit (k,α) = (2, 0.5) with initial conditions (21.28),
(21.28) |
Figure 21.10 is the solution of Equation 21.28; there is no limit cycle because a1(x); a0(x) do not satisfy the required conditions.
This chapter has presented two methods, frequency-domain analysis and time-domain analysis, to study circuits and systems with nonlinear nanodevices; the corresponding mathematical concepts are the DF and the Liénard equation, respectively.
The DF for a general v–i characteristic was calculated and the final formula, with particular parameters, was used to derive the DFs for a memristor and for DNA and CNTs-based nanodevices.
In the time domain, the Liénard equation was formulated for RLC series circuits with nanodevices; two types of nonlinearities were considered, v = φ(i) and ϕ = f(q). Conditions for the existence of a limit cycle were formulated and the theory was illustrated with simulations.
Future work could consider a time-varying DF to improve the modeling of frequency-dependent nanodevices such as the memristor [5]; also, the design of nanodevices to synthesize particular v–i characteristics is an open problem.
1. Khalil, H.K. Nonlinear Systems. Macmillan Publishing, New York, NY, 1992.
2. Delgado, A. The memristor as controller. IEEE Nanotechnology Materials and Devices Conference, Monterey, California, Oct. 12–15, 2010.
3. Strogatz, S.H. Nonlinear Dynamics and Chaos. Perseus Books Publishing, Cambridge, 1994.
4. Chua, L.O. Memristor—The missing circuit element. IEEE Transactions on Circuit Theory, 18(5), 507–519, 1971.
5. Strukov, D.B., Snider, G.S., Stewart, D.R., and Williams, R.S. The missing memristor found. Nature, 453, 80–83, 2008.
6. Jeong, H. J., Lee, J. Y., Ryu, M. K., and Choi, S. Y. Bipolar resistive switching in amorphous titanium oxide thin film. Physica Status Solidi RRL 4(1–2), 28–30, 2010.
7. Cohen, H., Nogues, C., Naaman, R., and Porath, D. Direct measurement of electrical transport through single DNA molecules of complex sequence. Proceedings of the National Academy of S ciences, USA, 102(33), 11589–11593, 2005.
8. Kaun, C.-C., Larade, B., Mehrez, H., Taylor, J., and Guo, H. Current—voltage characteristics of carbon nanotubes with substitutional nitrogen. Physical Review B, 65, 205416 (1–5), 2002.
9. Wolfram Alpha LLC. 2012. Wolfram|Alpha. http://www.wolframalpha.com/input/?i=integrate (access May 22, 2012). The use of Wolfram|Alpha does not imply in any way that Wolfram endorses the particular substance or general quality of this work.
10. Consortium Scilab–Digiteo 2011. Scilab: Free and Open Source software for numerical computation (OS, Version 5.XX) [Software]. Available from: http://www.scilab.org.