3410 Lecture Notes v1.0
3410 Lecture Notes v1.0
MICROELECTRONICS I
LECTURE NOTES
Spring 2018
Chris Winstead
Associate Professor
Electrical and Computer Engineering
[email protected]
Copyright © 2018
http://www.ece.usu.edu
List of Netlists . . . . . . . . . . . . . . . . . . . . . . . . . 5
List of Examples . . . . . . . . . . . . . . . . . . . . . . . . 7
List of EveryCircuit Demos . . . . . . . . . . . . . . . . . 10
Introduction 13
Signal sources . . . . . . . . . . . . . . . . . . . . . . . . . 13
Ideal Amplifier Models . . . . . . . . . . . . . . . . . . . . 17
Real Amplifiers . . . . . . . . . . . . . . . . . . . . . . . . 21
Equivalent small-signal resistance and impedance . . . . 26
Frequency Response of Amplifiers . . . . . . . . . . . . . 26
Harmonic distortion . . . . . . . . . . . . . . . . . . . . . 33
Introduction to Diodes 63
Ideal switch model . . . . . . . . . . . . . . . . . . . . . . 63
Exponential model . . . . . . . . . . . . . . . . . . . . . . 65
Constant voltage-drop model . . . . . . . . . . . . . . . . 65
Iterative Analysis . . . . . . . . . . . . . . . . . . . . . . . 66
Linearized Model . . . . . . . . . . . . . . . . . . . . . . . 69
Diode Circuits 71
Half-Wave Rectifier . . . . . . . . . . . . . . . . . . . . . . 71
Resistor-diode regulator . . . . . . . . . . . . . . . . . . . 73
Peak rectifier . . . . . . . . . . . . . . . . . . . . . . . . . . 74
Envelope detector . . . . . . . . . . . . . . . . . . . . . . . 75
Bridge Rectifier . . . . . . . . . . . . . . . . . . . . . . . . 77
4
Voltage Regulators . . . . . . . . . . . . . . . . . . . . . . 79
Super Diode, Precision Rectifier . . . . . . . . . . . . . . . 83
DC Restoration, Clamped Capacitor . . . . . . . . . . . . 85
Boost converter . . . . . . . . . . . . . . . . . . . . . . . . 86
Memristors 89
Axiomatic Circuit Theory . . . . . . . . . . . . . . . . . . 89
Simulating Memristors . . . . . . . . . . . . . . . . . . . . 94
Memristor Applications . . . . . . . . . . . . . . . . . . . 99
Exploring Memristor Controversies . . . . . . . . . . . . . 105
Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
Bibliography 225
List of Netlists
1 envelope_detector.sp . . . . . . . . . . . . . . . . . . 76
2 bridge_rectifier.sp . . . . . . . . . . . . . . . . . . . . 78
3 basic_regulator.sp . . . . . . . . . . . . . . . . . . . . 81
4 superdiode.sp . . . . . . . . . . . . . . . . . . . . . . 84
5 741.sp (top lines showing port order) . . . . . . . . 84
6 dc_restorer.sp . . . . . . . . . . . . . . . . . . . . . . 86
7 boost_converter.sp . . . . . . . . . . . . . . . . . . . 87
8 DC sweep of common-source degeneration resistances138
9 AC simulation of CS configuration . . . . . . . . . . 156
10 Gain/BW tradeoff in CS configuration . . . . . . . . 157
11 diode_dc.sp . . . . . . . . . . . . . . . . . . . . . . . 199
12 tunnel_diode_model.sp . . . . . . . . . . . . . . . . 208
13 tunnel_diode_oscillator.sp . . . . . . . . . . . . . . . 209
List of Examples
1 Low-pass RC circuit . . . . . . . . . . . . . . . . . . . . . 15
2 High-pass RC circuit . . . . . . . . . . . . . . . . . . . . . 15
3 Linearization of a sensor . . . . . . . . . . . . . . . . . . . 23
4 Linearization of a thermistor model . . . . . . . . . . . . 24
5 Bias current in inverting configuration . . . . . . . . . . . 46
6 Maximum resistance due to I bias . . . . . . . . . . . . . . 46
7 Inverting configuration with offset voltage . . . . . . . . 48
8 Closed-loop frequency response, low-gain . . . . . . . . 52
9 Closed-loop frequency response, high-gain . . . . . . . . 52
11 Max-Value Circuit . . . . . . . . . . . . . . . . . . . . . . . 64
12 Min-Value Circuit . . . . . . . . . . . . . . . . . . . . . . . 64
13 Min-value circuit with 0.7V drop model . . . . . . . . . . 65
14 Iterative analysis . . . . . . . . . . . . . . . . . . . . . . . 66
15 Iteration with linearized model . . . . . . . . . . . . . . . 70
16 Half-wave rectifier with vin < 0 . . . . . . . . . . . . . . . 71
17 Half-wave rectifier with vin = 1V . . . . . . . . . . . . . . 71
18 Four-diode voltage regulator design . . . . . . . . . . . . 79
19 Two Diode Regulator with Op Amp Buffer . . . . . . . . 82
20 Diode as a nonlinear resistor . . . . . . . . . . . . . . . . 92
21 Static CMOS logic design . . . . . . . . . . . . . . . . . . 125
22 Static CMOS XOR gate . . . . . . . . . . . . . . . . . . . . 126
23 Passive-biased CS amp with source degeneration . . . . 135
24 Common-gate configurations . . . . . . . . . . . . . . . . 142
25 Source Follower output resistance . . . . . . . . . . . . . 145
26 Current-mirror active bias . . . . . . . . . . . . . . . . . . 147
27 Output offset with current-mirror bias . . . . . . . . . . . 149
28 Voltage divider bias . . . . . . . . . . . . . . . . . . . . . . 161
29 Feedback bias . . . . . . . . . . . . . . . . . . . . . . . . . 162
30 Charge density in doped semiconductors . . . . . . . . . 172
31 Fermi levels in doped silicon . . . . . . . . . . . . . . . . 174
32 Drift current in a homogeneous material. . . . . . . . . . 177
33 Drift and diffusion currents . . . . . . . . . . . . . . . . . 180
34 Junction capacitance due to depletion . . . . . . . . . . . 188
36 Zener temperature coefficients. . . . . . . . . . . . . . . . 191
10
Rsig
Signal sources +
+
vsig −
Electronic circuits and systems can be loosely divided into two
classes: those that process signals, which are electrical repre- −
sentations of information, and those that convey or convert Figure 1: Thévenin equivalent voltage signal source.
electrical power. In this course we are primarily concerned with
circuits that process signals, in the broadest possible sense. A
signal may convey physical information, e.g. an audio signal +
information, either a current or voltage. Figure 2: Norton equivalent current signal source.
A transducer is a device that converts a physical signal into
an electrical one. From the circuit perspective, we usually model1
a transducer as either a voltage source or a current source. 1
A model is a useful approximation of a physical
Every signal source has an associated internal impedance, device or system. We use models to simplify our
understanding of complex electronic components.
represented by Thévenin or Norton equivalent circuits.
| X ( jω )|2 = X ( jω ) × X ∗ ( jω ).
VS (s) (dB)
in degrees (0◦ to 360◦ ).
The individual sinusoidal signal components are expressed 20
as
v a (t) = VA sin (ω0 t + φ) ,
0
where VA is the zero-to-peak amplitude, ωs is the signal fre- 101 102 103 104 105
quency in radians per second, t is the time in seconds, and φ is ω (rad/sec)
the phase-shift in radians. Figure 4: A sinusoid has a single Fourier component
A pure or single-tone sinusoid has a magnitude spectrum that appears as an impulse function on the spectral
representation. In this case the magnitude is 40 dB,
represented by a single impulse function which corresponds to a time-domain zero-to-peak
amplitude of 200 V.
1
VA (ω ) = V δ ( ω − ωs ) .
2 A
The impulse height in decibels is 20 log10 (VA /2). So a zero-to-
peak magnitude of 1 V corresponds to −6 dB, 10 V corresponds
to 14 dB, 100 V corresponds to 34 dB, and so on. Frequency Units:
When signals pass through electronic circuits, their magni- • ω = 2π f
tude and phase are altered. The circuit’s transfer function is the • f is in Hz (cycles per second)
ratio of the output spectrum to the input spectrum. If a circuit’s • ω is in radians per second.
input is X (s) and its output is Y (s), then the transfer function is • ω is “omega”, not w.
• A magnitude V in dB is 20 log10 (V )
Y (s)
H (s) = . • A power P in dB is 10 log10 ( P).
X (s)
| H (ω )| (dB) = 20 log10 | H (ω )|
= 10 log10 | H (ω )|2 .
tude is less than the input, and the circuit is said to attenuate
the signal.
For a linear circuit, the transfer function is obtained using
complex impedances for capacitors and inductors. A capacitor L Ls
with capacitance C has impedance 1/(sC ), and an inductor with
Figure 5: Passive linear components and their
equivalent Laplace-domain impedances.
introduction 15
1
H ( jω ) =
1 + jω/ω3dB
At low frequencies where ω ω3dB , the magnitude response is flat, approximately equal to one. At
higher frequencies where ω ω3dB , the magnitude drops rapidly. At these high frequencies, since
ω/ω3dB 1, we can make an approximation:
ω/ω3dB + 1 ≈ ω/ω3dB
⇒ | H (ω )| ≈ ω3dB /ω (high frequencies above ω3dB )
jω/ω3dB
H ( jω ) =
1 + jω/ω3dB
− jω/ω3dB
jω/ω3dB
| H (ω )|2 =
1 + jω/ω3dB 1 − jω/ω3dB
!
ω 2 /ω3dB
2
= 2
1 + ω 2 /ω3dB
At high frequencies where ω ω3dB , the magnitude response is flat, approximately equal to one. At
lower frequencies where ω ω3dB , the magnitude drops rapidly. Since ω/ω3dB 1, we can make
an approximation:
ω/ω3dB + 1 ≈ 1
⇒ | H (ω )| ≈ ω/ω3dB (low frequencies below ω3dB )
So as ω increases from very low frequencies, the magnitude increases by 20 dB per decade.
Rout
Ideal Amplifier Models + +
+
An ideal linear amplifier is a circuit which receives an input vin Rin Av vin − vout
signal X and produces an output signal Y = AX. In other
− −
words, the output is larger than the input by a constant multiple
Figure 8: Ideal linear voltage amplifier model at low
A, called the gain. The input/output signals can be either or mid-band frequencies.
current or voltage, which introduces four possible amplifier
configurations:
Rsig Rout
+ +
Rin
vIN = vSIG
+ Rin + Rsig
vsig +
− vin Rin Av vin − vout RL
RL
vOUT = Av vIN
Rout + R L
− −
vOUT
open-circuit gain: Avo , = Av
vIN
v Rin RL
loaded gain: AvL , OUT = Av .
vIN Rin + Rsig Rout + R L
This demonstration implements the ideal voltage amplifier model from Figure 9 with Rsig , Rin , Rout and
R L all equal to 1 kΩ, and a voltage gain Av = 10 V/V. The simulation traces show an attenuation by
half at each signal port due to the voltage-divider couplings.
Exercise: Increase Rin to 10 kΩ and then 100 kΩ, and observe what happens to the amplitude of vin
compared to vsig for these values. Then do the same for R L . You should notice that the coupling
effects disappear when R L Rout and Rin Rsig . Verify that your observations match the value of the
loaded gain predicted by our analysis in this section.
introduction 19
iin iout
The ideal linear current amplifier is very similar to the
voltage amplifier, except that we get current dividers instead of
voltage dividers at the input and output terminals. In a current vsig Rin Ai iin Rout load
divider, the opposite resistance appears in the numerator, so the
conditions for achieving maximum are reversed.
For current amplifiers, the most ideal gain is called the short-
circuit gain Ais , since we can eliminate the coupling ratios by Figure 10: Ideal linear current amplifier model at low
setting R L to zero, hence short-circuiting the output. The gain or mid-band frequencies.
expressions for a current amplifier are:
iOUT
short-circuit gain: Ais , = Ai
iIN
i Rsig Rout
loaded gain: AiL , OUT = Ai .
iIN Rin + Rsig Rout + R L
Maximum gain in current amp: Rout R L and
To maximize the current amplifier’s gain, we want to elimi- Rin Rsig . In the limit, a truly ideal current amp has
nate the coupling ratios by making them very close to one. This Rin → 0 and Rout → ∞.
iin iout
Rsig
iIN = vSIG
Rin + Rsig
isig Rsig Rin Ai iin Rout RL
Rout
iOUT = Ai iIN
Rout + R L
Rsig iout
+
Rin
vIN = vSIG
Rin + Rsig
vsig +
− vin Rin Gm vin Rout RL
Rout
iOUT = Gm vIN
Rout + R L
−
iin Rout
+
Rsig
iIN = vSIG
+ Rin + Rsig
isig Rsig Rin Rm iin − vout RL
RL
vOUT = Rm iIN
Rout + R L
−
Input Input
−15 −10 −5 5 10 15 −15 −10 −5 5 10 15
−100 −100
A temperature sensor provides a change of 2mV per ◦ C, connected to a load of 10kΩ. The output
changes by 10mV when T is changed by 10◦C. What is the source resistance of the sensor?
The sensor model is linearized:
dvS RS
vs = VS + ∆T +
dT T0
R (Ω)
1
d − B T1 − T1
R ≈ R0 + ∆T
R0 e 0
dT
T =T
0
1
−B T − T 1 d 1 1 0.99
= R0 + ∆T R0 e
0 −B −
dT T0 T
T = T0
! 280 290 300 310 320
B
= R0 − ∆TR0 T (Kelvin)
T02
Figure 21: Linearized approximation of thermistor
resistance for temperatures near 300 K
So if T0 = 300 K , R0 = 10 kΩ and B = 50 K−1 , then for
temperatures near 300 K we have
R ≈ 10 kΩ − ∆T × 5.5 Ω
Rsig Cf
Rout
+ +
+
vsig +
− vin Cin Rin Av vin − Cout vout RL
− −
∏k (1 − s/ωzk )
H (s) = K
∏m 1 − s/ω pm
where ωzk are the zeros, indexed by k and ω pm are the poles, + +
indexed by m. In this course we will concern ourselves almost
VY (s) vy
exclusively with “simple” poles and zeros in the left half-plane. vX H (s) = VX (s)
In other words, we’ll assume that all poles and zeros are real-
valued (not complex or imaginary), are all well separated (they − −
do not overlap in value), and are negative valued. If these
conditions are satisfied, then we can use a simplified “stick- Figure 23: General “black-box” model of a linearized
amplifier circuit. Zeros are roots of s in the numera-
figure” method to produce approximate magnitude and phase tor and poles are roots in the denominator.
response diagrams, which are called Bode plots.
28 ece3410 lecture notes
Low-pass systems
For every pole ω pm , the magnitude decreases by 20 dB per
decade at frequencies above the pole. The phase response
decreases by 90° between the frequencies 0.1ω pm and 10ω pm ,
and crosses −45° at ω pm .
100
ω p0
Gain Magnitude (dB)
-20dB/decade
50
100 101 102 103 104 105 106 107 108 109
0
45◦ Phase Loss
Phase (°)
−50
−100
100 101 102 103 104 105 106 107 108 109
Frequency
High-pass systems
For every zero ωzk , the transfer function increases by 20 dB
per decade at frequencies above the pole. The phase response
increases by 90° between the frequencies 0.1ω pm and 10ω pm ,
and crosses −45° at ω pm .
In high-pass systems, there is usually a zero at the origin.
In that case, there is no phase response associated with the
zero (it occurs at infinitely low frequency on the logarithmic
scale), and the zero must be canceled by one or more poles at
higher frequencies. On the Bode plot, the magnitude response
reaches its maximum value and becomes constant after the first
pole, ω p0 . For frequencies below ω p0 , the magnitude response
decreases by 20 dB per decade.
100
ω p0
+20dB/decade
Gain Magnitude (dB)
50
100 101 102 103 104 105 106 107 108 109
0
45◦ Phase Loss
Phase (°)
−50
−100
100 101 102 103 104 105 106 107 108 109
Frequency
Band-pass systems
− −
vout
H (s) =
vsig
RL sCC1 Rin
= Av
R L + Rout (1 + sCC1 ( Rin + Rsig )) (1 + sCout ( R L k Rout ))
In this case there is a zero at the origin and two poles located at
introduction 31
two frequencies:
1
ω p0 =
CC1 ( Rin + Rsig )
1
ω p1 =
Cout ( R L k Rout )
100
ω p0 ω p1
Gain Magnitude (dB)
+20dB/decade -20dB/decade
50
midband
100 101 102 103 104 105 106 107 108 109
0
45◦ Phase Loss
→ 90◦ Phase Loss
−50
Phase (°)
−100
135◦ Phase Loss
−150
→ 180◦ Phase Loss
100 101 102 103 104 105 106 107 108 109
Frequency
This demonstration shows an implementation of the band-pass model from Figure 26. Examine this
circuit and perform both AC simulations (frequency mode) and transient simulations (time mode).
Try increasing and decreasing both CC1 and Cout by 10× (for a total of four different cases), and
observe how the pole frequencies change. Verify that the observations match the predictions from our
analysis in this section.
This capacitive coupling demonstration shows how we can remove a signal’s DC offset and replace
it with a different offset. The circuit works through superposition of high-pass and low-pass signal
paths. The input AC signal has an offset of 10 V and a zero-to-peak amplitude of 1 V. At the out-
put, the original offset is rejected by the coupling capacitor. A new offset of 1 V is provided by an
independent DC voltage source, and is superimposed through the 1 kΩ resistor on the output side.
introduction 33
Harmonic distortion
40
When a pure sinusoid is input to perfectly linear amplifier, the
VS (s) (dB)
output is expected to be a pure sinusoid, and its magnitude
spectrum should have a single impulse. Real amplifiers are not
20
perfectly linear though, so the output is usually not a perfect
sinusoid.
As a result, unexpected features called harmonics appear
0
in the output magnitude spectrum. Harmonics are spuri- 103 104
ous impulses that appear at integer multiples of the original ω (rad/sec)
fundamental signal frequency. So if the input sinusoid has a Figure 28: Harmonic “spurs” appear at integer mul-
fundamental frequency component at f 0 , the distorted output tiples of the fundamental frequency, and represent
distortion.
sinusoid has harmonic components spaced at integer multiples
fk = k f0 .
Aliasing
Since the harmonic components can extend to very high fre-
quencies, they may contribute to aliasing effects in a digital
oscilliscope’s FFT display. Aliasing occurs when a signal vio-
lates the Shannon-Nyqvist Sampling Theorem, which states that
the sampling rate must be at least twice the highest frequency
present in the signal. On a typical digital oscilloscope, we must
be aware of the following considerations:
• If the signal frequency f > f S /2, then the scope will show an
image at f − f S /2. So if you increase f beyond f S /2, the signal
peak on the FFT will appear to move backwards.
v−
− −
Amplifiers with finite open-loop gain
vid vOUT
Inverting amplifier
vIN
i1 = .
R1
gain is
R2
Gi? = − .
R1
vout = A v+ − v−
⇒ vout = A 0 − v−
vout
⇒ v− = −
A
vout − v−
i2 = i1 =
R2
−
v − vin
=
R1
Then we have
v
1 out
R1 vout 1 + = R2 − − vin
A A
1 R2 R2
⇒ vout 1 + + = − vin
A R1 A R1
vout R2 A
⇒ Gi = = −
vin R1 A + 1 + R2 /R1
Gi = Gi? × e
R2
Gi? = −
R1
A
e=
A + 1 + R2 /R1
This circuit implements an inverting configuration where the op amp’s open-loop gain is A = 10 V/V
(i.e. 20 dB). The resistor values are R1 = 1 kΩ and R1 = 2 kΩ, so we expect an ideal closed-loop gain
of G ? = −2 V/V. The input signal has a zero-to-peak amplitude of 1 V, so the output amplitude
should be 2 V. Simulate this circuit and observe the output amplitude. It should be 1.54 V. To verify
that this matches the prediction from our theory, solve for e and G using the methods described in
this section. Then, try increasing the op amp’s open-loop gain to 20 V/V and repeat your calculations
to verify that the theory holds up.
38 ece3410 lecture notes
Non-inverting amplifier R2
+
vIN − vIN
i1 = −
R1
vOUT = vIN − i1 R2
Figure 31: Non-inverting amplifier configuration.
The “virtual short” effect causes the op-amp’s
R1 input terminals to have nearly equal potentials, so
= vIN + vIN v− ≈ v+ .
R2 Summary: This configuration’s characteristics are:
? R1
⇒ Gni = 1+
R2 R2
G? = 1 +
R1
vout = A vin − v−
vout
⇒ v− = vin −
A
vout − v−
i2 = i1 =
R2
v−
=
R1
Rearranging we get:
R1 vout − v− = R2 v −
vout vout
⇒ R1 vout − vin + = R2 vin −
A A
1 R2
⇒ vout 1 + (1 + R2 /R1 ) = vin 1 +
A R1
G?
?
⇒ vout 1 + ni = vin Gni
A
?
A + Gni ?
⇒ vout = vin Gni
A
vout ? A
⇒G= = Gni ?
vin A + Gni
? A
⇒ G = Gni
A + 1 + R2 /R1
operational amplifier circuits 39
? R2
Gni = 1+
R1
A
e=
A + 1 + R2 /R1
?
Gni = Gni ×e
Notice that the error coefficient, e, is the same for both the
inverting and non-inverting configurations.
Generalized Result
Since the error coefficient is the same in both configurations, the
closed-loop gain can be generally expressed as
G = G? × e
? A
=G
A + 1 + R2 /R1
Make a copy of the inverting configuration circuit and modify it to implement a non-inverting config-
uration. Keep the parameters from the original exampe, R1 = 1 kΩ, R2 = 2 kΩ and A = 10 V/V, and
set the input signal amplitude to 1 V. For these parameters, calculate the expected values of G ? , e and
G. Simulate the circuit and verify that the output amplitude agrees with your calculations.
40 ece3410 lecture notes
Voltage Follower
−
The voltage follower represents a slightly different case, since
vOUT
there are no resistors. vIN +
In this configuration, we have the following device equations:
Figure 32: Voltage follower configuration. Due to the
vOUT = A v+ − v−
“virtual short” effect, vOUT ≈ vIN .
Summary: This configuration’s characteristics are:
= A (vIN − vOUT )
v A
⇒ G = OUT = G? = 1
vIN A+1
G = eG ?
In this case, the gain can be expressed as A
e=
1+A
Gv? f = 1
A
ev f =
A+1
Gv f = Gv? f × ev f
Make a copy of the inverting configuration circuit and modify it to implement a voltage follower
configuration. Keep the same op amp gain from the original example, A = 10 V/V, and set the input
signal amplitude to 1 V. Calculate the expected values of G ? , e and G. Simulate the circuit and verify
that the output amplitude agrees with your calculations.
operational amplifier circuits 41
Difference amplifiers
? R2
Gni = 1+
R1
− −
R2
Gi? = − vIN R1 vOUT
R1 + +
R3
R3 R2 R1
=
R3 + R4 R1 R2 + R1
R2
=
R1 + R2
R3 R
1+ = 1+ 1
R4 R2
R3 R1
⇒ =
R4 R2
Ad
CMRR =
ACM
This figure is often specified in dB. Ideally it should be infinite.
operational amplifier circuits 43
A
e= R2
.
1+ A+ R1 + Rsig
44 ece3410 lecture notes
Instrumentation amplifiers
vi1 + R3 R4
vx
R20
Disadvantages:
• Needs three op amps.
− iy
R3 R4
• Higher power consumption. vy
vi2 +
R2 R
v x = vi1 1 + − vi2 2
R1 R1
R2 R
vy = vi12 1 + − vi1 2
R1 R1
The overall gain of the pre-amplifier stage is then
v x − vy
A D1 =
vi1 − vi2
2R
= 1+ 2
2R1
R
= 1+ 2.
R1
The difference amplifier contributes a gain of R4 /R3 , so the
total differential gain is
R2 R4
AD = 1 + .
R1 R3
i x − iy = 0
⇒ v x = vy
ACM1 = 0V/V.
• Offset voltage
Based on this example, we can see that the effect of Ibias Ibias
is to introduce a DC offset voltage on vout . This places Figure 35: Inverting configuration showing bias-
current sources.
a limitation on the size of R2 that can be used. Suppose,
for instance, that we have
Ibias = 10µA R2 = 1MΩ VR = 5V
and the op amp’s power rails are at ±VR . In this case, the bias current induces an output offset
voltage equal to
Ibias R2 = 10V,
which is greater than the rail of the op amp. As a result, the op amp will simply saturate.
operational amplifier circuits 47
Given Ibias , an input signal vIN and a desired closed-loop gain G, how can we determine the maxi-
mum allowable value for R2 ? Suppose vmax is the maximum value of vIN , and vmin is the minimum
(note that vmin can be a negative voltage). Then our circuit must satisfy
VR + Gvmin
Ibias R2 + Gvmin < VR ⇒ R2 < .
Ibias
So returning to our example where Ibias = 10 µA and VR = 5 V, and let G = −10 V/V and vmin =
−0.1 V, we find
(5 V) + (1 V)
R2 < = 600 kΩ.
10 µA
48 ece3410 lecture notes
Offset Voltage
vout = A v+ − v− + Vofs .
R2
R1
vin −
v− vout
+
Vofs
−
+
Suppose an inverting op amp configuration
Figure 36: Inverting configuration showing input
has supply rails equal to +5V and −5V, and offset voltage source.
is configured to have a closed loop gain Effect of Vofs
G = − R2 /R1 = −100V/V. The input signal
is a sinusoid with peak-to-peak amplitude 5 ∗
vout
Output Voltage [V]
This circuit implements models of both Ibias and VOFS in a non-inverting op amp configuration. Since
the EveryCircuit op amp model is very ideal, a slight circuit trick is used to model the bias current by
steering it into ground instead of into the op amp terminal. This trick doesn’t change anything at all
about the circuit’s behavior. The model uses typical values of 10 µA and −10 mV for the bias current
and offset voltage, respectively.
We see that the output waveform has a significant DC offset due to the bias and offset effects, and
part of the waveform is saturated. To get some experience with these effects, you can experiment
with larger and smaller values of each, and with positive and negative values of VOFS . Occasionally
the simulation will halt and complain that it can’t find a solution. Usually in these cases you can just
restart the simulation and will proceed without any problems.
Design question: how can the circuit be modified to minimize the undesirable offset and avoid
saturating the output waveform?
50 ece3410 lecture notes
100
Gain Magnitude (dB)
ωc
50
ωt
0
100 101 102 103 104 105 106 107 108
0
−50
Phase (°)
−100
−150 PM
PM = 180° + ∠ A ( jωt )
Unity-Gain Frequency
A typical op amp has a a very large DC open-loop gain, often
greater than 80 dB or 10 000 V/V. Then the magnitude response
can be approximated as
s
A20
| A (ω )| =
1 + ω 2 /ωc2
| A0 | ω c
≈
ω
ωCL = ωt / (1 − G ? )
= ωt / (1 + R2 /R1 ) .
100
Gain Magnitude (dB)
ωc
50
ωCL
ωt
0
100 101 102 103 104 105 106 107 108
ωt = 10MHz
R2 /R1 = 10
What are the closed-loop DC gain, the 3dB cutoff frequency, and the unity-gain frequency for these
parameters?
V
ACL = −10
V
ωt = 10MHz
ωc = 909kHz
operational amplifier circuits 53
ωt = 10MHz
R2 /R1 = 100
What are the closed-loop DC gain, the 3dB cutoff frequency, and the unity-gain frequency for these
parameters?
V
ACL = −100
V
ωt = 10MHz
ωc = 99kHz
Notice that as the desired gain G grows to be very large, the cutoff frequency approximates to ωt /G.
This circuit models an op amp with a single-pole transfer function connected in a non-inverting
configuration. Since EveryCircuit’s built-in op amp model is basically ideal, we have to insert extra
components to introduce a pole at the op amp’s output node. This is accomplished using an RC
low-pass network followed by a voltage buffer comprised of a dependent voltage source with a gain
of 1 V/V.
Perform a frequency simulation of the closed-loop system and observe the major parameters: the
DC gain ACL (in dB), the 3 dB cutoff frequency ωc , and the unity-gain frequency ωt . Next, increase
the value of R2 by 10× and then 100×, and observe how it changes ACL and ωc . Convert ACL to
V/V in order to test the predictions from the theory presented in this section. Verify that ωt remains
constant and is approximately equal to ACL ωc .
54 ece3410 lecture notes
Slewing
Slewing is very different from ordinary transfer-function based A signal affected by slewing
0 , i.e.
behavior. Slewing can be thought of as saturation of vout
a second-order saturation effect. As such, it is fundamentally 1
vout
SR = max vout given in V/µs. 0
dt
Slewing tends to turn the output signal into a triangle wave.
If the op amp’s input signal is a pure sinusoid, then we can
−1
determine if the output will be affected by slewing:
0 2 4 6 8 10
∗ Time
vout = | ACL | VA sin (2π f t)
Figure 40: Slew-rate distortion in an op amp circuit.
d ∗
⇒ v = 2π f | ACL | VA cos (2π f t)
dt out
where VA is the input signal amplitude. This tells us that slew-
ing may occur if VA is large, or if the frequency f is large, or if
the closed-loop gain | ACL | is large. The maximum rate of signal
change must be less than the slew-rate:
2π f | ACL | VA ≤ SR
Example 10.
SR = 1V/µs
ACL = 2 V/V
f = 100kHz
SR
VA max =
2π f | ACL |
= 0.796V.
Clearly the slew rate can present real limitations for a circuit.
SR
FPBW =
2πVR
SR
FPBW =
πVR
The circuit will process any frequency less than the FPBW
distortion-free. For higher frequencies, you may begin to see
spurious harmonics in the output spectrum.
Once the FPBW is known, the slewing limit can be predicted
as follows:
FPBW
VO max = VR
f
Hence if you know the FPBW and the rail voltage, you can
estimate the maximum allowable amplitude at a given high
frequency.
operational amplifier circuits 57
Z2
Op Amp Integrators and Differentiators
1
Z1 (s) =
sC1
⇒ ACL (s) = −sC1 R2 .
Integrator: Z2 is a capacitor R1
vin −
vn
If Z2 is a capacitor C2 , and Z1 is a resistor R1 , then vout
+
1
Z2 (s) = Figure 42: Ideal Miller integrator.
sC2
R
⇒ ACL (s) = − 1 .
sC2
In this transfer function, − R2 /C1 is just a scale constant.
The signal is multiplied by 1/s, resulting in integration in the C2
time-domain.
C1
Z1 and Z2 are both capacitors
vin −
vn
If both of the impedances are capacitors, then the behavior is vout
+
similar to an inverting configuration.
Practical Considerations RF
C2 C2
C1 C1
vin − vin −
vn vout vn vout
+ +
+ +
− VOFS − VOFS
Phase 1 Phase 2
Figure 46: The two switching phases of a capacitive
inverting configuration.
In Phase 2, the input is connected, and node vn is left float-
ing, so that no charge can be added or removed from vn . So
any charge added to the outside plate of C1 has to be balanced
by an opposite charge on the outside plate of C2 . The charge
Q1 = vIN /C1 must be balanced by Q2 = vOUT /C2 = − Q1 .
operational amplifier circuits 59
Therefore
vOUT v
= − IN
C2 C1
vOUT C2
⇒ =− .
vIN C1
Miller Integrator
Frequency Analysis
the magnitude:
1
ACL ( jω ) =
ωR1 C2
1
⇒ ωt = .
R1 C2
+ vD −
A diode is like a valve that lets current flow one direction
vA vB
but not the other. It is a nonlinear device, which means the
traditional linear analysis techniques cannot be directly applied.
iD
We begin with some simplified models that are useful for
Figure 49: Diode symbol and notation.
building intuition about diode circuits. After that, we’ll build up
to more accurate (but difficult) models and techniques. We’ll see
that accurate simulation using SPICE (or a similar software tool)
is essential for designing nonlinear circuits.
Analysis steps:
D1
1. Suppose both diodes are OFF. Then vC = 0. But then both vA
D1 and D2 have positive potentials across their terminals,
so they cannot both be off. v A = 4V
2. Observe that D2 has the larger forward potential across D2 v B = 7V
its terminals. Based on this, suppose that D2 is ON while vB vC
D1 remains OFF. In this case, vC = v B = 7V, hence the
potential across D1 is vC − v A = −3V, which is consistent R
with the hypothesis.
vC = max (v A , v B ) .
Analysis steps:
VDD
1. Suppose both diodes are OFF. Then vC = VDD .
Exponential model
A more accurate model of the diode is given by this expression: Diode physical device parameters:
IS = Scale current, typ. ≈ 1pA to 1nA
vD
i D = IS exp −1 n = Grading coefficient, typ. ≈ 1
nUT
kB T
UT = Thermal voltage, ≈ 26mV at room temp.
Notice that when v D = 0, the current is also zero. When q
v D > 0, the exponential part rapidly becomes much greater than k B = Boltzmann constant = 8.6173 × 10−5 eV/K
T = Temperature (K) ≈ 300K at room temp.
one. When v D < 0, the exponential part rapidly becomes much
q = Elementary charge = 1eV/V
smaller than one. This splits the diode into two different modes:
called forward bias and reverse bias, respectively.
10
vD
Forward bias i D ≈ IS exp nU T
.
8
Reverse bias i D ≈ − IS .
6
i D (mA)
Constant voltage-drop model 4
From the physical model, we can see that the i D curve becomes
2
very steep when v D ≥ 0.7V, so we can say this is approximately
the ON voltage of the diode. When a diode is turned ON, it 0
0 0.2 0.4 0.6 0.8 1
should have a nearly constant forward voltage drop equal to
vD
0.7V.
Figure 52: Diode transfer characteristic. The current
increases very rapidly when v D ≈ 0.7 V.
Analysis steps:
1. Complete the analysis using the ideal switch model. We find that D1 is ON and D2 is OFF.
2. Estimate a more accurate result by adding a 0.7V forward drop to every diode that is ON, hence
vC = v A + 0.7V = 4.7V
Based on our analysis, we can generalize the result and describe this circuit by the function
vC = min (v A , v B ) + 0.7V.
Important Note: If |v A − v B | < 0.7V, then this method does not yield any valid solution. In that case,
we must use the full physical model with iterative analysis to arrive at the correct solution.
66 ece3410 lecture notes
Iterative Analysis
Most simulators will allow you to adjust the abstol and reltol
parameters. Smaller values result in better accuracy, but will
take more time to finish.
introduction to diodes 67
68 ece3410 lecture notes
iD
Linearized Model
+
id
Yet another way of modeling the diode is to use a linear approxi- +
mation. vd rd
Recall the definitions of small-signal notation, ID and VD are
the operating point values, id and vd are small variations, and i D −
and v D are the actual physical signal values. Hence vD ID
i D = ID + i d
+
− VD
v D = VD + vd
VD = 0.7V
ID = 1mA ( k +1)
iD
rd = 26Ω +
id
+
Iteration with the linearized model
(k)
vd rd
Iterative analysis can be combined with small-signal analysis by
repeatedly recalculating rd . In this version, the diode’s circuit −
( k +1) (k)
model looks like this: vD iD
R i ( k +1)
vA
+ id
+
R vd (k)
rd
vA
−
≈ ( k +1)
vB i (k)
+ (k)
− vB
Figure 56: Iterative solution of the resistor-diode series configuration using the linearized model.
In this example, the iterative procedure yields the following table of results:
k i [mA] v B [V] rd [Ω] id [mA] vd [V]
0 1 0.7 26 2.25504 0.032943
1 3.5504 0.73294 7.3232 −1.0705 −0.009330
2 2.4799 0.72361 10.4843 −0.1934 −0.002112
3 2.2865 0.7215 11.3713 −0.0000079 −0.00009
4 2.2786 0.72141 11.4106 0 0
Diode Circuits
vin vout
Half-Wave Rectifier
R
The 1/2-wave rectifier circuit passes current only when vout >
0.7V. In this case, the diode’s forward voltage drop is close to
0.7V, regardless of the current that flows, so that vout ≈ vin − 0.7V. Figure 57: Half-wave rectifier circuit.
When the diode is OFF, no current flows, so vout ≈ 0V. This
behavior is approximately described by the expression
vin
1
vout (ideal)
Voltage
−1
Time
Figure 58: Behavior of the half-wave rectifier. The
ideal switch model is compared to the more accurate
constant-0.7 V drop model.
In this example, let vin = −1V and R = 1kΩ. We want to solve for vout . First, we assume the diode is
OFF and check for consistency. We find that vout = 0 and therefore the diode’s forward drop is v D =
−1V. Since v D is negative, the diode must be OFF, so vout = 0V.
72 ece3410 lecture notes
−
+
1 kΩ
vin − 0.7 26 Ω
1 1 0.7 V
⇒ vout + = 1mA +
26 1kΩ 26
vin − 0.7
⇒ vout = (1mA) (26Ω k 1kΩ) + (26Ω k 1kΩ)
26
Now notice that (26 k 1000) ≈ 26 (try it). Then we can simplify the approximation:
This provides a more accurate approximation when the resistor R is large. In the case where vout =
1V, we find that
vout ≈ 0.343V
i D ≈ 343µA
diode circuits 73
Resistor-diode regulator
The regulator circuit is similar to the 1/2-wave rectifier, only vin vout
it interchanges the positions of the diode and resistor. In this
circuit, when vin < 0.7V, the diode is either OFF or only weakly
ON, so the current is close to zero. In that case, the voltage drop
across R is nearly zero, so vout ≈ vin . When vin > 0.7V, the diode
is clearly ON. Using the constant voltage drop model, we find
Figure 60: Single-diode regulator circuit.
that vout ≈ 0.7V, so the waveform is “clipped” at 0.7V.
0.6
0.4
0.2
0
Time
Figure 61: Behavior of the single-diode regulator
circuit with R = 100. The results from SPICE
simulation are more accurate than hand analysis.
A more accurate analysis is obtained using the linearized
diode model. By applying the node-voltage method at vout , we
find that
1 mA
As the name implies, regulators are used to produce sta-
ble DC voltages. Ideally, a regulator should produce 0.7V re- +
gardless of vin (so long as vin > 0.7V). The preceding analysis 0.7V −
26Ω
∆vout = ∆vin . Figure 62: Linearized model of the single-diode
R + 26Ω regulator.
Peak rectifier
Z tF
1
vout = i D (t) dt.
C 0
When the diode is OFF, a small current still flows, and that
current is steadily accumulated by the capacitor’s integrating
behavior.
Consider the output from a SPICE simulation where C =
1nF, shown below. In this simulation, we can see that vout rises
initially to 0.263V, which is approximately vin − 0.437V. Clearly
the 0.7V model is not working.
vin
0.5 vout
Voltage
−0.5
Time
2
Current
0
Time
Figure 63: Current delivered into the capacitor in the
peak detector circuit.
Envelope detector
Envelope Detector
4
vin
vout
2
Voltage
−2
Time
Figure 65: Behavior of the envelope detector circuit.
When the didoe is OFF, the output waveform is described by
76 ece3410 lecture notes
Netlist 1: envelope_detector.sp
* envelope detector circuit
* Transient simulation:
.tran 1m 1.5
.end
diode circuits 77
Bridge Rectifier
D2
1
D
− load +
vIN +
0 3
D3
4
D
vOUT
−
Figure 66: Full wave bridge rectifier circuit.
Analysis:
Case 1: vIN > 0. In this case, we see that the most
positive potential appears at the anode of
D2 . Based on this, we may predict that D2 is
ON while D1 is OFF. Since the most negative
potential appears at the anode of D4 , we may
conclude that D4 is OFF.
Based on this reasoning, we infer that the
current flows in a zig-zag through D2 , then the
load, then D3 . The potential appearing across
the load is
vOUT ≈ vIN − 1.4V.
−5
−10
Time
Netlist 2: bridge_rectifier.sp
* bridge rectifier circuit
* Bridge rectifier:
D1 0 2 diode
D2 2 3 diode
D3 0 1 diode
D4 1 3 diode
* Load resistor:
Rload 3 0 1k
* Transient simulation:
.tran .1m 0.02
.end
diode circuits 79
VDD
Voltage Regulators
R
We previously considered a 0.7V regulator circuit. We can
vOUT
extend this concept to produce other regulated voltages by
connecting multiple diodes in series. For example, we may
connect four diodes in series to create a 2.8V regulator circuit:
Let
Basic analysis
Find R to get an average current of 1mA, resulting in vout = 2.8V.
vin − vout
1mA = I =
R
vin − vout
⇒R=
I
= 7.2kΩ
The behavior of this circuit is investigated using SPICE simulation. The results shown below include
a supply ripple with zero-to-peak amplitude of 0.5V at 120Hz. From SPICE simulations, we see
that the actual output voltage is 2.7863V, which is slightly less than the intended value. The ripple
amplitude is also found to be 14.158mV.
By using the small-signal model, we can obtain a reasonable estimate of vout , the small ripple wave-
form that is superimposed on the regulator’s output:
vout
LR =
vdd
7.12mV
=
500mV
= 0.014239
= 1.4239%
diode circuits 81
Netlist 3: basic_regulator.sp
* 2.8V regulator circuit
* Regulator circuit
* The output is at node 2
R1 1 2 300
D1 2 3 diode
D2 3 4 diode
D3 4 5 diode
D4 5 0 diode
* Transient simulation:
.tran .1m 0.02
.end
10 vin
vout
Voltage
0
Time
Figure 69: Behavior of the four-diode regulator from
SPICE simulation.
82 ece3410 lecture notes
2.79
Voltage
2.79
2.78
Time
Figure 70: Zoomed view of the ripple voltage on
vOUT .
General Analysis
NVT /ID
LR =
R + NVT /ID
NVT
=
RID + NVT
NVT
= .
VDD − VOUT + NVT
vin +
0
−2
−4
Time
84 ece3410 lecture notes
Netlist 4: superdiode.sp
* super-diode precision rectifier simulation
* 741 instance
* Pin order: v+ v- VR+ VR- vo
X1 1 2 10 0 3 uA741
D1 3 2 diode
RL 2 0 1k
.tran 1m 0.5
.end
C
vin vout
DC Restoration, Clamped Capacitor
q (t) = i D t
−vout
= IS exp t.
VT
vin
vout
2
Voltage
−2
Time
Figure 73: Behavior of DC restorer (clamped capaci-
tor) circuit simulated in SPICE.
86 ece3410 lecture notes
Netlist 6: dc_restorer.sp
* DC restoration circuit
* Transient simulation:
.tran 1m 2
.end
L
VIN VOUT
Boost converter iL
C
Diodes are frequently used in power conversion circuits. In this
appendix we look at one important step-up DC-to-DC converter
circuit, known as the boost converter. The boost converter con- Figure 74: Idealized boost converter circuit.
sists of an inductor, a diode, a switch and a load capacitance, as
in the schematic shown in Figure 74.
When the switch closes, the inductor is shorted to ground,
resulting in a large current. When the switch opens, the induc-
tor’s current cannot change instantly, so the current is forced
through the diode into the capacitor. This establishes a large
potential across the capacitor.
More precisely, suppose that the switch is initially open and
the current i L is zero. Then, at time t = 0, the switch closes
abruptly. The current is then
Z t
1 VIN
i L (t) = VIN dτ = × t.
L 0 L
Then if the switch opens again at some time t1 , the inductor will
possess a stored energy equal to
1 2
E ( t1 ) = LI .
2
Since the current must continue flowing through the inductor,
most of this energy is transferred into the capacitor, and some
of the energy is dissipated in the diode. For our purposes, we’ll
use the ideal switch model and ignore losses in the diode, hence
we will assume all the energy is transfered to the capacitor.
diode circuits 87
DT
∆i L (on) = VIN
L
(1 − D ) T
∆i L (off) = − (VIN − VOUT )
L
If we set the rise and fall equal to each other (as required for
steady-state operation), then we can solve for VOUT :
50
VIN
VOUT =
vOUT (V)
1−D 40
Note that this analysis only works if the switching is very fast,
30
so that the inductor current never drains completely to zero.
If the switching clock has a 50% duty cycle, then the circuit
20
acts as a voltage doubler. An example SPICE simulation follows
0 0.5 1 1.5
with VIN = 20 V and a switch frequency of 5 MHz with a duty
Time (s) ·10−3
cycle of 62%. The expected output is VOUT = 52.6 V. The simula-
Figure 75: Output voltage from the boost converter
tion output, shown in Figure 75, approaches the expected limit when initialized at zero.
after about 1 ms, which corresponds to five thousand switching
cycles in this example. Note that the simulation does account for
energy lost in the diode. The main effect of diode losses is that
the circuit takes longer to approach the steady state condition.
But just like with the peak detector circuit, the diode’s loss gets
smaller in each cycle, and some energy is delivered to the capac-
itor in each cycle, so the circuit asymptotically approaches the
ideal limit.
Netlist 7: boost_converter.sp
* Boost converter simulation
* switch model:
.model switch sw(Ron=5, Roff=100000, Vt=0.001, Vh=0.0001)
.end
Memristors
+ v (t) −
The “memristor” refers to a nonlinear two-terminal device for
which the resistance changes in response to an applied voltage
or current signal, and in which the resistance stays constant for Figure 76: Memristor symbol. The device acts like a
some period of time when the signal is removed. This behavior time-varyingR resistor, where the resistance changes in
response to v (t) dt.
is described as a “memory resistor,” hence the abbreviated
name, “memristor”. A device which exhibits resistance memory
is said to be a “memristive” device or system.
Memristive devices were first observed more than a century
ago, but they lacked a general theory to understand and apply
their behavior in the context of circuit engineering. The first
general theory was articulated by Leon Chua in two papers,
written in 1971 and 1976 2 . Chua’s theory began as a somewhat 2
Leon O Chua. Memristor-the missing circuit
speculative hypothesis. Interest surged in memristors after a element. Circuit Theory, IEEE Transactions on, 18(5):507–
519, 1971; and Leon O Chua and Sung Mo Kang.
group from HP Labs published a paper titled “The missing Memristive devices and systems. Proceedings of the
memristor found,” which showed that nano-scale resistance- IEEE, 64(2):209–223, 1976
switching devices have the characteristics predicted by Chua’s
theory 3 . 3
Dmitri B Strukov, Gregory S Snider, Duncan R
Stewart, and R Stanley Williams. The missing
memristor found. nature, 453(7191):80–83, 2008
Axiomatic Circuit Theory
• Voltage v (KVL)
• Current i (KCL)
capacitor
Constitutive Relations dq = C dv
voltage charge
In order to define the behavior of “black box” devices, we begin v q
with a standard format, called the constitutive relation, that
can be used to describe any passive two-terminal device. A
dv
=
constitutive relation is an expression that defines a conservative
φ
dt
relationship between two of the fundamental quantities. The resistor
dv = R di memristor
traditional linear components are defined by familiar relations,
dφ = M dq
expressed in both their usual and differential forms
t
id
=
dq
usual form differential form
current flux
v = Ri dv = R di
i φ
inductor
q = Cv dq = C dv dφ = L di
φ = Li dφ = L di Figure 77: An illustration of the constitutive relation-
ships corresponding to standard components and
In addition to the three standard passive elements, we may definitions. One relation is missing in the classical
define two more constitutive relations based on the definitions theory: the memristor.
i = IS ev/nUT .
dφ = M (φ, q) dq.
Current
0
reason behind the name “memory resistor.”
−2
Memristor Properties
−4
Chua authored several papers exploring the possible attributes
−4 −2 0 2 4
of a memristor device. The most important features can be
Applied Voltage
deduced from the constitutive relation itself. First of all, since
Figure 78: Pinched hysteresis in a memristor, driven
the memristor behaves like a time-varying resistance, we expect by a sinusoidal voltage across its terminals. In this
that when zero volts are applied, its current should also go model, we assume resistance increases when a rising
to zero, just as it would for any resistor. Second, if a positive positive voltage is applied (blue curve). As the
resistance increases, the blue curve’s slope becomes
voltage is applied, the flux and charge should begin to change, flatter. When the positive voltage begins to decrease
and we therefore expect the resistance to change over time. (red curve), the current is lower due to the increased
resistance in this portion of the curve. When the
Third, if a negative voltage is applied, it should reverse the voltage reaches zero, the current also goes to zero.
changes in the flux and charge, and so the resistance should When the applied voltage swings negative (green
change in the opposite direction. curve), the resistance decreases so the curve becomes
more steep. When the applied negative voltage
The above reasoning predicts what is known as pinched rises toward zero again (black curve), the current
hysteresis, which is the key attribute of any memristor device. is stronger (more negative) due to the decreased
resistance. Once the applied voltage returns to
Pinched hysteresis is observed on a type of plot called the Lis- zero, the memristor’s resistance should return to its
sajous figure, which displays the device’s current and voltage original value.
when being driven by a period source. For an ordinary resistor,
the Lissajous figure should be a straight line. For capacitors and
inductors, the Lissajous figure is a circle or ellipse. For a mem-
ristor, we see the “bow tie” pattern shown in Figure 78. This
pattern is not observed for any other device, and is considered
a “fingerprint” of memristance.
A second property of memristance is non-volatility: when A Real Memristor
integrator
erty can be useful, for example, in reconfigurable radio circuits,
where low-frequency signals can be used to “tune” the device’s C
Simulating Memristors
+
−
sensitive to simulation parameters and may crash the simulation
if conditions are not just right. Furthermore memristor models C
often become invalid outside of a limited range of voltage and
frequency, so they are not necessarily general-purpose models.
+
One of the simplest models uses an op amp integrator to − ro
model the flux in the memristor, as shown in Figure 80, where R +
B +
the integrator’s output voltage vφ should be proportional to vφ
v (t)
the flux φ. Then, to model the changing resistance, we have to co
utilize a special nonlinear dependent current source provided rs
memristor_integrator_model.sp
* Memristor integrator model
memristor_integrator_simulation.sp
* Memristor simulation based on integrator model
.include memristor_integrator_model.sp
.control
* Medium-frequency simulation:
tran 10n 50u uic
plot -i(v1) vs v(nplus)
plot v(X3.n3)
wrdata integrator_mid_freq v(nplus) i(v1) v(X3.n3)
* Low-frequency simulation:
alter @V1[sin]=[ 0 5 1k ]
tran .01m .0025 uic
plot -i(v1) vs v(nplus)
plot v(X3.n3)
wrdata integrator_low_freq v(nplus) i(v1) v(X3.n3)
* High-frequency simulation:
alter @V1[sin]= [ 0 5 1MEG ]
tran 1n 5u uic
plot -i(v1) vs v(nplus)
plot v(X3.n3)
wrdata integrator_high_freq v(nplus) i(v1) v(X3.n3)
.endc
.end
memristors 97
Lissajous figure for the mid-frequency case Integrator state (flux) for the mid-frequency case
100 0
50
Current (mA)
−5
vφ (V)
0
−10
−50
−15
−100
−6 −4 −2 0 2 4 6 0 10 20 30 40 50
Voltage (V) Time (µs)
Lissajous figure for the low-frequency case Integrator state (flux) for the low-frequency case
0
800
600 −50
Current (mA)
400
vφ (V)
−100
200
0 −150
−200
−200
60
0
40
−0.5
Current (mA)
20
vφ (V)
0
−1
−20
−40 −1.5
−60
−6 −4 −2 0 2 4 6 0 1 2 3 4 5
Voltage (V) Time (µs)
Figure 84: Simulation results for the integrator-based
memristor model at three different frequencies.
memristors 99
Memristor Applications
H
RRAM
memristor_integrator_model_with_threshold.sp
* Memristor integrator model with threshold
RRAM_threshold.sp
* Memristor resistive RAM demo
.endc
.end
104 ece3410 lecture notes
RRAM Simulation
5
READ
1 Vread
Voltage (V)
Vmemristor
high low
0.5 res res
0
−0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Time (s) ·10−3
0
vφ (V)
−100
−200
−0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2
Time (s) ·10−2
Claims of Credit
Pedantic Arguments
Philosophical Disputes
Pessimism
Conclusion
Introduction to MOSFETs +
Gate (G) v DS
+
MOSFET as switch
The MOSFET has three terminals, source, gate and drain. We iD
may first understand the MOSFET as a logic switch. In this Source (S)
+
model, the terminal potentials are interpreted as logic val- +
|vGS |
ues, i.e.the logic set {0, 1} is mapped to the potential values
{0, VDD }. Under this model we may consider all signals to be −
Gate (G) |v DS |
either HIGH or LOW. Then the behavior is as follows:
Device Type vG Device State v DS
NMOS HIGH ON small −
Drain (D)
LOW OFF large
iD
PMOS HIGH OFF large
LOW ON small
Figure 91: PMOS device symbol, which is comple-
Notice two things: mentary to the NMOS device. The drain and source
terminals are flipped (current is understood to flow
• The NMOS and PMOS have complementary behavior, i.e. they vertically downward from the source to the drain). A
“bubble” is commonly drawn at the gate to indicate
have opposite states in response to the gate voltage. that the device responds to the logical complement
of the gate signal. For PMOS devices, the substrate
• When the device is ON, the drain-source voltage v DS must be is usually connected to the circuit’s most positive
quite small. When the device is OFF, the drain-source voltage potential (usually VDD .
can be large. This is the behavior we expect from a switch.
112 ece3410 lecture notes
Electrical Characteristics
vOUT
5
the MOSFET must pull a constant current through the resistor
R in order to maintain a low output voltage. It is nevertheless
helpful to study the RTL inverter, because its properties are 4 I II
vOUT
three regions, representing the different modes of the NMOS
device: 2 VTh
vOUT
vGS = vIN and vDS = vOUT . If we imagine that vIN is initially
zero, and is slowly increased toward VDD , then we have three 2 VTh
subproblems:
I Cutoff— When vIN < VTh , verify that vOV < 0, therefore 1 III
iDS = 0. In that case, there is no current flowing through R, so
vOUT = VDD = 5 V. 0
0 1 2 3 4 5
vIN
II Saturation — When vIN > VTh , then vOV > 0 and the device’s
Figure 94: Comparison of analysis and simulation re-
current is given by the square-law equation. Then vout is sults in the three operating modes, cutoff, saturation,
determined by the voltage drop across R: and triode.
1
vOUT = VDD − Rkv2OV
2
= 5 V − 5 (vIN − 2 V)2
As vIN increases, vOUT will decrease until the MOSFET enters
the triode mode. That transition happens when vOUT = vOV ,
i.e.
1
vOV = 5 V − Rkv2OV
2
1
⇒ 0 = Rkv2OV + vOV − 5
2
Since the result is a quadratic equation, we can apply the
standard formula and solve:
r
−1 ± (1)2 − 4 12 Rk (−VDD )
vOV =
√ Rk
−1 ± 1 + 2RkVDD
=
r Rk
2VDD
if 2Rk 1 : vOV ≈
Rk
Note that we chose the positive result in the quadratic equa-
tion, since vOV has to be positive in both saturation and triode,
otherwise these equations wouldn’t apply. The results are:
III Triode — When vIN > 2.5188 V, the device should enter
triode, and the new device equation is
1 2
iDS = k vOV vOUT − vOUT
2
1 2
⇒ vOUT = VDD − Rk vOV vOUT − vOUT
2
1
0= Rkv2OUT − (1 + RkvOV ) vOUT + VDD
2 q
(1 + RkvOV ) ± (1 + RkvOV )2 − 2RkVDD
vOUT =
q Rk
(1 + RkvOV ) − 1 + ( RkvOV )2
=
Rk
Note that if RkvOV 1 the equation simplifies to vOUT = 0.
116 ece3410 lecture notes
II Saturation — When vOV > 0 and |vDS | > vOV , the device’s 3
vOUT
current is given by the square-law equation. This corresponds
to the case when: 2
5
And the corresponding values of vIN are:
vOUT
2
III Triode — When |vDS | < vOV , the device enters triode. This
corresponds to the case where
1 II I
VDD − vOUT < VDD − vIN − |VTh |
⇒ vOUT > vIN + |VTh | 0
0 1 2 3 4 5
vIN
or when vIN < vOUT − |VTh |
Figure 97: Comparison of analysis and simulation
results for the PMOS RTL inverter.
In this region, the device’s current and output voltage change
as follows:
1
iDS = k vOV |vDS | − v2DS
2
1
⇒ vOUT = Rk vOV |vDS | − v2DS
2
1 2
⇒ VDD − |vDS | = Rk vOV |vDS | − vDS
2
1
0= Rkv2DS − (1 + RkvOV ) |vDS | + VDD
2 q
(1 + RkvOV ) ± (1 + RkvOV )2 − 2RkVDD
|vDS | =
q Rk
(1 + RkvOV ) − 1 + ( RkvOV )2
=
Rk
Notice that this is the exact same result as before, only it’s
“upside down.” We can next get the solution for vOUT :
vOUT
tions have the same qualitative behavior. They both function
2 PMOS
as logic inverters. If we were to balance one of these circuits
right in the center of its saturation region (II), where the slope
is very steep, we could use it as an inverting amplifier. We will 0
0 2 4
soon introduce linearized amplifier models that apply in the vIN
saturation region; it will be important to recognize that both Figure 98: Overlay of the NMOS and PMOS RTL
NMOS and PMOS devices have the same linearized models in inverter transfer characteristics. Both devices behave
as an inverter. They differ slightly in the offset
saturation, just as they show the same behavior in the RTL in- voltage at which they “tip” from high to low.
verter configurations, even though they exhibit complementary
iDS
logical behavior.
Behavior in Saturation +
vDS −
+
We may now go one level deeper and examine the MOSFET’s −
behavior in the saturation mode. First, let’s understand why it’s
called “saturation.” Suppose we hold the gate potential fixed so
that vGS = 1 V and perform a DC sweep on vDS while measuring Figure 99: An experiment in which vGS is held
the current. constant while sweeping v DS .
According to the triode equation, the current should be a
parabola:
1 2
iDS = k (1)vDS − vDS .
2
But if that were true, the current would begin to decrease when
60
vDS > 1 (dashed curve below), and eventually the current would
swing negative, creating an impossible free-energy device. 40
Obviously this doesn’t happen, instead the device current rises
iDS (µA)
iDS (µA)
Channel Length Modulation 100
in i D that result from small variations in vGS and vDS . For hand
analysis, we can directly integrate the device equations:
d iD d 1
k (vGS − VTh )2
=
d vGS d vGS 2 DC
= k (VGS − VTh )
= kVOV
VGS = VD = VDD − ID R
1 2
ID = kV
2 OV
1
= k (VDD − ID R − VTh )2
2
√
To complete the solution, we define a variable x = ID , and The diode connection is guaranteed to always be in
then arrange the above equation into a quadratic equation: the saturation mode, since VDS = VGS it is always
assured that VDS > VOV .
r r
k 2 k
R x +x− (VDD − VTh ) = 0
2 2
Applying the quadratic formula:
p
−1 ± 1 + 2kR (VDD − VTh )
x= √
R 2k
p 2
−1 + 1 + 2kR (VDD − VTh )
⇒ ID = x 2 =
2kR2
A diode connected configuration is implemented with VDD = 5 V, R = 50 kΩ, and the MOSFET param-
eters are k = 500 µA/V2 , VTh = 0.5 V and λ = 0.05 V−1 . The analysis from this section predicts a bias
current of ID = 78.8 µA, which is quite close to the simulated value. We can also verify that the gate
voltage should be VG = 1.06 V, which is again quite close to the simulated result.
122 ece3410 lecture notes
Passive bias network: In this very general case, resistors are VDD
placed adjacent to both the source and drain terminals, and the
gate is biased at some constant voltage, such that the device is RD
held in its saturation mode. There are interactions at both the VD
drain and source terminal: ID
VG
VD = VDD − ID R D
VS = ID RS
VS
1
ID = k (VG − ID RS − VTh )2 RS
2
we solve this case in much the same way as the diode-connected
√
circuit. By defining x = ID , we can obtain a quadratic polyno-
mial:
r r
k 2 k
0 = RS x +x− (V − VTh )
2 2 G
p
−1 ± 1 + 2kRS (VG − VTh )
⇒x= √
RS 2k
and ID = x2 .
VDD
MOSFETs as Switches
because this is the voltage where vOV becomes zero, and device
switches OFF.
124 ece3410 lecture notes
This example shows the NMOS RTL inverter configuration. THe transient simulation shows that
when the NMOS device is ON, the active pull-down is very fast and effective. When the NMOS
device is OFF, however, the passive pull-up operation is very slow due to the RC delay. Note that this
example uses a manual switch to change the input signal. You have to click it to change the state.
This circuit shows a combination of the pull-up and pull-down configurations shown in Figure 105.
The two configurations are folded together, so there is no resistor current to overcome. The simula-
tion shows that the NMOS device pulls down very well, but the pull-up operation is both slow and
incomplete. A PMOS pull-down example is also available, showing the complementary behavior.
Suppose we are given a four-input logic function F = A( B + CD ). To synthesize a static CMOS imple-
mentation, we begin with a classic logic circuit, and then apply de Morgan’s laws to transform it:
C
D
B F C
A
1A. PMOS PUN: Using de Mor- D
B F
gan’s Laws, make gate transforma- A
tions to insert bubbles at the signal
inputs.
C
D
B F
A A B
1B. Cancel the double-bubbles and add an
inverter to eliminate the output bubble. D
Now there should be bubbles on all input
signals, but no interior bubbles. C
C
D
B F
A
2. NMOS PDN: Since there are already no interior bubbles,
simply place a bubble on the output signal and cancel it with
an inverter. Now both the PUN and PDN have inverters at A
the output; for the transistor implementation, we ignore the F
inverters. An inverter will be attached to the gate’s output after
we are finished.
The Exclusive OR (XOR) function is crucial for many logic and arithmetic circuits. The function is
F = A ⊕ B = AB + AB. To synthesize a static CMOS XOR gate, we begin with a classic logic circuit,
and then apply de Morgan’s laws to transform it as follows:
A
B 1. PMOS PUN: Just com-
F plement two of the input A
A signals to get bubbles on B
B all of them. F
A
B
Transmission Gates φ
Analog Switching
A
In addition to digital applications, transmission gates are useful
for switching analog signals. One common application is the
track and hold (T/H) circuit, which serves as the front-end for B
many sampling circuits, such as analog-to-digital converters.
The most basic T/H circuit contains only a capacitor and switch.
When the switch is ON, the capacitor “tracks” the voltage of Q
the input signal. When the switch is OFF, the capacitor is left
floating. In this condition, no charge can be added or removed
from the capacitor, so it “holds” whatever voltage it contained
at the moment when the switch turned OFF.
When a T/H circuit is used to sample a slow-changing ana-
log signal, then there will usually be a small signal difference
between each sampling event. As a result, |vDS | will tend to be Figure 111: XOR gate based on transmission gates.
small each time the transmission gate is switched ON, so the The logic function is Q = AB + AB. This circuit
devices start out in triode rather than saturation. This scenario requires a total of eight MOSFETs.
This circuit implements a passive T/H circuit based on the transmission gate switch. The circuit has
supply voltage VDD = 5 V and a 100 nF hold capacitor. The device parameters are as follows:
NMOS PMOS
kn = 500 µA/V2 kp = 250 µA/V2
VThN = 0.5 V VThP = 0.5 V
λn = 0.05 V−1 λp = 0.1 V−1
In parallel, the total ON resistance for this switch should be rON = 296 Ω. When the switch is turned
ON, the output signal’s rise time is determined by the time constant formed by rON and C, which
works out to be τ = 29.63 µs. Then the 10–90% rise time is tr = 2.2τ = 65.2 µs. This gives an indica-
tion as to the minimum switching period for this T/H circuit. In practice, 10–90% is not sufficient to
obtain a high-precision sample, so a switching period five to ten times slower may be required.
The simulation example shows an input signal at 100 Hz with a switching period of 500 µs, which is
almost 10× higher than the calculated rise time. By zooming in on the waveform, you can see that
the tracking accuracy improves gradually with more time in the tracking phase. The slow conver-
gence is a direct consequence of the ON resistance in the transmission gate.
130 ece3410 lecture notes
MOSFETs as Amplifiers
Common-Source Configuration
+
VIN −
As a first example, we consider the RTL inverter circuit, only
now we will try and balance the circuit at the point where
its transfer characteristic is steepest. We refer to this as the vout
quiescent point, Q point, bias point, or DC operating point.
vin
We then superimpose a small AC signal on top of the DC +
operating point. Amplifier design is therefore divided into two
vgs = vin gm vgs ro R
tasks: biasing and small-signal analysis. We’ll consider biasing
strategies later. In this section we focus on basic small-signal
−
analysis techniques, as they dictate amplifier behavior and
potential applications.
To begin with, we consider the common-source configura-
tion and assume it is appropriately biased at a suitable DC Figure 112: NMOS common-source amplifier config-
uration and its small-signal equivalent model. Since
operating point. To analyze the small-signal behavior, we re- VDD is shorted out in the small-signal model, the
place the MOSFET with its small-signal equivalent model (the bias resistor R appears in parallel with the device’s
internal resistance ro .
transconductance amplifier model). Second, we zero-out any
DC independent sources. This means that the VDD node gets
shorted to ground, so any devices connected to it are “folded
over” onto the ground node. Summary: CS amp
To analyze the amplifier characteristics, we use the small-
signal equivalent circuit to solve for the gain and output resis- • Inverting amplifier
tance. From the model in Figure 113, we see that the amplifier • Output resistance: ROUT = ro k R
consists of a current source and two resistors. Since the two • Gain: Av = − gm ROUT
This example shows a basic common-source configuration for an NMOS device with k n = 500 µA/V2 ,
λn = 0.05 V−1 and VThN = 0.5 V. The supply voltage is 5 V. The gate is biased with a DC operating
voltage of VIN = 0.9 V, and the bias resistor is R = 50 kΩ. Capacitive coupling is used at the gate to
separate the DC bias voltage from the AC small-signal input. Capacitive coupling is also used at the
drain to remove the DC offset from the output signal.
Based on these parameters, we can calculate the device’s small-signal characteristics, and then obtain
the gain and output resistance as follows:
Run the transient simulation and verify that the predicted gain and output offset are correct. You will
probably notice that the simulated output offset is 2.727 V. Can you explain this discrepancy? (Hint:
consider the effect of CLM with VDS = 2.727 V, then calculate new values for ID and VOUT ).
132 ece3410 lecture notes
voltage: when vIN rises, i D falls. Since the resistor is positioned vin
R
between the drain and ground, a smaller current means a
+
smaller output voltage at the drain. The result is that the small- VIN −
signal behavior is the same for both the NMOS and PMOS
versions.
vout
To obtain the small-signal equivalent circuit, we zero-out
VDD and VIN , so that the PMOS source terminal is connected to vin
+
small-signal ground. Even though the PMOS device current has
an inverse response to the gate voltage, we can flip the device vgs = vin gm vgs ro R
upside down so that the source terminal is folded back onto the
ground node. We then obtain the exact same model as we had −
for the NMOS version. What this means is that every NMOS
circuit configuration should have a complementary PMOS
version with the exact same behavior. The only differences will Figure 113: PMOS common-source amplifier configu-
be in the device’s k, VTh and λ parameters. ration. Its small-signal equivalent model is the same
as the NMOS version.
This example shows a PMOS version of the common-source amplifier. The parameters very similar to
the NMOS case: k p = 250 µA/V2 , λ p = 0.1 V−1 and VThP = 0.5 V. The supply voltage is 5 V. The gate
is biased with a DC operating voltage of VIN = VDD − 0.9 V, and the bias resistor is R = 50 kΩ. Then:
Run the transient simulation and verify that the predicted gain and output offset are correct. You will
probably notice that the simulated output offset is 1.364 V. Can you explain this discrepancy? (Hint:
consider the effect of CLM with |VDS | = VDD − 1.364 V, then calculate new values for ID and VOUT ).
introduction to mosfets 133
Avo = − gm ro
vout
Rout = ro
vin
+
The gain magnitude of this configuration, gm ro , is commonly
referred to as the intrinsic gain of the MOSFET, since it is the vgs = vin gm vgs ro
highest gain achievable with a single MOSFET device. When
−
the circuit is analyzed with no load attached, it is referred to as
the “open-circuit gain” and the subscript letter ‘o’ is added in
Avo to signify this.
In practice, a nearly-ideal current source can be implemented Figure 114: NMOS active-bias common-source am-
using a MOSFET device with a constant gate voltage. For plifier configuration and its small-signal equivalent
model. The current source directly forces a DC bias
example, a PMOS device can be substituted in place of the current of ID in the NMOS device. Since the bias
current source. The PMOS gate voltage, VGP , should be chosen current is forced by a DC independent source, it
is zeroed out in the small-signal model, leaving an
so that the device is biased in its saturation mode. In that open-circuit at the output node.
configuration, the PMOS device is insensitive to the voltage
at its drain terminal, so its constant gate voltage maintains a
constant bias current ID .
Since the PMOS device is not perfectly ideal, it contributes a VGP
load effect due to its intrinsic resistance ro . In the small-signal
ro,p
model, the NMOS and PMOS ro ’s will appear in parallel, so the vOUT
output resistance and gain are slightly modified: ro,n
vIN
Rout = ro,n k ro,p
Av = − gm Rout vin
By using a PMOS device the circuit’s gain is roughly cut in half VIN
+
−
due to the interaction of ro ’s. In general, an amplifier’s output
node is connected to two branches, one “going up” toward VDD
and another “going down” toward ground. The total output Figure 115: NMOS active-bias common-source
resistance is taken as the parallel combination of equivalent amplifier configuration with PMOS bias device. The
PMOS device acts as a current source.
resistance looking up with the resistance looking down, i.e.
Rout = Rup k Rdown .
134 ece3410 lecture notes
that case, the current flowing into the output branch must be
+ RS
zero. If a portion of the circuit is enclosed by the dashed box VIN −
shown in Figure 116, then the total current flowing into the box
has to equal the total current flowing out of the box (this is a
version of Kirchoff’s current law). The MOSFET does not allow id = 0 A
vout
any current at its gate terminal, so the gate current is zero. The
output terminal is open-circuited, so the drain current is also vin
+
zero. The only remaining branch is the source terminal, which ig = 0 A
must be zero since there is no other route for current to flow vgs = vin − vs gm vgs ro
into the box. Since is = 0, there is no voltage drop across RS , so
the source voltage is also zero. −
As a result of this analysis, the model for solving the gain of vs
this circuit is identical to the model in Figure 114, so the gain is
must be exactly the same, Av = − gm ro . Where the models RS
differ is in the output resistance. To find Rout for this circuit, we
zero out the input signal and apply a test voltage at vout . Then
we solve for the current that flows through the output branch. Figure 116: NMOS active-bias common-source
Since the output node is no longer open-circuited, a non-zero amplifier with source degeneration resistor RS . The
current flows through the drain and source terminals, with effect of RS is to reduce the amplifier’s gain while
improving error tolerance in the bias point.
id = is = iout . Also, since vin = 0, the gate-source voltage is
vgs = −vs = −is RS . Based on these considerations, we obtain
the circuit shown in Figure 117. is
To solve for the output resistance, we consider the voltage
drop across ro . Two downward currents are superimposed on
+
ro :
vgs = 0 − vs gm v s ro
vout = vs + ro ( gm vs + is ) +
− vout
= i s R S + r o ( gm i s R S + i s ) −
vout vs
⇒ Rout = = R S + r o + gm r o R S
is is
RS
So although the active-bias open-circuit gain is the same when
source degeneration is present, the output resistance is much
higher. This should result in a more significant coupling effect
when a load is connected. Figure 117: Finding the output resistance for the
degenerated amplifier.
introduction to mosfets 135
RD
AvL = (− gm ro )
R D + Rout vIN
− gm r o R D
=
R D + R S + r o + gm r o R S
vin
Another way of looking at it is that the resistance R D summa-
RS
rizes the circuit’s “up” branch, and the open-circuit amplifier
+
VIN −
summarizes the circuit’s “down” branch. The two branches can
be analyzed separately, and then joined together via a coupling
analysis. After coupling, the new overall output resistance is Figure 118: NMOS CS amplifier with passive bias
0 = R k R . and source degeneration.
Rout out D
In both the passive and PMOS biased circuits, we make use of vin
the idealized amplifier model shown below.
+ RS
VIN −
+
Rin → ∞ − − gm ro,n vin Rup
Consider a passive-biased common-source amplifier like the one shown in Figure 119. The NMOS
device has parameters k = 500 µA/V2 , VTh = 0.5 V and λ = 0.05 V−1 . The bias resistor is R D = 50 kΩ
and the degeneration resistor is RS = 20 kΩ. If the input offset voltage is VIN = 1.25 V, what is the
circuit’s gain?
To solve this problem, we first solve the DC operating point and then calculate the small-signal pa-
rameters. Referring back to the bias configurations studied earlier in the chapter, we see that this
circuit is already covered by the passive bias network analysis. In a previous example, we found that
√
the bias current should be ID = 22.5 µA. Then gm = 2kID = 150 µA/V2 , ro = (λID )−1 = 889 kΩ,
gm ro = 133 V/V, and Rdown = 3.58 MΩ. The high intrinsic gain looks pretty promising, but the
resistive coupling effect is going to ruin it. Putting all this together, the amplifier’s loaded gain is
RD
AvL = (− gm ro )
R D + Rdown
= −1.83 V/V.
With the large source degeneration resistance, the circuit does not make a very good amplifier.
If the degeneration resistance could be removed (while keeping the bias current the same), then the
output resistance would be a much smaller value of ro , and in that case the loaded gain would be
much better:
RD
AvL = (− gm ro ) = −7 V/V.
R D + ro
introduction to mosfets 137
VDD ndd 0 DC 5V
VIN ndc 0 DC 1V
vsig nsig 0 SIN(0 0.1 1k)
RL nout 0 1Meg
Cout nd nout 10uF
.control
* Foreach loop to scan through RS values:
foreach RSval 0 500 1000 5000 10000 20000
alter RS = $RSval
dc VIN 0 5 0.05
end
This example implements the NMOS common-source amplifier described in the SPICE netlist above,
with a passive bias resistor R D = 50 kΩ, a degeneration resistor RS = 20 kΩ, and a bypass capacitor
CB = 10 µF. With the bypass capacitor in place, the gain is close to −7 V/V as we predicted in exam-
ple 23. If you remove the bypass capacitor, you should notice that the gain drops to about −1.8 V/V
as predicted.
The circuit is highly tolerant of different DC input bias voltages. Try adjusting the DC gate voltage
source between 1 V and 1.8 V. The circuit continues to function throughout this range, while main-
taining a gain close to the −7 V/V target. Try repeating the simulations with RS and CB removed (i.e.
shorted out). The results will not be as robust for different gate offset voltages.
introduction to mosfets 139
RD
−
Rdrain = RS + ro + gm ro RS
vs
is +
This result covers all possible cases. When the source terminal vs −
is connected directly to ground, RS = 0, then Rdrain = ro . If
there is an ideal current source under the source terminal, then
RS → ∞, in which case Rdrain → ∞.
Figure 126: Model for solving the resistance looking
Resistance into the source: To find the resistance looking into the source terminal.
into the source terminal, we summarize any circuitry present
above the drain as an equivalent resistance, R D . We then apply
a test voltage at the source and solve for the current that flows
140 ece3410 lecture notes
vin
Common-Gate amplifier configuration
+
VIN −
In the common-gate (CG) configuration, the input signal is
applied to the source terminal, the output is sampled from the
drain terminal, while the gate terminal is held at a constant
bias voltage. In the small-signal equivalent model, the gate vout
Rout = Rdrain = RS + ro + gm ro RS .
Figure 127: NMOS Common-Gate amplifier config-
uration with ideal active bias, and its small-signal
The input resistance is the resistance seen looking into the
equivalent circuit model. The signal source is as-
source terminal. Since there is an ideal current source connected sumed to have a series resistance of RS .
above the drain, the effective resistance above the drain is
Rd → ∞, so for this configuration,
Rin = Rsource = ∞.
Consider a CG amplifier with passive-bias where R D = 50 kΩ, RS = 10 kΩ, VG = 2.5 V and VIN =
1.6 V. The NMOS device parameters are k = 500 µA/V2 , VTh = 0.5 V and λ = 0.05 V−1 , and the supply
voltage is VDD = 5 V. What gain and output resistance will be achieved in the series configuration (??)
and the bypass configuration (??)?
To begin with, we solve the DC operating point, then the small-signal parameters, and the open-
circuit characteristics of this amplifier:
VOV = VG − VS − VTh
VS = VIN + ID RS
1 2 1
ID = kV = k (VG − VTh − VIN − ID RS )
2 OV 2
√
By defining x = ID , we can express a quadratic equation:
r r
k 2 k
R x +x− (V − VTh − VIN ) = 0
2 2 G
then solve using the quadratic formula:
p !2
−1 + 1 + 2kR (VG − VTh − VIN )
ID = √ = 15.3 µA
R 2k
From this point, the two circuits will diverge in the value of Rdrain at higher frequencies. In the by-
pass configuration, RS is masked for AC signals, so Rdrain = ro . For the series configuration, however,
RS has a big effect:
This demonstration shows a basic common-gate amplifier with a passive bias resistor R D = 50 kΩ.
The NMOS device has the familiar characteristics: k = 500 µA/V2 , VTh = 0.5 V and λ = 0.05 V−1 . The
supply voltage is VDD = 5 V, the NMOS gate is biased at a constant VG = 2 V, and the input signal has
a DC offset of VIN = 1.1 V, and the input AC small signal vin has an amplitude of 100 mV. The input
signal has zero series resistance.
To determine the DC operating point and small-signal characteristics, we can start by directly calculat-
ing VOV since there is no resistor below the source terminal:
VOV = VG − VS − VTh
= VG − VIN − VTh = 0.4 V
Then the DC bias current, output offset, and small-signal parameters are:
1 2
ID = kV = 40 µA
2 OV
VD = VDD − ID R D = 3 V
VDS = VD − VIN = 1.9 V > VOV
gm = k VOV = 200 µA/V
ro = (λID )−1 = 500 kΩ
Now measure the amplitude of the output signal in EveryCircuit, and verify that the gain is a little
over 9 V/V, as predicted by our analysis.
As an exercise, try removing R D and replace it with an ideal current source (pointing down) sup-
plying 40 µA. Predict the effect this will have on the circuit, and verify your prediction in the Ev-
eryCircuit simulation (note: for this exercise you will need to reduce the amplitude of vin to 1 mV, and
carefully increase the DC offset to VIN = 1.11 V).
144 ece3410 lecture notes
This circuit implements the passive-bias CG configuration described in example 24. The bypass
version is shown. Run the transient simulation and verify that the results align with the predictions
from example 24. As an exercise, remove the bypass capacitor and reposition the AC input in series
between the DC offset and RS . Repeat the simulation, and verify that the gain decreases to a value
close to what was predicted.
introduction to mosfets 145
If the input signal is applied to the gate while the output is sam- vIN
gm r o
Avo = . vout
1 + gm r o
The output resistance is a little more subtle. Since the resistance vin − vout ro
gm r o gm R S
AvL =
1 + gm r o 1 + gm R S
gm R S
≈
1 + gm R S
Suppose a CS configuration has an open-circuit gain of Avo = 20 V/V an output resistance Rout =
500 kΩ, and needs to drive a load R L = 10 kΩ. If the CS amplifier is connected directly to the load,
the gain will be attenuated so that AvL = 20 × 10/(10 + 500) = 0.392 V/V. Now suppose a SF
configuration is inserted in between the CS output and the load resistor, and the SF circuit has gm =
200 µA/V, ro = 500 kΩ, and the load itself acts as a passive bias resistance of RS = 10 kΩ. Then the
overall gain of this two-stage circuit will be
gm r o gm R L
AvL = 20 V/V = 19.99 V/V,
1 + gm r o 1 + gm R L
vIN RL vIN
RL
introduction to mosfets 147
load
passive-biased CS amplifier may deliver an output amplitude ID
no greater than 1 V.
amplifier
device. The diode connection regulates vGS to support the forced
current ID . On the output side, the device can be used as a
current source that delivers iout = ID to a load connected at its + +
vGS 1 vGS 1
drain terminal. − −
Current mirrors can be used to generate bias currents, and
to source multiple copies of a reference current. This can be Figure 133: Example showing two current mirrors. In
the NMOS mirror, the reference current is setup by
quite useful for amplifier biasing. In the structure shown in the interaction between R and the diode-connected
Figure 133, the reference current (i.e. the input) is initially NMOS device. The PMOS mirror then generates a
determined by a resistor R. The current in this configuration new copy of the reference current, which is used to
bias an amplifier.
was previously found to be
p 2
−1 + 1 + 2kR (VDD − VTh )
ID = .
2kR2
148 ece3410 lecture notes
Consider the CS-SF configuration from the previous example. By using current mirrors, we can
convert this to an active-bias configuration, replacing the resistors with MOSFETs as shown below.
In the expanded circuit below, current mirroring is used to generate both the VIN offset voltage and
the current sources for ID . There are two current mirrors, an NMOS mirror producing three copies of
ID via vGS 1 , and a PMOS mirror producing one copy of ID via |vGS 2 |.
vin
introduction to mosfets 149
how can we control this? The answer lies with symmetry. Rref
vY vX
In the two-mirror bias structure shown in Figure 134, the
NMOS mirror is used to set the gate offset voltage, and the
PMOS mirror is used to source the active bias current. For this
setup to work, all devices must be perfectly matched, i.e. they
must have the same physical parameters (k, VTh , λ, etc) and the MCS
+ +
same geometry (W and L). They must furthermore all operate vGS 1 vGS 1
− −
in saturation, so that the sensitivity to their drain voltages is
minimized. Then they should all have the same device current,
Figure 134: Two current mirrors used to bias a
ID . common-source amplifier. By symmetry, we can infer
Using the square-law device equation, we can solve for that v X = vY .
Suppose the circuit of Figure 134 is constructed with the following parameters: k n = 5 mA/V2 , k p =
3 mA/V2 , VThN = VThP = 0.5 V, λn = 0.01 V−1 , λ p = 0.05 V−1 , Rref = 50 kΩ and VDD = 5 V. What is the
circuit’s complete DC operating point?
Using the previous analysis of the diode-connected MOSFET, we calculate ID and then vGS 1 and vGS 2
p 2
−1 + 1 + 2kR (VDD − VTh )
ID = = 86.3 µA
2kR2 s
s 2ID
as 2ID VOV 2 = = 0.24 V
VOV 1 = = 0.186 V kp
kn
|VGS 2 | = VOV 2 + VThP = 0.74 V
VGS 1 = VOV 1 + VThN = 0.686 V
In the result from Example 27, notice that the output offset
is very close to the maximum output voltage. It is near the top
of its range. That means the positive leg of an output signal
will be clipped. To resolve this problem, we need to break the
symmetry by a small amount. One option is to slightly increase
the k of device MCS (by increasing its width), so that its device
current is slightly greater than ID . That will pull down the
output offset. A second option is to slightly decrease the gate
bias voltage at MCS . Both of these methods are risky, since it
can be challenging to calculate the exact variation required. Real
MOSFETs may deviate slightly from our model equations, and
manufacturing variations can result in physical parameters that
are slightly different from the ones on the data sheet.
This demonstration implements the current-mirror bias network from Example 26. In order to correct
for the output bias problem discussed in Example 27, the width of device MCS is slightly increased.
Based on the parameters from Example 27, we found that the bias current is ID = 86.3 µA. Continu-
ing this analysis, we find the small-signal parameters are
p
gm = 2k n ID = 929 µA/V
ro,n = (λn ID )−1 = 231.7 kΩ
ro,p = (λ p ID )−1 = 115.9 kΩ
Feedback biasing
In order to achieve a more reliable bias solution that is highly −
LPF
each other, so
1 2 1
k p VDD − Vreg − VThP = k p (v Z − v X − VThP )2
2 2 s !
2ID
⇒ VDD − Vreg − VThP = VDD − VThP − − v X − VThP .
kp
This example modifies the design from Example 27 to use feedback bias with an ideal op amp as the
error amplifier, like the solution shown in Figure 135. In this configuration, all devices are matched,
so there is no need to manipulate the width of MCS .
This example modifies the design from Example 27 to use feedback bias with an error amplifier like
the one from Figure 136. In this design, the control voltage is Vreg = 3.0 V, corresponding to VX? + VThP .
All NMOS devices are matched, and the error amplifier feedback is achieved without using an op
amp.
introduction to mosfets 153
It will often be the case that one pole is much larger than the
other. In that case may consider the smaller pole to be domi-
nant, and treat the circuit as a one-pole system governed by the
dominant pole.
154 ece3410 lecture notes
| H ( f )| (dB)
312.5 kΩ
tradeoff is clear:
RL
AvL = Avo
R L + Rout
RL
≈ Avo ID
Rout
f c = [2π ( Rout k R L )Cout ]−1 vOUT
−1
≈ [2πR L Cout ]
vIN Cout
We see that there is a one-to-one exchange between gain and
bandwidth. In other words, the gain-bandwidth product is
constant: GBW = Avo /(2πRout Cout ).
For a larger load, or when no resistive load is present, we Figure 141: Active-bias CS configuration, when the
output pole is dominant,
√ shows a strong gain/BW
also see the tradeoff as a property of the bias current. The −1
tradeoff: AV ∝ ID , whereas f c ∝ ID .
amplifier’s open-circuit gain is
Avo = − gm ro
p 1
=− 2kID
λID
s
2k
=−
λ2 ID R
| H ( f )| (dB)
10
The table below summarizes the results. From the data, we
see that the prediction is more accurate for larger R. The netlists 5 kΩ
used to simulate the input-dominant frequency response and 0
the gain-bandwidth tradeoff are shown after the data table.
−10
101 103 105 107 109
Predicted Simulated Frequency (Hz)
R (kΩ) ID (µA) Av (dB) f c (MHz) ID (µA) Av (dB) f c (MHz)
Figure 143: Simulated transfer functions for the CS
5 400 7 31.8 223 8 32
configuration from Figure 142. The NMOS device
10 200 10 15.9 131 11.6 15.7 has k = 500 µA/V2 , VTh = 0.5 V, VDD = 5 V and
25 80 14 6.4 61 16 12.7 VIN = 2.5 V. R is varied from 5 kΩ up to 100 kΩ.
50 40 17 3.18 33 19.5 3.7
100 20 20 1.59 17.5 22.7 1.86
VDD ndd 0 DC 5V
VIN ndc 0 DC 2.5V
vsig nsig 0 DC=0 AC=1 SIN(0 0.01 1k)
RL nout 0 1Meg
CC2 nd nout 10uF
Cout nd 0 1fF
VDD ndd 0 DC 5V
VIN ndc 0 DC 2.5V
vsig nsig 0 DC=0 AC=1 SIN(0 0.01 1k)
RL nout 0 1Meg
CC2 nd nout 10uF
Cout nd 0 1pF
vB
iB +
v BE i E = i B + iC
The bipolar junction transistor was the first practical transistor −
vE
device for mass production, and defined the semiconductor
Emitter
industry from the 1950s into the 1980s. Today, BJTs are not as
widespread as MOSFETs, but are still very important for niche Figure 144: The NPN BJT device, showing its be-
havior in the active mode (when v BE ≈ 0.7 V and
applications. Some areas where BJTs excel include high-voltage vCE > 0.3 V). In this mode, the device can be de-
applications, and radio-frequency power amplifiers, where scribed as a current amplifier with constant gain
BJTs are able to drive antennas and transmission lines with β.
very good linearity and, hence, low distortion. BJTs also have a
Emitter
high transconductance, and it is usually easier to make a good vE
+
discrete BJT amplifier, whereas MOSFETs may need several i E = i B + iC
v BE
devices in order to achieve a good bias configuration, making
them more suited for integrated circuit designs. −
The chief drawbacks to BJT devices are high power consump- iB vB
Base
tion (BJT bias currents must usually be on the order of 1 mA
iC = βi B = αi E
to 100 mA), and comparatively high voltage overhead (unlike
vC
MOSFETs, the BJT’s overdrive voltage is not adjustable and
Collector
cannot be reduced for low-voltage applications). In addition,
BJT devices cannot be miniaturized to nano-scale dimensions, Figure 145: The PNP BJT device, showing its be-
havior in the active mode (when v EB ≈ 0.7 V and
so they cannot achieve the same performance enhancements or v EC > 0.3 V). The PNP device’s behavior is comple-
cost improvements that come with MOSFET scaling. Lastly, BJT mentary to the NPN.
devices pass current through their base terminals (comparable
to the MOSFET’s gate), which makes them inefficient for logic
circuits, and complicates amplifier analysis and design.
BJTs are built out of PN junctions (diodes), and normally
have three operating modes corresponding to the diode states: Active Mode Summary:
• cuttoff: both junctions are not forward biased, i.e. v BE < 0.4 V
v BE
i B = IS exp
and vCB > −0.4. Note that the junctions do not have to nUT
Note that these conditions imply that vCE > 0.3 V. In this
mode, the base-current is defined by the forward-bias diode
equation, i B = IS exp (v BE /nUT ), and the collector current
is iC = βi B , where β is the device’s current gain, n is the
forward emission coefficient (usually close to 1.0), UT is
the thermal voltage (26 mV at room temperature), and IS
is a scale current on the order of pA. This mode roughly
corresponds to the MOSFET’s saturation mode, and is the
appropriate DC bias mode for amplifier circuits. Collector Emitter
VCC = IE R E + v BE + IB R B + IE RC
IE R B
⇒ VCC = IE R E + 0.7 V + + IE R C
β+1
β+1
⇒ RB = (VCC − 0.7 V − IE ( R E + RC ))
IE
162 ece3410 lecture notes
A voltage-divider network like the one in Figure 147 has VCC = 5 V, RC = 2 kΩ, R E = 0.5 kΩ,
R B1 = 10 kΩ, and the device has β = 100 A/A. If the desired emitter current is 1 mA, what is the cor-
rect value for R B2 ?
Based on the given parameters, we can directly calculate VE = IE R E = 0.5 V, and VB = VE + 0.7 V =
1.2 V. Then
R B2
VCC = 1.2 V
R B1 + R B2
⇒ VCC R B2 = 1.2 V ( R B1 + R B2 )
1.2 V
⇒ R B2 = R B1
VCC − 1.2 V
= 3.16 kΩ
A voltage-divider network like the one in Figure 148 has VCC = 5 V, RC = 3 kΩ, R E = 0.5 kΩ, and the
device has β = 100 A/A. If the desired emitter current is 1 mA, what is the correct value for R B ?
As in Example 28, we can directly calculate VE = IE R E = 0.5 V, and VB = VE + 0.7 V = 1.2 V. Then
β+1
RB = (VCC − 0.7 V − IE ( R E + RC ))
IE
= 80.8 kΩ.
vout R E
vbe ≈ vin − ie R E = vin + RC Rcoll.
RC
vout R E ie ≈ ic
⇒ vout ≈ − gm RC vin + RE
RC
vout − gm R C
⇒ ≈
vin 1 + gm R E
R Figure 150: Common-Emitter configuration based on
≈ − C. the voltage-divider bias network, and its simplified
RE
small-signal model.
In the feedback-biased case, we may reach a similar conclu-
sion if the value of R B is much larger than RC and also much VCC
larger than Rsig . In that case, the negative feedback loop created
by R B will have little impact on the small-signal analysis, and
we arrive at the same result for the gain. RC
CC2
RB
vOUT
VC
VB
vIN
CC1
VE
RE
Using the bias network from Example 29, we introduce input and output signals using the capacitive-
coupled connections shown in Figure 151. A bypass capacitor is used to eliminate the AC influence of
R E , so the AC gain should be − gm R E . Since IC ≈ 1 mA, the transconductance and gain should be
1 mA
gm = = 38.46 mA/V
26 mV
⇒ − gm RC = −115 V/V
The simulation verifies a gain of 104.5 V/V, which is close to our prediction. The small discrepancy is
due to the assumptions and approximations made in our analysis.
Basic Electronic Device Theory
Si C
distribution equation:
1
f ( E) = .
1 + exp Ek− ETF
B
This equation is also called the Fermi function for short. The
Fermi function gives the probability that an energy state S at
energy E is occupied by an electron. We could say the Fermi
function has units of electrons per state. Some common bandgap energies (eV):
The Fermi function is completely general, and applies to Material Eg
Silicon (Si) 1.12
any solid material in thermal equilibrium. To understand the Gallium Arsenide (GaAs) 1.424
distribution of charges in a specific material, the Fermi function Germanium (Ge) 0.664
Indium Antimonide (InSb) 0.230
must be joined together with a function for the material’s
Indium Phosphide (InP) 1.344
Density of States. In a semiconductor or insulator, not every Indium Arsenide (InAs) 0.354
energy level is allowed. For example, there are no allowed states Zinc Oxide (ZnO) 3.3
Zinc selenide (ZnSe) 2.822
within the forbidden gap, as illustrated in the figure below.
Energy states within the gap are excluded; they are “chopped
out” from the distribution. Then the amount of mobile charge is
determined by the distribution’s “tail” , the amount surviving
in the conduction band.
E E E
−
−
− − −
− − −
EC
− −
− − − − − − EF − − − − − − EF EF
− − − − − − − − − −
− − − − − − − − − − − − − − − − − − EV
− − − − − − − − − − − − − − − − − −
− − − − − − − − − − − − − − − − − −
− − − − − − − − − − − − − − − − − −
− − − − − − − − − − − − − − − −
f ( E) f ( E) f ( E)
T = 0 K (absolute zero) T ≈ 300 K (room temperature) T ≈ 300 K (room temperature)
known to be
3/2
NC = 2 2π η me k B T/h2
where η is a material-dependent constant, me is the electron rest Physical constants for electrons (e− ):
mass in free space, and h is Planck’s constant. The effective mass
η me is an adjustment that accounts for a variety of material- me = 9.1 × 10−31 kg
specific effects, such as bond lattice geometry and random h = 6.626 070 04 × 10−34 m2 kg/s
collisions with other particles. In general, a lower effective mass k B = 1.380 648 52 × 10−23 m2 kgs−2 K−1
indicates that the material can be used to make faster devices.
Material η
If we consider “number of electrons” as a physical unit, then
Silicon (Si) 1.09
we can interpret the Fermi function as “number of electrons Gallium Arsenide (GaAs) 0.067
per state” in the material. The density of states reveals the Germanium (Ge) 0.55
Indium Antimonide (InSb) 0.013
“number of states per energy level per volume.” Hence when Indium Phosphide (InP) 0.08
we multiply N ( E) and f ( E) together, the resulting units are Indium Arsenide (InAs) 0.023
Zinc Oxide (ZnO) 0.29
Zinc selenide (ZnSe)
e− states e− 0.17
× = ,
state energy level × volume energy level × volume
If EC − EF > 4k B T, then we may further simplify the result: Physical constants for holes (h+ ):
Material ηh
E − EF Silicon (Si)
0.55
n ≈ NC exp − C . Gallium Arsenide (GaAs) 0.45
kB T Germanium (Ge) 0.37
Indium Antimonide (InSb) 0.4
When an e− transitions into the conduction band, it leaves Indium Phosphide (InP) 0.64
behind a hole (h+ ) in the valence band. Holes are able to move Indium Arsenide (InAs) 0.4
majority minority
n2i
N-type (donors) n = ND p= ND
n2i
P-type (acceptors) p = NA n= NA
Now let’s take a closer look at the P-type material illustrated positive mobile carriers h+
in Figure 161. The sum of all charges must balance to zero at negative mobile carriers e−
equilibrium (i.e. when no external energy is applied). Athough negative fixed charges
the P-type material has fewer mobile e− , there are correspond-
Figure 161: Extrinsic P-type single-crystal Si doped
ingly more fixed negative charges at the dopant sites, so the with acceptors, yielding a surplus of mobile h+
material has zero net charge. The situation is the same for balanced by fixed negative charges at the acceptor
sites.
N-type material, where the charge polarities are swapped.
basic electronic device theory 173
A typical N-type silicon material might have a donor density of ND = 1016 dopants per cm3 . In this
material, at thermal equilibrium (with no externally applied potential and no net current flow), the
equilibrium charge concentrations will be
dopants
nn ≈ ND = 1016
cm3
n2i
pn ≈
nn
2
1.5 × 1010 carriers per cm3
=
ND = 1016 carriers per cm3
= 2.25 × 104 carriers per cm3 .
174 ece3410 lecture notes
EF − Ei
n = ni exp
kT −
−
Ei − EF
−
p = ni exp − −
EC
kT EF
− − − − − − EV
− − − − − −
− − − − − −
− − − − − −
EF
− − − EV
− − − − −
− − − − − −
ND − − − − − −
(N-type) EF = Ei + kT ln − − − − − −
ni f ( E)
NA P-Type Material with Acceptors
(P-type) EF = Ei − kT ln
ni
Figure 162: How doping affects the Fermi level. Note
the presence of conduction-band electrons in the
N-type case, and the presence of valence-band holes
(empty spots) in the P-type case.
basic electronic device theory 175
Part A: A piece of single-crystal Si is doped with Phosphorous (a donor species) with a dopant den-
sity of ND = 1016 dopants per cm3 . Calculate the concentration of majority carriers (electrons) and
minority carriers (holes) in this material at room temperature. Also calculate the change in Fermi
level due to doping, and compare it to the bandgap.
Solution A: For moderately high doping concentration (greater than 1010 dopants per cm3 ), the den-
sity of majority carriers approximately equals the doping concentration. The mass-action law dictates
that n × p = n2i , and for Si at room temperature, ni ≈ 1.45 × 1010 carriers per cm3 . Then the densities
are
n ≈ ND = 1 × 1016 carriers/cm3
n2i
p= = 2.1 × 104 carriers per cm3
n
Compared to intrinsic Si, the Fermi level is moved by this amount:
ND
∆EF = EF − Ei = kT ln
ni
= (25.7 meV) ln 6.5 × 105
= 345.5 meV
Since Si has a bandgap of Eg = 1.12 eV, the Fermi level is shifted toward the conduction band edge
by 0.3455
1.12 , or about 30% of the band gap. Part B: A piece of single-crystal Si is doped with Boron (an
acceptor species) with a dopant density of NA = 1016 dopants per cm3 . Calculate the concentration of
majority carriers (holes) and minority carriers (electrons) in this material at room temperature. Also
calculate the change in Fermi level due to doping, and compare it to the bandgap.
Solution B: This problem is completely symmetric to Part A, only the carrier’s polarity is reversed.
The majority carriers are now holes, and the Fermi level will be moved downward toward the valence
band. The calculations are otherwise identical:
p ≈ NA
n2i
n= = 2.1 × 104 carriers per cm3
p
NA
∆EF = EF − Ei = −kT ln
ni
= − (25.7 meV) ln 6.5 × 105
= −345.5 meV
So in this case the Fermi level is shifted downward by about 30% of the band gap.
176 ece3410 lecture notes
Resistivity
The drift current mechanism supplies the electronic theory
underlying the classical lineary property of resistivity. If the
material is homogeneous, meaning it has the same composition
+
−
Idrift = Acs qE µ p p + µn n
V x (position)
= Acs qµ p p + µn n
L
V A 1 Figure 166: Energy bands are tilted in the presence
⇒R= = CS of an E -field. Electrons tend to “roll” down the
Idrift L q µ p p + µn n
slope from higher energies toward lower energies.
Finally we can see that the expression has the same form as Note that the same process will cause holes to be
transferred up the hill. Since a current is flowing, the
classical resistivity, R = ρACS /L, where the ρ is the material’s
system is not in equilibrium. The Fermi level is only
resistivity defined as well-defined for systems in equilibrium, so there is
no EF here.
1
ρ= .
q µ p p + µn n
Solution: First note that the electric field has units of V/cm. Since 1 V is applied across the “long
dimension” of 1 cm, the electric field should be E = 1 V/cm. Since the material is P-type, the majority
carriers are h+ . We can ignore the negative carriers with little consequence. Then the h+ velocity and
drift current density are
Lastly, the cross-sectional area is 1 mm × 1 mm. Since material units are most often given in terms of
centimeters, we must take care to convert all dimensions to centimeters and ensure that the units
properly balance. Then the cross-sectional area and total drift current are:
∆N ∆x
φdiff = × (1)
2 ∆t
!
∆N (∆x )2
= × (2)
∆x 2∆t
= ∇ N (x)D (3)
∂N ( x )
∇ N (x) , .
∂x
This tells us that the diffusion current is proportional to the
concentration gradient; particles flow from regions of high
concentration toward regions of lower concentration.
The second new concept is the diffusivity D, with units
cm2 /s. Diffusivity indicates the speed with which particles
move under thermal excitation. It is a property of Brownian
basic electronic device theory 179
A piece of N-type Si has a “built-in” uniform E -field of 1 V cm−1 and a cross-sectional area of
0.01 cm2 . At a point x = 0 along the direction of current flow, the electron concentration is
n (0) = 16 − 3
10 e /cm . If the material is in thermal equilibrium, solve an expression for n( x ). As-
sume that the electron mobility is constant at µn = 1200 cm2 /Vs, and the thermal voltage is
UT = k B T/q = 27 mV. Furthermore assume that the minority carriers (holes) are negligible. Lastly,
calculate the drift and diffusion currents at a position of x = 1 µm, and state the electron concen-
tration at that position. How can the material be in equilibrium if an electric field and currents are
present?
Solution: Since the material is in equilibrium, the drift and diffusion currents must cancel each other
out. Therefore:
qE µn n = q∇nDn
µn
⇒ ∇n = nE
Dn
E
=n
UT
A differential equation of this form has only one non-trivial solution: n( x ) must be exponential with
respect to x:
Ex
n( x ) = n(0) exp .
UT
At the position x = 1 µm = 0.001 cm, the electron concentration is n( x ) = 1.037 73 × 1016 e− /cm3 .
Then the drift current is 19.95 mA and the diffusion current is −19.95 mA. In spite of all this activity,
the material is in equilibtrium since the net current from all contributors balances out to zero.
basic electronic device theory 181
Neutral bias
P xp xn N
E
− +
∆E = kT ln ( NA /ni ) + kT ln ( ND /ni )
!
NA ND
= kT ln [eV]
n2i
W = xn + x p ,
ND
x p = xn
N
A
ND
⇒ W = xn 1 +
NA
W NA
⇒ xn =
NA + ND
Now we can substitute for xn into the solution for V0 :
qW ND W NA
V0 =
2e ND + NA
2 1 NA ND
=W
2e NA + ND
s
2e 1 1
⇒W= + V0
q NA ND
W ≈ x p (one sided, ND NA ).
Reverse Bias
1
W Emax
V0 + VR =
2
s
2e 1 1
⇒W= + (V0 + VR )
q NA ND
++++
2 (V0 + VR )
Emax =
W − − − − EC
From this we see that the E -field strength increases with VR1/2 ,
which enhances the opposition to current flow. The energy band
interpretation is illustrated in Figure 171.
EV
Leakage Current
Then the junction capacitance is Cj0 = ACS e/W = 86.4 fF. When a 2 V reverse bias is applied:
Cj0
Cj = 1/2 = 47.2 fF,
VR
1+ V0
Reverse Breakdown
Avalanche Breakdown
2 (V0 + VBR )
Ecrit =
s W
2e 1 1
W= + (V0 + VBR )
q NA ND
2 (V0 + VBR )
⇒ Ecrit = r
2e 1 1
q NA + ND (V0 + VBR )
2
qEcrit
1 1
⇒ VBR = + − V0 .
2e NA ND
Example 35.
Avalanche breakdown voltage. Suppose a juntion is doped with NA = 1015 dopants/cm3 on the P
side and ND = 1017 dopants/cm3 on the N side. The reverse breakdown potential is found to be
VBR = 300 V. Now suppose another junction is fabricated with NA increased to 101 6 dopants/cm3 .
How will VBR change in the new junction?
Solution: In the first case, the built-in potential is V0 = UT ln NAnN
2
D
= 0.725 97 V. Let V00 and VBR
0 be
i
the built-in potential and breakdown potential for the second case, respectively. Since the doping on
the P side is now equal to NA0 = 1016 dopants/cm3 , we find that V 0 = 0.788 14 V Then, since V + V
0 0 BR
1 1
is proportional to N + ND , we can say that
A
1 1
V00 + VBR
0
NA0 + ND
= 1 1
V0 + VBR NA + ND
1 1
NA0 + ND
0 − V00
⇒ VBR = (V0 + VBR ) 1 1
NA + ND
= 31.96 V
P ←− E N
Tunneling Breakdown
When the doping concentration is especially high, the built-in
potential increases and the depletion width decreases. This re-
sults in a steep energy band transition, as shown in Figure 173.
At the level where Ec aligns with Ev , separated by a small
distance L, it is possible for electrons to undergo quantum me- −−−−
chanical tunneling across the barrier. Quantum theory reveals ++++
that the tunneling current density is given by EC
!
4L 2m∗ Eg
p
Jtunnel = qv R n exp − ,
3h̄
·10−4
In practical terms, the main difference between tunneling and
avalanche mechanisms is in their response to temperature. In
5
each case, temperature effects are measured via the temperature
coefficient parameter, defined as
dVBR
TC , ,
dT 0
A Zener diode with a nominal breakdown voltage of 4.5 V has a temperature coefficient TC =
−0.3 mV/K. If the diode’s temperature rises to 400 K, how will that alter the breakdown voltage?
Solution: Since the nominal breakdown voltage is defined for room temperature, 300 K, the change
in temperature from nominal conditions is ∆T = 100 K. Then the change in breakdown voltage is
∆VBR = TC × ∆T = −40 mV. Then the actual breakdown voltage should shift to 4.46 V.
192 ece3410 lecture notes
Forward Bias
of the Junction:
VF
np = n2i exp .
UT
Since the majority concentration is set by the dopant concen-
tration, the increased concentration appears mainly in the
minority carriers. As a very rough approximation, we suppose
those new carriers appear as an impulse function at the very
edge of depletion, so the minority carrier gradient is equal to
the impulse height. Minority carriers will diffuse due to this
sharp gradient in a process called minority carrier injection.
Then the diffusion current should be
VF
I = IS exp ,
UT
where IS is a scale constant. It so happens that IS is equal to the
reverse bias leakage current. A full analysis of generation and
recombination reveals the scale current to be
Dn n p0 D p pn0
IS = qACS + ,
Ln Lp
where Ln and L p are the mean diffusion length parameters for
electrons and holes, respectively. The diffusion length refers
to how far a particle diffuses, on average, before it recombines.
The diffusion lengths are obtainable from the diffusivity and the
average carrier lifetime τ:
√
Ln = Dn τn
q
L p = D p τp ,
dQ
Cdiff =
dV
d
= Iτ
dV
I
= D τ,
UT
hence the capacitance is directly proportional to the forward-
bias current ID . Note that in a one-sided junction, where the
depletion region extends mostly into the weakly doped side
(usually that’s the P side), the characteristics are dominated by
the weaker doping. In that case, the average lifetime is τ = τp .
Recombination Current
In strong forward bias, the diode’s current is primarily driven
by minority carrier diffusion. Parameters like IS and n are
defined to accurately model the forward bias current, but they
don’t always provide the best model for weak forward bias
or reverse bias. In these cases, a much smaller drift current is
present, caused by thermal generation of mobile charge pairs
in the depletion region. A small amount of mobile pairs appear
the pn junction 195
High-Injection Effects
When the junction is placed in a sufficiently strong foward bias,
the minority carrier concentrations exceed the material’s doping
concentration. In order to preserve net charge neutrality in the
quasi-neutral regions, the concentration of majority carriers
must also increase. This sets up a process of diffusion and
recombination that partially cancels out the minority diffusion
current, hence the diode’s current starts to flatten out with
increasing voltage. This effect is well modeled by setting the
ideality factor to n = 2 when the forward current exceeds a
critical threshold called the high-injection knee current, IKF .
The modified diode expression is
vD
i D ≈ IS exp (when i D > IKF ).
2UT
196 ece3410 lecture notes
Temperature Dependence
The ideal diode model contains UT in its exponent, giving it a
strong temperature dependence. We have also seen that mobil-
ity, diffusivity and the minority carrier concentrations are all
sensitive to temperature, and they all contribute to the diode’s
scale current. All of these combined effects can overpower the
contribution of UT in the exponent. In most datasheets, the
scale current IS is measured at room temperature T0 . If the tem-
perature (in K) is changed by ∆T, the resulting change in IS is
approximately given by
IS ( T0 + ∆T ) ≈ IS ( T0 ) exp ( xti ∆T ) ,
Suppose a diode has ideal characteristics IS = 6 nA and n = 1.9 measured at room temperature, T0 =
300 K. In addition, it has non-ideal characteristics ISR = 11 nA, n R = 5, RS = 0.33 Ω, IKF = 43 mA and
xti = 4.0. Consider the following cases.
VF= 0.8 V, calculate the diode’s current using both the ideal forward bias
(A) In forward bias, with
vF
model, i D ≈ IS exp nU T
, and again using the complete model including recombination current.
How do they compare? Solution: Computing the two expressions yields the approximate result of
i D ≈ 1.112 98 mA and the more accurate model gives i D = 1.113 08 mA. In forward bias, the difference
is about 0.09%. We may conclude that the simplified exponential formula is sufficiently accurate for
“normal” forward bias conditions.
(B) In what range of voltages does the recombination current become a significant influence? Solution:
Recombination is significant when
vD vD
IS exp ≤ ISR exp
nUT n R UT
1 1 I
⇒ vD − ≤ ln SR
nUT n R UT IS
nn R UT ln IISR
S
⇒ vF ≤
nR − n
= 0.048 V
It’s easily verified that at this voltage, both the diffusion current and recombination current are
around 16 nA. For voltages less than this and for reverse bias, recombination will play a significant
role.
(C) What will be the reverse leakage current? Solution: In reverse bias, the total leakage current
should be sum of diffusion and recombination currents, IS + ISR = 17 nA.
(D) What forward voltage is required to exceed the high-injection knee current? Solution: This is a
more challenging question. The knee current is IKF = 43 mA, and using the exponential forward-
bias model we can solve for the forward voltage as VKF ≈ nUT ln ( IKF /IS ) = 0.779 78 V. But if
the current is 43 mA then we should see a combined voltage drop in the quasi-neutral regions of
2RS IKF = 0.0258 V, which reduces the forward voltage and drops the current down to 25.5 mA. But
the lesser current also reduces the voltage drop across RS ... To get the real answer we need to iterate.
After iterating a few times, the solution converges to 29.9 mA when VF = 0.779 78 V, which does not
actually surpass the knee current. Using the iterative method, we can search by increasing v F until
the iterated solution exceeds IKF . This occurs when v F ≈ 0.805 V.
198 ece3410 lecture notes
* DC Simulation:
.control
DC Vin -1.0 1.25 0.001
plot ylog ylimit 1e-10 10 abs(i(Vin))
.endc
.end
200 ece3410 lecture notes
100
Diffusion only
10−1 Diffusion and Recombination High-Injection
10−2
10−3
as
Bi
10−4
| i D | [A]
d
ar
w
n)
r
io
Fo
10−5
us
iff
(d
Recombination
10−6
10−7
Reverse Bias
10−8
10−9
−1.2 −1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4
v D [V]
P depletion N
Photons in the Depletion Region
When a photon is absorbed in the quasi-neutral bulk regions of
T
a PN diode, it adds to the pool of mobile carriers and slightly al- W
ters the conductivity, which is a fairly weak interaction. When a Figure 178: Photons pass through the top surface of
photon is absorbed in the depletion region, the result is more in- this junction. They are able to be usefully collected in
teresting: a new pair of mobile carriers is generated, which are the depletion region, which has width W. The junc-
tion has thickness T, so the total area for collecting
immediately swept apart by the electric field, setting up a net photons is W × T.
reverse current through the junction. This current is called the
photo current, and is directly proportional to the rate of photon
absorption in the depletion region. The photon flux Φ is de-
fined as the rate of photons of a particular wavelength passing
through a material section, usually with units ofphotons/scm2 .
As illustrated in Figure 178, photons incident on the surface
plane are “eligible” to be absorbed as they pass through the
bulk. If the junction thickness is T and the depletion width is W,
then the photocurrent should be
Iph = qηΦTW,
A/W, given by
qλη
Rλ = .
hc
the pn junction 203
Photodetector Circuits
A photo detector can be made from a single photodiode and a Iph
biasing circuit. In Figure 179, a resistor is used to achieve the
reverse bias. Since the photocurrent is directly proportional
to the incident optical power, and the resistor’s voltage drop vOUT
is proportional to the photocurrent, overall this is a linear
photodetector. Linear detectors are not necessarily the best
option for all applications. Considering that the difference in
irradiance between outdoor and indoor lighting can span four
orders of magnitude, it is not feasible to design a passive linear
photodetector for general use. Figure 179: A linear resistor-diode photodetector.
The resistor and power supply are used to setup a
An alternative design is the active pixel sensor (APS) shown
reverse bias condition in the diode, with an output
in Figure 180. In this design, a MOSFET switch (M1) is first voltage close to VDD . When light is present the
switched ON to drive node v x up to VDD − VTh , placing the diode photocurrent induces a drop across R, giving an
output vOUT = VDD − Iph R.
in reverse bias. Then M1 is switched OFF, leaving v x floating on
the capacitor CP . The photocurrent discharges CP , causing v x to
drop in proportion to the time integral of Iph . The output signal
is isolated by a source-follower (M2). Suppose M1 is turned M1
rst
OFF at time 0, and the signal at vOUT is sampled after a short
time ∆t. Then the change in output voltage should be
Z ∆t M2
∆vOUT = − Iph dt. Iph
0 CP vOUT
out circuit similar to the one shown in Figure 181. In this circuit, CS
This example circuit models an active pixel sensor with correlated dobule sampling. Since EveryCir-
cuit doesn’t have any models for photonic signals, we simulate the photocurrent by superimposing
a parallel current source with the reverse-biased diode. There are three different MOSFET configu-
rations here: switching mode (in red), amplifiers (blue) and bias transistors (green). The reset switch
is used to initialize the charge on capacitor CP . We don’t need to know the precise initial value, just
that it is the same after each reset. The sampler switches are used with capacitors as track-and-hold
(T/H) circuits: when the switches are ON, the capacitors are connected to their respective signals.
When the switches are OFF, the capacitors are floating, and the remember whatever signal value was
present at the time of switch-off. The amplifier configurations are Common-Source (CS) gain stages
and Source-Follower (SF) buffers. We previously studied all of these configurations, but this is one of
our first examples showing a system with multiple configurations working together.
CS Amp. CS Amp.
Diff. Sampler
Vbias M7 M9
Reset
+ ∆vOUT −
SF Buffer φ2
S/E Sampler
M10 M11
φ1 M1
VP CO
M6
M2 M4 φ4 φ4
CS
Iph
M3
VSF M5
M8
Photodiode Model
CE
φ3
the pn junction 205
∆t
M1 ON
φ1
M4 ON
φ2
M5 ON
φ3
M10 , M11
φ4
ON
/ R
Reset Int ( )
VP
Track / Hold
VS
VE /
Track Hold
P-I-N Diodes
Since the photocurrent only responds to photons collected
within the depletion region, to make a sensitive photodetector
we should make W as large as possible. Recall the expression
for W in reverse bias: P I N
s
2e 1 1
W= + (V0 + VR ). Figure 184: Structure of a PIN diode. The large
q NA ND intrinsic region in the middle is fully depleted. In
√ some versions, the P-type surface is exposed to
We see that W ∝ 1/NA + 1/ND . Suppose we let NA or ND go photons. If the P region is kept very thin, the photons
to zero. In that extreme limit, W goes to infinity, so the weakly will pass through it to be absorbed in the thick I
region below.
doped side of the junction will be fully depleted throughout
its volume. By placing an undoped region of intrinsic semicon-
ductor in between the P and N regions of a normal diode, a P N
very wide depletion zone is created, resulting in a highly sensi-
tive photodetector. The wide depletion width also reduces the
junction capacitance and gives it a high-bandwidth in reverse I
Avalanche Photodiodes
Tunnel Diodes EV
VD = VP
When a junction has very high doping, so that the depletion P ←− E N
region becomes very thin, there can be quantum tunneling
−−−− EC
during forward bias operation as well as in reverse breakdown.
The concept is illustrated by the energy band diagrams in Fig- ++++
EV
ure 186. The forward tunnel current comes to a maximum at
the peak voltage VP , where the conduction and valence bands are VD > VP
aligned. This is because the allowed energy states are mostly
squeezed in very close to EC and EV . Most of the mobile elec- Figure 186: Energy bands in a weakly forward
biased tunnel diode. When VD is slightly above zero,
trons are sitting at EC , and they need valid energy states near electrons can tunnel across the thin depletion region.
EV in order to tunnel. The tunnel current is maximum at the peak voltage
VP , where the valence band on one side aligns with
When the diode’s foward voltage is increased higher than
the conduction band on the other side, since the
VP , the bands are moved out of alignment so that the tunneling available energy states are concentrated near the
current actually decreases with increased forward bias. This band edges. When VD > VP , the available states are
poorly aligned so the tunnel current decreases.
causes a behavior called negative differential resistance – the
slope di/dv swings negative. This means that in the small-signal
10−1
equivalent circuit, the tunnel diode behaves like a negative resis-
tor. Negative resistance has a variety of exotic applications, but
the most important use is in high-frequency RF and microwave 10−2
Nobel prizes for the tunnel diode’s inventor, Leo Esaki (tunnel 10−3
diodes are often called “Esaki diodes” in his honor). “peak”
The tunnel diode’s current consists of a primary tunneling
current, superimposed with secondary or “excess” tunneling 10−4
VD − VV
VD
Iex = exp ,
RV Vex
+
to minimize error around the negative differential resistance − R2 C L
(NDR) portion of the transfer characteristic.
A typical tunnel diode application is in RF oscillators like the
one shown in Figure 188. A voltage source and resistor pair are
used to bias the tunnel diode in the middle of the NDR region. Figure 188: Tunnel diode oscillator circuit. Node
numbers corresponding to the SPICE netlist are
In a classical RLC circuit, the resistor dissipates energy with indicated in blue.
each period of the oscillation. If the resistor is negative, however,
then the dissipation is also negative. This means that energy
is added with each period of the oscillation. The tunnel diode
can therefore be used to deliver energy into a resonant LC tank
oscillator, so long as the signal oscillations occur within its NDR
region of operation. The resulting frequency is approximately 0.1
that of the LC tank,
1
ω= √ .
v(2) [V]
0
LC
In the example SPICE netlist shown below, the parameters −0.1
are L = 100 µH and C = 100 pH, yielding a resonant frequency
of about 1.5 MHz. The diode is biased so that its differential
−0.2
resistance is about rd = −1 kΩ. Simulation results are shown
in Figure 189. The waveform is distorted due to the diode’s
0 0.2 0.4 0.6 0.8 1
nonlinear transfer characteristic, but the waveform can be made
Time [s]
·10−5
more sinusoidal by adding resonant filters to the circuit.
Figure 189: Transient SPICE simulation results for a
tunnel diode oscillator.
V1 1 0 DC 1.5V
X1 2 3 MBD1057
R1 1 2 200.0
R2 2 0 20
C1 3 0 100p
L1 3 0 100u
.nodeset v(3)=0.1
.control
tran 1n 10us uic
plot v(3)
.endc
.end
210 ece3410 lecture notes
n+ source n+ drain
L
The Field Effect Transistor (FET) was first described by Julius p− substrate
Edgar Lilienfeld in a 1926 patent. The basic concept is the same
as modern FET devices, but the devices could not be manufac- Figure 191: Cross section of a standard bulk MOSFET.
tured in Lilienfeld’s day since they require highly perfect crystal Devices are most often fabricated on a P-type
substrate with weak doping, indicated as p− . The
structure, precise doping control, and microscopic dimensions. Source and Drain regions are implanted with strong
The first modern Metal Oxide Semiconductor (MOS) FET device N-type doping, indicated as n+ . A layer of insulation,
was demonstrated at Bell Labs in 1960 by Atalla and Kahng. usually comprised of silicon dioxide (SiO2 , often
called the “glass” layer), is grown or deposited on
One of the key innovations in fabricating successful MOSFETs the surface. The region between the Source and
was to grow a layer of insulating Silicon Dioxide on the device’s Drain terminals is called the channel with length L.
Suspended above the channel is the Gate terminal,
surface, thereby controlling surface imperfections that had es- which is a deposited layer of poly-crystaline Si
sentially created a short-circuit between the switching terminals (“poly” for short).
in older devices. The basic structure of the 1960 MOSFET is
still used up to the present day. Variations on this structure
Gate
have been at the cutting edge until very recently. Today we see
competition from alternative nano-wire and “fin” based FET
structures which will be discussed later; first we’ll study the
classic bulk MOSFET device that is used in nearly all integrated
Source Drain
circuits produced since the 1980s.
W
Physical Structure
eox
Cox0 = [F/µm2 ],
tox
214 ece3410 lecture notes
and the total capacitance depends on the area under the gate,
which can be adjusted by an integrated circuit designer. One
of the main benefits of a capacitive gate is that it prevents any
current from flowing into the device’s gate terminal. This makes
the MOSFET nearly ideal for use as a voltage amplifier, and is
a crucial property for energy efficient large-scale logic circuits.
For some very small devices, with channel lengths below 90 nm,
tox is small enough to permit quantum tunneling through the
oxide. This phenomenon, often called “gate leakage,” is one of
many problems that has led the industry to modify or abandon
the classic bulk MOSFET structure. More complex structures,
like the FinFET, are increasingly favored in new technologies.
We will examine the FinFET structure in a later chapter.
On either side of the gate, strong N-type doping (indicated
as N + ) is implanted into the wafer to form the source and drain
terminals of the device. Terminal contacts are created by etching
small holes in the oxide where metal can be deposited. By
alternately depositing oxide, etching oxide, depositing metal,
and so forth, complex circuits can be constructed across the
wafer’s surface. A device with a P-type substrate and N-type
source/drain terminals is referred to as an N-type MOSFET.
We will discuss the complementary P-type MOSFET in a later
section.
Manufacturers and process designers determine many of
the device’s properties, like the doping concentrations, oxide
thickness below the gate, minimum dimensions and spacings
of objects, depth of the source and drain implants, and so on.
But chip-level designers have the freedom to manipulate two
important dimensions: the device’s length L, defined as the
spacing between the source and drain terminals, and the width
W. These dimensions are determined by the device’s physical
layout. An example layout is shown in Figure 192. A layout
designer has considerable freedom to manipulate the device
geometry so long as the basic process design rules are obeyed.
field effect transistors 215
When VG is slightly greater than zero, but less than the device’s
threshold voltage, majority carriers (holes) are repelled from
gate
the region beneath the gate. Minority carriers (electrons) are + + +
attracted from the bulk. The minority carriers aggregate at the
Si surface, just beneath the oxide layer, resulting in an inversion SiO2
n+ source − − − n+ drain
of mobile charge polarity in the MOSFET channel between the
source and drain. In this mode, the potential along the channel
is nearly constant, and the voltage drop between VD and VS − − −
occurs close to their junctions. Hence there is no net E -field in
p− substrate
the channel, and no drift current flows. The excess minority
charges tend to diffuse into the Source and Drain junctions,
Figure 194: MOSFET cross section with weakly
in a process similar to a weakly forward biased diode. Since inverted channel.
the weak inversion current flows primarily by diffusion across
junctions, the device has an exponential current similar to a
diode. The Source and Drain diffusion currents are
φ − VS
Isource = I0 exp
UT gate
φ − VD
Idrain = I0 exp
UT Cox
n+ source − − − n+ drain
where I0 is a scale current and φ is the channel potential.
As Figure 195 shows, a capacitive divider forms between the Cdep
κVG − VS κVG − VD
I = I0 −
UT UT
κVG − VS −VDS
= I0 exp 1 − exp
UT UT
For a bulk MOSFET device, the scale current I0 can be shown to
be !
0 U2
−κVTh
2µn Cox T W
I0 = exp .
κ L UT
When VDS is greater than about 100 mV, the current becomes
insensitive to VDS and the device is said to be in substhreshold
saturation. In this mode we can use the approximate expression
κVG − VS
ID ≈ I0 exp .
UT
field effect transistors 217
A MOSFET device has parameters κ = 0.7, µn = 700 cm2 /Vs, Cox 0 = 5 fF/µm2 , V
Th = 0.5 V, W =
0.4 µm and L = 0.06 µm. What is the device’s current when VS = 0 V, VD = 1 V and VG = 0.2 V?
Solution: First, the device’s scale current is found to be I0 = 228 nA (note: remember to reconcile the
units cm2 in µn with µm2 in Cox 0 !). Since V
DS 100 mV, we can use the approximate expression for
subthreshold saturation. In that case, the current is ID ≈ 40.7 µA.
One of the most important consequences of subthreshold operation is that it sets up a DC leakage
current in digital circuits, which creates static power dissipation in CMOS integrated circuits and
systems. Suppose a MOSFET has scale current I0 = 100 nA, and is sitting idle in its OFF state with
VG = VS = 0 and VD > 100 mV. Will any current flow in the device? From the weak inversion model,
what should we expect the current to be? Will the current have any dependence on VD ? If the system
contains 106 such devices and operates with a power supply of 2 V, what will be the total static power
consumption?
Solution: Since VD > 100 mV, the device is in the subthreshold saturation mode (this is typical of
switched-OFF MOSFETs in digital circuits). In this mode, the current is given by
κVG − VS
Ileak = I0 exp ,
UT
but since VG and VS are both zero, the leakage current is just I0 or 100 nA. As long as VD stays
above 100 mV or so, the leakage current will stay roughly the same. If VD is reduced below this
level, the current will exponentially decrease. If VD → 0 V then the leakage current can be fully
6
suppressed (according to this model). If there are 10 leaking devices then the total static power is
2 V × Ileak × 106 = 200 mW. This may sound small, but it would be enough to drain a high-capacity
cell phone battery, holding a charge of about 2900 mA h, in just over one day, assuming the device
is sitting idle (and in airplane mode) the entire time. Static power is really just the standby power
consumed while the system is idle; active usage will consume dynamic power, which is usually much
greater. Hence it is desirable to save static power so that the energy can be used for actual tasks. As
we saw here, one way of saving static power is to dynamically reduce the supply voltage while the
system is idle, so that the VD for every device is made small. This practice is called dynamic voltage
scaling.
218 ece3410 lecture notes
resulting in a drift current between the drain and source. This SiO2
transition occurs when VG ≈ VTh , and the channel is said to be n+ source − − − − −− n+ drain
strongly inverted. In this mode, the potential along the channel
changes gradually from VS to VD . Let φs ( x ) be the surface − − −
potential at some point x along the channel, so that φ(0) = VS
and φ( L) = VD . Now consider a small “slice” of width dx, the p− substrate
total capacitance of the oxide under the gate is Figure 196: MOSFET cross section with strongly
inverted channel.
Cslice = Cox0 Wdx.
and then we can simultaneously integrate the left side over the
length of the channel from x = 0 up to x = L, and on the right
side over the potential from φs (0) = VS up to φs ( L) = VD :
Z L Z VD
ID dx = µn Cox0 W (VG − φs ( x ) − VTh ) dφs
0 VS
VD
1
⇒ ID L = µn Cox0 W (VG − φs − VTh )2
2 VS
1 W h i
⇒ ID = µn Cox0 (VGS − VTh )2 − (VGD − VTh )2
2 L
field effect transistors 219
1 Wh 2 i
ID = µn Cox0 2
VGS − VGD − 2VGS VTh + VTh2 + 2VGD VTh − VTh2
2 L
1 0 W
h i
= µn Cox 2VG VDS − VD2 − VS2 + 2VTh (VGD − VGS )
2 L
1 W h i
= µn Cox0 2VG VDS −2VS VDS − 2VTh VDS − VD2 − VS2 +2VS VDS
2 L
1 W h i
= µn Cox0 2 (VGS − VTh ) VDS − VD2 − 2VD VS + VS2
2 L
0 W 1 2
= µn Cox (VGS − VTh ) VDS − VDS
L 2
220 ece3410 lecture notes
d 1 W
−1
ROUT = 0
µn Cox (VGS − VTh )2
dvDS 2 L
1 d 1
= µn Cox0 W (VGS − VTh )2
2 dvDS L
1 1 dL
= µn Cox0 W (VGS − VTh )2 − 2
2 L dvDS
1 W 1 dL
= µn Cox0 (VGS − VTh )2 −
2 L L dvDS
? 1 dL
= ID −
L dvDS
E [eV]
When metal, oxide and semiconductor are brought into contact,
the energy bands are altered around the material interfaces.
−qVFB
This so-called “band bending” effect contributes to forming
Φs
the device’s channel. To understand band-bending, we need
to introduce the concept of a material’s work function, which
Ec
is the amount of energy needed to eject an electron from the
Φm
material. This concept was first postulated by Einstein as part of
EF
his work on the photoelectric effect, for which he was awarded
Ev
the Nobel Prize. For metals, the work function energy can be
in the range of visible light, so if a visible photon is absorbed
in metal, it can excite an electron to fly out of the material.
Semiconductors tend to have a much larger work function.
In the MOS structure, shown for the equilibrium case in Fig-
M O S
ure 199, all materials share the same Fermi level, so the vacuum d [nm]
energy level must bend to accommodate differences in the metal Figure 199: MOS band-bending at zero gate bias. The
and semiconductor work functions Φm and Φs , respectively. In bands are bent due to differing metal and semicon-
the bulk (on the far right side of the diagram), the Fermi level ductor work functions, Φm and Φs , respectively.
VTh [V]
⇒ φb = UT ln .
ni
To truly invert the channel, the bands at the surface must be 0.3
bent so that the Fermi level is φb above the intrinsic level instead
0.2
of below it. Since the amount of band bending is equal to
ψs − VFB , the threshold voltage is where
0.1
ψs − VFB = 2φb 1014 1015 1016 1017 1018
NA [dopants/cm3 ]
⇒ ψs = VFB + 2φb
Now to complete the solution for VTh , we need to know how the
surface potential ψs is related to vG . We mentioned before that
there is a capacitive divider formed between the gate oxide and
the depletion layer beneath the channel. The analysis is rather
complex, but it basically reduces to the depletion width under
the gate and has a form similar to our previous solutions for
depletion width. The answer can be shown to be
p
4es qNA φb
vG = ψs + ,
Cox0
so the threshold voltage at the gate is equal to
p
4es qNA φb
VTh = VFB + 2φb + .
Cox0
From this result we can see that the threshold voltage is heav-
ily dependent on doping concentration and oxide capacitance,
and secondarily dependent on temperature. Since most MOS-
FETs are made using poly-Si gates over Si substrates, we can
focus our attention on those materials. For a typical poly-Si
gate, the work function is Φm ≈ 4.05 V. Since the work func-
tion is defined as E0 − EF , and in the Si substrate EF depends
on the doping concentration, the Si work function is doping
dependent:
Φs = E0 − Ei + ( Ei − EF )
= E0 − Ei + φb = 4.15 V + φb
⇒ VFB = 4.05 V − 4.15 V − φb
p
4es qNA φb
⇒ VTh = −0.1 V + φb + .
Cox0
This expression is evaluated for different doping concentrations
and temperatures in ??.
224 ece3410 lecture notes