-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathtest014_nmosSubthresholdLeakageTemp.gnuplot
152 lines (125 loc) · 4.44 KB
/
test014_nmosSubthresholdLeakageTemp.gnuplot
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
# TEST 014 - gnuplot
#
# TSMC 180nm Subthreshold Leakage Model - Temperature dependency
#
# This test allows to simulate NMOS transistor subthreshold leakage
#
# Author: Jan Belohoubek, 01/2021
#
# https://users.fit.cvut.cz/~belohja4/
#
#
#
# Parameters were extracted from BSIM3v3.3 and from TSMC model file
#
# For details about model implementation, please refer:
# WESTE, Neil HE; HARRIS, David. CMOS VLSI design: a circuits and systems perspective. Pearson Education India, 2015.
# BSIM3v3.3 MOSFET Model: http://ngspice.sourceforge.net/external-documents/models/bsim330_manual.pdf
#
# --- Model parameters ---
# Particular model parameters were extracted from simulation "test014", where:
# NMOS channel length is 0.2u
# NMOS channel width is 2.0u
#
W=2
L=0.2
# The subthreshold slope, S reaches the common value 100mV/dec; the value below is extracted from simmulation
S = 0.0945
# Off Leakage; the value below is extracted from simmulation
Ids0 = 90e-9 # ~ 90n
# Threshold Voltage
VT = 0.354505
# DIBL coefficient in subthreshold region;
ETA = 0.0231564
# body effect coefficient
K1 = 0.5733393
# Thermal voltage in at room temperature
Tv = 0.026
# Boltzman's constant
k = 1.38e-23
# electron charge
q = 1.602e-19
# subthreshold swing -- an estimated value
n = 1.3
# Offset voltage in the subthresh-old region at large W and L
Voff = -0.0652968
#Mobility at Temp = Tnom
U0 = 128.7704538
#Supply Voltage
Vsupp = 1.8
# Oxide thickness
Tox = 4.1e-9
#
beta = U0 * Tox * W/L
# Ids0 = beta* (((k * TEMP) / q)**2) * exp(1.8)
# Body-source voltage
Vsb0 = -0.05
Vsb1 = -0.04
Vsb2 = -0.03
Vsb3 = -0.02
Vsb4 = -0.01
Vsb5 = -0.0
# --- Graph settings ---
# print to PDF -- coment to enable live 3D view
set terminal pdf size 24cm, 12cm font "Sans,20"
set output "output/test014_Vgs_Vsb_Temp.pdf"
#set title "Subthreshold Leakage in TSMC 180nm NMOS"
set ticslevel 0.01
set tics
set xrange [0:0.1]
set yrange [250:400]
set xlabel "V_{gs} [V]"
set ylabel "T [K]"
set zlabel "I_{ds} [nA/um^2]" rotate by 90
# set log z
#set pm3d
set palette
#unset surface
set samples 70,45
set isosamples 70,45
set hidden3d
set angle degree
# --- Data Sets ---
IDS0(x,y) = ((beta* (((k * y) / q)**2) * exp(1.8)) * (exp( ( x - VT + (ETA * Vsupp) - K1 * Vsb0 ) / ( n * ((k * y) / q)) * (1 - exp(-1*Vsupp/(k*y/q))) )) ) * 2.5 * 10**9 # nA/um2
set table $Data00
splot IDS0(x,y)
unset table
IDS1(x,y) = ((beta* (((k * y) / q)**2) * exp(1.8)) * (exp( ( x - VT + (ETA * Vsupp) - K1 * Vsb1 ) / ( n * ((k * y) / q)) * (1 - exp(-1*Vsupp/(k*y/q))) )) ) * 2.5 * 10**9 # nA/um2
set table $Data01
splot IDS1(x,y)
unset table
IDS2(x,y) = ((beta* (((k * y) / q)**2) * exp(1.8)) * (exp( ( x - VT + (ETA * Vsupp) - K1 * Vsb2 ) / ( n * ((k * y) / q)) * (1 - exp(-1*Vsupp/(k*y/q))) )) ) * 2.5 * 10**9 # nA/um2
set table $Data02
splot IDS2(x,y)
unset table
IDS3(x,y) = ((beta* (((k * y) / q)**2) * exp(1.8)) * (exp( ( x - VT + (ETA * Vsupp) - K1 * Vsb3 ) / ( n * ((k * y) / q)) * (1 - exp(-1*Vsupp/(k*y/q))) )) ) * 2.5 * 10**9 # nA/um2
set table $Data03
splot IDS3(x,y)
unset table
IDS4(x,y) = ((beta* (((k * y) / q)**2) * exp(1.8)) * (exp( ( x - VT + (ETA * Vsupp) - K1 * Vsb4 ) / ( n * ((k * y) / q)) * (1 - exp(-1*Vsupp/(k*y/q))) )) ) * 2.5 * 10**9 # nA/um2
set table $Data04
splot IDS4(x,y)
unset table
IDS5(x,y) = ((beta* (((k * x) / q)**2) * exp(1.8)) * (exp( ( x - VT + (ETA * Vsupp) - K1 * Vsb5 ) / ( n * ((k * y) / q)) * (1 - exp(-1*Vsupp/(k*y/q))) )) ) * 2.5 * 10**9 # nA/um2
set table $Data05
splot IDS5(x,y)
unset table
# --- Graph ---
# Color palete
Zmin = 0
Zmax = 0.35
Frac(z) = (z-Zmin)/(Zmax-Zmin)
# MyPalette
Red(z) = 65536 * ( Frac(z) > 0.75 ? 255 : int(255*abs(2*Frac(z)-0.5)))
Green(z) = 256 * (Frac(z) > 0.50 ? 128 : int(255*sin(180*Frac(z))))
Blue(z) = (Frac(z) > 0.25 ? 64 : int(255*cos(150*Frac(z))))
MyPalette(z) = Red(z) + Green(z) + Blue(z)
set view 44,316
splot $Data00 u 1:2:3:(MyPalette($3)) w l lc rgb var title sprintf("V_{sb} = %0.2f", Vsb0), \
$Data01 u 1:2:3:(MyPalette($3)) w l lc rgb var title sprintf("V_{sb} = %0.2f", Vsb1), \
$Data02 u 1:2:3:(MyPalette($3)) w l lc rgb var title sprintf("V_{sb} = %0.2f", Vsb2), \
$Data03 u 1:2:3:(MyPalette($3)) w l lc rgb var title sprintf("V_{sb} = %0.2f", Vsb3), \
$Data04 u 1:2:3:(MyPalette($3)) w l lc rgb var title sprintf("V_{sb} = %0.2f", Vsb4), \
$Data05 u 1:2:3:(MyPalette($3)) w l lc rgb var title sprintf("V_{sb} = %0.2f", Vsb5)
pause -1 "Hit any key to continue"