Active and Passive Electronic Components

Active and Passive Electronic Components / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 830182 |

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.

Modelling of an Esaki Tunnel Diode in a Circuit Simulator

Academic Editor: G. Ghibaudo
Received15 Jan 2011
Accepted26 Feb 2011
Published17 Apr 2011


A method for circuit-level 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 single-valued.

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 semiconductor-based tunnel diodes, molecular diodes, and in current leakage through the thin oxides of advanced metal-oxide-metal (MOS) transistors, see Olivo et al. [1]. Tunnel diodes, see Hall [2], have been used in fast-switching 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 flux-potential 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 current-voltage (I-V) 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 sub-65 nm MOS transistors where the gate oxide is necessarily thin. The typical current-voltage (I-V) 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 forward-biased p-n 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 quantum-physical 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 I-V 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 quantum-physical 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 curve-fitting to a relatively smooth function, and modifications to the underlying algorithms that solve a circuit containing a tunnel diode. The macro-modelling 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 voltage-controlled current sources arranged such that they will approximately mimic the tunnel diode I-V characteristic. Curve-fitting has also been used extensively wherein the I-V characteristic was approximated by a piecewise-linear 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 curve-fitting techniques, Chang et al. [8], have also been used. Macro-modelling and curve fitting are essentially nonphysical techniques introducing well-behaved functions to model a device whose operation is quite complex.

Physics-based modelling is essential when a strong correlation between simulation and measurement is desired—particularly with steadily decreasing device geometries. A physics-based 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 Newton-Raphson iterative method and implemented in a circuit simulator. Based on the onset of certain conditions during simulation, the algorithm switches from voltage-based iterations to current-based iterations. Since the I-V 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 homotopy-based 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 well-established 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 I-V 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,𝐯(𝑡)=𝐯𝐱(𝑡),𝑑𝐱𝑑𝑑𝑡,…,𝑚𝐱𝑑𝑡𝑚,𝐱𝐷,(𝑡)𝐢(𝑡)=𝐢𝐱(𝑡),𝑑𝐱𝑑𝑑𝑡,…,𝑚𝐱𝑑𝑡𝑚,𝐱𝐷,(𝑡)(1) 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 time-delayed 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 ad-hoc modifications of the iterative algorithm required, enables state-of-the-art, off-the-shelf numerical libraries to be employed.

