Analytical Modeling of Output Conductanc PDF
Analytical Modeling of Output Conductanc PDF
Analytical Modeling of Output Conductanc PDF
Fig. 2. Ids versus Vds curves at Vgs = 0.6 V for the halo-doped device and
the uniformly doped device. The Vdsat for the halo-doped device is less than
the uniformly doped device. The parameters used for the device simulation are: Fig. 3. Channel potential along the surface of a uniformly doped device
W = 1 m, L = 1 m, and tox = 1.2 nm. and nonuniformly doped halo device at Vgs = 0.6 V, Vbs = 0.0 V, Vds =
1.2 V, 0.5 V, 0.2 V.
doping where the center region is inverted before the source and
drain ends are inverted. It is also distinctly clear that the output observe that once the transistor reaches saturation, there is no
resistance is much lower for device B. Another interesting point further drop in the significant portion of the channel, and as ex-
to observe is that although both devices were designed to have pected, the additional drain bias is dropped in the small pinch-
the same threshold voltage, the device B pinches off sooner than off region close to the drain. Thus, in saturation, the current
the uniformly doped device A. Both devices are long enough does not increase for the uniformly doped device. For the non-
that the drain saturation occurs because of channel pinch off. uniformly doped device B, however, it is interesting to see that
Under pinch-off conditions, the drain saturation voltage Vdsat additional potential drop continues to occur along the entire
is approximately given by the expression [16] channel length even in saturation. Thus, transistor B behaves
almost as if it were still a resistor in the saturation region.
Vgs Vt This anomalous behavior can be readily explained using a
Vdsat = (1)
two-transistor model. As previously observed [10], [14], only
dQb the halo-doped region on the drain side influences the Rout
=1 + . (2) degradation. Hence, the channel can be split into two transistors
ds
in series, namely the lower doped transistor of length LLhalo
In the above equations, Vt stands for the threshold voltage of on the source side and halo-doped transistor of length Lhalo
the device, alpha represents the depletion charge sensitivity to on the drain side as illustrated in Fig. 4(a). As the drain bias
surface potential, Qb is the depletion charge normalized to the is increased, the transistor operates in the linear region until
gate capacitance, and s is the surface potential. the heavily doped halo transistor saturates. However, the lower
It should be noted that the Vt used in (1) is the threshold doped transistor continues to operate in the linear region. In
voltage corresponding to the doping at the drain end of the this region of operation, two mechanisms contribute to the
channel. Therefore, even though the overall threshold voltage anomalous behavior: 1) DIBL, which is modeled as long-
of the devices A and B are the same, the threshold voltage channel DIBL and 2) movement of the pinch-off point toward
corresponding to the halo region of device B is much higher the source occurring in the halo transistor. Thus, when the
than the uniformly doped regions of device A. Thus, device B halo-doped transistor is operating in saturation, the following
pinches off sooner than device A because the threshold voltage drops in potential occur in the following different regions of the
corresponding to the heavily doped halo region on the drain side transistor.
is higher than that of the threshold voltage corresponding to the
uniformly doping in device A. In other words, as the drain bias 1) Vds Vdsathalo is dropped in the pinch-off region of the
is increased, the drain end of the channel for device B pinches halo transistor [Region III in Fig. 4(b)].
off earlier because it requires a higher vertical field to support 2) Vdsathalo Vd is dropped in the linear region of the halo
the inversion charges. Furthermore, for the heavily doped transistor [Region II in Fig. 4(b)].
halo region is greater than for the uniformly doped region. 3) Vd is dropped in low-doped transistor [Region I in
Thus, the impact of is to further reduce the Vdsat for device B. Fig. 4(b)].
The potential along the channel is shown for both devices at With increasing Vds , the pinch-off point moves closer to the
multiple drain biases in Fig. 3. At the lowest drain bias (0.2 V), lower doped region. This causes a reduction in the length
the transistor is still in the linear region, while at the other two of Region II and thus reducing the effective resistance of
drain biases (0.5 and 1.2 V), the transistor is in saturation. Com- Region II. Since Vdsathalo is dropped across both Regions I and
paring the different curves of the uniformly doped device A, we II, a larger value of Vd must now be dropped across Region I
MUDANAI et al.: ANALYTICAL MODELING OF OUTPUT CONDUCTANCE IN MOSFETs 2093
Fig. 4. (a) Illustration of the two-transistor model. The nonuniform transistor can be modeled as two transistors in series. The source-side transistor is uniformly
doped with Cbulk , while the drain-side transistor is uniformly doped with Chalo . (b) Low-doped Region I operates in the linear region; the halo-doped region
behaves like a transistor in saturation. Region III is the pinch-off region. (c) Equivalent-circuit model to approximate the equivalent Rout for the two-transistor
system.
From [17]
Lhalo = Lhalo l log y + 1 + y 2 (10)
Vdd Vdsat
halo
y= . (11)
l Esat
In (10), l is a characteristic length, which is a function of
the oxide thickness, junction depth, and dielectric constants,
and Esat is the lateral field needed for velocity saturation.The
saturation voltage for the entire device can be written as
halo
Vdsat = Vd s + Vdsat . (12)
halo
Using (9), (10), and (12), the Rout can be derived as
2 2
halo L Vdsat Vds Vdsat
Rout = halo + . (13)
Id L l
Bulk
C. Rout
Fig. 5. (a) Id versus Vds fits for a 2-m drawn gate length and 10-m-wide
As discussed in Section II, the bulk-doped transistor on the NMOSFET at Vbs = 0.0 V and Vgs = 0.5, 0.6, 0.7, 1.1, 1.2 V. (b) Gds versus
source side operates in the linear region. Hence, the channel Vds fits with and without the Rout model for a 2-m drawn gate length and
10-m-wide NMOSFET at Vbs = 0.0 V and Vgs = 0.5, 0.6, 0.7, 1.1, 1.2 V.
current can be expressed as The plot without the Rout model includes long-channel DIBL model.
Bulk
Fig. 7. Comparison of Gds versus length between the experimental data and
the compact model with and without the analytic Rout model at Vds = 0.8 V
and Vgs = 0.6 V.
Fig. 6. (a) Id versus Vds fits for 2.0-m drawn gate length and 10.0-m-
wide PMOSFET at Vbs = 0.0 V and Vgs = 0.5, 0.6, 0.7, 1.1, 1.2 V.
(b) Gds versus Vds fits with and without the Rout model for 2-m
drawn gate length and 10-m-wide PMOSFET at Vbs = 0.0 V and Vgs =
0.5, 0.6, 0.7, 1.1, 1.2 V. The plot without the Rout model includes
long-channel DIBL model.
Id = l ln + 1 + 2
Lhalo 1 + 2Id LBulk
Lhalo (Vt )
(22)
where
Vds Vdsat L
= . (23)
Vdsat l
The analytic Rout model was fit to the NMOS and PMOS IV
data from Intels 90-nm technology [18]. The Ids and Gds as a Fig. 8. (a) Third-derivative fits for a 2-m drawn gate length and 10-m-wide
NMOSFET at Vbs = 0.0 V and Vgs = 0.5, 0.7, 1.2 V. (b) Third-derivative fits
function of the drain bias for the NMOSFET at L = 2.0 m for 2-m drawn gate length and 10-m-wide PMOSFET at Vbs = 0.0 V and
are compared well with the experimental data, as shown in Vgs = 0.5, 0.7, 1.2 V.
Fig. 5(a) and (b). A similar fit is shown for the PMOSFET
IV data at L = 2.0 m at Vbs = 0.0 in Fig. 6(a) and (b). It model with and without the analytic Rout model. The impact of
should be noted that in Figs. 5 and 6, the curves without the the Rout model is lucid in this comparison. As the gate length is
Rout model do include the long-channel DIBL model that was reduced, the halo implants from the source and the drain sides
developed in house. The long-channel model was calibrated begin to merge, as a result, the Rout degradation effect slowly
based on the observed DIBL in the subthreshold region. begins to reduce. This region can be seen at lengths somewhat
Thus, the comparison in Figs. 5(b) and 6(b) shows the effect greater than 0.1 m as a mild hump. Another important result
of CLM only in the halo transistor and not the long-channel to be noted is that at lengths below 0.1 m, where the two halo
DIBL. In Fig. 7, a comparison of Gds versus L is shown for the regions are already merged, the Rout degradation model has no
NMOSFET device between the experimental and the compact effect on the characteristics of the very short devices. This is
2096 IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. 53, NO. 9, SEPTEMBER 2006
achieved through the Vt term in (22), which becomes zero for [14] M. Miyamoto et al., Asymmetrically-doped buried layer (ADB) structure
very short devices. Finally, for use in low-distortion analog and CMOS for low-voltage mixed analog-digital applications, in VLSI Symp.
Tech. Dig., 1996, p. 102.
RF applications, the compact model must be able to predict the [15] P. R. Gray and R. G. Meyer, Analysis and Design of Analog Integrated
harmonic distortion. It is well known [19] that when a sinusoid Circuits. Hoboken, NJ: Wiley, 1993.
is applied to the input terminals of a MOSFET, the output signal [16] Y. Tsividis, Operation and Modeling of the MOS Transistor. London,
U.K.: Oxford Univ. Press, 1999.
contains higher order harmonics in addition to the fundamental [17] N. Arora et al., PCIM: A physically based continuous short-channel
frequency, and that harmonic distortion is related to the higher IGFET model for circuit simulations, IEEE Trans. Electron Devices,
derivatives of the draincurrent. Thus, the third derivative of vol. 41, no. 6, pp. 988997, Jun. 1994.
[18] K. Kuhn et al., A comparison of State-of-the-Art NMOS and SiGe HBT
the measured draincurrent with respect to the drain bias is devices for Analog/Mixed-signal/RF circuit applications, in VLSI Symp.
compared against the model to ensure that the model meets the Tech. Dig., 2004, p. 224.
requirements of distortion simulations in Fig. 8(a) and (b). The [19] R. Van Langevelde et al., Accurate drain conductance modeling
for distortion analysis in MOSFETs, in IEDM Tech. Dig., 1997,
harmonic distortion of the draincurrent in the saturation region pp. 313316.
originates from the channel-length-modulation model used to
describe the halo transistor, and it can be seen that the compact
model captures the third derivative in this region accurately. Sivakumar Mudanai (M04) received the B.Tech.
degree from the Indian Institute of Technology,
Chennai, Madras, India, and the M.S. and Ph.D. de-
grees from the University of Texas at Austin, Austin.
V. CONCLUSION He is currently with the Compact Device-
Modeling Group, Intel Corporation, Hillsboro, OR.
The degradation of Rout for long-channel devices has been His research interests include MOSFET modeling
explained using the 2-D numerical device simulation. A phys- for digital, analog, and RF applications.
ical argument was presented to help understand how the pres-
ence of a potential barrier on the drain side impacts the channel
resistance using a simple two-transistor model. An analytic
approach to model the Rout degradation was developed and
implemented in conjunction with our previous three-transistor Wei-Kai Shih received the B.S. degree in physics
from National Taiwan University, Taipei, Taiwan,
threshold-voltage model. The model was validated on several in 1989, the M.A. degree in physics and the Ph.D.
Intels technologies, but the model results are shown against degree in electrical engineering from the Univer-
Intels 90-nm technology data. The model was also shown to sity of Texas at Austin, Austin, in 1994 and 1997,
respectively.
capture the hump in Gds observed as a function of the length in In 1997, he joined TCAD, Intel Corporation,
technologies that use a strong halo. Hillsboro, OR, to develop electrothermal modeling
and simulation tools for full-chip reliability verifi-
cation. He later moved into the compact transistor
modeling area, focusing on gate leakage modeling
R EFERENCES and analog-RF modeling. His current research interests include modeling front-
[1] D. Buss, Device issues in the integration of analog/RF functions in deep end, back-end, and process variations for aggressively scaled technology, as
submicron digital CMOS, in IEDM Tech. Dig., 1999, pp. 423425. well as improving computer-aided design methodology to facilitate design and
[2] T. Hori, A 0.1 m CMOS technology with Tilt-Implanted Punch through process convergence in all phases of a product development cycle.
Stopper(TIPS), in IEDM Tech. Dig., 1994, pp. 7578.
[3] Y. Okumura et al., A novel source-to-drain non-uniformly doped chan-
nel (NUDC) MOSFET for high current drivability and threshold voltage
controllability, in IEDM Tech. Dig., 1990, pp. 391394. Rafael Rios (S84M85) studied engineering sci-
[4] Y. Taur et al., CMOS devices below 0.1 m: How high will performance ences in Montevideo, Uruguay. He received the M.S.
go?, in IEDM Tech Dig., 1997, pp. 215218. and Ph.D. degrees in electrical engineering from
[5] S. Saha, Design considerations for 25 nm MOSFET devices, Solid State Drexel University, Philadelphia, PA.
Electron., vol. 45, no. 10, pp. 18511857, Oct. 2001. He is currently with Intel Corporation, Hillsboro,
[6] Y. Cheng et al., Modeling reverse short channel and narrow width effects OR, where he works on the development and valida-
in small size MOSFETs for circuit simulation, in Proc. SISPAD, 1997, tion of models and tools for circuit design.
pp. 249252.
[7] H. Ueno et al., Impurity-profile-based threshold-voltage model of
pocket-implanted MOSFETs for circuit simulation, IEEE Trans.
Electron Devices, vol. 49, no. 10, pp. 17831789, Oct. 2002.
[8] Y.-S. Pang et al., Analytical subthreshold surface potential model for
pocket n-MOSFETs, IEEE Trans. Electron Devices, vol. 49, no. 12,
pp. 22092216, Dec. 2002. Xuemei (Jane) Xi (M02) received the B.S., M.S.,
[9] R. F. M. Roes et al., Implications of pocket optimisation on analog and Ph.D. degrees in microelectronics from Peking
performance in deep sub-micron CMOS, in Proc. ESSDERC, 1999, University, Beijing, China.
pp. 176179. In 1995, she joined the Institute of Microelec-
[10] K. M. Cao et al., Modeling of pocket implanted MOSFETs for anom- tronics, Peking University, where she became an
alous analog behavior, in IEDM Tech. Dig., 1999, pp. 171174. Associate Professor in 1997. From 1999 to 2000, she
[11] A. Chatterjee et al., Transistor design issues in integrating analog func- was a Senior Engineer with Motorola China. From
tions with high performance digital MOS, in VLSI Symp. Tech. Dig., 2000 to 2005, she was a Research Staff with the
1999, p. 147. Department of Electrical Engineering and Computer
[12] R. Rios et al., A three-transistor threshold voltage model for halo Science, University of California, Berkeley, where
processes, in IEDM Tech. Dig., 2002, pp. 113116. she worked on the Berkeley Short-Channel IGFET
[13] J. B. Jacobs et al., Channel profile engineering for MOSFETs with Model 3 (BSIM3)/BSIM4/BSIMSOI model development and offered technical
100 nm channel lengths, IEEE Trans. Electron Devices, vol. 42, no. 5, support to BSIM users from both the industry and the academics. In 2006, she
pp. 870875, May 1995. joined Intel Corporation, Hillsboro, OR.
MUDANAI et al.: ANALYTICAL MODELING OF OUTPUT CONDUCTANCE IN MOSFETs 2097
Jung-Hoon Rhew received the B.S. degree from Po- Paul Packan received the B.S. degree in electri-
hang University of Science and Technology, Pohang, cal engineering from the University of Washington,
Korea, in 1994, the M.S. degree from the University Seattle, in 1984, and the M.S. and Ph.D. degrees
of Washington, Seattle, in 1999, and the Ph.D. de- in electrical engineering from Stanford University,
gree from Purdue University, West Lafayette, IN, in Stanford, CA, in 1985 and 1991, respectively.
2003, all in electrical engineering. His study and re- After graduation, he joined Siemens AG, Munich,
search focused on the physics and modeling of semi- Germany, as a Guest Scientist, working on high-
conductor devices. speed bipolar transistor design. In 1992, he joined
In 2003, he joined the Compact Device-Modeling Intel Corporation, Hillsboro, OR, working on MOS
Group, Intel Corporation, Hillsboro, OR, where transistor technology development. After working on
he works on device modeling and associated four generations of Intel technologies and managing
automation. the Process and Device-Modeling Group, he moved into circuit-level device
modeling, where he managed the Compact Device-Modeling Group. He is
currently managing the Device Development Group for the 32-nm technology
node. He is the author and coauthor of more than 30 papers and is also the
Kelin Kuhn received the B.S. degree in electri- holder of multiple patents in the area of transistor design and architecture.
cal engineering from the University of Washington,
Seattle, in 1980 and the M.S. and Ph.D. degrees
in electrical engineering from Stanford University,
Stanford, CA, in 1985.
Prior to joining Intel, she was a Tenured Fac-
ulty Member with the Department of Electrical and
Computer Engineering, University of Washington. In
1997, she joined Intel Corporation, Hillsboro, OR,
where she has held a variety of technical positions
on the 0.35-$\mu\hbox{m}$, 130-nm, and 90-nm
nodes and, most recently, was the Integration Manager for the RF/analog deriv-
ative technology on 90 nm. She is currently the Device Manager for the Intel
45-nm CMOS technology, which is being developed at the Portland Technology
Development Center, Intel Corporation. She is the author of numerous technical
publications on optical and electronic properties of semiconductor devices.