Research Article  Open Access
Nikhil M. Kriplani, Stephen Bowyer, Jennifer Huckaby, Michael B. Steer, "Modelling of an Esaki Tunnel Diode in a Circuit Simulator", Active and Passive Electronic Components, vol. 2011, Article ID 830182, 8 pages, 2011. https://doi.org/10.1155/2011/830182
Modelling of an Esaki Tunnel Diode in a Circuit Simulator
Abstract
A method for circuitlevel modelling a physically realistic Esaki tunnel diode model is presented. A paramaterisation technique that transforms the strongly nonlinear characteristic of a tunnel diode into two relatively modest nonlinear characteristics is demonstrated. The introduction of an intermediate state variable results in a physically realistic mathematical model that is not only moderately nonlinear and therefore robust, but also singlevalued.
1. Introduction
Tunnelling here refers to the phenomenon of movement of charge across a potential barrier with the result that there is a nonunique relationship between current and applied voltage. The effect occurs in semiconductorbased tunnel diodes, molecular diodes, and in current leakage through the thin oxides of advanced metaloxidemetal (MOS) transistors, see Olivo et al. [1]. Tunnel diodes, see Hall [2], have been used in fastswitching electronic circuits ever since their discovery, Esaki [3]. While tunnelling is a critical effect, it can be difficult to model in a circuit simulator because of the multivalued fluxpotential relationship which is also strongly nonlinear. While this paper focuses on the modelling of tunnel diodes, the technique presented can be used to model many types of tunnelling effects.
One unique feature of an Esaki tunnel diode is the Negative Differential Resistance (NDR) region of its currentvoltage (IV) characteristic. Under certain circuit conditions, this creates a bistable or multistable circuit, and switching between the two stable states can be rapid. In other situations, tunnelling is undesirable and sets limits on performance. Determining the onset of tunnelling is important in designing integrated circuits using sub65โnm MOS transistors where the gate oxide is necessarily thin. The typical currentvoltage (IV) characteristic of an Esaki tunnel diode can be divided into three regions in the direction of increasing voltage, see Figure 1. The first region involves an increase in the output current until it reaches a peak value and represents a Positive Differential Resistance (PDR). The region following this involves a decrease in current with increasing voltage and this corresponds to the NDR region. Here the current decreases until it reaches a minimum or valley value. The third region corresponds to conduction of a forwardbiased pn junction diode and is the second PDR region. All three regions are represented by exponential functions and are therefore strongly nonlinear. These exponential functions are derived from a quantumphysical description of the diode, Sze [4], and the total current is a sum of these exponential functions. In order to simulate circuits consisting of Esaki tunnel diodes, it is essential to be able to accurately model the nonlinear behaviour of the IV characteristic as well as the negative resistance. Ideally, one would like to model the Esaki tunnel diode in a circuit simulator with minimal changes to the quantumphysical representation of the device.
An Esaki tunnel diode is not directly modelled in commercial circuit simulators. Three strategies have been used to mitigate the bistable and/or strongly nonlinear physical description of the diode. Techniques used include: macromodelling approximation by curvefitting to a relatively smooth function, and modifications to the underlying algorithms that solve a circuit containing a tunnel diode. The macromodelling was used in Kuo et al. [5] wherein a tunnel diode was represented as a set of linear current sources which were appropriately switched on or off depending on bias conditions. This approach uses combinations of switches, resistors, capacitors, voltage sources, and voltagecontrolled current sources arranged such that they will approximately mimic the tunnel diode IV characteristic. Curvefitting has also been used extensively wherein the IV characteristic was approximated by a piecewiselinear model, see Mohan et al. [6], or by a piecewise nonlinear model utilising a diode for each of the PDR regions, see Neculoiu and Tebeanu [7]. Similarly, polynomial curvefitting techniques, Chang et al. [8], have also been used. Macromodelling and curve fitting are essentially nonphysical techniques introducing wellbehaved functions to model a device whose operation is quite complex.
Physicsbased modelling is essential when a strong correlation between simulation and measurement is desiredโparticularly with steadily decreasing device geometries. A physicsbased description of the diode was used in Bhattacharya and Mazumder [9] where the iterative technique used to solve the nonlinear set of differential equations describing a circuit was modified. In particular, algorithmic modifications were made to the NewtonRaphson iterative method and implemented in a circuit simulator. Based on the onset of certain conditions during simulation, the algorithm switches from voltagebased iterations to currentbased iterations. Since the IV characteristic of the tunnel diode can have multiple values of voltage for a single value of current, the authors proposed a technique to approximately assign a single value of voltage and resolve this uncertainty. In order to calculate this voltage value, the authors introduced additional parameters in the diode model. Other techniques proposed to model circuits with tunnel diodes can be thought of as variations to homotopybased methods as described in Melville et al. [10], and in general require a steep learning curve on the part of the user to be able to understand and implement the model of a nonlinear device. Also, the generality of the modified iterative techniques and the homotopy methods is limited.
In this paper, a technique for modelling a physically based tunnel diode is demonstrated that makes no changes to the basic equations describing the operation of a device and requires no modification to the underlying algorithms of a circuit simulator. A similar approach can be used for other strongly nonlinear devices.
2. Parameterised Device Modelling with State Variables
Simulating circuits with strongly nonlinear devices can often be plagued by nonconvergence. In the wellestablished circuit simulator tradition, modifications are made at the individual element level to obtain convergence. This philosophical procedure is used in the SPICE simulator to obtain convergence for circuits with junction diodes that have exponential characteristics. In SPICE, heuristics are used including limiting the fraction by which current and voltage can change from one iterative step to the next. This cannot be extended to handling tunnel diodes. In essence, the technique described in this paper introduces an intermediate state variable between the current and voltage variables in order to mitigate the nonlinear diode IV characteristics. The iterative solver then acts on this intermediate state variable while trying to minimize the error function associated with the circuit under consideration. As outlined in Christofferssen et al. [11], when an element is modelled using state variables, the current and voltage are expressed as a function of these state variables. This is, where and are vectors of voltages and currents at the ports of the nonlinear device, is a vector of state variables, and is a vector of timedelayed state variables, that is, . Equation (1) defines the modelling scope. This can be a problem if the describing equations are strongly nonlinear as local convergence control cannot be used. However, avoiding the use of local convergence control, and the adhoc modifications of the iterative algorithm required, enables stateoftheart, offtheshelf numerical libraries to be employed.
Consider the equation for a microwave diode whose IV relationship is described by an exponential function. Once the diode is forward biased, a small change in input voltage can produce a large change in output current. This nonlinear behaviour can be a source of serious simulation errors that can manifest themselves as an incomplete or nonconvergent solution, or worse still, an incorrect solution. In order to continue to use the physicsbased equations for a nonlinear device and overcome the problem of strong nonlinear behaviour, the technique of parameterisation was introduced in Rizzoli and Neri [12]. This is illustrated by considering the example of a microwave diode.
The conventional current equation for a diode is expressed as where is the voltage across the diode, is the reverse saturation current, and is a parameter representing the slope of the conductance curve of the diode characteristic. Exponential relationships of this type are responsible for most of the numerical instabilities in circuit simulation. To circumvent this, the state variable in (2) is changed from voltage to a fictitious nonconventional variable, . This state variable is identical to the junction voltage below some nonzero threshold value, say , and is defined as a linear function of the current in (2) above . Care is taken to ensure that the current and its derivatives are continuous at . The equations for diode voltage and current are now expressed as where (3) is the exact parametric representations of the diode IV characteristic. The highly nonlinear IV characteristic is now transformed into currentstatevariable (IX) and voltagestatevariable (VX) characteristics that are not as strongly nonlinear, without any changes to the essential mathematical model of the device. The solution is well behaved and local convergence or heuristics within the diode model are not required. After the transformations of (3), the reduced nonlinear nature of the relationship is illustrated in Figure 2, the IX relationship, and in Figure 3, the VX relationship, with โaA, , and โV.
3. Modelling the Esaki Tunnel Diode
The intrinsic Esaki tunnel diode model consists of an ideal tunnel diode a nonlinear capacitor and a nonlinear resistor, see Figure 4. There are two controlling variables in this model, which are an as yet unspecified timedependent state variable, and its time derivative. The ideal tunnel diode is modelled using current and voltage as a function of the primary state variable . The capacitance is modelled by using the evaluated voltage and the time derivative of the state variable, , to calculate the derivative of charge as a function of time or . From this, the current through the capacitor can be determined. The resistor is modelled using the evaluated voltage and the total current. An outline of the model calculation of the total current and voltage is shown in Figure 5.
The tunnel diode current density is typically described as the sum of three exponential functions derived from quantum mechanical considerations. This formulation appears in Sze [4], although here the physics is limited only to the forwardbias direction. Referring to Figure 1, this is expressed as where The first term is a closedform expression of the tunnelling current density which describes the behaviour particular to the tunnel diode. This includes the negative resistance region which captures the core functionality of the tunnel diode. The second term describes the excess tunnelling current density while the third term is the normal diode characteristic. In (5), is the peak current density and is the corresponding peak voltage. In (6), is the peak current density and is the corresponding peak voltage. The parameter represents an excess current prefactor. Finally, in (7), is the saturation current density, is the charge of an electron, is Boltzmann's constant, and is the temperature in degrees Kelvin.
As mentioned in the previous section, the exponential functions can lead to a rapid change in current for a relatively small change in voltage which can create convergence issues for the simulator's nonlinear solver. The tunnel diode equation cannot be solved for current as a single and unique function of voltage. Recognising that the three major components of current density can be thought of as three diodes with exponential IV characteristics, they are parameterised as previously described. After individual parameterisation of the three regions based on equating current density and first derivatives at threshold points, the tunnel diode equations become As can be seen from the expressions above, each current is a function of a different . Since the voltage across each current source is the same, the voltage equations can be set equal to each other to obtain a relationship between , , and . The linear version of is used for while the linear versions of and are respectively used for and . Therefore, when combining the terms as a function of a single state variable , they should be a function of for and a function of or for and , respectively.
Careful study of the various possible relationships between , , and reveals that the transition between the various regions will be continuous as long as is chosen to be less than or equal to and . Using this constraint, continuous equations for the current densities and voltage as a function of a single are derived. Note that has been used as a criteria to determine whether the or should be used for large where and are calculated based on the first derivative of the current. The expressions for the diode current density and voltage as a function of a single state variable are expressed as which include a number of constant terms. The expressions for the voltage thresholds, , , and , define the boundary between the current regions. These are derived by taking the first derivative of the current equation and solving for voltage as a function of slope while maintaining current continuity: When the other constants are, and when ,
The tunnel diode model is completed by incorporating the diode's junction capacitance . Referring to the parameter values for the diode shown in Table 1, the model of depends on whether the zerobias depletion capacitance parameter CT0 is set. If it is set, then where , and is the diode junction voltage. If the depletion capacitance parameter, CD0, is set, then is incremented by the amount .

