[ Chapter start ] [ Previous page ] [ Next page ]
13.10 TransistorLevel Simulation
Sometimes we need to simulate a logic circuit with more accuracy than provided by switchlevel simulation. In this case we turn to simulators that can solve circuit equations exactly, given models for the nonlinear transistors, and predict the analog behavior of the node voltages and currents in continuous time. This type of transistorlevel simulation or circuitlevel simulation is costly in computer time. It is impossible to simulate more than a few hundred logic cells using a circuitlevel simulator. Virtually all circuitlevel simulators used for ASIC design are commercial versions of the SPICE (or Spice , Simulation Program with Integrated Circuit Emphasis ) developed at UC Berkeley.
FIGURE 13.2 Output buffer (OB.IN) schematic (created using Capilano’s DesignWorks) 

13.10.1 A PSpice Example
Figure 13.2 shows the schematic for the output section of a CMOS I/O buffer driving a 10 pF output capacitor representing an offchip load. The PSpice input file that follows is called a deck (from the days of punched cards):
VIN input Ground PWL(0us 5V 10ns 5V 12ns 0V 20ns 0V)
m1 output input Ground Ground NMOS W=100u L=2u
m2 output input +5V +5V PMOS W=200u L=2u
.model nmos nmos level=2 vto=0.78 tox=400e10 nsub=8.0e15 xj=0.15e6
+ ld=0.20e6 uo=650 ucrit=0.62e5 uexp=0.125 vmax=5.1e4 neff=4.0
+ delta=1.4 rsh=37 cgso=2.95e10 cgdo=2.95e10 cj=195e6 cjsw=500e12
.model pmos pmos level=2 vto=0.8 tox=400e10 nsub=6.0e15 xj=0.05e6
+ ld=0.20e6 uo=255 ucrit=0.86e5 uexp=0.29 vmax=3.0e4 neff=2.65
+ delta=1 rsh=125 cgso=2.65e10 cgdo=2.65e10 cj=250e6 cjsw=350e12
Figure 13.3 shows the input and output waveforms as well as the current flowing in the devices.We can quickly check our circuit simulation results as follows. The total charge transferred to the 10 pF load capacitor as it charges from 0 V to 5 V is 50 pC (equal to 5 V ¥ 10 pF). This total charge should be very nearly equal to the integral of the drain current of the pullup ( p channel) transistor I _{ L} ( m2 ). We can get a quick estimate of the integral of the current by approximating the area under the waveform for id(m2) in Figure 13.3 as a triangle—half the base (about 12 ns) multiplied by the height (about 8 mA), so that
Notice that the two estimates for the transferred charge are equal.

FIGURE 13.3 Output Buffer (OB.IN). (Top) The input and output voltage waveforms. (Bottom) The current flowing in the drains of the output devices. 
Next, we can check the time derivative of the pullup current. (We can also do this by using the Probe program and requesting a plot of did(m2) ; the symbol dn represents the time derivative of quantity n for Probe. The symbol id(m2) requests Probe to plot the drain current of m2 .) The maximum derivative should be roughly equal to the maximum change of the drain current ( D I _{ L} ( m 2) = 8 mA) divided by the time taken for that change (about D t = 2 ns from Figure 13.3 ) or
The large time derivative of the device current, here 4 MAs^{ –1} , causes problems in highspeed CMOS I/O. This sharp change in current must flow in the supply leads to the chip, and through the inductance associated with the bonding wires to the chip which may be of the order of 10 nanohenrys. An electromotive force (emf ), V _{ P } , will be generated in the inductance as follows,
The result is a glitch in the power supply voltage during the buffer output transient. This is known as supply bounce or ground bounce . To limit the amount of bounce we may do one of two things:
We can work on the first solution by careful design of the packages and by using parallel bonding wires (inductors add in series, reduce in parallel).
13.10.2 SPICE Models
Table 13.14 shows the SPICE parameters for the typical 0.5 m m CMOS process (0.6 m m drawn gate length), G5, that we used in Section 2.1 . These LEVEL = 3 parameters may be used with Spice3, PSpice, and HSPICE (see also Table 2.1 and Figure 2.4 ).
TABLE 13.14 SPICE transistor model parameters ( LEVEL = 3 ). 

SPICE 
Units 2 

There are several levels of the SPICE MOSFET models, the following is a simplified overview (a huge number of confusing variations, fixes, and options have been added to these models—see Meta Software’s HSPICE User’s Manual, Vol. II, for a comprehensive description [ 1996]):
 LEVEL = 1 ( Schichman–Hodges model ) uses the simple squarelaw I _{ DS} – V _{ DS} relation we derived in Section 2.1 (Eqs. 2.9 and 2.12).
 LEVEL = 2 ( Grove–Frohman model ) uses the 3/2 power equations that result if we include the variation of threshold voltage across the channel.
 LEVEL = 3 ( empirical model ) uses empirical equations.
 The UCB BSIM1 model (~1984, PSpice LEVEL = 4 , HSPICE LEVEL = 13 ) focuses on modeling observed device data rather than on device physics. A commercial derivative (HSPICE LEVEL = 28 ) is widely used by ASIC vendors.
 The UCB BSIM2 model (~1991, the commercial derivative is HSPICE LEVEL = 39 ) improves modeling of subthreshold conduction.
 The UCB BSIM3 model (~1995, the commercial derivative is HSPICE LEVEL = 49 ) corrects potential nonphysical behavior of earlier models.