Consider the equation for a microwave diode whose I-V 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 non-convergent solution, or worse still, an incorrect solution. In order to continue to use the physics-based 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𝑖(𝑡)=𝐼𝑠exp(𝛼𝑣(𝑡))−1,(2) 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 𝑉0, and is defined as a linear function of the current in (2) above 𝑉0. Care is taken to ensure that the current and its derivatives are continuous at 𝑥=𝑉0. The equations for diode voltage and current are now expressed as𝑉𝑣(𝑡)=0+1𝛼ln1+𝛼𝑥(𝑡)−𝑉0,𝑉0≤𝑥(𝑡),𝑥(𝑡),𝑥(𝑡)≤𝑉0,𝑖𝐼(𝑡)=𝑠exp𝛼𝑉01+𝛼𝑥(𝑡)−𝑉0−𝐼𝑠,𝑉0𝐼≤𝑥(𝑡),𝑠exp(𝛼𝑣(𝑡))−1,𝑥(𝑡)≤𝑉0,(3) where (3) is the exact parametric representations of the diode I-V characteristic. The highly nonlinear I-V characteristic is now transformed into current-state-variable (I-X) and voltage-state-variable (V-X) 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 I-X relationship, and in Figure 3, the V-X relationship, with 𝐼𝑠=100 aA, 𝛼=38.686, and 𝑉0=0.8 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 time-dependent 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 forward-bias direction. Referring to Figure 1, this is expressed as𝐽total=𝑉(𝑡)𝑉𝑃𝐽𝑇+𝐽𝑋+𝐽TH,(4) where𝐽𝑇=𝐽𝑃exp1−𝑉(𝑡)𝑉𝑃𝐽,(5)𝑋=𝐽𝑉𝐴exp2𝑉(𝑡)−𝑉𝑉𝐽,(6)TH=𝐽𝑆expğ‘žğ‘‰(𝑡)𝑘𝑇−1.(7) The first term is a closed-form 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 𝐴2 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 I-V 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 ğ½ğ‘‡âŽ§âŽªâŽªâŽ¨âŽªâŽªâŽ©ğ½(𝑡)=𝑃𝑉exp1−𝑇𝑉𝑃×11−𝑉𝑃𝑥𝑇(𝑡)−𝑉𝑇,𝑥𝑇(𝑡)≤𝑉𝑇,𝐽𝑃𝑥exp1−𝑇(𝑡)𝑉𝑃,𝑥𝑇(𝑡)>𝑉𝑇,ğ‘‰âŽ§âŽªâŽªâŽ¨âŽªâŽªâŽ©ğ‘‰(𝑡)=𝑇−𝑉𝑃1×ln1−𝑉𝑃𝑥𝑇(𝑡)−𝑉𝑇,𝑥𝑇(𝑡)≤𝑉𝑇,𝑥𝑇(𝑡),𝑥(𝑡)>𝑉𝑇,ğ½ğ‘‹âŽ§âŽªâŽ¨âŽªâŽ©ğ½(𝑡)=𝑉𝐴exp2𝑥X(𝑡)−𝑉𝑉,𝑥𝑋(𝑡)≤𝑉𝑋,𝐽𝑉𝐴exp2𝑉𝑋(𝑡)−𝑉𝑉×1+𝐴2𝑥𝑋(𝑡)−𝑉𝑋,𝑥𝑋(𝑡)>𝑉𝑋,âŽ§âŽªâŽ¨âŽªâŽ©ğ‘¥ğ‘‰(𝑡)=𝑋(𝑡),𝑥𝑋(𝑡)≤𝑉𝑋,𝑉𝑋+1𝐴2×ln1+𝐴2𝑥𝑋(𝑡)−𝑉𝑋,𝑥𝑋(𝑡)>𝑉𝑋,𝐽TH=âŽ§âŽªâŽªâŽ¨âŽªâŽªâŽ©ğ½ğ‘†î‚€expğ‘žğ‘¥TH𝑘𝑇−𝐽𝑆,𝑥TH≤𝑉TH,𝐽𝑆expğ‘žğ‘‰THî‚¶Ã—î‚€ğ‘žğ‘˜ğ‘‡1+𝑥𝑘𝑇TH−𝑉TH−𝐽𝑆,𝑥TH(𝑡)>𝑉TH,âŽ§âŽªâŽ¨âŽªâŽ©ğ‘¥ğ‘‰(𝑡)=TH(𝑡),𝑥TH≤𝑉TH,𝑉TH+ğ‘˜ğ‘‡ğ‘žî‚€ğ‘žln1+𝑥𝑘𝑇TH(𝑡)−𝑉TH,𝑥TH>𝑉TH.(8) 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 𝑥TH(𝑡). The linear version of 𝐽𝑇 is used for 𝑥𝑇(𝑡)≤𝑉𝑇 while the linear versions of 𝐽𝑋 and 𝐽TH are respectively used for 𝑥𝑋(𝑡)>𝑉𝑋 and 𝑥TH(𝑡)>𝑉TH. 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 𝑥TH(𝑡) for 𝑥𝑋(𝑡)>𝑉𝑋 and 𝑥TH(𝑡)>𝑉TH, respectively.

