Introduction to the
PSP
MOSFET model
- R. van Langevelde, G.D.J. Smit, A.J. Scholten and D.B.M. Klaassen
- G. Gildenblat, X. Li, H. Wang and W. Wu
MOS-AK, Grenoble, September 16th, 2005
PSP MOSFET model R. van Langevelde, G.D.J. Smit , A.J. Scholten and - - PowerPoint PPT Presentation
Introduction to the PSP MOSFET model R. van Langevelde, G.D.J. Smit , A.J. Scholten and D.B.M. Klaassen G. Gildenblat, X. Li, H. Wang and W. Wu MOS-AK, Grenoble, September 16 th , 2005 outline introduction & history introduction
MOS-AK, Grenoble, September 16th, 2005
PSP, MOS-AK, Grenoble, September 16th, 2005 2/50
PSP, MOS-AK, Grenoble, September 16th, 2005 3/50
introduction: merger SP and MM11 MOS Model 11
(Philips)
SP
(Penn State)
PSP, MOS-AK, Grenoble, September 16th, 2005 4/50
introduction: motivation
industrial standard model for compact CMOS modeling
Research) were both worthy candidates
competition
Why PSP?
PSP, MOS-AK, Grenoble, September 16th, 2005 5/50
introduction: similarities MM11 and SP
global (maxiset) model parameters
MOS Model 11 (MM11) and SP are both compact MOSFET models that:
PSP, MOS-AK, Grenoble, September 16th, 2005 6/50
introduction: CMC requirements
– distortion – symmetry
– including accumulation region for varactors
requirements for new standard model:
PSP, MOS-AK, Grenoble, September 16th, 2005 7/50
introduction: CMC selection procedure
07/2004 candidate models:
selection procedure for next generation model:
UC Berkeley Qinv iterative
EPFL Qinv iterative
Hiroshima Univ. ψs iterative
PSU/Philips ψs explicit
model type solution from
10/2005 final selection
PSP, MOS-AK, Grenoble, September 16th, 2005 8/50
PSP, MOS-AK, Grenoble, September 16th, 2005 9/50
general: structure of PSP temperature scaling model equations
T
DS GS SB
V V V local model
currents
charges noise local parameter set stress scaling W, L geometry scaling SA, SB global par set
PSP, MOS-AK, Grenoble, September 16th, 2005 10/50
general: physical effects included
PSP, MOS-AK, Grenoble, September 16th, 2005 11/50
general: physical effects included (ii)
PSP, MOS-AK, Grenoble, September 16th, 2005 12/50
general: PSP’s unique features
– reproduces details in long channel gds
– most complete and robust – correctly including velocity-saturation
– verified against channel segmentation
unavailable outside PSP context; verified against device simulations and several modern CMOS processes
PSP, MOS-AK, Grenoble, September 16th, 2005 13/50
PSP, MOS-AK, Grenoble, September 16th, 2005 14/50
PSP, MOS-AK, Grenoble, September 16th, 2005 15/50
physics: DC: surface-potential-based
Idrift = f(VGB ,ψs0 ,ψsL) Idiff = g(VGB ,ψs0 ,ψsL) IDS = Idrift + Idiff ψs-based model:
single equation for whole operation range:
1 2 10-11 10-10 10-9 10-8 10-7 10-6 10-5 10-4
Idiff I drift
IDS= Idrift+ Idiff
VSB = 0 V VDS = 1 V VGS (V) IDS (A)
PSP is surface-potential (ψs) based
PSP, MOS-AK, Grenoble, September 16th, 2005 16/50
physics: DC: ψs-calculation
− ⋅ + − ⋅ ⋅ + = − −
− − −
1 1
T s T s T B
T T s 2 s FB GB ϕ ψ ϕ ψ ϕ ϕ
ϕ ϕ ψ γ ψ e e e V V
V
0.4 0.8 1.2
1 2
V GB - V FB (V) ψs (V)
V = 0 V
[R. Rios et al., IEDM2004]: approximation of ψs becomes inaccurate for high positive back bias
ψs obtained from explicit
analytical approximation, accurate within 1nV under all relevant conditions
PSP, MOS-AK, Grenoble, September 16th, 2005 17/50
physics: DC: ψs–based model ψs-calculation: description of ideal long-channel MOSFET
in ‘real’ device:
PSP, MOS-AK, Grenoble, September 16th, 2005 18/50
physics: DC: mobility reduction mobility reduction due to phonon, surface- roughness, and Coulomb scattering
E eff µ
b b inv C
Q Q Q + ∝ µ
3 1 eff ph −
∝ E µ
2 eff sr −
∝ E µ
Si inv b eff
ε η Q Q E ⋅ + − =
PSP:
( )
THEMU eff
MUE E ⋅
BETN, MUE, THEMU, CS, XCOR
PSP, MOS-AK, Grenoble, September 16th, 2005 19/50
0.0 0.4 0.8 1.2 50 100 V
GS (V)
g
m (
µ A/V)
physics: DC: mobility reduction (ii) measurement
gm versus VGS (100µm/10µm n-MOS, VDS=30mV)
VGS (V) gm (µA/V)
VSB=0V VSB=1.2V
model
CS switched
0.0 0.4 0.8 1.2 50 100 V
GS (V)
g
m (
µ A/V)
VGS (V)
CS switched
gm (µA/V)
PSP, MOS-AK, Grenoble, September 16th, 2005 20/50
physics: DC: velocity saturation
in PSP:
along channel
⋅ ⋅ ⋅ − = x v Q L W I d
inv D
E x v
sat
v v =
|| sat ||
1 E v E v ⋅ + ⋅ = µ µ
2 || sat ||
1 ⋅ + ⋅ = E v E v µ µ
v = µ · E||
THESAT, (THESATG, THESATB)
PSP, MOS-AK, Grenoble, September 16th, 2005 21/50
0.0 0.4 0.8 1.2 2 4 V
GS (V)
I
D (mA)
0.0 0.4 0.8 1.2 2 4 V
DS (V)
I
D (mA)
physics: DC: velocity saturation (ii) PSP VDS (V) VGS (V) ID (mA) ID (mA) measurement ID vs VGS ID vs VDS
10/0.12 n-MOS VGS=1.2V VGS=0.8V VGS=0.4V VDS=1.2V 0.6V 50mV
PSP, MOS-AK, Grenoble, September 16th, 2005 22/50
0.0 0.4 0.8 1.2 10
10
10
10
10
V
GS (V)
g
m
i (A/V i )
0.0 0.3 0.6 0.9 10
10
10
V
DS (V)
g
DSi (A/V i )
gmi (= ∂iID/ ∂VGS
i) vs VGS
VGS=1.2V
physics: DC: velocity saturation (iii) VDS (V) VGS (V) gDSi (A/Vi) gmi (A/Vi) measurement PSP
10/0.12 NMOS gDSi (= ∂iID/ ∂VDS
i) vs VDS
VDS=1.2V
gm1 gm2 gm3 gDS1 gDS2 gDS3
PSP, MOS-AK, Grenoble, September 16th, 2005 23/50
physics: DC: short-channel effects
Si 2 2 2 2
Poisson equation:
Si A 2 2 2 2
(in subthreshold)
2 2 Si A 2 2
Si eff 2 2
non-ideal subthreshold slope and DIBL (due to 2D-effects): lateral gradient factor f:
SB DS GS A eff
n+ p n+
+++
n+
x - -
n+ p n+
+++
n+
x
PSP-parameters: F0, AF, BF, CF, (CFB)
PSP, MOS-AK, Grenoble, September 16th, 2005 24/50
physics: DC: channel length modulation
VG VD Nsub uniformly doped device VD VG Npck Npck Nsub pocket implanted device VDS electron potential
in PSP:
conductance included
threshold
PSP, MOS-AK, Grenoble, September 16th, 2005 25/50
physics: DC: CLM (ii)
0.0 0.4 0.8 1.2 0.00 0.05 0.10 0.15 0.20 V
DS (V)
I
D (mA)
ID-VDS and gDS-VDS for VSB=0V and T=25°C 10µm/10µm NMOS (90nm process) ID (mA)
VDS (V)
0.0 0.4 0.8 1.2 10
10
10
10
V
DS (V)
g
DS (A/V)
gDS (A/V)
VDS (V)
PSP, MOS-AK, Grenoble, September 16th, 2005 26/50
physics: DC: CLM (iii)
ID-VDS and gDS-VDS for VSB=0V and T=25°C 10µm/0.1µm NMOS (90nm process)
0.0 0.4 0.8 1.2 2 4 6 V
DS (V)
I
D (mA)
ID (mA)
VDS (V)
0.0 0.4 0.8 1.2 10
10
10
10
V
DS (V)
g
DS (A/V)
gDS (A/V)
VDS (V)
PSP, MOS-AK, Grenoble, September 16th, 2005 27/50
physics: DC: gate leakage current
Gate EV Oxide Oxide EC Ei EF Substrate
V q ⋅
q χ ⋅
Gate EV Oxide Oxide EC Ei EF Substrate
V q ⋅
q χ ⋅ gate current density: tunneling probability
( )
G
V P F J ⋅ ∝
S
Esaki-Tsu supply function
PSP-parameters: IGINV, IGOV, GC0, GC2, GC3
PSP, MOS-AK, Grenoble, September 16th, 2005 28/50
1.E-14 1.E-13 1.E-12 1.E-11 1.E-10 1.E-09
0.5 1 1.5
V GS (V) I G (A)
Measurements MM11
I GB I GS+ I GD I GOV
physics: DC: gate leakage current (ii)
gate leakage model includes 3 components:
S S D D I IG
G
IGD IGS IGB IGOV
model measurements
PSP, MOS-AK, Grenoble, September 16th, 2005 29/50
physics: DC: gate leakage current (iii)
IG-VGS curves for VSB = 0 V and T = 25 °C
0.36µm/0.1µm NMOS (CMOS-90)
IG (A) VGS (V)
VDS=25mV VDS=1.0V
model measurements
PSP, MOS-AK, Grenoble, September 16th, 2005 30/50
PSP, MOS-AK, Grenoble, September 16th, 2005 31/50
physics: AC: intrinsic charges
n+ p n+
+ + +
n+
⋅ ⋅ =
L
x Q W Q
gate G
d
∫
⋅ ⋅ ⋅ =
L
x Q L x W Q
inv D
d
∫
⋅ ⋅ − ⋅ =
L
x Q L x W Q
inv S
d 1
G S D B
Q Q Q Q − − − =
intrinsic capacitances: :
≠ ∂ ∂ − = ∂ ∂ = j i V Q j i V Q C for for
j i j i ij
where i, j =G, S, D or B
PSP, MOS-AK, Grenoble, September 16th, 2005 32/50
physics: AC: intrinsic charges (ii)
2
L
s s m
ψ ψ ψ + =
linearization at mid-point bias: which allows for simple charge expressions
Vds = 2V Vbs= 0 V Vfb=-1V
1 2 3 4 0.0 0.3 0.6 0.9 Linearized CSM Original CSM Cdg Cbg Csg Cgg
Normalized Transcapacitances Vgs (V)
PSP, MOS-AK, Grenoble, September 16th, 2005 33/50
physics: AC: intrinsic charges (iii)
0.4 0.8 1.2
1 2 3
V GS (V) C GG (nF)
Measurements MOS Model 11
tox=3.6nm
0.4 0.8 1.2
1 2 3
V GS (V) C GG (nF)
Measurements MOS Model 11
tox=3.6nm
effects
0.4 0.8 1.2
1 2 3
V GS (V) C GG (nF)
Measurements MOS Model 11
tox=3.2nm
charge model includes:
PMOS, VDS=0 V, W/L=80*612/2.5µm
model
PSP, MOS-AK, Grenoble, September 16th, 2005 34/50
physics: noise
noise model in PSP includes: + correlation HF noise model in PSP has the same accuracy as the segmentation approach
PSP, MOS-AK, Grenoble, September 16th, 2005 35/50
results for 90nm CMOS gate-source bias [V] Sid , Sig [A2/Hz]
NMOS VDS=1.0V f=10GHz L=0.1µm
physics: noise: thermal/induced gate noise
measurement PSP model
PSP, MOS-AK, Grenoble, September 16th, 2005 36/50
physics: NQS model
continuity equation
terminal currents, and all operation regions
modulation, poly-depletion, and QM effects
NQS model in PSP
PSP, MOS-AK, Grenoble, September 16th, 2005 37/50
physics: NQS: subcircuit implemention
⋅ ⋅ − = ⋅ = x V q x x I W t q d d d d d d 1 d d
i D i
µ
+ spline approximation system of (coupled) differential equations
( )
N 1 k k
, , u u f dt du K − =
implemented as sub-circuits, solved by circuit simulator:
( ) ( )
2 1 2 2 2 1 1 1
, , , V V f f V V f f = =
C R C R Cf1 Cf2 V1=u1 V2=u2
R is large carries negligible current
Example: (N = 2) current continuity equation
PSP, MOS-AK, Grenoble, September 16th, 2005 38/50
physics: NQS : results
1E9 1E10 1E11 0.6 0.8 1.0 1.2 1.4 1E7 1E8 1E9 1E10 1E11
QS NQS with N=1 NQS with N=2 NQS with N=3 NQS with N=5 NQS with N=9
L=1µm Tox=4nm Nsub=2E17cm
|ym|/gm
Frequency (Hz)
L=1µm Tox=4nm Nsub=2E17cm
QS NQS with N=1 NQS with N=2 NQS with N=3 NQS with N=5 NQS with N=9
Phase (degree) Frequency (Hz)
normalized transconductance Ym versus f
magnitude of Ym phase of Ym
PSP, MOS-AK, Grenoble, September 16th, 2005 39/50
physics: NQS : results (ii)
VDS=1.5 V VGS=0.5 V VGS=1.0 V VGS=1.5 V PSP, SWNQS=5 MM11, 5 segments
Re(Y11) versus f
VDS=1.5 V VGS=0.5 V VGS=1.0 V VGS=1.5 V measurement QS model NQS model
PSP, MOS-AK, Grenoble, September 16th, 2005 40/50
physics: NQS : resistor subcircuit
D RG G S B Rjun,S Rjun,D Rbulk Cjun,S Rwell Cjun,D
PSP NQS
JUNCAP2 JUNCAP2
PSP, MOS-AK, Grenoble, September 16th, 2005 41/50
– depletion capacitance
– ideal current – Shockley-Read-Hall – trap-assisted tunneling – band-to-band tunneling – avalanche & breakdown – shot noise
physics: diode: JUNCAP2 model JUNCAP2: A new junction diode model
forward & reverse
PSP, MOS-AK, Grenoble, September 16th, 2005 42/50
bottom STI-edge gate-edge
physics: diode: I-V curves
T=-400C T=2000C
0.12µm CMOS, n+/p-well
ideal BBT SRH+TAT
PSP, MOS-AK, Grenoble, September 16th, 2005 43/50
PSP, MOS-AK, Grenoble, September 16th, 2005 44/50
parameters: structure of PSP temperature scaling model equations
T
DS GS SB
V V V local model
currents
charges noise local parameter set stress scaling W, L geometry scaling SA, SB global par set
PSP, MOS-AK, Grenoble, September 16th, 2005 45/50
parameters: parameter extraction
local parameters for each dut geometry scaling global parameter set
temperature scaling measurements
0.1 0.15 0.2 0.25 5 10 15
local global
1/LE (1/µm) P
PSP, MOS-AK, Grenoble, September 16th, 2005 46/50
availability of PSP (in public domain)
–April 2005: first release –August 2005: optimized version
–QS model –extended code incl. NQS and gate/bulk resistors
using ADMS)
–interface to Pstar, Spectre, ADS, and UltraSim –interface to Aplac in progress
PSP, MOS-AK, Grenoble, September 16th, 2005 47/50
availability of PSP (ii) circuit simulator PSP available
Spectre (Cadence) yes ADS (Agilent) yes Hspice (Synopsis) 09/2005 (β version: 07/2005) HSIM (Synopsis) next release NanoSiM (Synopsis) 03/2006 Eldo (Mentor Graphics) 10/2005 (β version: 08/2005) Spice (Berkeley) 09/2005 (from PSP-homepage)
PSP, MOS-AK, Grenoble, September 16th, 2005 48/50
acknowledgements
from Freescale
we would like to thank: compact modeling work at Penn State is supported in part by SRC, Freescale and IBM
PSP, MOS-AK, Grenoble, September 16th, 2005 49/50
selected references
– A.J. Scholten et al., IEDM 2005
– H. Wang et al., TED 2003 – H. Wang et al., CICC 2005
– R. van Langevelde et al., IEDM 2003 – J. Paasschens, TED 2005 (accepted).
PSP, MOS-AK, Grenoble, September 16th, 2005 50/50
MM11 and SP
characteristics over complete bias, temperature and geometry range
characteristics including noise and NQS effects
conclusions
http://pspmodel.ee.psu.edu