Abstract

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(𝙰𝙡𝙰𝙲𝑉𝑗).

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 (http://www.freeda.org/). 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.

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.