LEVEL 3 IDS: Empirical Model
This sections provides the LEVEL 3 IDS: Empirical model parameters and equations.
LEVEL 3 Model Parameters
The LEVEL 3 model parameters follow.
Basic Model Parameters
Name (Alias)
|
Units
|
Default
|
Description
|
LEVEL
|
|
1.0
|
DC model selector. LEVEL=3 is an empirical model.
|
COX
|
F/m
2
|
3.453e-4
|
Oxide capacitance per unit gate area. If this parameter is not specified, it is calculated from TOX.
|
DERIV
|
|
1
|
Derivative method selector
DERIV=0: analytic
DERIV=1: finite difference
|
KAPPA
|
V
-1
|
0.2
|
Saturation field factor. This parameter is used in the channel length modulation equation.
|
KP (BET, BETA)
|
A/V
2
|
2.0e-5
|
Intrinsic transconductance parameter. If this parameter is not specified and UO and TOX are entered, KP is calculated from
KP = UO
·
COX
|
TOX
|
m
|
1e-7
|
Gate oxide thickness
|
VMAX (VMX)
|
m/s
|
0.0
|
Maximum drift velocity of carriers. Use zero to indicate an infinite value.
|
Effective Width and Length Parameters
Name (Alias)
|
Units
|
Default
|
Description
|
DEL
|
m
|
0.0
|
Channel length reduction on each side
DELscaled = DEL
·
SCALM
|
LD (DLAT, (LATD)
|
m
|
|
Lateral diffusion into channel from source and drain diffusion,
If LD and XJ are unspecified, LD Default= 0.0.
If LD is unspecified but XJ is specified, LD is calculated from XJ as LD = 0.75
·
XJ.
|
LDAC
|
m
|
|
This parameter is the same as LD, but if LDAC is included in the .MODEL statement, it replaces LD in the Leff calculation for AC gate capacitance.
|
LREF
|
m
|
0.0
|
Channel length reference
LREFscaled = LREF
·
SCALM
|
LMLT
|
|
1.0
|
Length shrink factor
|
WD
|
m
|
0.0
|
Lateral diffusion into channel width from bulk
WDscaled = WD
·
SCALM
|
WDAC
|
m
|
|
This parameter is the same as WD, but if WDAC is included in the .MODEL statement, it replaces WD in the Weff calculation for AC gate capacitance.
|
WMLT
|
|
1.0
|
Diffusion layer and width shrink factor
|
WREF
|
m
|
0.0
|
Channel width reference
WREFscaled = WREF
·
SCALM
|
XJ
|
m
|
0.0
|
Metallurgical junction depth
XJscaled = XJ
·
SCALM
|
XL (DL, LDEL)
|
m
|
0.0
|
Length bias accounts for masking and etching effects
XLscaled = XL
·
SCALM
|
XW (DW, WDEL)
|
m
|
0.0
|
Width bias accounts for masking and etching effects
XWscaled = XW
·
SCALM
|
Threshold Voltage Parameters
Name (Alias)
|
Units
|
Default
|
Description
|
DELTA
|
|
0.0
|
Narrow width factor for adjusting threshold
|
ETA
|
|
0.0
|
Static feedback factor for adjusting threshold
|
GAMMA
|
V
1/2
|
0.5276
|
Body effect factor. This parameter is calculated from NSUB if not specified (See Common Threshold Voltage Parameters).
|
LND
|
µm/V
|
0.0
|
ND length sensitivity
|
LN0
|
µm
|
0.0
|
N0 length sensitivity
|
ND
|
V
-1
|
0.0
|
Drain subthreshold factor
|
N0
|
|
0.0
|
Gate subthreshold factor (typical value=1)
|
NFS (DFS,NF, DNF)
|
cm
-2
·
V
-1
|
0.0
|
Fast surface state density
|
NSUB (DNB, NB)
|
cm
-3
|
1e15
|
Bulk surface doping. This parameter is calculated from GAMMA if not specified.
|
PHI
|
V
|
0.576
|
Surface inversion potential. This parameter is calculated from NSUB if not specified (see Common Threshold Voltage Parameters).
|
VTO (VT)
|
V
|
|
Zero-bias threshold voltage. This parameter is calculated if not specified (see Common Threshold Voltage Parameters).
|
WIC
|
|
0.0
|
Sub-threshold model selector
|
WND
|
µm/V
|
0.0
|
ND width sensitivity
|
WN0
|
µm
|
0.0
|
N0 width sensitivity
|
Mobility Parameters
Name (Alias)
|
Units
|
Default
|
Description
|
THETA
|
V
-1
|
0.0
|
Mobility degradation factor
|
UO (UB,UBO)
|
cm
2
/
(V
·
s)
|
600(N) 250(P)
|
Low field bulk mobility. This parameter is calculated from KP if KP is specified.
|
LEVEL 3 Model Equations
The LEVEL 3 model equations follow.
IDS Equations
The following describes the way the LEVEL 3 MOSFET model calculates the drain current, I
ds
.
Cutoff Region,
(See subthreshold current)
On Region, vgs > vth
where:
and:
NOTE: In the above equation the factor 4 should be 2, but since SPICE uses a factor of 4, Star-Hspice uses factor of 4 as well.
The narrow width effect is included through the f
n
parameter:
The term f
s
expresses the effect of the short channel and is determined as:
Effective Channel Length and Width
The model determines effective channel length and width in the LEVEL 3 model as follows:
Threshold Voltage,
vth
The effective threshold voltage, including the device size and terminal voltage effects, is calculated by:
where:
or
The VTO is the extrapolated zero-bias threshold voltage of a large device. If VTO, GAMMA, and PHI are not specified, Star-Hspice computes them (see Common Threshold Voltage Parameters).
Saturation Voltage,
vdsat
For the LEVEL 3 model, Star-Hspice determines saturation voltage due to channel pinch-off at the drain side. The model uses the parameter VMAX to include the reduction of the saturation voltage due to carrier velocity saturation effect.
where:
The surface mobility parameter "us" is defined in the next section. If the model parameter VMAX is not specified, then:
Effective Mobility,
ueff
The model defines the carrier mobility reduction due to the normal field as the effective surface mobility (us).
vgs>vth
The model determines the degradation of mobility due to the lateral field and the carrier velocity saturation if you specify the VMAX model parameter.
VMAX>0
otherwise,
Channel Length Modulation
For vds>v
dsat
, the channel length modulation factor is computed. The model determines the channel length reduction
differently, depending on the VMAX model parameter value.
VMAX = 0
VMAX>0
where E
p
is the lateral electric field at the pinch off point. Its value is approximated by:
The current I
ds
in the saturation region is computed as:
In order to prevent the denominator from going to zero, Star-Hspice limits the
value as follows:
If
then
Subthreshold Current, I
ds
This region of operation is characterized by the model parameter for fast surface state (NFS). The modified threshold voltage (von) is determined as follows:
NFS>0
where:
The current I
ds
is given by:
vgs<von
NOTE: The model does not use the modified threshold voltage in strong inversion.
If WIC=3, the model calculates subthreshold current differently. In this case, the I
ds
current is:
Subthreshold current isub for LEVEL=3 is the same as for LEVEL=13 (see Subthreshold Current ids).
N0eff and NDeff are functions of effective device width and length.
Compatibility Notes
This section describes compatibility issues.
Star-Hspice versus SPICE3
Differences between Star-Hspice and Berkeley SPICE3 can arise in the following situations:
Small XJ
Star-Hspice and SPICE3 differ for small values of XJ, typically less than 0.05 microns. Such small values for XJ are physically unreasonable and should be avoided. XJ is used to calculate the short-channel reduction of the GAMMA effect,
f
s
is normally less than or equal to 1. For very small values of XJ, f
s
can be greater than one. Star-Hspice imposes the limit f
s
<=
1.0, while SPICE3 allows f
s
>1.0.
ETA
Star-Hspice uses 8.14 as the constant in the ETA equation, which provides the variation in threshold with v
ds
. Berkeley SPICE3 uses 8.15.
Solution:
To convert a SPICE3 model to Star-Hspice, multiply ETA by 815/814.
NSUB Missing
When NSUB is missing in SPICE3, the KAPPA equation becomes inactive. In Star-Hspice, a default NSUB is generated from GAMMA, and the KAPPA equation is active.
Solution:
If NSUB is missing in the SPICE3 model, set KAPPA=0 in the Star-Hspice model.
LD Missing
If LD is missing, Star-Hspice uses the default 0.75·XJ. SPICE3 defaults LD to zero.
Solution:
If LD is missing in the SPICE3 model, set LD=0 in the HSPICE model.
Constants
Boltzmann constant
|
k
|
= 1.3806226e-23J
·
K
-1
|
Electron charge
|
e
|
= 1.6021918e-19C
|
Permittivity of silicon dioxide
|
|
= 3.45314379969e-11F/m
|
Permittivity of silicon
|
|
= 1.035943139907e-10F/m
|
Temperature Compensation
The example below verifies temperature dependence for LEVEL 3.
Input file
$ test of temp dependence for LEVEL=3 Tlevc=0 Tlev=1
.option ingold=2 numdgt=6
.temp 25 100
vd d 0 5
vg g 0 2
m1 d g 0 0 nch w=10u L=1u
.op
.print id=lx4(m1) vdsat=lv10(m1)
.model nch nmos LEVEL=3 tlev=1 tlevc=0 acm=3
+ uo=600 tox=172.6572
+ vto=0.8 gamma=0.8 phi=0.64
+ kappa=0 xj=0
+ nsub=1e16 rsh=0
+ tcv=1.5e-3 bex=-1.5
.end
This simple model, with XJ=0 and KAPPA=0, has a saturation current:
Using the model parameters in the input file and the equations from the previous page:
At room temperature:
At T=100:
Star-Hspice results:
T=25, id=6.91200e-04
T=100, id=5.72451e-04
These results are in agreement with the hand calculations.
Star-Hspice Manual - Release 2001.2 - June 2001