BSIM4 Gate Leakage Model Including Source-Drain Partition
BSIM4 Gate Leakage Model Including Source-Drain Partition
BSIM4 Gate Leakage Model Including Source-Drain Partition
ECB,SiO2
= 3.1 eV
J
HVB
J
ECB
J
EVB
ECB,SiO2
= 3.1 eV
J
HVB
Figure 1 Three mechanisms of gate dielectric direct
tunneling leakage.
Table 1 Major mechanisms (as shown in Fig.1)
contributing to current components illustrated in Fig. 2.
The mechanism governing I
gso
and I
gdo
under all bias
conditions is ECB.
Current
Component
Igc
Igb
Region of
Operation
inversion
V
g
>0
V
g
<0
PMOS
HVB
ECB
EVB
NMOS
ECB
ECB
EVB
0-7803-6441-4/00/$10.00 (C) 2000 IEEE
mechanism, the leakage current density can be modeled
by:
( )( )
ox ox ox
t V V B
ox
aux g
ntox
ox
oxref
g
e
t
V V
t
T
A J
+
,
_
1
2
where
b
h q A 8
2
, h qm B
b ox
3 2 8
2 / 3
, m
ox
is the
effective carrier mass in oxide,
b
is the tunneling
barrier height (shown in Fig. 1),
ox
t is the oxide
thickness, T
oxref
is the reference oxide thickness at
which all the parameters are extracted, ntox is a fitting
parameter that defaults to 1, and
aux
V is an auxiliary
function which approximates the density of tunneling
carriers as well as available states and tabulated in
Table 2.
MOSFET Gate Leakage Model
As illustrated in Fig. 2, the gate tunneling current is
composed of several components. I
gb
is the gate-to-
substrate leakage current; I
gso
and I
gdo
are parasitic
leakage currents through gate-to-S/D extension overlap
region; and I
gc
is the gate-to-inverted channel tunneling
current. Part of I
gc
is collected by the source (I
gcs
) while
the rest goes to the drain (I
gcd
). I
gb
, I
gso
, I
gdo
, and I
gc
can
be easily determined from the MOS capacitor model
introduced in the previous section. The remaining
challenge is modeling I
gc
for nonzero V
ds
and
partitioning I
gc
into I
gcs
and I
gcd
.
To derive I
gcs
and I
gcd
, current continuity equation is
solved for the voltage along the channel as a function
of x, which is 0 at the source and L at the drain.
From drift equation, the channel current density in
source-drain direction is
( ) dx dV V V V C
th gs ox
.
Current continuity requires:
( ) ( )
( ) 0 +
ox g
th gs
C x J
dx
dx dV V V V d
where V(x) is the channel potential relative to the
source. The gate tunneling current density at x can
be approximated with
( ) V B
g
V V BT
oxs
E B
ox g
e J e AE e AE J
oxs ox ox
*
0
2 2
where J
go
is the gate tunneling current density with
V
ds
=0 (can be modelled by Eq. 1),
2
gcd
*
oxs ox i
V BT P B and
gs oxs
V V is the voltage
across the gate oxide at x=0. P
igcd
is a fitting
parameter added for flexibility with a default value
of 1.
Assuming the gate leakage current is small
compared to the drain current, the perturbation of
J
g
on V(x) is small,
( ) ( ) ( ) x V x V x V
1 0
+
where V
0
(x) is the solution of Eq. (2) with J
g
(x)=0
and V
0
(x)>>V
1
(x). V
0
(x) can be solved and
approximated to be:
( ) ( ) ( )
x K
L x V V V V V V V V x V
ds ds th gs th gs th gs
2 2
2
where
( )L V V V V V V K
th gs ds ds th gs
) 2 (
. This approximation
is not very accurate when x and V
ds
are both large
compared to L and V
gs
-V
th
respectively. However,
in that case J
g
is large only at small x, where the
approximation is accurate. Therefore, the
simplification is reasonable.
(1)
(2)
(3)
(4)
(5)
Table 2 The V
aux
function, which models the density of
tunneling carriers and available states.
Region of
Operation
Current V
aux
acc. I
gb
,
_
,
_
+
tm
fb gb
tm
V nigbacc
V V
V nigbacc exp 1 log
I
gb
,
_
,
_
+
tm
oxdepinv
tm
V nigbinv
eigbinv V
V nigbinv exp 1 log
inv.+dep.
I
gc
,
_
,
_
+
tm
th eff gs
tm
V nigc
V V
V nigc
0 _
exp 1 log
both I
gso
gs
V
both I
gdo
gd
V
Igso Igdo
Igb
Igcs Igcd
Igc
Figure 2 Components of tunneling current. I
gso
and I
gdo
are parasitic tunneling currents through gate-
source/drain overlap regions. I
gb
flows between gate and
substrate. I
gc
is the gate-to-channel tunneling current. It
is partitioned into I
gcs
, which flows to the source, and
I
gcd
, which flows to the drain.
0-7803-6441-4/00/$10.00 (C) 2000 IEEE
Eq. (2) can be further simplified considering
V
0
(x)>>V
1
(x).
( )
0 ) ( ) 0 (
, 0 2
1 1
0 1 0
2
0
2
1
2
1
2
0
0
*
+
L V V
C e J dx dV dx dV dx V d V dx V d V V
ox
V B
g gst
V
1
can then be solved as
( ) ( )
( ) Kx V K B C
e
L
x
e J
V
gst ox
KL B Kx B
g
'
+
2
2
*
0
1
* *
1 1
and I
gcs
and I
gcd
can be derived as:
( )
( )
( )
2 2
2
*
*
0 1
gcd
2 2
2
*
*
0
0
1
1
1
* *
*
L K B
e KLe B WL J
dx
dV
W KL V C I
L K B
e KL B WL J
dx
dV
W V C I
KL B KL B
g
L x
gst ox
KL B
g
x
gst ox gcs
+
+
Using the same derivation, better accuracy but more
complex equation can be obtained by including the
second order term in the exponent of Eq. (3).
From Eq. (8), the I
gc
-V
ds
function can then be easily
calculated as:
( )
KL B
e WL J
I
KL B
g
gc
*
0
*
1
Parameter extraction can follow two steps. The first
step is to extract the parameters in Eq. 1 from the
tunneling currents I
gc
(from shorted source and drain)
and I
gb
. The second step is to extract P
igcd
from the
gate-to-channel tunneling current as a function of drain
and gate bias using Eq. (9). Afterwards, the gate-to-
channel current partition can be calculated using Eq.
(8) without any parameter. The proposed partition
model is based on physical derivation and requires no
empirical parameters. This is especially important
because the partition ratio is extremely difficult or
impossible to measure.
Verification of the Model
We found that the proposed tunneling model (Eq. 1) for
MOS capacitors works well for all major mechanisms
(Fig. 3) and for a wide range of oxide thickness (Fig.
8). Our simple I
gc
-V
ds
model (Eq. 9) also agrees with
simulation results and experimental data very well as
shown in Fig. 4,5, and 11. The proposed partition
model (Eq. 8) is verified by simulation data as a
function of V
gs
(Fig. 6) and oxide thickness (Fig. 7). Fig.
9 shows the model agrees well with I
gc
-V
ds
data for
different V
gs
.
Summary
We have developed a simple gate dielectric leakage
current model suitable for compact MOSFET
models. This model considers all three major
leakage components and is accurate in direct-
tunneling regime with excellent T
ox
scaling
capability. Physical gate-to-channel current
partition models are developed for the first time.
With its great accuracy and satisfying test results
from the member companies of the Compact
Modeling Council (CMC), the proposed models
have been implemented into the BSIM4 and can
already help to address the impacts of gate leakage
current on circuit performance.
Acknowledgement: The work is supported by SRC
under Contract 98-SJ-417 and TI, Mentor
Graphics, Xilinx under MICRO program.
References
[1] 1999 International Technology Roadmap for Semiconductors.
[2] H. S. Momose et al., 1.5nm direct-tunneling gate oxide Si
MOSFETs, IEEE Trans. Elec. Dev., vol.43(8), pp. 1233, Aug. 1996.
[3] P. J. Wright and K. C. Saraswat, Thickness limitations of SiO2
gate dielectrics for MOS ULSI, IEEE Trans. Elec. Dev., vol.37(8),
pp. 1884, Aug. 1990.
[4]K. F. Schuegraf and C. Hu, Hole injection SiO2 breakdown model
for very low voltage lifetime extrapolation, IEEE Trans. Elec. Dev.,
vol. 41(5), pp.761, May 1994.
[5] W.-C. Lee and C. Hu, Modeling gate and substrate currents due
to conduction- and valence-band electron and hole tunneling, 2000
Symp. on VLSI Tech., pp.198, 2000.
[6] C.-H. Choi et al., Direct tunneling current model for circuit
simulation, IEDM Technical Digest pp.735, 1999.
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
10
-8
10
-7
10
-6
10
-5
10
-4
10
-3
10
-2
10
-1
10
0
10
1
T
ox
=24.8 A
W/L=10/10
NMOS I
gc
(ECB)
NMOS I
gb
(EVB)
PMOS I
gc
(HVB)
Model
J
g
(
A
/
c
m
2
)
|V
g
| (V)
Figure 3 The proposed MOS model agrees with
measured data for tunneling mechanisms ECB, EVB,
and HVB.
(6)
(7)
(8)
(9)
0-7803-6441-4/00/$10.00 (C) 2000 IEEE
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4
0.0
0.2
0.4
0.6
0.8
1.0
V
gs
=1.5V
V
gs
=2.0V
Model
I
g
c
/
I
g
c
(
V
d
s
=
0
)
V
ds
(V)
Figure 4 The model and normalized I
gc
-V
ds
2D
simulation data at different V
gs
.
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4
0.0
0.2
0.4
0.6
0.8
1.0
T
ox
=20A
T
ox
=15A
Model
I
g
c
/
I
g
c
(
V
d
s
=
0
)
V
ds
(V)
Figure 5 The model and normalized I
gc
-V
ds
2D
Simulation data at different oxide thickness.
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4
0.0
0.1
0.2
0.3
0.4
0.5
V
gs
=1.5V
V
gs
=2.0V
Model
I
g
c
d
/
I
g
c
V
ds
(V)
Figure 6 The current partition model agrees with
simulation data at different V
gs
bias.
0.0 0.2 0.4 0.6 0.8 1.0 1.2
1.0x10
-9
1.5x10
-9
2.0x10
-9
2.5x10
-9
3.0x10
-9
V
gs
=0.6V x10
V
gs
=0.9V x3
V
gs
=1.2V x1
Experimental Data
Model
T
ox
=20A W/L=10m/2.7m
I
g
(
A
)
V
ds
(V)
Figure 9 Agreement of the model with experimental I
g
-
V
ds
data at different V
gs
bias.
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4
0.0
0.1
0.2
0.3
0.4
0.5
T
ox
=20A
T
ox
=15A
Model
I
g
c
d
/
I
g
c
V
ds
(V)
Figure 7 The current partition model agrees with
simulation data at different oxide thickness.
0.0 0.2 0.4 0.6 0.8 1.0 1.2
10
-4
10
-3
10
-2
10
-1
10
0
10
1
10
2
10
3
10
4
20A
18A
16A
12A
Data
Model
I
g
c
(
J
/
c
m
2
)
V
gs
(V)
Figure 8 The proposed MOS tunneling model is
scalable with oxide thickness (NMOS).
0-7803-6441-4/00/$10.00 (C) 2000 IEEE