Models Analog Design 9
Models Analog Design 9
Models Analog Design 9
1)
Calucalate and simulate the values of ID and VGS in the Following circuit
µF
KPn = 120
V
VTHN = .8V
KPn w
VDS = 2 − 100 k * I D ⇒ I D = * * (VGS − VTHN ) 2
2 l
VGS = VDS
VDS = 2 − 100 k * ( ID = KPn / 2 * w / l * (V DS − VTHN ) 2 )
VDS = 2 − 60 * V DS + 95 * V DS _ 36 .4
2
VDS = .934V
I D = 10 .8 µA
100K
2V
10/1
Circuit: *** Example 9.5 CMOS: Circuit Design, Layout, and Simulation ***
.control
destroy all
run
print all
.endc
.option scale=1u
.op
VDD VDD 0 DC 2
VDR VDD VDR DC 0
R1 VDR VR 100K
M1 VR VR 0 0 NMOS L=1 W=10
.MODEL NMOS NMOS LEVEL = 3
+ TOX = 200E-10 NSUB = 1E17 GAMMA = 0.5
+ PHI = 0.7 VTO = 0.8 DELTA = 3.0
+ UO = 650 ETA = 3.0E-6 THETA = 0.1
+ KP = 120E-6 VMAX = 1E5 KAPPA = 0.3
+ RSH = 0 NFS = 1E12 TPG = 1
+ XJ = 500E-9 LD = 100E-9
+ CGDO = 200E-12 CGSO = 200E-12 CGBO = 1E-10
+ CJ = 400E-6 PB = 1 MJ = 0.5
+ CJSW = 300E-12 MJSW = 0.5
*
.MODEL PMOS PMOS LEVEL = 3
+ TOX = 200E-10 NSUB = 1E17 GAMMA = 0.6
+ PHI = 0.7 VTO = -0.9 DELTA = 0.1
+ UO = 250 ETA = 0 THETA = 0.1
+ KP = 40E-6 VMAX = 5E4 KAPPA = 1
+ RSH = 0 NFS = 1E12 TPG = -1
+ XJ = 500E-9 LD = 100E-9
+ CGDO = 200E-12 CGSO = 200E-12 CGBO = 1E-10
+ CJ = 400E-6 PB = 1 MJ = 0.5
+ CJSW = 300E-12 MJSW = 0.5
.end
Problem9.2
A MOSFET with its drain connected to gate and a current flowing through it is always in
saturation. The current flowing through the MOSFET can be written as:
1 W 2
ID = * KPn * * (Vd −Vthn) ……………..(1) Also,
2 L
2 − Vd
ID = …………………………...(2)
100 K
Substituting the values of KPn =120uA/V, Vthn =0.8V(from table 9.1) and W=L=10um,
and equating equations (1) and (2), the following quadratic equation is obtained.
1 10 2 − Vd
* (120uA / V ) * * (Vd − 0.8) 2 =
2 10 100 K
R1 vdd d 100k
Vdd vdd 0 DC 2
. op
Problem 9.3:-
Calculate and simulate the values of ID and VSG in the following circuit (use long-channel
process information given in table 9.1).
Solution:-
The above PMOS is gate-drain connected. It will be in saturation if VSG ≥ VTHP.
We know for saturation, I D =
KPp W
2 L
bVSG − VTHP g
2
ID =
40 µA 10
2 V 12 b g
VSG − 0.9
2
(1)
(5 − VSG )
Also by KVL, I D = (2)
100k
Equating (1) and (2),
(5 − VSG ) 40 µA 10
100k
=
2 V 22 b g
VSG − 0.9
2
By solving the above quadratic equation, we get two values for VSG (1.33V or 0 .42v). As
the MOSFET is in saturation region, current ID flows in the circuit if VSG ≥ VTHP.
Therefore, VSG = 1.33v & I D = 36.7µA
Spice Netlist:-
**Problem#9.3
.control
destroy all
run
print all
.endc
.option scale=1u
.op
VDD VDD 0 DC 5
R1 VD 0 100K
M1 VD VD VDD VDD PMOS L=1 W=10
Simulation Results:-
DC Operating Point ... 100%
vd = 3.713491e+00
vdd = 5.000000e+00
vdd#branch = -3.71349e-05
Vsg = 1.286509e+00
vds -1.2865
vbs 0.
vth -834.8123m
vdsat -345.5119m
beta 396.4747u
gam eff 521.2973m
gm 147.9632u
gds 2.3812u
gmb 0.
cdtot 2.0237f
cgtot 13.4667f
cstot 11.2084f
cbtot 234.5729a
cgs 11.2084f
cgd 2.0237f
Problem 9.04
Since the above figure is a gate drain connected mosfet and current is flowing its
operating in saturation region.
Using KVL
5-ID.100k =VSG
5-KPP.W/2L.(VSG-VTHP)2 . 100K =VSG ;where KPP=40uA/V2, VTHP=0.9V
solving for VSG we have VSG =2.1v and –0.8v,since mosfet is in saturation therefore
VSG=2.1v
ID= KPP.W/2L.(VSG-VTHP)2
Spice simulations
**Problem 9.04
.control
destroy all
run
let vsg=vdd-vd1
print vd1
print vsg
print mag(vdd#branch)
.endc
.option scale=1u
M1 vd1 vd1 vdd vdd PMOS W=10 L=10
Vdd vdd 0 DC 5
R1 vd1 0 100k
*.dc Vdd 0 5 1m
.op
Results:
DC Operating Point ... 100%
vd1 = 2.732722e+00
vsg = 2.267278e+00
mag(vdd#branch) = 2.732722e-05
Problem 9.5:-
Calculate ID , VDS , and estimate the small-signal resistance looking into the drain of the
MOSFET in the following circuit.
Solution:-
From the above figure, the gate overdrive voltage is high and the resistor value is also
large. It is a good indication that the transistor may be in triode (because even if a small
current flows in the circuit, the voltage drop across the resistor is high). So let’s start with
the assumption that the MOSFET is in triode.
I D = KPn
W
L
MNLb g V 2
VGS − VTHN VDS − DS
2
PQO
I = 120
µA 10 LMb3 − 0.8gV − V DS
2
OP →(1)
V 2 N Q
D 2 DS
2
But, I =
D
b5 − V g →(2)
DS
200k
By solving the above quadratic equation, we get two values for VDS (4.39V or 19mv). As
the MOSFET is in triode region, VDS = 19mv and I D = 24.9 µA .
As the MOSFET is in triode region, the small signal resistance looking into the drain of
the MOSFET will be its channel resistance.
1
Rch ≈
KPn
W
L
b
VGS − VTHN g
⇒ R ch = 757.5Ω
9.5 & 9.6 Simulations
.control
destroy all
run
print all
* plot VD1#branch
.endc
* ======================================
M11 d1 g1 0 0 nmos W=10 L=2
R1 d1 Vdd 200k
Cbig d1 Vt 1
Vt Vt 0 DC 0 AC 1m
VG1 g1 0 DC 3.0
* ======================================
* .AC DEC 10 1 10K
.options scale=1u
* ABSTOL=1u VNTOL=1mv RELTOL=0.01
.op
* Level 3 models
*
* 1 um models created by RJB. These models are for educational purposes only! They are
*not*
* extracted from actual silicon.
*
* Don't forget the .options scale=1u if using an Lmin of 1
* 1<Ldrawn<200 10<Wdrawn<10000 Vdd=5V
.end
Problem 9.6
To determine the small signal resistance we can run an AC analysis in spice and then plot
the ratio of the drain voltage divided by the current through the 10mV AC source “vt”.
We sweep the frequency from 10Hz to 1MHz. The netlist for this circuit is shown below.
Problem 9.6
.opt scale=1u
.control
destroy all
run
plot test/I(Vt)
plot I(Vt)
plot I(Vin)
.endc
.ac DEC 10 100 1MEG
Vin in 0 DC=5
Vg gate 0 DC=3
Vt test 0 DC=0 AC=10m
C1 drain test 1
R1 in drain 200k
M1 drain gate 00 nmos L=2.00 W=10.00
.include C:\ 1u_models.txt
The first plot is of the vd/id and is shown in the plot below.
Looking at the plot we can see that the small signal resistance is 847 Ω.
To determine how the 200kΩ resistor effects the circuit, we can see that the 200kΩ
resistor is in parallel with the small signal resistance when we do AC analysis.
Lets look at how much current is flowing through the capacitor branch. This is shown in
the plot below.
There is about 11.8 uA flowing through the capacitor. Now lets look at how much current
is flowing through the 200kΩ branch
There is about 50nA flowing through the 200kΩ resistor. The current through the 200kΩ
resistor is much smaller than the current flowing through the capacitor. Thus we can
neglect the effect of the 200kΩ resistor to the current through the capacitor.
Problem #9.7
Figure 9.42
When Ibias is increased, VS3 and VGS3 increase. When VGS3 is increased, VGS4 also
increases because their gates are tied together. When VGS increases, the Drain to Source
Voltage decreases, as illustrated below in Figure 9.42.A. Therefore VDS4 decreases.
ID
VGS2
VGS1
VDS
Ibias2 Ibias1
Figure 9.42.A.
9.8.) Describe qualitatively what happens if we steal or inject current at the point
indicated in Fig 9.43. How does this affect the operation of M1 and M2? Verify your
answer with SPICE. (See Fig 9.43 for circuit schematic).
Since the gates of M2 and M4 are tied together and the sources of M2 and M4 are tied
together, VSG2=VSG4. Additionally, since we will assume a current of 1µA for Ibias at
the drain of M1, that current must flow through both M1 and M2, as well as M4 and M3.
Hence VSD2=VSG2=VSG4=VSD4 and VSD1=VSG1=VSG3 since all the PMOS
devices are assumed to be sized the same (we will assume W/L=30/2 for the simulations).
When current is stolen from the drain of M4, it causes the drain voltage node (here called
VD4M, which is also the source voltage of M3, VS3) to decrease which results in a
decreased VSD3 (i.e. M3 begins to lose its ability to pass current, VSG3 also goes down).
If 1 µA of current is stolen, M3 completely shuts off and all of the 1 µA flows through
M4 on that side of the circuit. (Note that the x-coordinate should be labeled in units of
µA as it is a sweep of -1µA of injected current to +1µA of injected current).
In the plot below VSG3 decreases due to VS3 decreasing when current is stolen from the
node at VD4M.
VSD3 is actually the same as VD4M since VSD3=VS3-VD3=VS3-0=VD4M
Conversely, if current is injected into the node mentioned above, it will cause the voltage
VD4M to increase to the point where the current flowing through M3 will be 2µA when
the injected current has reached 1µA (1µA provided by Iinject and 1µA by M4).
The simulation results below show what happens to VSG1 and VSG2 due to Iinject, note
that they do not change.
Similar results are obtained for VSD1 and VSD2, since current and VSG is constant,
VSD must also be constant.
The SPICE NetList for this simulation is provided below with the Long Channel
Simulation Model that was used:
let VSG1=VG2-VG1
let VSG2=VDD-VG2
let VSG3=VD4M-VG1
let VSG4=VDD-VG2
let VSD1=VG2-VG1
let VSD2=VDD-VG2
let VSD3=VD4M-V3S
let VSD4=VDD-VD4
plot vmeas#branch
plot v3test#branch
plot VD4M
plot VSG1
plot VSG2
plot VSG3
plot VSG4
plot VSD1
plot VSD2
plot VSD3
plot VSD4
plot VG1
plot VG2
*print all
.endc
.option scale=1u
*.OP
.dc Iinject -1uA 1uA 10nA
VDD VDD 0 DC 5
Ibias VG1 0 DC 1u
Iinject 0 VD4M DC 1u
VMEAS VD4 VD4M DC 0
VTEST2 VG2 VG2M DC 0
V3TEST V3S 0 DC 0
.end
Problem 9.9:-
Using simulations, generate the plot seen in Fig. 9.12 for both NMOS and PMOS devices.
Solution:-
For the MOSFETs shown in the above figure, a current source of 1uA is connected at the
drain so that the gate overdrive voltage is less. (Then VGS or VSG will be almost equal to
the threshold voltage. The plots below show the variation of the threshold voltage with
respect to the variation of source to bulk potential for both NMOS and PMOS devices.
Spice Netlists:-
***NMOS ***PMOS
.control .control
destroy all destroy all
run run
plot VGS let VSG=VDD-VG
.endc plot VSG
.option scale=1u .endc
.option scale=1u
.dc VSB 0 5 .01
.dc VBS 0 5 .01
VDD VDD 0 DC 5
VSB 0 VSB DC 0 VDD VDD 0 DC 5
IBIASP VDD VGS DC 1u VBS VBS VDD DC 0
IBIASP VG 0 DC 1u
M1 VGS VGS 0 VSB NMOS L=1 W=10 M1 VG VG VDD VBS PMOS L=1 W=10
9.10) Design a circuit that will linearly convert an input voltage that ranges from 0 to
4V into a current that ranges from roughly 50µA to 0. Simulate the operation of
the design showing the linearity of the voltage to current conversion. How does
the MOSFET’s length affect the linearity?
Solution: For voltage to current conversion we use an inactive, (resistor tied to the
source) source follower. With this configuration the source of our PMOS transistor is not
held to a fixed voltage so to reduce the source to body potential difference we tie the
body of the transistor to the source of the PMOS.
For a quick sizing of the resistor we assume that our Vthp, transistor must stay in the
saturation region, is approximately 1V which gives us a resistor value of 80k if we are to
source 50µA.
For the sizing of the PMOS let’s start with the values in the book and run a simulation to
see what our linearity is. We will plot the ID vs VSG and the derivative of ID to see the
linearity. W=30µ L=2µ
From the ID plot we see that we are almost at 50uA with the 80k resistor, a value of 74k
would have got us to 50uA but we’ll leave the 80k for comparison as we change the
sizing of the PMOS to improve linearity. From the plot of the slope of ID we see that we
have a variation of around 24%.
To improve the voltage to current linearity we’ll push out the point at which the PMOS
starts to enter the subthreshold region, we want to keep the PMOS in the saturation or
linear region as long as possible. We do this by increasing our W and decreasing our L.
We can clearly see from the log of ID plot that the linearity is almost ideal until around
3.5V, where the MOSFET starts to turn off and enter the subthreshold region. Plotting
the voltage at the source of the PMOS confirms this at the swept gate voltage of 3.5V the
source voltage is around 4.4V which is right at the expected point where the PMOS
enters the subthreshold region.
Problem 9.11:-
Using a PMOS device, discuss and show with simulations how it can be used to
implement a 10k resistor. Are there any limitations to the voltage across the PMOS
resistor? Explain.
Solution:-
The MOSFET looks like a resistor when we are using in a triode region. The equation for
the channel resistance is given as below
1
Rch = . If VSG – VTHP >>VSD , this equation can be written as,
KPp
W
L
b
VSG − VTHP − VSD g
1
Rch ≈ . Use VSG =5V so that the bias is fixed and also gives a
KPp
W
L
b
VSG − VTHP g
maximum range of VSD,sat. Given Rch=10K.
1 W
Rch = = 10k ⇒ = 0.61. If L=10, then W=61.
µA W
40 2
V L
b g
5 − 0.9 L
Limitations: - One end of the resistor should be VDD and the other end should be
VD>VTHP. Though we VD>VTHP, the simulation shows that VD>2.5V for better linearity of
the resistor.
Simulation:-
.control
destroy all
run
plot VD#branch
.endc
.option scale=1u
.dc VD 0 5 .1
VDD VDD 0 DC 5
VD VD 0 DC 0
Problem 9.12
Using SPICE (and ensuring the MOSFET is operating in the saturation region with sufficient
VDS) generate the iD vs. vGS curve seen in Fig 9.15. Using SPICE take the derivative of iD to plot
the device’s gm (versus VGS). How does the result compare to Eq. 9.22? Does the level 3 model
used in the simulation show a continuous change form sub-threshold to strong inversion?
Solution:
The first step to completing this task is to create a netlist that enables will allow us to supply a
specific VDS to the part and allow us to vary VGS. VGS will be our main x-axis variable because
gm is defined as the derivate of ID with respect to VGS. Therefore, in our netlist we will supply a
constant VDS and using a DC sweep vary VGS. VDS should be well above VDS,sat so we will
arbitrarily choose 250mV. The .DC statement varies VGS from 0 to 2.5V in 10mV increments.
Netlist
.control
destroy all
run
*** DC Analysis ***
let ID=-VDS#branch
let gm=deriv(ID)
let gmcalc=(120E-6*(10/2)*(VGS-0.8))
plot ID gm gmcalc
*plot gm
.endc
.option scale=1u
.DC VGS 0 2.5 10m
1
Chapter 9: CMOS Circuit Design, Layout, and Simulation
Next, three definitions are issued with ‘let’ statements to define ID, gm, and calculated gm. gm,calc
is merely equation 9.22 using the parameters from the level 3 NMOS model for βn and VTHN.
Figure 1 shows the results when this netlist is run.
The red line is ID. ID behaves as expected because there is almost no current until VGS > VTHN
(~0.8V) and then there is a rapid increase in the ID for an increase in VGS. The green line is gm
and also behaves as expected. gm is defined as the rate of change of id with respect to vgs. This
definition is only value for the area of the ID vs VGS curve around the VGS bias point. For this
particular case the VGS is approximately 1.05V.
The blue line is the calculated value of gm. It is obvious that the calculated value of gm correlated
very well with the simulated version for the immediate region surrounding the VGS bias point.
2
Chapter 9: CMOS Circuit Design, Layout, and Simulation
This model shows a very continuous change from the subthreshold to the strong inversion
regions. This can be seen in Figure 2 by observing ID and gm with no discontinuities below for
VGS < VTHN.
3
9.13)
VDD VDD=5V
1.3V
M1
21/3
1mv id
To find weather the MOSFET is in saturation or triode we need to check the following
two conditions
VSG > VTHP => 1.3V > 0.9V
VSD > VSG – VTHP => 5V > 1.3V – 0.9V. (Were VSG= 1.3 V; VTHP= 0.9V)
As both the conditions are satisfied the MOSFET is in Saturation.
The relation between the AC gate voltage to the AC drain current is given by
id = gm x vgs.
And the value of gm is given by
gm = Kpp x W/L (VSG – VTHP)
gm = 40µA/V x 21/3 x (1.3-0.9) = 112 µA/V.
id = gm x vgs = 112 µA/V x 1mV = 112 nA.
The AC drain current id = 112 nA Sin2Πf.
And ID = Kpp/2 x W/L (VSG – VTHP)2 = 40/2 x 21/3 x (1.3-0.9)2 = 22.4 µA.
SPICE SIMULATIONS
*** Problem 9.13 CMOS: Circuit Design, Layout, and Simulation ***
* AC ANALYSIS
.control
destroy all
run
plot -VDD#BRANCH
.endc
.option scale=1u
.ac dec 100 1 10k
VDD VDD 0 DC 5
VG1 VDD VG1 DC 1.3 AC 1m SIN 2.5 1m
M1 0 VG1 VDD VDD PMOS L=3 W=21
Simulation result:
*** Problem 9.13 CMOS: Circuit Design, Layout, and Simulation ***
VDD VDD 0 DC 5
VG1 VDD VG1 DC 1.3 AC 1m SIN 2.5 1m
M1 0 VG1 VDD VDD PMOS L=3 W=21
solving the above quadratic equation we get VSD=8.19v and 0.00813v. The former value
is not a valid value as the source voltage applied is only 5V. So taking the latter value
VD3=Vs3-VSD3=4.991V
As the gate to source voltage of M3 and M4 are same and same ID flows
VSD3= VSD4 thus VD3 = VD4=4.991V
Consider M3
VSG3=5-0=5V and VSD3=0.00813V
therefore VSG3>Vthn and VSD3<VSG3-Vthp thus our assumption that M3 is triode is correct.
Consider M1
VGS1=2.4-1.44=1.06V and VDS1=4.967V
therefore VGS1> Vthn and VDS1> VGS1- Vthn thus our assumption is correct.
Similarly M2 and M4 are operating in saturation and triode regions respectively.
AC currents
DC Operating Point
vd13 = 4.989699e+00
vd24 = 4.989699e+00
vdd = 5.000000e+00
vdd#branch = -4.00000e-05
vg1 = 2.500000e+00
vg1#branch = 0.000000e+00
vg2 = 2.500000e+00
vg2#branch = 0.000000e+00
vs12 = 1.171839e+00
.
AC voltages
from simulations mag(vd13)=74.2uV,Vg1=1mv The hand calculated values and simulations results
are almost close
9.15 To calculate the AC, DC voltages and currents in the circuit below:
To find the operating points of the circuit do the DC analysis. Hence short all the AC voltage sources.
DC Analysis:
Assume both M1 and M2 are in saturation, we will verify our assumption shortly.
Since the sources of both M1 and M2 are tied together and gates are connected to 2.5V
VGS1=VGS2
2 IL
VGS1 or VGs2 = + Vthn
KPnW
2.10.2
= + 0.8 = 0.9825V
120.10
Coming to M3 and M4, since both the gates are tied to ground and source is at VDD by intuition they
might be in triode. We will verify this too shortly.
In triode:
ID = KPp.W/L. [(VSG-Vthp).VSD3- (VSD3)2 /2)]
Since ID = 10 µA and other parameters are known except VSD3, solve for VSD3
30 Vsd 2
10 = 40. ((5 − 0.9)Vsd − )
2 2
Solving for VSD3 we get VSD3 = 0.00406 V or 8.195 V. The later value for VSD3 doesnot make any sense.
∴ VD3 = 5 -VSD3=4.995V
Since VDS > VGS - Vthn ,our assumption that M1 is in saturation is correct.
Now for M3:
VSD = 5 - 4.995 = 5mV
VSG -Vthp = 5 - 0.90 = 4.1V
Since VSD < VSG - Vthp ,our assumption that M3 is in triode is verified.
Since M3 and M4 are in triode they behave as a resistors whose resistance is given by:
W
gm1 = gm1 = KPn (Vgs − Vthn )
L
= 120.10/2(0.9825-0.8) = 109.5 µA/V
Having calculated all the DC parameters lets check our calculations with simulations:
The discrepancy in vs12 is because of the fact that we did not take body effect into account.
AC Analysis:
vgs1-vgs2=2mV
id
and since = vgs it follows
gm
id 1 id 2
− = 2mV
gm1 gm 2
and we know that id1= -id2 = id3= -id4 and gm1=gm2
Therefore it follows that vgs1= -vgs2 = 1mV
iD1 = 10u+0.1095sin2πf
iD2 = 10u-0.1095sin2πf
AC input voltage
Drain voltages of
M1,M2,M3,M4
SPICE NETLIST:
.control
destroy all
run
** for the operating point analysis
*print all
*
** for the AC analysis
plot mag(vd13) mag(vg1)
.endc
.option scale=1u
*.op
.ac dec 100 1 10k
VDD VDD 0 DC 5
VG1 VG1 0 DC 2.5 AC 1m SIN 2.5 1m 10k
VG2 VG2 0 DC 2.5 AC -1m SIN 2.5 1m 10k
Ibias VS12 0 DC 20u
*Prob 9.16*
.control
destroy all
run
v1 vd13 va1 0v
v2 vd24 va2 0v
vdd vdd 0 dc 5
vg1 vg1 0 dc 2.5 ac 1m sin 2.5 1m 10k
vg2 vg2 0 dc 2.5
ibias vs12 0 dc 40u
*Level 3 models
.MODEL NMOS NMOS LEVEL = 3
-----------------------------
.MODEL PMOS PMOS LEVEL = 3
--------------------------------
.end
AC Analysis
Transient Analysis
Problem 9.17
When the potential of the source of a MOSFET is increasing, a point is reached where the gate to
source potential goes below the threshold voltage. The MOSFET then shuts off. That’s the point
where the graph cuts the x-axis. When the source voltage (or the source to bulk potential, since
the body is grounded) is 0, the value of ID is where the graph cuts the y-axis.
SPICE Netlist:
*** PROBLEM 9.17 CMOS: Circuit Design, Layout, and Simulation ***
.control
destroy all
run
LET ID = -VMETER#branch
PLOT ID
.endc
.option scale=1u
VDD VDD 0 DC 5
VG VG 0 DC 5
VSB VSB 0 DC 5 AC 1m
VMETER VD VDD DC 0
.MODEL NMOS NMOS LEVEL = 3
+ TOX = 200E-10 NSUB = 1E17 GAMMA = 0.5
+ PHI = 0.7 VTO = 0.8 DELTA = 3.0
+ UO = 650 ETA = 3.0E-6 THETA = 0.1
+ KP = 120E-6 VMAX = 1E5 KAPPA = 0.3
+ RSH = 0 NFS = 1E12 TPG = 1
+ XJ = 500E-9 LD = 100E-9
+ CGDO = 200E-12 CGSO = 200E-12 CGBO = 1E-10
+ CJ = 400E-6 PB = 1 MJ = 0.5
+ CJSW = 300E-12 MJSW = 0.5
*
.MODEL PMOS PMOS LEVEL = 3
+ TOX = 200E-10 NSUB = 1E17 GAMMA = 0.6
+ PHI = 0.7 VTO = -0.9 DELTA = 0.1
+ UO = 250 ETA = 0 THETA = 0.1
+ KP = 40E-6 VMAX = 5E4 KAPPA = 1
+ RSH = 0 NFS = 1E12 TPG = -1
+ XJ = 500E-9 LD = 100E-9
+ CGDO = 200E-12 CGSO = 200E-12 CGBO = 1E-10
+ CJ = 400E-6 PB = 1 MJ = 0.5
+ CJSW = 300E-12 MJSW = 0.5
.end
Problem 9.18
For n-MOSFET W/L=100/20, p-MOSFET W/L=300/20
For a long channel n-MOSFET in saturation the equation for drain current is
KPn W
ID = . .(VGS −VTHP )2
2 L
when W and L are both multiplied by 10 as seen in the equation the drain current does not
change.
2I D L 2.20 20
VSG = . +VTHP = . +0.9 = 1.158 V
KPP W 40 300
The VGS (NMOS) & VSG (PMOS) of the MOSFET’s are 1.058 V and 1.158 V respectively they
remain the same as in Ex. 9.5 as the drain current (roughly 20µ bias current) is almost the same.
Both the MOSFET’s have a VDS.SAT of 250mV
Figure shows the IV plots and the output resistance (the reciprocal of the derivative of the drain
current) for the MOSFET’s.
The output resistances (from the plots) of NMOS & PMOS are 500 MEG & 300 MEG respectively
as against 5 MEG & 4 MEG in Ex.9.6.
As λα 1/L, r0 α L2/V2 DS,SAT
Since L increases by 10 the output resistance increased by a factor of 100 and λ decreases.
Netlist
*** Problem 9.18
NMOS VDD
.control
destroy all
run VGN 100/20
let ID=-VDD#branch M1
plot ID
VGN=1.05 V
let ro=1/deriv(ID)
plot ro
.endc
.option scale=1u
.dc VDD 0 5 1m
VDD VDD 0 DC 5
VGN VGN 0 DC 1.05
.option scale=1u
.dc VDD 0 5 1m
VDD VDD 0 DC 5
VGP VDD VGP DC 1.15
M1 0 VGP VDD VDD PMOS L=20 W=300
NMOS NMOS
PMOS PMOS
Problem 9.19
ft1/ft2=( L2 /L1)2
ft1/ft2=( 20 /2 )2 = 100
ft2 = ft1/100
From this relationship we can see that ft depends only on channel length. In other words
ft decreases with channel length increase (quadratic dependence on channel length
relationships).
To illustrate change, WinSpice simulations were performed and obtained following
graphs:
From figures 1 and 2 we can see that ft1 = 2 x 109 Hz (ft1 intersect 0-axes) and
ft2 = 2 x 107 Hz (ft2 intersect 0-axes) approximately. So,
This problem asks us to estimate the threshold voltage of a PMOS and a NMOS transistor which
are fabricated using a short channel process. The dimensions of both MOSFETS are 50/5.
We will start with the NMOS. The threshold voltage may be estimated by examining two
different plots. First, if we plot ID versus VGS, the threshold voltage may be estimated by linearly
extrapolating back to the x-axis. As seen in figure 1, the extrapolation yields an approximate
value of 280mV for VTHN.
The second method of threshold estimation involves plotting the derivative of ID versus VGS.
Once again, we will linearly extrapolate back to the x-axis to find the approximate value of VTHN.
Figure 2 shows that the threshold voltage is approximately 210mV. This measurement is
considered more accurate than the first. For more explanation, please see Figure 9.27 and the
accompanying paragraph in chapter 9 of the text.
*** Problem 9.21, from Figure 9.27 CMOS: Circuit Design, Layout, and Simulation ***
.control
destroy all
run
let ID=-VDS#branch
let gm=deriv(ID)
plot gm
plot ID
.endc
.option scale=50n
.DC VGS 0 1 1m
VDS VDS 0 DC .1
VGS VGS 0 DC 0
The situation is the same for the PMOS. Figure 4 shows that VTHP is about 250mV. Figure 5
shows the VTHP is about 190mV. Finally the netlist used for the PMOS simulations is shown in
figure 6.
Figure 4: ID verses VSG, showing the Threshold Voltage at 250mV
Figure 5: The Derivative of ID (gm) verses VSG, showing the Threshold Voltage at ~190mV
*** Problem 9.21, from Figure 9.27 (PMOS) CMOS: Circuit Design, Layout, and Simulation **
.control
destroy all
run
let ID=-VSD#branch
let gm=deriv(ID)
plot gm
plot ID
.endc
.option scale=50n
.DC VSG 0 1 1m
VSD 0 VSD DC .1
VSG 0 VSG DC 0
Show the details leading to Eq. (9.43). Show, as an example, that the approximation is
'
valid if Q bo ' = 30mv . (Remember: temperature is in Kelvin.)
C ox
Solution:-
From chapter 6, Eq.(6.17),
Q ' bo − Q ' ss
VTHN = −Vms − 2V fp + , where Q ' ss is a constant.
C ' ox
kT N A kT N D , poly
Q ' bo = 2qN Aε si / − 2V fp / , V fp = − ln and Vms = ln − V fp
q ni q ni
∂VT ∂ kT k
N D , poly = 10 20 , N A = 1015 , = = = 0.085mv / K .
∂T ∂T q q
∂VTHN k 10 20
Plugging in the numbers we get = − ln 15 = −0.9786 mV / K .
∂T q 10
Q ' bo
With the contributions from = 30mv included at T=300K we get
C ' ox
∂VTHN k N D , poly Q 'bo
= − ln + = -0.9786 + 0.05 =-0.9286 mV/K.
∂T q NA 2.T .C ' ox
PROBLEM 9.23
(1) The NMOS is in saturation as a constant current is flowing through a gate- drain
connected MOSFET.The expression for drain current for the short channel process is:
Substituting the values for υsatn , W, Cox ' and VDSsat from table9.2.
10*10-6=(9*109um/s)*(100um)*(25fF/um2)*(1-VG-0.28-0.05)
Deriving the equation for the variation of VSG with temperature for PMOS just like
derived for NMOS results in the following equation:
*Ibias d 0 DC 10u
Ibiasn vdd d DC 10u
*Vdd s 0 DC 1
Vddn vdd 0 DC 1
.op
9.24. Prove using simulation outputs that the GFT product is constant in the saturation
region for different biasing conditions on a long channel process. Since GFT=gm*r0*fT,
we will generate these three parameter for three different biasing using our 1u models.
The three biasing currents will be set by forcing VGS to 1.5V, 2V, and 2.5V. The
following file was used to simulate gm and ID with Winspice:
vdd vdd 0 dc 5V
vg vg 0 dc 1V
m1 vdd vg 0 0 nmos w=10 l=2
.option scale=1u
** Include Model File Here ****
.end
vdd vdd 0 dc 3V
vg vg 0 dc 1.5V
m1 vdd vg 0 0 nmos w=10 l=2
.options scale=1um
The plot for ID versus VDS at VGS=1.5V and r0 at VGS=1.5V, 2V, and 2.5V are shown
below:
Therefore: r0 @ VDS=5V and VGS=1.5V, 2V, 2.5V = 980K, 370K, and 200K since R1
is r0 at VGS=1.5V, R2 is r0 at VGS=2V, and R3 is r0 at VGS=2.5V. IDSSAT at
VGS=1.5V, 2V, 2.5V = 110uA, 280uA, and 500uA. Since VDSSAT = VGS – VTHN,
VDSSAT=0.7V at VGS=1.5V, VDSSAT=1.2V at VGS=2V, and VDSSAT=1.7V at
VGS=2.5V.
To calculate the GFT product, we need to find the fT at our three bias points. FT is found
by running an AC analysis of a transistor in the saturation region with the AC signal
applied to the gate terminal for the three VGS bias values. The value of fT is the
frequency where |id| / |ig| = 0db. Plotting the logrithum of |id| / |ig| versus frequency and
fT is the frequency where this quanity = 0. The following file was used to calculate fT for
each bias condition and must be simulated three times with the different vg voltages:
vdd vdd 0 dc 5V
vg vg 0 dc 1.5V ac 1mV
m1 vdd vg 0 0 nmos w=10 l=2
.options scale=1um
Therefore: fT @ VGS=1.5V, 2V, 2.5V = 1.9GHz, 2.7GHz, and 3.2GHz since FT1 is fT
at VGS=1.5V, FT2 is fT at VGS=2V, and FT3 is fT at VGS=2.5V.
run
Let ft=mag(VDS#branch)/mag(VGS#branch)
plot 20*log(ft)
.endc
.option scale=50n
.ac dec 100 100M 10G
VDS VDS 0 DC 1
VGS VGS 0 DC 350m AC 1m
M1 VDS VGS 0 0 NMOS L=2 W=50
100/2
M3
50/2 50/2
M1 M2
1mV
+ 0.5V
0.5V + -
-
20 uA
DC operating points:
The gates of both M1 and M2 are at 0.5v and the sources of M1 and M2 are tied together VGS1=VGS2 and
ID1 =ID2=10uA
VGS1=VGS2=350mv
ID
VSG 3 = + VTHP + VSDsat
v sat .C OX '.W
Substituting the values from the Table 9.2 for short channel devices we get
VSG3=350mv
VD1=650mv
Ac Analysis:
id1 id 2
1mv= vgs1 - vgs2 = −
gm1 gm2
Here gm1=gm2 and id1=-id2
Therefore vgs1=-vgs2=0.5mv
AC drain voltages of M2
vd2=id2. ro =12.45mv
Since vD2 is 1+ac voltage there fore vD2 will stay at vdd since it can’t go beyond vdd.
Simulations
*** Problem 9.26 CMOS: Circuit Design, Layout, and Simulation ***
.control
destroy all
run
** for the operating point analysis
*print all
*let vgs1=vg1-vs12
*let vgs2=vg2-vs12
*print vgs1 vgs2 vd13
let vd2=vdd
** for the AC analysis
*plot mag(vd13) mag(vdd)
** for the transient analysis
plot vd13
plot vd2
.endc
.option scale=50n
*.op
*.ac dec 100 1 10k
.tran 1u 300u
VDD VDD 0 DC 1
VG1 VG1 0 DC 0.5 AC 1m SIN 2.5 1m 10k
VG2 VG2 0 Dc 0.5
Ibias VS12 0 DC 20u
M1 VD13 VG1 VS12 VS12 NMOS L=2 W=50
M2 VDD VG2 VS12 VS12 NMOS L=2 W=50
M3 VD13 VD13 VDD VDD PMOS L=2 W=100
DC Operating Point
vgs1 = 3.434905e-01
vgs2 = 3.434905e-01
vd13 = 6.525805e-01
Problem 9.27
The thermal noise of a MOSFET’s drain current when operating in saturation is given by Equation
(9.63).
4 KT
IR (f) =
2
3 1
*
2 gm
3 1
where the resistance is *
2 gm
When the MOSFET is in DEEP TRIODE region, resistance is given by
1
Rch = from equation (9.16)
gm
In between the deep triode and saturation regions, Thermal noise can be modeled with
1
γ×
gm
where 1 < γ < 3/2.
Problem 9.28
If we look at equation 9.66 we see that the total PSD of the mosfet drain noise current is
proportional to the gm
IM2 ∝ gm
It will also decrease the signal to noise ratio (SNR). We would like to have as large of a
signal to noise ratio as possible. Lets take a closer look.
For the current through the transistor, let us take its square so we have;
. ids2 = gm2*vgs2
gm 2 • v gs
2
SNR =
KF • I D
AF
8kT
+ • gm
f • Cox • LW
2
3
V gs
SNR =
KF •
AF
ID 8kT
+
3 • gm
f •COX • LWgm 2
Since the gm2 is proportional to ID then the first term of the denominator wont change
much as gm decreases. When we decrease gm the second term of the denominator will
increase which will cause the SNR to decrease. For an amplifier we would like to have a
large SNR and decreasing gm will decrease the SNR so although at first it appears to be a
good way to decrease the noise, for an amplifier it negatively affects the SNR and is not a
good idea.