Careful study of the various possible relationships between 𝑉𝑇, 𝑉𝑋, and 𝑉TH reveals that the transition between the various regions will be continuous as long as 𝑉𝑇 is chosen to be less than or equal to 𝑉𝑋 and 𝑉TH. Using this constraint, continuous equations for the current densities and voltage as a function of a single 𝑥(𝑡) are derived. Note that min{𝑉𝑋,𝑉TH} has been used as a criteria to determine whether the 𝑥𝑋(𝑡) or 𝑥TH(𝑡) should be used for large 𝑥(𝑡) where 𝑉X and 𝑉TH 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ğ½ğ‘‡âŽ§âŽªâŽªâŽªâŽ¨âŽªâŽªâŽªâŽ©ğ½(𝑡)=𝑃𝑉exp1−𝑇𝑉𝑃×11−𝑉𝑃𝑥(𝑡)−𝑉𝑇,𝑥(𝑡)≤𝑉𝑇,𝐽𝑃exp1−𝑥(𝑡)𝑉𝑃,𝑉𝑇<𝑥(𝑡)≤𝑉XTH,𝐴𝑇2×1+𝐵𝑥(𝑡)−𝑉XTH𝐶𝑇,𝑥(𝑡)>𝑉XTH,ğ½ğ‘‹âŽ§âŽªâŽªâŽªâŽ¨âŽªâŽªâŽªâŽ©ğ½(𝑡)=𝑉𝐴exp2𝑉𝑇−𝑉𝑉×11−𝑉𝑃𝑥(𝑡)−𝑉𝑇−𝐴2𝑉𝑃,𝑥(𝑡)≤𝑉𝑇,𝐽𝑉𝐴exp2𝑥(𝑡)−𝑉𝑉,𝑉𝑇<𝑥(𝑡)≤𝑉XTH,𝐴𝑋21+𝐵𝑥(𝑡)−𝑉XTH𝐶𝑋,𝑥(𝑡)>𝑉XTH,𝐽TH⎧⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎩𝐽(𝑡)=𝑆expğ‘žğ‘‰ğ‘‡î‚¶Ã—î‚µ1𝑘𝑇1−𝑉𝑃(𝑥(𝑡)−𝑉𝑇)î‚¶âˆ’ğ‘žğ‘‰ğ‘ƒ/𝑘𝑇−𝐽𝑆,𝑥(𝑡)≤𝑉𝑇,𝐽𝑆expğ‘žğ‘¥(𝑡)𝑘𝑇−𝐽𝑆,𝑉𝑇<𝑥(𝑡)≤𝑉XTH,𝐴TH2×1+𝐵𝑥(𝑡)−𝑉XTH𝐶TH−𝐽𝑆,𝑥(𝑡)>𝑉XTH,âŽ§âŽªâŽªâŽªâŽ¨âŽªâŽªâŽªâŽ©ğ‘‰ğ‘‰(𝑡)=𝑇−𝑉𝑃1×ln1−𝑉𝑃𝑥(𝑡)−𝑉𝑇,𝑥(𝑡)≤𝑉𝑇,𝑥(𝑡),𝑉𝑇<𝑥(𝑡)<𝑉XTH,𝑉XTH+1𝐵×ln1+𝐵𝑥(𝑡)−𝑉XTH,𝑥(𝑡)>𝑉XTH,(9) which include a number of constant terms. The expressions for the voltage thresholds, 𝑉𝑇, 𝑉𝑋, and 𝑉TH, 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:𝑉𝑇=𝑉𝑃−𝑉𝑃ln−𝑚𝑇𝑉𝑃𝐽𝑃,slope𝑚𝑇𝑉<0,𝑋=𝑉𝑉+1𝐴2𝑚ln𝑋𝐴2𝐽𝑉,slope𝑚𝑋𝑉>0,TH=ğ‘˜ğ‘‡ğ‘žî‚µğ‘šlnTHğ‘˜ğ‘‡ğ½ğ‘†ğ‘žî‚¶,slope𝑚TH>0.(10) When 𝑉𝑋≤𝑉TH the other constants are,𝑉XTH=𝑉𝑋,𝐴T2=𝐽𝑃𝑉exp1−𝑋𝑉𝑃,𝐴X2=𝐽𝑉𝐴exp2𝑉𝑋−𝑉𝑉,𝐴TH2=𝐽𝑆expğ‘žğ‘‰ğ‘‹î‚¶,𝑘𝑇𝐵=𝐴2,𝐶𝑇1=−𝐴2𝑉𝑃,𝐶𝑋𝐶=1,TH=ğ‘žğ‘˜ğ‘‡ğ´2,(11) and when 𝑉𝑋>𝑉TH,𝑉XTH=𝑉TH,𝐴T2=𝐽𝑃𝑉exp1−TH𝑉𝑃,𝐴X2=𝐽𝑉𝐴exp2𝑉TH−𝑉𝑉,𝐴TH2=𝐽𝑆expğ‘žğ‘‰TH,ğ‘žğ‘˜ğ‘‡ğµ=,𝐶𝑘𝑇𝑇=âˆ’ğ‘˜ğ‘‡ğ‘žğ‘‰ğ‘ƒ,𝐶𝑋=𝐴2ğ‘˜ğ‘‡ğ‘ž,𝐶TH=1.(12)

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 zero-bias depletion capacitance parameter CT0 is set. If it is set, then𝐶𝑗=âŽ§âŽªâŽªâŽ¨âŽªâŽªâŽ©ğ™²ğšƒğŸ¶î€·1−𝑉𝑗/𝙵𝙸𝙶𝙰𝙼𝙰,−𝑉𝑗>0,𝙲𝚃𝟶+𝑉𝑗𝑘1+𝙲𝚃𝟶𝑉𝑗1+𝑘3+𝑘2𝑘31+𝑘3,−𝑉𝑗<0,(13) where 𝑘1=𝙲𝚃𝟶×𝙶𝙰𝙼𝙰/𝙵𝙸,𝑘2=𝙲𝚃𝟶/𝙶𝙰𝙼𝙰0.2,𝑘3=exp(10(𝑉𝑗−0.8𝙵𝙸)), and 𝑉𝑗 is the diode junction voltage. If the depletion capacitance parameter, CD0, is set, then 𝐶𝑗 is incremented by the amount 𝙲𝙳𝟶×exp(𝙰𝙵𝙰𝙲𝑉𝑗).