The parameterised IV relationship using the constants listed in Table 1 are compared with the original IV equation and in Figure 6. For the sake of brevity, only the total current equation is plotted instead of the individual current components. As can be seen, the results are identical and the parameterisation process has not affected the accuracy of the underlying equations for the diode.
The benefits of parameterisation are demonstrated by plotting total current and voltage with respect to the suitably introduced state variable . As can be seen in Figure 7, which plots tunnel diode voltage with respect to a state variable, and Figure 8, which plots tunnel diode current with respect to , the highly nonlinear relationship between current and voltage is replaced by two milder nonlinearities.
4. Simulation and Results
The parameterised Esaki tunnel diode model described above was implemented in an opensource circuit simulator (http://www.freeda.org/). Simulated IV characteristics of the tunnel diode using the parameter values in Table 1 are generated by sweeping the voltage from โ0.05 volts to 0.6 volts. This voltage range is the typical range of operation for this device. The simulated curve is compared with measured IV characteristics and this is displayed in Figure 9. The diode used in measurements is a generalpurpose germanium diode with a nominal peak voltage of 105โmV, valley voltage of 380โmV peaktovalley current ratio of 8. A match obtained between simulation and measurement thereby validates the effectiveness of this modelling technique. In other words, it is possible to easily and efficiently implement a highly nonlinear device keeping the physics of the device intact while requiring no change to the underlying algorithm that solves the nonlinear equations associated with the circuit.
A transient simulation was also performed on a canonical oscillator circuit containing the tunnel diode to demonstrate the switching characteristics of the device. The circuit used along with the values of the components is shown in Figure 10. A snapshot of the simulated output voltage versus time in Figure 11 shows the circuit oscillating at approximately 150 MHz. The values of the diode parameters used for this simulation that differ from those in Table 1 are indicated in Table 2.

5. Conclusions
Tunnelling is an important aspect of charge transport in semiconductor and molecular devices. While a semiconductor tunnel diode with characteristics described by sums of exponential functions was considered, the parameterisation technique introduced can be used with any analytic expression. The central result is transcribing a circuit simulation problem from one dealing with strong nonlinearities to one working with wellbehaved, moderately nonlinear element characteristics. This is achieved without an increase in problem size. The introduction of an alternative state (or state variable) enables the strongly nonlinear character of tunnelling to be modelled, with full accuracy, by intermediate equations which are relatively well behaved and moderately nonlinear. Thus, offtheshelf numerics can be used in a circuit simulator and there is no need for local heuristics, homotopy, or functional approximation to obtain convergence. The fundamental requirement is that the circuit simulator support state variables. The state variables replace the (usual) nodal voltages as the unknowns and the error function becomes the energy norm rather than solely based on Kirchoff's Current Law while the total number of unknowns remains unchanged.
References
 P. Olivo, T. N. Nguyen, and B. Ricco, โHighfieldinduced degradation in ultrathin SiO_{2} films,โ IEEE Transactions on Electron Devices, vol. 35, no. 12, pp. 2259โ2267, 1988. View at: Publisher Site  Google Scholar
 R. N. Hall, โTunnel diodes,โ IRE Transactions on Electron Devices, vol. 7, no. 1, pp. 1โ9, 1960. View at: Google Scholar
 L. Esaki, โNew phenomenon in narrow germanium pn junctions,โ Physical Review, vol. 109, no. 2, pp. 603โ604, 1958. View at: Publisher Site  Google Scholar
 S. M. Sze, Physics of Semiconductor Devices, John Wiley & Sons, New York, NY, USA, 2nd edition, 1981.
 T. H. Kuo, H. C. Lin, U. Anandakrishnan, R. C. Potter, and D. Shupe, โLargesignal resonant tunneling diode model for SPICE3 simulation,โ in Proceedings of the International Electron Devices Meeting, pp. 567โ570, Washington, DC, USA, December 1989. View at: Google Scholar
 S. Mohan, J. P. Sun, P. Mazumder, and G. I. Haddad, โDevice and circuit simulation of quantum electronic devices,โ IEEE Transactions on ComputerAided Design, vol. 14, no. 6, pp. 653โ662, 1995. View at: Google Scholar
 D. Neculoiu and T. Tebeanu, โSPICE implementation of double barrier resonant tunnel diode model,โ in Proceedings of the International Semiconductor Conference (CAS '96), pp. 181โ184, Bucharest, Romania, October 1996. View at: Google Scholar
 C. E. Chang, P. M. Asbeck, K. C. Wang, and E. R. Brown, โAnalysis of heterojunction bipolar transistor/resonant tunneling diode logic for lowpower and highspeed digital applications,โ IEEE Transactions on Electron Devices, vol. 40, no. 4, pp. 685โ691, 1993. View at: Publisher Site  Google Scholar
 M. Bhattacharya and P. Mazumder, โAugmentation of SPICE for simulation of circuits containing resonant tunneling diodes,โ IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems, vol. 20, no. 1, pp. 39โ50, 2001. View at: Publisher Site  Google Scholar
 R. C. Melville, L. Trajkovic, S.C. Fang, and L. T. Watson, โArtificial parameter homotopy methods for the DC operating point problem,โ IEEE Transactions on ComputerAided Design of Integrated Circuits and Systems, vol. 12, no. 6, pp. 861โ877, 1993. View at: Publisher Site  Google Scholar
 C. E. Christoffersen, U. A. Mughal, and M. B. Steer, โObject oriented microwave circuit simulation,โ International Journal of RF and Microwave ComputerAided Engineering, vol. 10, no. 3, pp. 164โ182, 2000. View at: Publisher Site  Google Scholar
 V. Rizzoli and A. Neri, โExpanding the powerhandling capabilities of harmonicbalance analysis by a parametric formulation of the MESFET model,โ Electronics Letters, vol. 26, no. 17, pp. 1359โ1361, 1990. View at: Google Scholar
Copyright
Copyright © 2011 Nikhil M. Kriplani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.