TABLE 13.15 PSpice parameters for process G5 (PSpice LEVEL = 4 ). 3 

+ VFB=0.7, LVFB=4E2, WVFB=5E2 + K1=0.78, LK1=8E4, WK1=5E2 + K2=2.7E2, LK2=5E2, WK2=3E2 + ETA=2E3, LETA=2E02, WETA=5E3 + U1=3.3E2, LU1=3E2, WU1=1E2 + X2E=4.4E4, LX2E=3E3, WX2E=9E4 + X3E=5E5, LX3E=2E3, WX3E=1E3 + X2U0=1E2, LX2U0=1E3, WX2U0=5E3 + X2U1=1E3, LX2U1=1E3, WX2U1=7E4 + X2MS=6E2, LX2MS=1, WX2MS=4 + X3U1=9E3, LX3U1=2E4, WX3U1=5E3 + CGDO=3E10, CGSO=3E10, CGBO=4E10 
+ VFB=0.2, LVFB=4E2, WVFB=0.1 + K1=0.35, LK1=7E02, WK1=0.2 + K2=4.5E2, LK2=9E3, WK2=4E2 + ETA=1E2, LETA=2E2, WETA=4E4 + X2E= 5E5, LX2E=1E3, WX2E=2E4 + X3E=8E4, LX3E=2E4, WX3E=1E3 + X2U0=9E3, LX2U0=2E3, WX2U0=2E3 + X2U1=6E4, LX2U1=5E4, WX2U1=3E4 + X3MS=1E2, LX3MS=2, WX3MS=1 + X3U1=1E3, LX3U1=5E4, WX3U1=1E3 + CGDO=2.4E10, CGSO=2.4E10, CGBO=3.8E10 
Table 13.15 shows the BSIM1 parameters (in the PSpice LEVEL = 4 format) for the G5 process. The Berkeley shortchannel IGFET model ( BSIM ) family models capacitance in terms of charge. In Sections 2.1 and 3.2 we treated the gate–drain capacitance, C _{ GD} , for example, as if it were a reciprocal capacitance , and could be written assuming there was charge associated with the gate, Q _{ G} , and the drain, Q _{ D} , as follows:
Equation 13.31 (the Meyer model ) would be true if the gate and drain formed a parallel plate capacitor and Q _{ G} = – Q _{ D} , but they do not. In general, Q _{ G} ≠ – Q _{ D} and Eq. 13.31 is not true. In an MOS transistor we have four regions of charge: Q _{ G} (gate), Q _{ D} (channel charge associated with the drain), Q _{ S} (channel charge associated with the drain), and Q _{ B} (charge in the bulk depletion region). These charges are not independent, since
We can form a 4 ¥ 4 matrix, M , whose entries are ∂ Q _{ i} /∂ V _{ j} , where V _{ j} = V _{ G} , V _{ S} , V _{ D} , and V _{ B} . Then C _{ ii} = M _{ ii} are the terminal capacitances; and C _{ ij} = – M _{ ij} , where i ≠ j , is a transcapacitance . Equation 13.32 forces the sum of each column of M to be zero. Since the charges depend on voltage differences, there are only three independent voltages ( V _{ GB} , V _{ DB} , and V _{ SB} , for example) and each row of M must sum to zero. Thus, we have nine (= 16 – 7) independent entries in the matrix M . In general, C _{ ij} is not necessarily equal to C _{ ji} . For example, using PSpice and a LEVEL = 4 BSIM model, there are nine independent partial derivatives, printed as follows:
Derivatives of gate (dQg/dVxy) and bulk (dQb/dVxy) charges
From these derivatives we may compute six nonreciprocal capacitances :
∂ Q _{ G} /∂ V _{ GB} + ∂ Q _{ G} /∂ V _{ DB} + ∂ Q _{ G} /∂ V _{ SB} 

∂ Q _{ G} /∂ V _{ GB} + ∂ Q _{ B} /∂ V _{ GB} + ∂ Q _{ D} /∂ V _{ GB} 

and three terminal capacitances:
–(∂ Q _{ G} /∂ V _{ SB} + ∂ Q _{ B} /∂ V _{ SB} + ∂ Q _{ D} /∂ V _{ SB} ) 
Nonreciprocal transistor capacitances cast a cloud over our analysis of gate capacitance in Section 3.2, but the error we made in neglecting this effect is small compared to the approximations we made in the sections that followed. Even though we now find the theoretical analysis was simplified, the conclusions in our treatment of logical effort and delay modeling are still sound. Sections 7.3 and 9.2 in the book on transistor modeling by Tsividis [ 1987] describe nonreciprocal capacitance in detail. Pages 1542 to 1544 in Vol. II of Meta Software’s HSPICE User Manual [ 1996] also gives an explanation of transcapacitance.
1. Meta Software’s HSPICE User’s Manual [ 1996], p. 1536 and pp.1613 to 1615, explains these parameters.
[ Chapter start ] [ Previous page ] [ Next page ]