Parameter Description Default Units

JS Saturation current 1 × 1 0 − 1 6 A
CT0 Zero bias depletion capacitance 0 F
FI Built-in barrier potential 0.8 V
GAMA Capacitance power law parameter 0.5 —
CD0 Zero bias diffusion capacitance 0 F
AFAC Slope factor of diffusion capacitance 38.696 1/V
R0 Series resistance in forward bias 0 Ω
TAU Intrinsic depletion layer time constant 0 sec
AREA Area multiplier 1 —
JV Valley current 1 × 1 0 − 4 A
JP Peak current 1 × 1 0 − 3 A
VV Valley voltage 0.5 V
VPK Peak voltage 0.1 V
A2 Excess current prefactor 30 —
MT Tunnel current slope factor −1 —
MX Excess current slope factor 1 —
MTH Thermal current slope factor 1 —
TEMP Temperature 3 0 0 ∘ K

The parameterised I-V relationship using the constants listed in Table 1 are compared with the original I-V 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 open-source circuit simulator ( Simulated I-V 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 I-V characteristics and this is displayed in Figure 9. The diode used in measurements is a general-purpose germanium diode with a nominal peak voltage 𝑉𝑝 of 105 mV, valley voltage 𝑉𝑣 of 380 mV peak-to-valley 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.

Parameter Value

JS 2 × 1 0 − 1 2
JV 1 . 1 × 1 0 − 4
VV 0.35
VPK 0.065
A2 10
CT0 0 . 0 3 × 1 0 − 1 2
GAMA 0.19
FI 0.0368
R0 1.5

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 well-behaved, 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, off-the-shelf 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.


  1. P. Olivo, T. N. Nguyen, and B. Ricco, “High-field-induced degradation in ultra-thin SiO2 films,” IEEE Transactions on Electron Devices, vol. 35, no. 12, pp. 2259–2267, 1988. View at: Publisher Site | Google Scholar
  2. R. N. Hall, “Tunnel diodes,” IRE Transactions on Electron Devices, vol. 7, no. 1, pp. 1–9, 1960. View at: Google Scholar
  3. L. Esaki, “New phenomenon in narrow germanium p-n junctions,” Physical Review, vol. 109, no. 2, pp. 603–604, 1958. View at: Publisher Site | Google Scholar
  4. S. M. Sze, Physics of Semiconductor Devices, John Wiley & Sons, New York, NY, USA, 2nd edition, 1981.
  5. T. H. Kuo, H. C. Lin, U. Anandakrishnan, R. C. Potter, and D. Shupe, “Large-signal 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
  6. S. Mohan, J. P. Sun, P. Mazumder, and G. I. Haddad, “Device and circuit simulation of quantum electronic devices,” IEEE Transactions on Computer-Aided Design, vol. 14, no. 6, pp. 653–662, 1995. View at: Google Scholar
  7. 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
  8. C. E. Chang, P. M. Asbeck, K. C. Wang, and E. R. Brown, “Analysis of heterojunction bipolar transistor/resonant tunneling diode logic for low-power and high-speed digital applications,” IEEE Transactions on Electron Devices, vol. 40, no. 4, pp. 685–691, 1993. View at: Publisher Site | Google Scholar
  9. M. Bhattacharya and P. Mazumder, “Augmentation of SPICE for simulation of circuits containing resonant tunneling diodes,” IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, vol. 20, no. 1, pp. 39–50, 2001. View at: Publisher Site | Google Scholar
  10. 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 Computer-Aided Design of Integrated Circuits and Systems, vol. 12, no. 6, pp. 861–877, 1993. View at: Publisher Site | Google Scholar
  11. C. E. Christoffersen, U. A. Mughal, and M. B. Steer, “Object oriented microwave circuit simulation,” International Journal of RF and Microwave Computer-Aided Engineering, vol. 10, no. 3, pp. 164–182, 2000. View at: Publisher Site | Google Scholar
  12. V. Rizzoli and A. Neri, “Expanding the power-handling capabilities of harmonic-balance analysis by a parametric formulation of the MESFET model,” Electronics Letters, vol. 26, no. 17, pp. 1359–1361, 1990. View at: Google Scholar

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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles