HSPICE® Reference Manual Elements and Device Models
HSPICE® Reference Manual Elements and Device Models
SystemC is a trademark of the Open SystemC Initiative and is used under license.
ARM and AMBA are registered trademarks of ARM Limited.
Saber is a registered trademark of SabreMark Limited Partnership and is used under license.
All other product or company names may be trademarks of their respective owners.
1. Overview of Models . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
Using Models to Define Netlist Elements. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Supported Models for Specific Simulators . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
Selecting Models . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
Subcircuits . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
Safe Operating Area Warnings . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Model LEVEL Parameter Must be a Constant Value . . . . . . . . . . . . . . . . . . . . 6
Use of Example Syntax . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
iii
Contents
3. Diodes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Diode Types . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
Using Diode Model Statements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
Setting Control Options . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
Bypassing Latent Devices (HSPICE only) . . . . . . . . . . . . . . . . . . . . 51
Setting Scaling Options . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
Using the Capacitor Equation Selector Option — DCAP . . . . . . . . . 51
Using Control Options for Convergence . . . . . . . . . . . . . . . . . . . . . . 52
Specifying Junction Diode Models . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
Using the Junction Model Statement . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
Using Junction Model Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
Geometric Scaling for Diode Models . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
LEVEL=1 Scaling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
LEVEL=3 Scaling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
Defining Diode Models . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
Diode Current . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
Using Diode Equivalent Circuits . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
Determining Temperature Effects on Junction Diodes . . . . . . . . . . . . . . . 66
Using Junction Diode Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
iv
Contents
v
Contents
Scaling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
JFET Current Conventions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
JFET Equivalent Circuits . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
Transconductance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
Output Conductance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
JFET and MESFET Model Statements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
JFET and MESFET Model Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . 120
ACM (Area Calculation Method) Parameter Equations . . . . . . . . . . 128
JFET and MESFET Capacitances . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
Gate Capacitance CAPOP=0 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
Gate Capacitance CAPOP=1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
Gate Capacitance CAPOP=2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
Capacitance Comparison (CAPOP=1 and CAPOP=2) . . . . . . . . . . . . . . 134
JFET and MESFET DC Equations. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135
DC Model LEVEL=1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135
DC Model LEVEL=2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 136
DC Model LEVEL=3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 136
JFET and MESFET Noise Models . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
Noise Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
JFET and MESFET Temperature Equations . . . . . . . . . . . . . . . . . . . . . . . . . . 140
Temperature Compensation Equations . . . . . . . . . . . . . . . . . . . . . . . . . . 142
Energy Gap Temperature Equations . . . . . . . . . . . . . . . . . . . . . . . . 142
Saturation Current Temperature Equations . . . . . . . . . . . . . . . . . . . 143
Gate Capacitance Temperature Equations . . . . . . . . . . . . . . . . . . . 143
Threshold Voltage Temperature Equation . . . . . . . . . . . . . . . . . . . . 144
Mobility Temperature Equation. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
Parasitic Resistor Temperature Equations . . . . . . . . . . . . . . . . . . . . 145
TriQuint (TOM) Extensions to LEVEL=3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 145
LEVEL=7 TOM3 (TriQuint’s Own Model III) . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
Using the TOM3 Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
Model Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 148
DC Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 148
Capacitance Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 149
LEVEL=8 Materka Model. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 152
Using the Materka Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 153
DC Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 153
Gate Capacitance Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 154
Noise Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 155
References. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 157
vi
Contents
vii
Contents
viii
Contents
ix
Contents
Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 331
x
About This Manual
This manual describes standard models that you can use when simulating your
circuit designs in HSPICE or HSPICE RF:
■
Passive devices
■ Diodes
■
JFET and MESFET devices
■
BJT devices
Chapter Description
Chapter 1, Overview of Describes the elements and models you can use
Models to create a netlist in HSPICE.
Appendix A, Finding Device Lists device libraries you can use in HSPICE.
Libraries
Manual Description
HSPICE User Guide: RF Describes how to use special set of analysis and
Analysis design capabilities added to HSPICE to support RF
and high-speed circuit design.
AvanWaves User Guide Describes the AvanWaves tool, which you can use to
display waveforms generated during HSPICE circuit
design simulation.
Note: To use this feature, the HSPICE documentation files, the Index
directory, and the index.pdx file must reside in the same
directory. (This is the default installation for Synopsys
documentation.) Also, Adobe Acrobat must be invoked as a
standalone application rather than as a plug-in to your web
browser.
Conventions
The following typographical conventions are used in Synopsys HSPICE
documentation.
Convention Description
Bold Indicates user input — text you type verbatim — in syntax and
examples. Bold indicates a GUI element.
Convention Description
Edit > Copy Indicates a path to a menu command, such as opening the
Edit menu and choosing Copy.
Customer Support
Customer support is available through SolvNet online customer support and
through contacting the Synopsys Technical Support Center.
Accessing SolvNet
SolvNet includes an electronic knowledge base of technical articles and
answers to frequently asked questions about Synopsys tools. SolvNet also
gives you access to a wide range of Synopsys online services, which include
downloading software, viewing Documentation on the Web, and entering a call
to the Support Center.
To access SolvNet:
1. Go to the SolvNet Web page at http://solvnet.synopsys.com.
2. If prompted, enter your user name and password. (If you do not have a
Synopsys user name and password, follow the instructions to register with
SolvNet.)
If you need help using SolvNet, click Help on the SolvNet menu bar.
The link to any recorded training is
https://solvnet.synopsys.com/trainingcenter/view.faces
Acknowledgments
Portions Copyright (c) 1985-90 by Kenneth S. Kundert and the University of
California.
Portions Copyright (c) 1988-90 Regents of the University of California.
Describes the elements and models you can use to create a netlist in HSPICE
and HSPICE RF.
• T element
• U element
• W element
■
Metal Oxide Semiconductor Field Effect Transistors (MOSFETs), see the
HSPICE MOSFET Models Manual.
This chapter discusses these topics:
■
Using Models to Define Netlist Elements
■
Supported Models for Specific Simulators
■
Subcircuits
■
Safe Operating Area Warnings
■
Model LEVEL Parameter Must be a Constant Value
■
Use of Example Syntax
Within your netlist, each element that refers to a model is known as an instance
of that model. When your netlist refers to predefined device models, you reduce
both the time required to create and simulate a netlist, and the risk of errors,
compared to fully defining each element within your netlist.
Selecting Models
To specify a device in your netlist, use both an element and a model statement.
The element statement uses the simulation device model name to reference
the model statement. The following example uses the MOD1 name to refer to a
BJT model. The example uses an NPN model type to describe an NPN
transistor.
Q3 3 2 5 MOD1 <parameters>
.MODEL MOD1 NPN <parameters>
You can specify parameters in both element and model statements. If you
specify the same parameter in both an element and a model, then the element
parameter (local to the specific instance of the model) always overrides the
model parameter (global default for all instances of the model, if you do not
define the parameter locally).
The model statement specifies the type of device — for example, a BJT, the
device type might be NPN or PNP.
Subcircuits
X<subcircuit_name> adds an instance of a subcircuit to your netlist. You must
already have defined that subcircuit in your netlist by using a .MACRO
or .SUBCKT command.
Argument Definition
parnam A parameter name set to a value (val) for use only in the
subcircuit. It overrides a parameter value in the subcircuit
definition, but is overridden by a value set in a .PARAM statement.
Example 1
The following example calls a subcircuit model named MULTI. It assigns the
WN = 100 and LN = 5 parameters in the .SUBCKT statement (not shown).
The subcircuit name is X1. All subcircuit names must begin with X.
X1 2 4 17 31 MULTI WN = 100 LN = 5
Example 2
This example defines a subcircuit named YYY. The subcircuit consists of two 1
ohm resistors in series. The .IC statement uses the VCC passed parameter to
initialize the NODEX subcircuit node.
.SUBCKT YYY NODE1 NODE2 VCC = 5V
.IC NODEX = VCC
R1 NODE1 NODEX 1
R2 NODEX NODE2 1
.EOM
XYYY 5 6 YYY VCC = 3V
HSPICE ships hundreds of examples for your use; see Listing of Input
Demonstration Files for paths to demo files.
You can use the set of passive model definitions in conjunction with element
definitions to construct a wide range of board and integrated circuit-level
designs. Passive device models let you include the following in any analysis:
■ Transformers
■
PC board trace interconnects
■
Coaxial cables
■ Transmission lines
The wire element model is specifically designed to model the RC delay and RC
transmission line effects of interconnects, at both the IC level and the PC board
level.
To aid in designing power supplies, a mutual-inductor model includes switching
regulators, and several other magnetic circuits, including a magnetic-core
model and element. To specify precision modeling of passive elements, you
can use the following types of model parameters:
■
Geometric
■ Temperature
■
Parasitic
These topics are covered in the following sections:
■ Resistor Device Model and Equations
■
Resistor Temperature Equations
■
LEVEL2 CMC R2 Resistor Model
■
LEVEL 5 CMC R3 Resistor Model
■
Capacitor Device Model and Equations
■
Inductor Device Model and Equations
Wire RC Model
You can use the .MODEL statement to include a Wire RC model in your
HSPICE netlist and evaluate both thermal noise and flicker noise in HSPICE.
For a general description of the .MODEL statement, see .MODEL in the
HSPICE and RF Command Reference.
Syntax
.MODEL MNAME R KEYWORD=value [CRATIO=val]
The wire element RC model is a CRC (pi) model. Use the CRATIO wire model
parameter to allocate the parasitic capacitance of the wire element (between
the input capacitor and the output capacitor of the model). This allows for
symmetric node impedance for bidirectional circuits, such as buses.
Parameter Description
mname Model name. Elements use this name to reference the model.
Parameter Description
CRATIO Ratio to allocate the total wire element parasitic capacitance. This is the
capacitance between the capacitor connected to the input node, and the
capacitor connected to the output node of the wire element pi model.
You can assign a value between 0 and 1 to CRATIO. Default=0.5
0.5 Assigns half of the parasitic capacitance to the input node, and half
to the output node.
If you set a CRATIO value outside the 0 to 1.0 range, simulation shows a
warning, sets CRATIO to 0.5, and continues the analysis.
in out
name. In the following example, R1 assumes that REXX refers to the model, and
not to the parameter.
Resistor Syntax
Rxxx n1 n2 mname [R=]resistance [TC1=val]
+ [TC2=val] [VC1=val] [VC2=val] [SCALE=val] [M=val] [AC=val]
+ [DTEMP=val] [L=val] [W=val] [C=val] [NOISE=val]
You can then use the standard resistor model call to map the models to an
element declaration:
R1 1 2 REXX L=0.6 W=0.5
■
Leff represents the effective length of the resistor, from physical edge to
physical edge.
■
DWeff is the distance from the drawn edge of the resistor to the physical
edge of the resistor.
The effective width is the same as the width used in the resistor calculation.
Leff = Lscaled – 2 ⇒DLeff
Weff = Wscaled – 2 ⇒DWeff
If you specify the element capacitance, C:
CAPeff = C ⋅ SCALE ( element ) ⋅ M
Otherwise, the equation calculates the capacitance from the Leff, Weff, and COX
values:
CAPeff = M ⋅ SCALE ( element ) ⋅ [ Leff ⋅ Weff ⋅ COX ]
+ 2 ⋅ ( Leff + Weff ) ⋅ CAPSW ]
Computing the bottom-wall capacitance, COX, is based on a hierarchy of
defaults and specified values, involving:
■
dielectric thickness (THICK)
■
relative dielectric constant (DI)
Whether you specify a COX value affects how HSPICE uses the equation:
■
If you specify COX=value, then the equation uses the value.
■
If you do not specify COX, but you do specify a value other than zero for
THICK (the dielectric thickness):
• If you specify a non-zero value for DI=value, then:
DI ⋅ ε o
COX = --------------------
THICK
• If you do not specify a DI value or if the value is zero, then:
ε ox
COX = -------------------
THICK
The following values apply to the preceding equation:
ε o = 8.8542149e-12 F/meter
ε ox = 3.453148e-11 F/meter
■
If you do not specify COX, and THICK = 0, this is an error.
• If you specify only the model capacitance (CAP), then:
CAPeff = CAP ⋅ SCALE ( element ) ⋅ M
• If you specify the capacitance, but you do not specify the bulk node, then
the resistor model does not evaluate the capacitance, and issues a
warning message.
4kT 1 / 2
inr = ⎛⎝ NOISE ⋅ ------------⎞⎠
Rval
In the preceding equation, NOISE is a model parameter (default=1). To
eliminate the contribution of resistor noise, specify the NOISE parameter in a
resistor model.
Parameter Description
RX Transfers the function of thermal noise to the output. This is not noise,
but is a transfer coefficient, which reflects the contribution of thermal
noise to the output. For example:
V(output) = I(local)*rx(from local to output)
Where V(output) is the noise voltage at the output port, I(local) is the
local noise current in the specific noise element.
Parameter Description
Δt t - tnom
AF
KF × ⎛⎝ I ⎞⎠
Noise ( sid ) = -----------------------------------------------------------------------
⎛ Leff Lf⎞ × ⎛ Weff Wf⎞ × ⎛ f Ef⎞
⎝ ⎠ ⎝ ⎠ ⎝ ⎠
Where
I Current A
f Frequency HZ
The following lists the parameters and descriptions for the flicker noise model
KF flicker noise ( 2 – AF ) ( LF + WF ) ( EF – 1 ) 0 no
A × m × Hz
coefficient( ≥ 0
)
AF exponent of none 2 no
current (>0)
LF exponent of none 1 no
effective length
(>0)
WF exponent of none 1 no
effective
width(>0)
EF exponent of none 1 no
frequency(>0)
If these parameters are not specified, HSPICE uses their default values, and
the flicker noise will be 0. Observe that these parameters all have limitations:
KF ≥ 0 , AF, LF, WF and EF >0. If their values exceed the limitations, HSPICE
will issue warning messages.
The parameters l, w, dlr, dw must also be specified in the model card for
HSPICE to evaluate the effective length and width. If these parameters are not
specified, their default value is 0, so the effective length and width are 0, and
there will be a divide-by-zero error in the evaluation of the flicker noise
equation. HSPICE will abort and issue an error message.
Controlling the Evaluation of Noise
HSPICE uses the NOISE parameter in a model as a switch to control
evaluation of thermal and flicker noise. Its value is 0 or 1. For example:
.model Res1 R noise=1 kf=6.0e-28 af=2 lf=1 wf=1 ef=1 l=1u
dlr=0.01u w=10u dw=0.01u
.model Res2 R noise=0 kf=6.0e-28 af=2 lf=1 wf=1 ef=1 l=1u
dlr=0.01u w=10u dw=0.01u
R1 1 2 Res1 50
R2 2 3 Res2 50
In the above example, resistor R1 uses model Res1 which specifies noise=1,
and resistor R2 uses model Res2 which specifies noise=0. The noise analysis,
R1 has noise and R2 is noise-free. The default value of model parameter
NOISE is 1.
Instance of Parameter NOISE: Noiseless Resistor
Sometimes a noiseless resistor is required. The resistor can be noiseless even
if there are noise parameters set in the model. This is because the parameter
NOISE can be both an instance parameter and a model parameter. If NOISE is
set in both the instance and the model, the instance parameter will override the
model parameter. For example:
In the example, resistors R1 and R3 use model Res1 and resistors R2 and R4
use model Res2. R1 will have noise and R2 will be noiseless since these
resistors will use the noise parameter defined in the model. Resistor R3 is
noiseless and R4 will have the noise because the noise instance parameter is
specified for each.
Test Case for Both Thermal and Flicker Noise
The following example uses the .noise analysis syntax to test for both thermal
and flicker noise in a resistor.
v1 1 0 10 ac=1
R1 1 2 Res1 50 $noise=1
R2 2 3 Res2 50 $noise=1
R3 3 4 Res1 50 noise=0
R4 4 0 Res2 50 noise=1
.model Res1 R noise=1 kf=6.0e-28 af=2 lf=1 wf=1 ef=1
+ l=1u dlr=0.01u w=10u dw=0.01u
.model Res2 R noise=0 kf=6.0e-28 af=2 lf=1 wf=1 ef=1
+ l=1u dlr=0.01u w=10u dw=0.01u
.model Res3 R $noise=1
.model Res4 R $noise=1
.options post list
.ac dec 10 10k 100k
.noise v(2,3) v1 10
.print noise inoise onoise
.op
.end
In the .list file below, the output shows rs is thermal noise, 1/f is flicker noise:
**** resistor squared noise voltages (sq v/hz)
element 0:r1 0:r2 0:r3 0:r4
rs 5.145e-20 0. 0. 5.145e-20
1/f 5.924e-20 0. 0. 5.924e-20
total 1.107e-19 0. 0. 1.107e-19
rx 12.5000 37.5000 12.5000 12.5000
tc2l 0.0 μm ⁄ K
2 Resistance quadratic TC length coefficient
tc2w 0.0 μm ⁄ K
2 Resistance quadratic TC width coefficient
Usage
With Model Card:
rinstanceName (node1 node2) modelName instanceParameters
.model modelName r2_cmc modelParameters
Example
r137 (n1 n2) rnpoly1 w=1u l=10u
.model rnpoly1 r2_cmc
+ rsh=100.0 xl=0.2u xw=-0.05u
+ p3=0.12 q3=1.63 p2=0.014 q2=3.79
r138 (n2 n3) r=110.0 tc1=1.0e-3
Re1 Re2
n1 i1 i2 n2
i2n
dt
nc
R TH CTH
ITH=VI
c1 0 0 ∞ # contacts at n1 terminal
c2 0 0 ∞ # contacts at n2 terminal
fdxwinf 0.0 μm Finite doping width offset width value for wide
devices
ca 0.0 0.0 ∞ F ⁄ μm
2 Fixed capacitance per unit area
tc2l 0.0 μm ⁄ K
2 Resistance quadratic TC length coefficient
tc2w 0.0 μm ⁄ K
2 Resistance quadratic TC width coefficient
tc2rc 0.0 μm ⁄ K
2 Contact resistance quadratic TC
tc2vbv 0.0 ⁄K
2 Breakdown voltage quadratic TC
CMC R3 Usage
With Model card:
rinstanceName ([node1 node2 node3) [mname]
instanceParameters
.MODEL modelName r modelParameters
Without Model card
rname (node1 node2) r=resistanceValue [tc1=tc1Value]
[tc2=tc2Value]
Example
r137 n1 n2 n3 rnpoly1 w=1u l=10u
.model rnpoly1 r
+ level=5 rsh=100.0
r138 n2 n3 r=110.0 tc1=1.0e-3
Capacitance Model
You can use the .MODEL statement to include a capacitance model in your
HSPICE netlist. For a general description of the .MODEL statement, see
.MODEL in the HSPICE Reference Manual: Commands and Control Options.
Syntax
.MODEL mname C parameter=value
Parameter Description
Otherwise, the equation calculates the capacitance from the Leff, Weff, and
COX values:
CAPeff = M ⋅ SCALE ( element ) ⋅
[ Leff ⋅ Weff ⋅ COX + 2 ⋅ ( Leff + Weff ) ⋅ CAPSW ]
If you do not specify COX, but THICK is not zero, then:
DI ⋅ ε o
COX = -------------------- if DI not zero
THICK
ε ox
COX = ------------------- if DI=0
THICK
The following values apply to the preceding equation:
ε o = 8.8542149e-12
F
---------------
meter
ε ox = 3.453148e-11
F -
--------------
meter
If you specify only the model capacitance (CAP), then:
CAPeff = CAP ⋅ SCALE ( element ) ⋅ M
Parameter Description
Δt t - tnom
Parameter Description
tnom+273.15 + TNOM
To use this example, obtain the core model parameters from the manufacturer’s
data. Figure 4 on page 38 illustrates the required b-h loop parameters for the
model.
The model includes:
■
Core area
■
Length
■
Gap size
■
Core growth time constant
Example 2
This example is located in the following directory:
$installdir/demo/hspice/mag/bhloop.sp
AREA, (AC) cm2 1 Mean of the cross section for the magnetic
core. AC is an alias of AREA.
LX5
dm an
Slope of the anhysteretic magnetization curve,
dh
Parameter Description
Δt t - tnom
tnom=273.15 + TNOM
Parameter Description
V2 = ( j ⋅ ω ⋅ M ) ⋅ I1 + ( j ⋅ ω ⋅ L2 ) ⋅ I2
· he
m an = MS ⋅ ⎛⎝ coth ⎛ -----⎞ – -----⎞⎠
A
⎝ A ⎠ he
h e = h + ALPHA ⋅ m an
Parameter Description
h Magnetic field.
ALPHA This model parameter represents the coupling between the magnetic
domains.
dm an 1
= --------------------------------------------
dh A
3 ⋅ -------- – ALPHA
MS
The slope must be positive; therefore, the denominator of the above equation
must be positive. If the slope becomes negative, an error message appears.
Anhysteretic magnetization represents the global energy state of the material, if
the domain walls could move freely, but the walls are displaced and bent in the
material.
If you express the bulk magnetization (m) as the sum of an irreversible
component (due to wall displacement), and a reversible component (due to
domain wall bending), then:
( m an – m ) dm an dm⎞
= ------------------------ + C ⋅ ⎛⎝
dm
– ⎠
dh K dh dh
-or-
dm ( m an – m ) C dm an
= ----------------------------- + ------------- ⋅
dh (1 + C) ⋅ K 1 + C dh
Solving the above differential equation obtains the bulk magnetization value, m.
The following equation uses m to compute the flux density (b):
b = μ0 ⋅ ( h + m )
The following values apply to the preceding equation:
■
μ0 The permeability of free space, is 4π ⋅ 10 –7
■
The units of h and m are in amp/meter.
■
The units of b are in Tesla (Wb/meter2).
Example
This example demonstrates the effects of varying the ALPHA, A, and K model
parameters, on the b-h curve.
■
Figure 5 on page 44 shows b-h curves for three values of ALPHA.
■
Figure 6 on page 44 shows b-h curves for three values of A.
■
Figure 7 on page 45 shows b-h curves for three values of K.
Input File
This input file is located in the following directory:
$installdir/demo/hspice/mag/jiles.sp
Figure 9 Output Curves, Before Optimization (top), and After Optimization (bottom)
Diode model types include the junction diode model, and the Fowler-Nordheim
model. Junction diode models have two variations: geometric and
nongeometric.
Diode Types
Use the geometric junction diode to model:
■
IC-based, standard silicon-diffused diodes
■
Schottky barrier diodes
■
Zener diodes.
Use the geometric parameter to specify dimensions for pn junction poly and
metal capacitance for a particular IC process technology.
Use the non-geometric junction diode to model discrete diode devices, such as
standard and Zener diodes. You can use the non-geometric model to scale
currents, resistances, and capacitances by using dimensionless area
parameters.
The Fowler-Nordheim diode defines a tunneling current-flow, through
insulators. The model simulates diode effects in nonvolatile EEPROM memory.
Control options, related to the analysis of diode circuits and other models,
include:
■ DCAP
■
DCCAP
■
GMIN
■
GMINDC
■
SCALE
■
SCALM
Include DCAP=<val> in the .MODEL statement for the diode to override the
global depletion capacitance equation that the .OPTION DCAP=<val>
statement selects.
■
Geometric junction
■
Capacitance parameters
Both LEVEL=1 and LEVEL=3 junction diode models share the same element
parameter set. Poly and metal capacitor parameters of LM, LP, WM and WP, do
not share the same element parameter.
Element parameters have precedence over model parameters, if you repeat
them as model parameters in the .MODEL statement.
Parameters common to both element and model statements are:
AREA, PJ, M, LM, LP, WM, WP, W, and L.
Table 12 Junction Diode Element Parameters
Function Parameters
Temperature DTEMP
Parameter Description
mname Model name. The diode element uses this name to refer to the
model.
Parameter Description
Example
.MODEL D D (CO=2PF, RS=1, IS=1P)
.MODEL DFOWLER D (LEVEL=2, TOX=100, JF=1E-10, EF=1E8)
.MODEL DGEO D (LEVEL=3, JS=1E-4, JSW=1E-8)
.MODEL d1n750a D
+ LEVEL=1 XP =0.0 EG =1.1
+ XOI =0.0 XOM =0.0 XM =0.0
+ WP =0.0 WM =0.0 LP =0.0
+ LM =0.0 AF =1.0 JSW =0.0
+ PB =0.65 PHP =0.8 M =0.2994
+ FC =0.95 FCS =0.4 MJSW=0.5
+ TT =2.446e-9 BV =4.65 RS =19
+ IS =1.485e-11 CJO =1.09e-9 CJP =0.0
+ PJ =0.0 N =1.615 IK =0.0
+ IKR =1.100e-2 IBV =2.00e-2
Function Parameters
geometric capacitance L, LM, LP, SHRINK, W, WM, WP, XM, XOJ, XOM, XP, XW
(LEVEL=3 only)
noise AK, KF
Name (Alias)
Units Default Description
Name (Alias)
Units Default Description
JSW (ISP) amp/ 0.0 Sidewall saturation current, per unit junction
PJeff periphery.
■
For LEVEL=1:
JSWeff = PJeff⋅ JSW
■ For LEVEL=3:
JSWeff = PJeff⋅ JSW/SCALM
Name (Alias)
Units Default Description
IS (JS) amp/ LEVEL Saturation current per unit area. If the IS value
AREAeff 1= 1.0e- is less than EPSMIN, the program resets the
14 value of IS to EPSMIN, and shows a warning
message.
LEVEL EPSMIN default=1.0e-28
3= 0.0
If the value of IS is too large, the program
displays a warning.
■ For LEVEL=1:
ISeff = AREAeff⋅ IS
■ For LEVEL=3:
ISeff = AREAeff ⋅ IS/SCALM2
Name (Alias)
Units Default Description
Note: If you use a diode model than does not specify an AREA, then
AREA defaults to 1, and RS is in units of ohms. If you specify the
AREA in square meters (m2) in the netlist, then the units of RS are
ohms/m2.
Table 15 Junction Capacitance Parameters
SHRINK
LEVEL=1 Scaling
LEVEL=1 uses the AREA and M Element parameters to scale the following
element and model parameters: IK, IKR, JS, CJO, and RS.
For AREA and M, default=1.
This element is not geometric, because it uses dimensionless values to
measure both the area (AREA) and the periphery (PJ). The .OPTION SCALE
and .OPTION SCALM statements do not affect these parameters.
Here is the process:
• Diode models multiply the periphery junction parameter by M (the multiplier
parameter) to scale a dimensionless periphery junction:
PJeff = PJ ⋅ M
The diode models then use PJeff to scale CJP (the zero-bias junction
capacitance), and the sidewall saturation current (JSW).
CJPeff = PJeff ⋅ CJP
JSWeff = PJeff ⋅ JSW
JTUNSWeff = PJeff ⋅ JTUNSW
The models use the AREA and M values to obtain AREAeff.
AREAeff = AREA ⋅ M
Models multiply CJO, IK, IKR, IBV, and IS by AREAeff to obtain their
effective scaled values. However, diode models divide RS by AREAeff.
IKeff = AREAeff ⋅ IK
IKReff = AREAeff ⋅ IKR
ISeff = AREAeff ⋅ IS
RSeff = RS/AREAeff
CJOeff = CJO ⋅ AREAeff
JTUNeff = JTUN ⋅ AREAeff
LEVEL=3 Scaling
The SCALM, SCALE, SHRINK, and M parameters affect LEVEL=3 scaling.
■
SCALE affects the following LEVEL=3 element parameters:
AREA, LM, LP, PJ, WM, WP, W, L
■
SCALM affects the following model parameters:
AREA, IBV, IK, IKR, IS, PJ, JSW, RS, CJO, CJP, LM, LP,
WP, XM, XP, W, L, XW, JTUN, JTUNSW
Diode Current
Figure 10 shows the direction of current flow through the diode. Use either
I(D1) or I1(D1) syntax to print the diode current.
If the voltage on node1 is 0.6V greater than the voltage on node2, then the
diode is forward biased or turned on. The anode is the p-doped side of a diode,
and the cathode is the n-doped side.
Anode
rs
+
vd id cd
-
Cathode
Figure 11 Equivalent Circuit for Diode in Transient Analysis
Anode
rs
gd cd
Cathode
Anode
rs
Inrs
id cd
Inrd
Cathode
Variable Parameter
Variable Parameter
GAP1 eV/× 7.02e-4 First bandgap correction factor. From Sze, alpha
term.
TPB (TVJ) V/× 0.0 PB temperature coefficient. If you set the TLEVC
parameter to 1 or 2, TPB overrides default
temperature compensation.
Variable Definition
f frequency
gd diode conductance
id diode DC current
Quantity Definition
t temperature in °Kelvin
Δt t - tnom
Forward Bias: vd ≥ – 10 ⋅ vt
vd
---------------
id = ISeff ⋅ ⎛⎝ e N ⋅ vt – 1⎞⎠
vd = v node1 – v node2
In reverse-bias, the anode (node1) is more negative than the cathode. The
diode is turned off and conducts a small leakage current.
Breakdown: vd < – BV
– vd
– ⎛⎝ ---------------------------⎞⎠
vd + BVeff
N ⋅ vt ⎛ ---------------------------
NTUN ⋅ vt ⎞
-
id = – ISeff ⋅ e – JTUNeff ⋅ ⎜ e – 1⎟
⎝ ⎠
The device model adjusts the BV parameter to obtain BVeff :
– BV
---------------
ibreak = – ISeff ⋅ ⎛ e N ⋅ vt – 1⎞
⎝ ⎠
Otherwise:
IBVeff = ibreak
Most diodes do not behave as ideal diodes. The IK and IKR parameters are
called high-level injection parameters. They tend to limit the exponential
increase in current.
Note: Diode models use the exponential equation in both the forward
and reverse regions.
Forward Bias
id1
id = ----------------------------------
1/2
-
⎛ id1
1 + ------------- ⎞
⎝ IKeff⎠
Reverse Bias
id1
id = --------------------------------------
1/2
-
1 + ⎛ -----------------⎞
id1
⎝ IKReff⎠
For vd ≥ – BVeff :
vd – vd
⎛ ---------------
N ⋅ vt ⎞ ⎛ ---------------------------
NTUN ⋅ vt ⎞
-
id1 = ISeff ⋅ ⎜ e – 1⎟ – JTUNeff ⋅ ⎜ e – 1⎟
⎝ ⎠ ⎝ ⎠
Otherwise:
vd ⎛ vd + BVeff⎞
– ⎝ ---------------------------⎠ – vd
⎛ ---------------
N ⋅ vt ⎞ NBV ⋅ vt ⎛ ---------------------------
NTUN ⋅ vt ⎞
-
id1 = – ISeff ⋅ ⎜ e – 1⎟ – ISeff ⋅ e – JTUNeff ⋅ ⎜ e – 1⎟
⎝ ⎠ ⎝ ⎠
From DC measurements of the forward-biased diode characteristics, you can
estimate:
■
Reverse-saturation current (IS)
■
Emission coefficient (N)
■ Model parameter (RS)
You can determine N from the slope of the diode characteristic in the ideal
region. In most cases, the emission coefficient is the value of the unit, but is
closer to 2 for MOS diodes.
At higher bias levels, the diode current deviates from the ideal exponential
characteristic, due to ohmic resistance in the diode, and the effects of high-level
injection. The deviation of the actual diode voltage (from the ideal exponential
characteristic), at a specific current, determines the value of RS. In practice,
simulation of diode device models estimates RS at several values of id, and
averages them, because the value of RS depends upon diode current.
vd – MJ
cdepa = CJeff ⋅ ⎛ 1 – -------⎞
⎝ PB⎠
vd ≥ FC ⋅ PB
vd
1 – FC ⇒1 ( + MJ ) + MJ ⋅ -------
PB
cdepa = CJeff ⋅ --------------------------------------------------------------------------
-
( 1 – FC ) ( 1 + MJ )
The capacitance formula for the junction periphery is: vd < FCS ⋅ PHP
vd – MJSW
cdepp = CJPeff ⋅ ⎛ 1 – ------------⎞
⎝ PHP⎠
vd ≥ FCS ⋅ PHP
vd
1 – FCS Þ ( 1 + MJSW ) + MJSW ⋅ ------------
PHP
cdepp = CJPeff ⋅ ---------------------------------------------------------------------------------------------------
( 1 – FCS ) ( 1 + MJSW )
cdep = cdepa + cdepp
DCAP = 2 (default)
The capacitance formula for the total depletion is:
vd < 0
vd – MJ vd – MJSW
cdep = CJeff ⋅ ⎛ 1 – -------⎞ + CJPeff ⋅ ⎛ 1 – ------------⎞
⎝ PB⎠ ⎝ PHP⎠
vd ≥ 0
ε ox
cpoly = ⎛ -------------⎞ ⋅ ( WMeff + XMeff ) ⋅ ( LMeff + XMeff ) ⋅ M
⎝ XOM⎠
4 ⋅ k ⋅ t 1/2
inrs = ⎛ ---------------------⎞
⎝ RSeff ⎠
The ind current source models the shot and flicker noise of the diode. The
following equation defines ind:
KF ⋅ id AF 1 / 2
ind = ⎛⎝ 2 ⋅ q ⋅ id + --------------------------⎞⎠
f
t2
eg ( t ) = 1.16 – 7.02e-4 ⇒------------------------
t + 1108.0
TLEV = 2
tnom 2
egnom = EG – GAP1 ⇒-----------------------------------
tnom + GAP2
t2
eg ( t ) = EG – GAP1 ⇒-----------------------
t + GAP2
Leakage Current Temperature Equations
facln
-------------
JS ( t ) = JS ⋅ e N
facln
-------------
JSW ( t ) = JSW ⋅ e N
TLEV = 0 or 1
TLEV = 2
eg ( t )
facln = ---------------------- – ------------ + XTI ⋅ ln ⎛ -------------⎞
egnom t
vt ( tnom ) vt ( t ) ⎝ tnom⎠
TLEV = 2
eg ( t )
faclnt = KEG ⋅ ⎛---------------------- – ------------ ⎞ + XTITUN ⋅ ln ⎛ -------------⎞
egnom t
⎝ vt ( tnom ) vt ( t ) ⎠ ⎝ tnom⎠
egnom eg ( t )
PB ( t ) = PB ⋅ ⎛ -------------⎞ – vt ( t ) ⇒ 3 ⋅ ln ⎛ -------------⎞ + ---------------------- – ------------
t t
⎝ tnom⎠ ⎝ tnom⎠ vt ( tnom ) vt ( t )
eg ( t )
PHP ( t ) = PHP ⋅ ------------- – vt ( t ) ⇒ 3 ⋅ ln ⎛ -------------⎞ + ---------------------- – ------------
t t egnom
tnom ⎝ tnom ⎠ vt ( tnom ) vt ( t )
TLEVC = 1 or 2
PB ( t ) = PB – TPB ⇒Δt
PHP ( t ) = PHP – TPHP ⇒Δt
TLEVC = 3
PB ( t ) = PB + dpbdt ⋅ Δt
PHP ( t ) = PHP + dphpdt ⋅ Δt
PB ( t )
CJ ( t ) = CJ ⋅ 1 + MJ ⋅ ⎛⎝ 4.0e-4 ⇒Δt – -------------- + 1⎞⎠
PB
PHP ( t )
CJSW ( t ) = CJSW ⋅ 1 + MJSW ⋅ ⎛ 4.0e-4 ⇒Δt – ------------------- + 1⎞
⎝ PHP ⎠
TLEVC = 1
CJ ( t ) = CJ ⋅ ( 1 + CTA ⋅ Δt )
CJSW ( t ) = CJSW ⋅ ( 1 + CTP ⋅ Δt )
TLEVC = 2
PB MJ
CJ ( t ) = CJ ⋅ ⎛ --------------⎞
⎝ PB ( t )⎠
PHP MJSW
CJSW ( t ) = CJSW ⋅ ⎛ -------------------⎞
⎝ PHP ( t )⎠
TLEVC = 3
Δt
CJ ( t ) = CJ ⋅ ⎛ 1 – 0.5 ⇒dpbdt Þ -------⎞
⎝ PB⎠
Δt
CJSW ( t ) = CJSW ⋅ ⎛ 1 – 0.5 ⋅ dphpdt ⋅ ------------⎞
⎝ PHP⎠
Grading Coefficient Temperature Equation
MJ ( t ) = MJ ⋅ ( 1 + TM1 ⋅ Δt + TM2 ⋅ Δt 2 )
Resistance Temperature Equations
RS ( t ) = RS ⋅ ( 1 + TRS ⋅ Δt )
Parameter Description
Parameter Description
na Positive terminal (anode) node name. The series resistor for the
equivalent circuit is attached to this terminal.
area Diode area. In the model card, AB can use this value.
pgate Length of the side-wall in the AB diffusion area, which is under the
gate. In the model card, LG uses this value.
dtemp The difference between the element temperature and the circuit
temperature in degrees celsius. The default is DTA.
ic Initial voltage across a diode element. Use this value when you
specify the UIC option in the .TRAN statement. The .IC statement
overrides this value.
.OPTION list Prints the updated temperature parameters for the JUNCAP diode
model.
JUNCAP1 Model
To use this model, specify:
.MODEL mname D LEVEL=4 [keyword=val]
Parameter Description
mname Model name. The diode element uses this name to refer to the
model.
Parameter Description
Example
.model MD D LEVEL=4
+ AB=2E-12 LS=2E-6 LG=1.3E-6 DTA=0 TR=30 VR=0.3 JSGBR=1.2e-3
+ JSDBR=1.3e-3 JSGSR=1.1e-3 JSDSR=1.3e-3 JSGGR=1.4e-3
+ JSDGR=1.4e-3 NB=1.6 NS=1.3 NG=1.3 VB=0.9 CJBR=1.2e-12
+ CJSR=1.2e-12 CJGR=1.3e-12 VDBR=1.6 VDSR=1.3 VGDR=1.2 PB=0.5
+ PS=0.6 PG=0.4
Model Parameters
The JUNCAP1 model parameters are listed in Table 22.
Table 22 JUNCAP1 Model Parameters
Theory
This section summarizes the elementary physics of a junction diode. Refer to
semiconductor textbooks for additional information.
You can represent the current voltage characteristics as follows:
exp ⎛ ------⎞ – 1
2 qv
J = {J d ( n i + ( J g ( n i, V ) ) }⋅
⎝ kT⎠
3---
–Eg
ni T ⋅ exp ⎛ ---------⎞
2
⎝ 2kT⎠
Eg J Energy gap.
T K Temperature.
For V<VD, this equation describes the charge of the junction capacitance:
V 1–P
Q = Q j 1 – ⎛ 1 – -------⎞
⎝ V D⎠
Nomenclature
Table 23 lists the electrical variable parameters:
Table 23 JUNCAP Model Electrical Variable Parameters
Note: The model card lists the parameters. See JUNCAP model syntax
earlier in this chapter.
Variable or Program
No. Parameter Name Units Description
Variable or Program
No. Parameter Name Units Description
(I = IA = -IK)
ON/OFF Condition
Solving a circuit involves successive calculations. The calculations start from a
set of initial guesses for the electrical quantities of the non-linear elements. The
devices start in the default state.
T KD = T 0 + T A + DT A
T KR
V TR = k ⋅ ---------
q
T KD
V TD = K ⋅ ----------
q
( 7.02 ⋅ 10e – 4 ⋅ T KR ⋅ T KR )
V gR = 1.16 – --------------------------------------------------------------------------
( 1108.0 + T KR )
( 7.02 ⋅ 10e – 4 ⋅ T KR ⋅ T KD )
V gD = 1.16 – --------------------------------------------------------------------------
( 1108.0 + T KD )
T KD 1.5 V gR V gD
F TD = ⎛ ----------⎞ ⋅ exp ⎛⎝ ------------------------ – ------------------------- ⎞⎠
⎝ T KR ⎠ ( 2 ⋅ V TR ) ( 2 ⋅ V TD )
Internal Reference
The following equations specify the internal reference parameters for the
bottom component:
V DBR ⋅ T KD
⎛ V = ----------------------------------------------------------
-⎞
⎝ DB T KR – 2 ⋅ V TD ⋅ 1nF TD⎠
( V DBR – V R ) PB
C JB = C JBR ⋅ A B ⋅ ⎛ ------------------------------⎞
⎝ V DB ⎠
V DB PB
I SGB = J SGBR ⋅ F TD ⋅ A B ⋅ ⎛⎝ ------------------------------⎞⎠
( V DBR – V R )
I SDB = J SDBR ⋅ F TD ⋅ F TD ⋅ A B
Locos-edge and gate-edge components use similar formulations:
■ Replace the B (bottom) index with S (locos-edge) or G (gate-edge).
■
Replace the AB (bottom) area with LS (locos-edge) or LG (gate-edge).
For the locos-edge:
V DSR ⋅ T KR
V DS = -----------------------------------------------------------
T KR – 2 ⋅ V TD ⋅ 1nF TD
( V DSR – V R ) PS
C JS = C JSR ⋅ L S ⋅ ⎛ ------------------------------⎞
⎝ V DS ⎠
V DS PS
I SGS = J SGSR ⋅ F TD ⋅ L S ⋅ ⎛ ------------------------------⎞
⎝ ( V DSR – V R )⎠
I SDS = J SDSR ⋅ F TD ⋅ F TD ⋅ L S
For the gate-edge:
V DGR ⋅ T KD
V DG = -----------------------------------------------------------
T KR – 2 ⋅ V TD ⋅ 1nF TD
( V DGR – V R ) PG
C JG = C JGR ⋅ L G ⋅ ⎛ -------------------------------⎞
⎝ V DG ⎠
V DG PG
I SGS = J SGGR ⋅ F TD ⋅ L G ⋅ ⎛ -------------------------------⎞
⎝ ( V DGR – V R )⎠
I SDS = J SDGR ⋅ F TD ⋅ F TD ⋅ L G
Note: The remainder of this section shows the equations only for the
bottom component.
C LB = C JB ( 1 – F CB ) –P B
Q LB = Q JDB ( 1 – ( 1 – F CB ) ( 1 – P B ))
( 1 – ( 1 – V ) ) ( 1 – PB )
Q JBV = Q JCB ⋅ ⎛ -------------------------------⎞ V< V LB
⎝ V DB ⎠
1 + ( P B ( V – V LB ) )
Q LB + C LB ( V – V LB ) ⋅ ⎛ ---------------------------------------------------⎞ V >= V LB ( 12.63 )
⎝ 2 ⋅ V DB ⋅ ( 1 – F CB )⎠
Use similar expressions for the locos-edge (QJSV) and gate-edge (QJGV)
charges.
The following equation describes the total charge characteristic:
Q = Q JBV + Q JSV + Q JGV
From Equation 12.63 (above), you can use elementary mathematics to derive
simple equations for the capacitance of the bottom area:
1 – v PB
C JBV = C JB ⋅ ⎛ 1 ⁄ ⎛ -----------⎞ ⎞ V < V LB
⎝ ⎝ V DB ⎠ ⎠
( V – V LB )
C LB + C LB ⋅ P B ⋅ ⎛ -----------------------------------------⎞ V ≥ V LB
⎝ V DB ⋅ ( 1 – F CB )⎠
( V DB – V ) PB
I GB = I SGB ⋅ ⎛ ------------------------⎞ ⋅ ⎛ exp ⎛ ⎛ ------------------------⎞ – 1⎞ ⎞
V
V ≤V DB
⎝ V DB ⎠ ⎝ ⎝ ⎝ N B ⋅ V TD⎠ ⎠⎠
0 V > V DB
■
The first relation, concerning the diffusion component, is valid over the
whole operating range.
■ The second relation, describing the generation current, shows an unlimited
increase in the derivative of this function, at V=VDB.
Therefore, the power function merges at V=0.0 with a hyperbolic function in the
forward-bias range. Simulating the model then divides the exponential part by
exp ( V ⁄ ( N B ⋅ V TD ) ) . This enables a gradual decrease in the generation-current
component.
This calculation uses the hyperbolic function:
I HYP = F SB ( V + V AB ) –B
The B parameter controls the decrease in current for voltages V>0.0 for all
generation components. The model sets B to a fixed value of 2. The continuity
constraints of the function and derivative in the merge point lead to the
following relations for FSB and VAB:
V DB
V AB = B ⋅ ----------
PB
F SB = I SGB ⋅ V AB B
The generation current voltage characteristic in the forward region, becomes:
I GB = F SB ⁄ ( ( V + V AB ) B ) ⋅ ( 1 – exp ( – v ) ⁄ ( N B ⋅ V TD ) )
V DB – V PB
I GB = I SGB ⋅ ⎛ -------------------⎞ ⋅ exp ⎛ -----------------------------⎞ – 1
V
V ≤0
⎝ V DB ⎠ ⎝ ( N B ⋅ V TD )⎠
B
Isgb ⋅ ⎛ -------------------------⎞ ⋅ ⎛ 1 – exp ⎛ ----------------------------⎞ ⎞
Vab –V
V > 0.0
⎝ ( V + Vab )⎠ ⎝ ⎝ ( Nb ⋅ Vtd )⎠ ⎠
JUNCAP2 Model
The JUNCAP2 model is a compact MOS model intended to describe the
behavior of the diodes that are formed by the source, drain, or well-to-bulk
junctions in MOSFETs. It is the successor to the JUNCAP1 model, and was
developed by Koninklijke Philips Electronics N.V.
Physical Effects
The following physical effects have been included in the JUNCAP2 model:
■
Geometrical scaling
■
Depletion capacitance
■
Ideal current
■
Shockley-Read-Hall current
■
Trap-assisted tunneling current
■
Band-to-band tunneling current
■
Avalanche breakdown
■
Noise
For a full description of the JUNCAP2 model, see
http://www.semiconductors.philips.com/Philips_Models/.
Usage in HSPICE
The JUNCAP2 model is LEVEL=6 in the Synopsys diode models. Each version
can be identified with model parameter VERSION.
To use this model, specify:
.MODEL mname D LEVEL=6 [keyword=val]
Parameter Description
mname Model name. The diode element uses this name to refer to the model.
Example
.MODEL NDIO D LEVEL=6 VERSION=200.33
+ TYPE=1.000E+00 TRJ=21.000E+00 DTA=0.000E+00 IMAX=1.000E+03
+ CJORBOT=1.000E-03 CJORSTI=1.000E-09 CJORGAT=1.000E-09
+ VBIRBOT=1.000E+00 VBIRSTI=1.000E+00 VBIRGAT=1.000E+00
+ PBOT=500.000E-03 PSTI=500.000E-03 PGAT=500.000E-03
+ PHIGBOT=1.160E+00 PHIGSTI=1.160E+00 PHIGGAT=1.160E+00
+ IDSATRBOT=1.000E-12 IDSATRSTI=1.000E-18 IDSATRGAT=1.000E-18
+ CSRHBOT=100.000E+00 CSRHSTI=100.000E-06 CSRHGAT=100.000E-06
+ XJUNSTI=100.000E-09 XJUNGAT=100.000E-09 CTATBOT=100.000E+00
+ CTATSTI=100.000E-06 CTATGAT=100.000E-06 MEFFTATBOT=250.000E-03
+ MEFFTATSTI=250.000E-03 MEFFTATGAT=250.000E-03
+ CBBTBOT=1.000E-12 CBBTSTI=1.000E-18 CBBTGAT=1.000E-18
+ FBBTRBOT=1.000E+09 FBBTRSTI=1.000E+09 FBBTRGAT=1.000E+09
+ STFBBTBOT=-1.000E-03 STFBBTSTI=-1.000E-03 STFBBTGAT=-1.000E-03
+ VBRBOT=10.000E+00 VBRSTI=10.000E+00 VBRGAT=10.000E+00
+ PBRBOT=4.000E+00 PBRSTI=4.000E+00 PBRGAT=4.000E+00
Model Parameters
Table 25 lists the JUNCAP2 model parameters.
Table 25 JUNCAP2 Model Parameters
Capacitance parameters
Ideal-current parameters
Shockley-Read-Hall parameters
MEFFTATBOT - 0:25 Effective mass (in units of m0) for trap-assisted tunneling
of bottom component
MEFFTATSTI - 0:25 Effective mass (in units of m0) for trap-assisted tunneling
of STI-edge component
MEFFTATGAT - 0:25 Effective mass (in units of m0) for trap-assisted tunneling
of gate-edge component
Reverse Bias
vd < 0
ER ⋅ TOX
vd 2 ---------------------------
id = – AREAeff ⇒JR ⇒⎛ ------------⎞ ⇒e vd
⎝ TOX⎠
ID
RS
Qd
A K
Qr
Temperature Effects
The actual simulation temperature is denoted by TEMP (in °C). The
temperature at which the parameters are determined is specified by TREF
(in °C.)
■
Conversion to Kelvins
TK = TEMP + 273.15 + DTA
TRK = TREF + 273.15
■
Thermal Voltages
– 23 –1
k = 1.3806226 ⋅ 10 JK
– 19
q = 1.6021918 ⋅ 10 C
V T = ⎛ ---⎞ ⋅ T K
k
⎝ q⎠
V TR = ⎛⎝ ---⎞⎠ ⋅ T RK
k
q
■
Depletion Capacitances
TK 3
F = ⎛ ---------⎞ ⋅ exp VG ⋅ ⎛ --------- – ------⎞
1 1
⎝ T RK⎠ ⎝ V TR V T⎠
VD
VD T = --------- – 1n ( F ) ⋅ V T
V TR
VD P
CJ T = CJ ⋅ ⎛ -----------⎞
⎝ VD T⎠
■
Transit Times
T K 1.8
TAU T = TAU ⎛ ---------⎞
⎝ T RK⎠
■
Saturation Current
T K 1.8
IS T = IS ⋅ ⎛ ---------⎞ ⋅ exp -------- ⋅ ⎛ --------- – ------⎞
VG 1 1
⎝ T RK⎠ N ⎝ V TR V T⎠
■
Shockley-Read-Hall generation and trap assisted tunneling
TK 3⁄ 2
T up = ⎛ ---------⎞ ⋅ exp -------------------------- ⋅ ⎛ --------- – ------⎞
VG + VLC 1 1
⎝ T RK⎠ 2 ⎝ V TR V T⎠
CSRH T = CSRH ⋅ T up
CTAT T = CTAT ⋅ T up
3⁄ 2
ETAT T = 70.8 ⋅ T K
■
Band-to-Band Tunneling
CBBT T = CBBT (temperature independent)
–4 2
⎛ 7 4.21 ⋅ 10 ⋅ T K ⎞
F0 = 1.9 ⋅ 10 ⋅ ⎜ 1.04 – ------------------------------------------------⎟
⎝ 636 + T K ⎠
■
Avalanche Multiplication
dT = TEMP + DTA – 25° C
6
Bn = 1.23 ⋅ 10
–4 –6 2
Bn T = Bn ⋅ 1 + 7.2 ⋅ 10 ⋅ dT – 1.6 ⋅ 10 ⋅ dT
■
Breakdown
T K 0.1
VBR T = VBR ⋅ ⎛ ---------⎞
⎝ T RK⎠
VD T + VBR T 1 – P
EMVBRT = EMVBR ⋅ ⎛⎝ -------------------------------⎞⎠
VD + VBR
■
Resistance
T K PTRS
RS T = RS ⋅ ⎛⎝ ---------⎞⎠
T RK
Model Constants and Parameter-Related Constants
K = 0.01
KET = 0.1
ETM = 3
■
Maximum Electric Field and Depletion Layer width at zero bias:
EMVBR T
E 0 = --------------------------------------
1–P
-
⎛ 1 + VBR T⎞
--------------
⎝ VD T ⎠
VD T
W 0 = ------------------------------
E0 ⋅ ( 1 – P )
Diode Currents
First, the maximum reverse junction voltage is defined.
■
Above this voltage the current will be extrapolated on a logarithmic scale.
– 0.99 ⋅ VBR T, V AK1 < – 0.99 ⋅ VBR T
Vj =
V AK1 , V AK1 ≥ – 0.99 ⋅ VBR T
■
Ideal Forward Current
⎧ Vj ⎫
Id f = IS T ⎨ exp ⎛ -----------------⎞ – 1 ⎬
⎝ N ⋅ V T⎠
⎩ ⎭
■
Maximum Electric Field and Depletion Layer Width
⎧⎛ Vj ⎞ 2 ⎛ Vj ⎫ Vj ⎞
⎨ ⎝ 1 – ---------- - + ----------- ⋅ K⎞ ⎬ ⋅ ⎛ 1 – ---------- -
VD T ⎠ ⎝ VD T ⎠ ⎝ VD T⎠
⎩ ⎭
VD j = ----------------------------------------------------------------------------------------------------------
2
(1 – P)
E m = E 0 ⋅ VD j
P
W d = W 0 ⋅ VD j
■
Shockley-Read-Hall Generation
I srh = CSRH T ⋅ ( W d – W 0 )
■
Trap-Assisted Tunneling
E0 E0 2
----------------- + ETM – ⎛ ----------------- – ETM⎞ + KET
ETAT T ⎝ ETAT T ⎠
ET 0 = --------------------------------------------------------------------------------------------------------------
2
Em Em 2
----------------- + ETM – ⎛ ----------------- – ETM⎞ + KET
ETAT T ⎝ ETAT T ⎠
ET = --------------------------------------------------------------------------------------------------------------
2
⎧ 2 2 ⎫
⎪ exp ( ET ) – exp ( ET 0 ) ⎪
I tat = CTAT T ⋅ W d ⋅ ⎨ ----------------------------------------------------------- ⎬
⎪ Em ⎪
-----------------
⎩ ETAT T ⎭
■
Non-Ideal Forward Current including Tunneling
⎧ Em ⎫ VT
Is lf = CSRH T ⋅ ⎨ 6.28 + 38.58 ⋅ ⎛ -----------------⎞ ⋅ exp ( ET 0 ) ⎬ ⋅ -------
2
⎝ ETAT T⎠
⎩ ⎭ Em
Vj
exp ⎛ -----------------⎞ – 1
⎝ N ⋅ V T⎠
I lf = Is lf ⋅ ------------------------------------------------------------------------------------------------------ ⋅ exp ⎛⎝ ---------------------------⎞⎠
VLC
Vj 2 ⋅ N ⋅ VT
4 ⋅ exp ⎛ ---------------------------⎞ + exp ⎛ ---------------------------⎞
VLC
⎝ 2 ⋅ N ⋅ V T⎠ ⎝ 2 ⋅ N ⋅ V T⎠
■
Band-to Band-Tunneling
– CBBT T ⋅ V j
I bbt = -----------------------------------------------
F0⎞ 1.5
⎛ ------ ⋅ exp ⎛ -------⎞
F0
-
⎝ E m⎠ ⎝ E m⎠
■
Avalanche Multiplication
Em 2 Bn T Bn T
μ = 0.3295 ⋅ ⎛⎝ -----------------------⎞⎠ ⋅ exp ----------------------- – ---------
EMVBR T EMVBR T E m
■
Total Diode Current
1 + exp ( – 2 ⋅ μ)
( Id f + I lf – I srh ) ⋅ --------------------------------------- – I bbt + I tat ⋅ exp ( – μ)
2
I d = ------------------------------------------------------------------------------------------------------------------------------------------
1 – 2 ⋅ μ ⋅ {1 + exp ( – 2 ⋅ μ) }
■
Extrapolation of the Reverse Current
I dBR = I d at V j = – 0.99 VBR T
dI d
G dBR = -------- at V j = – 0.99 VBR T
dV j
Transient Model
Transient behavior is modeled using the DC equations.
■
Diffusion charge
Q D = TAU T ⋅ Id f
■ Depletion charge
⎛ --1-⎞
⎝ P⎠
FC = 1 – ⎛ -------------⎞
1+P
⎝ 3 ⎠
VD T
Q AT = CJ T ⋅ ⎛ ------------⎞
⎝ 1 – P⎠
V L = FC ⋅ VD T
–P
C L = CJ T ⋅ ( 1 – FC )
(1 – P)
Q L = Q AT ⋅ {1 – ( 1 – FC ) }
Or, if V AK1 ≥ V L
⎧ P ⋅ ( V AK'1 – V L ) ⎫
Q T = V L + C L ⋅ ( V AK'1 – V L ) – ⎨ 1 + -------------------------------------------------- ⎬
⎩ 2 ⋅ VD T ⋅ ( 1 – FC ) ⎭
AC Linearized model
Using the appropriate definitions for the various circuit elements leads to the
following equations:
1
R D = ---------------------------------
dID ⁄ dV AK1
Where dID ⁄ dV AK1 is the first derivative of the total diode current with respect
to the Internal voltage Internal voltage V AK1 .
The capacitances are defined as:
–P
⎧ V AK1 ⎫
C T = CJ T ⋅ ⎨ 1 – ⎛ ------------⎞ ⎬ for V AK1 < V L
⎝ VD T ⎠
⎩ ⎭
⎧ P ⋅ ( V AK1 – V L ) ⎫
C T = C L ⋅ ⎨ 1 + ---------------------------------------- ⎬ for V AK1 ≥ V L
⎩ VD T ⋅ ( 1 – FC ) ⎭
Id f + IS T
C 1 = C T + TAU T ⋅ ⎛ ---------------------⎞
⎝ N ⋅ VT ⎠
RD iNRS
RS
C1
A K1 K
iN
Noise Model
For noise analysis, noise sources are added to the small signal model as
shown in Figure 15. In these equations f represents the operation frequency of
the transistor and is the Δf
bandwidth. When Δf is taken as 1 Hz, a noise density is obtained.
■
Thermal noise
2 4 ⋅ K ⋅ T K ⋅ Δf
iN RS = -----------------------------------------
RS T
■
Current noise (shot noise and 1/f noise)
2 Id f AF Δf
iN = 2 ⋅ q ⋅ Id f ⋅ Δf + KF + MULT ⋅ ----------------
- ⋅ -----
MULT f
You must define the values for all parameters used in this model in either
a .PARAM statement or the .SUBCKT call. Circuit simulation then replaces the
diode statements with the call to the SDIODE subcircuit; for example,
XDS 14 1048 SDIODE SIS=67.32 SCJA=67.32 SRS=1.2285E-2
Describes how to use JFET and MESFET models in HSPICE circuit simulations.
Overview of JFETs
JFETs form by diffusing a gate diode between the source and drain. MESFETs
form by applying a metal layer over the gate region, and creating a Schottky
diode. To control the flow of carriers, both technologies modulate the gate diode
depletion region. These field effect devices are called bulk semiconductors and
are in the same category as bipolar transistors. Compared to surface effect
devices such as MOSFETs, bulk semiconductors have higher gain, because
bulk semiconductor mobility is always higher than surface mobility.
Enhanced characteristics of JFETs and MESFETs, relative to surface effect
devices, include lower noise generation rates and higher immunity to radiation.
These advantages have created the need for newer and more advanced
models.
Features for JFET and MESFET modeling include:
■
Charge-conserving gate capacitors
■ Backgating substrate node
■
Mobility degradation due to gate field
■
Computationally efficient DC model (Curtice and Statz)
■
Subthreshold equation
■
Physically correct width and length (ACM)
GaAs model LEVEL=3[1] assumes that GaAs device velocity saturates at very
low drain voltages. This model includes drain voltage induced threshold
modulation and user-selectable materials constants. These features let you use
the model for other materials, such as silicon, indium phosphide, and gallium
aluminum arsenide. The models that have been provided include a revised
Curtice model[2], and a TriQuint model (TOM) that extends the earlier Statz
model.
Specifying a Model
To specify a JFET or MESFET model, use a JFET element statement and a
JFET model statement. The model parameter LEVEL selects either the JFET or
MESFET model. LEVEL=1 and LEVEL=2 select the JFET, and LEVEL=3
selects the MESFET. Different submodels for the MESFET LEVEL=3 equations
are selected using the parameter SAT.
Parameter Description
Parameter Description
SAT=2 Cubic approximation of Curtice model with gate field degradation (Statz
model)
Parameter Description
You can use CAPOP=0, 1, 2 for any model level. CAPOP=1 and 2 are most often
used for the MESFET LEVEL=3 model.
The ACM model parameter selects the area calculation method:
Parameter Description
Example 1
The following example selects the n channel MESFET model, LEVEL=3. It
uses the SAT, ALPHA, and CAPOP=1 parameter:
J1 7 2 3 GAASFET
.MODEL GAASFET NJF LEVEL=3 CAPOP=1 SAT=1 VTO=-2.5
+ BETA=2.8E-3 LAMBDA=2.2M RS=70 RD=70 IS=1.7E-14
+ CGS=14P CGD=5P UCRIT=1.5 ALPHA=2
Example 2
The following example selects an n-channel JFET:
J2 7 1 4 JM1
.MODEL JM1 NJF (VTO=-1.5, BETA=5E-3, CGS=5P, CGD=1P,
+ CAPOP=1 ALPHA=2)
Example 3
The following example selects a p-channel JFET:
J3 8 3 5 JX
.MODEL JX PJF (VTO=-1.2, BETA=.179M, LAMBDA=2.2M
+ CGS=100P CGD=20P CAPOP=1 ALPHA=2)
Model Applications
Use MESFETs to model GaAs transistors for high speed applications. Using
MESFET models, transimpedance amplifiers for fiber optic transmitters up to
50 GHz can be designed and simulated.
scaling SCALM
Convergence
Enhance convergence for JFET and MESFET by using the GEAR method of
computation (.OPTION METHOD = GEAR) when you include the transit time
model parameter. Use the GMIN, GMINDC, and GRAMP options to increase the
parasitic conductance value in parallel with pn junctions of the device.
Capacitor Equations
The DCAP option selects the equation used to calculate the gate-to-source and
gate-to-drain capacitance for CAPOP=0. DCAP can be set to 1, 2 or 3. Default is
2.
Scaling
The AREA and M Element parameters, together with the SCALE and SCALM
control options, control scaling. For all three model levels, the model
parameters IS, CGD, CGS, RD, RS, BETA, LDEL, and WDEL, are scaled using the
same equations.
The SCALM option affects A, L, W, LDEL, and WDEL scaled parameters. SCALM
defaults to 1.0. For example, to enter the W parameter with micron units, set
SCALM to 1e-6 then enter W=5. The default setting is W=5e-6 meters or 5
microns.
To override global scaling that uses the .OPTION SCALM = <val> statement
in a JFET or MESFET model, include SCALM=<val> in the .MODEL statement.
nd
(drain node)
I1 (Jxxx)
ng
(gate node) nb
I2 (Jxxx) (bulk node)
ns
(source node)
I3 (Jxxx)
Transconductance
∂( ids )
gm = ----------------
∂( vgs ) vds = const.
Output Conductance
∂( ids )
gds = ----------------
∂( vds ) vgs = const.
The ids equation accounts for all DC currents of the JFET. Gate capacitances
are assumed to account for transient currents of the JFET equations. The two
diodes shown in Figure 17 are modeled by these ideal diode equations:
vgd
---------------
igd = ISeff ⋅ ⎛ e N ⋅ vt – 1⎞
⎝ ⎠
vgd > – 10 ⋅ N ⋅ vt
igd = – ISeff
vgd ≤– 10 ⋅ N ⋅ vt
vgs
---------------
igs = ISeff ⋅ ⎛ e N ⋅ vt – 1⎞
⎝ ⎠
vgs > – 10 ⋅ N ⋅ vt
igs = – ISeff
vgs ≤– 10 ⋅ N ⋅ vt
Gate
+ +
igs vgs cgs vgd igd
Source rs - cgd - rd Drain
ids
Note: For DC analysis, the capacitances are not part of the model.
Gate
gm(vgs vbs)
gds
Gate
gm(vgs vbs)
ind
Variable/
Quantity Definitions
Variable/
Quantity Definitions
rd Drain resistance
rs Source resistance
f Frequency
t Temperature in ° K
Dt t - tnom
Variable/
Quantity Definitions
vt(t) k ⋅ t/q
vt(tnom) k ⋅ tnom/q
Quantities Definitions
gm transconductance
Syntax
.MODEL <mname> NJF <LEVEL=val> <pname1=val1> ...
.MODEL mname PJF <LEVEL=val> <pname1=val1> ...
Parameter Description
Geometric ACM, ALIGN, AREA, HDIF, L, LDEL, LDIF, RD, RG, RS, RSH,
RSHG, RSHL, W, WDEL
Subthreshold ND, NG
Noise AF, KF
Name
(Alias) Units Default Description
ACM=0 AREA=Weff/Leff
ACM=1 AREA=Weff ⋅ Leff
AREAeff=M ⋅ AREA
ISeff = IS ⋅ AREAeff
Name
(Alias) Units Default Description
Name
(Alias) Units Default Description
CAPDSeff=CAPDS ⋅ Weff
------------ ⋅ M
Leff
Name
(Alias) Units Default Description
Name
(Alias) Units Default Description
Name
(Alias) Units Default Description
Name
(Alias) Units Default Description
Name Description
(Alias) Units Default
Name Description
(Alias) Units Default
ACM=1
ASPEC model, parameters function of element width.
AREAeff = Weff ⋅ Leff ⋅ m
RD
RDeff = --------
m
Or, if RD=0:
HDIF LDIF + ALIGN
RDeff = RSH ⋅ ----------------------- + RSHL ⋅ ---------------------------------------
Weff ⋅ m Weff ⋅ m
RG
RGeff = --------
m
Or, if RG=0:
Weff
RGeff = RSHG ⋅ ----------------------
Leff ⋅ m
RS
RSeff = -------
m
Or, if RS=0:
HDIF LDIF – ALIGN
RSeff = RSH ⋅ ----------------------- + RSHL ⋅ ---------------------------------------
Weff ⋅ m Weff ⋅ m
ACM=2
HSPICE model, combination of ACM=0,1 and provisions for lightly doped drain
technology.
ACM=3
Extends ACM=2 model to deal with stacked devices (shared source/drains) and
source/drain periphery capacitance along a gate edge.
Resulting calculations:
ISeff = IS ⋅ AREAeff
CGSeff = CGS ⋅ AREAeff
CGDeff = CGD ⋅ AREAeff
Weff
BETAeff = BETA ⋅ ------------ ⋅ m
Leff
Note: The model parameter units for IS, CGS, and CGD are unitless in
ACM=0 and per square meter for ACM=1.
Example
j1 10 20 0 40 nj_acm0 w=10u l=1u
j2a 10 20 0 41 nj_acm1 w=10u l=1u
vgd –M
cgd = CGDeff ⋅ ⎛ 1 – ---------⎞
⎝ PB ⎠
vgs < FC ⋅ PB
vgs – M
cgs = CGSeff ⋅ ⎛ 1 – --------⎞
⎝ PB ⎠
Forward Bias:
vgd >= FC ⋅ PB
vgd
1 – FC ⇒1 ( + M ) + M ⋅ ---------
∂ igd PB
cgd = TT ⋅ ------------ + CGDeff ⋅ ----------------------------------------------------------------------
∂vgd ( 1 – FC ) M + 1
vgs >= FC ⋅ PB
vgs
1 – FC ⇒1 ( + M ) + M ⋅ --------
∂igs PB
cgs = TT ⋅ ----------- + CGSeff ⋅ ---------------------------------------------------------------------
-
∂vgs ( 1 – FC ) M + 1
DCAP=2 (Default)
Reverse Bias:
vgd < 0
vgd –M
cgd = CGDeff ⋅ ⎛⎝ 1 – ---------⎞⎠
PB
vgs < 0
vgs – M
cgs = CGSeff ⋅ ⎛ 1 – --------⎞
⎝ PB ⎠
Forward Bias:
vgd ≥ 0
∂igd vgd
cgd = TT ⋅ ------------ + CGDeff ⋅ ⎛⎝ 1 + M ⋅ ---------⎞⎠
∂vgd PB
vgs ≥ 0
∂igs vgs
cgs = TT ⋅ ----------- + CGSeff ⋅ ⎛ 1 + M ⋅ --------⎞
∂vgs ⎝ PB ⎠
DCAP=3
Limits peak depletion capacitance to FC ⋅ CGDeff or FC ⋅ CGSeff with proper
fall-off when forward bias exceeds PB (FC > 1).
⎛ ⎞
CGD ⎜ vds ⎟
------------- ⋅ ⎜ 1 – --------------------------------------------------⎟
2 ⎜ 2⎟
vds + ⎛ --------------------⎞ ⎠
2 1
⎝ ⎝ ALPHA⎠
⎛ ⎞
⎜ vds ⎟
CGS veff – vte
Cgd = ⎜ ------------------------------ ⋅ 1 + -------------------------------------------------------- ⋅ 1 – -------------------------------------------------- ⎟+
⎜ vnew 2 2 2 ⎟
vds + ⎛ --------------------⎞
⎜ 4 1 – ------------- ( veff – vte ) + ( 0.2 ) 2 1 ⎟
⎝ PB ⎝ ALPHA⎠ ⎠
⎛ ⎞
⎜ ⎟
⎜ CGD
------------- ⋅ vds
1 + -------------------------------------------------- ⎟
⎜ 2 2 ⎟
vds + ⎛ --------------------⎞ ⎟
⎜ 2 1
⎝ ⎝ ALPHA⎠ ⎠
1 2 2
vnew = --- [ veff + vte + ( veff – vte ) + ( 0.2 ) ]
2
CGD = High -vds Cgd at vgs = 0
CGS = High -vds Cgs at vgs = 0
CGD - CGDeff
CGS - CGSeff
Figure 20 CAPOP=1 vs. CAPOP=2. Cgs, Cgd vs. Vgs for Vds=0, 1, 2, 3, 4
In Figure 21 on page 135, the Cgs-Cgd characteristic curve “flips over” below
the threshold for CAPOP=1, whereas for CAPOP=2, it is does not.
Figure 21 CAPOP=1 vs. CAPOP=2. Cgs, Cgd vs. Vds for Vgs=-1.5, -1.0, -0.5, 0
DC Model LEVEL=1
JFET DC characteristics are represented by the nonlinear current source
(ids). The value of ids is determined by the equations:
vgst = vgs – VTO
vgst< 0 Channel pinched off
ids = 0
0<vgst<vds Saturated region
ids = BETAeff ⋅ vgst 2 ⋅ ( 1 + LAMBDA ⋅ vds )
0<vds<vgst Linear region
LAMBDA = ⎛⎝ -----------------------------------------------------------------⎞⎠
ids2 – ids1
ids1 ⋅ vds2 – ids2 Þvds1
DC Model LEVEL=2
The DC characteristics of the JFET LEVEL=2 model are represented by the
nonlinear current source (ids). The value of ids is determined by the
equations:
vgst = vgs – VTO
vgst<0 Channel pinched off
ids = 0
0<vgst<vds, vgs<0 Saturated region, forward bias
ids = BETAeff ⋅ vgst 2 ⋅ [ 1 + LAMBDA ⋅ ( vds – vgst ) ⋅ ( 1 + LAM1 ⋅ vgs ) ]
0<vgst<vds, vgs<0 Saturated region, reverse bias
vgst
ids = BETAeff ⋅ vgst 2 ⋅ 1 – LAMBDA ⇒vds
( – vgst ) ⇒------------
VTO
0<vds<vgst Linear region
ids = BETAeff ⋅ vds ( 2 ⋅ vgst – vds )
DC Model LEVEL=3
The DC characteristics of the MESFET LEVEL=3 model are represented by the
nonlinear hyperbolic tangent current source (ids). The value of ids is
determined by the equations:
vds>0 Forward region
If model parameters VP and VTO are not specified they are calculated as:
q ⋅ NCHAN ⋅ Aeff 2
VP = – ----------------------------------------------------
2 ⋅ D ⋅ εo
VTO = VP + VBI
then:
vgst = vgs – [ VTO + GAMDS ⋅ vds + K1 ( vbs ) ]
BETAeff
beteff = ---------------------------------------------------
( 1 + UCRIT ⋅ vgst )
vgst<0 Channel pinched off
ids = idsubthreshold ( N0, ND, vds, vgs )
vgst>0, SAT=0 On region
ids = beteff ⋅ ( vgst VGEXP ) ⋅ ( 1 + LAMBDA ⋅ vds ) ⋅ tanh ( ALPHA ⋅ vds )
idsubthreshold ( N0, ND, vds, vgs )
vgst>0, SAT=1 On region
vds 3
ids = beteff ⋅ vgst 2 ⋅ ( 1 + LAMBDA ⋅ vds ) ⋅ 1 – ⎛ 1 – ALPHA ⇒--------⎞
⎝ 3 ⎠
Noise Equations
Figure 19 on page 117 shows the JFET noise model. Thermal noise generation
in the drain and source regions (RD and RS resistances) is modeled by the two
current sources, inrd and inrs. The inrd and inrs units are:
4 ⋅ k ⋅ t 1⁄ 2
inrd = ⎛ ---------------------⎞
⎝ rd ⎠
4⋅ k ⋅ t 1⁄ 2
inrs = ⎛ ---------------------⎞
⎝ rs ⎠
Channel thermal and flicker noise are modeled by the current source ind and
defined by the equation:
ind = channelthermalnoise + flickernoise
If the model parameter NLEV is less than 3, then:
8 ⋅ k ⋅ t ⋅ gm 1 ⁄ 2
channelthermalnoise = ⎛⎝ ------------------------------------⎞⎠
3
The previous formula, used in both saturation and linear regions, might lead to
wrong results in the linear region. For example, at VDS=0, channel thermal
noise is 0, because gm=0. This is physically impossible. If you set the NLEV
parameter to 3, simulation uses an equation that is valid in both linear and
saturation regions.[4]
For NLEV = 3
2
channelthermalnoise = ⎛ -------- ⋅ BETAeff ⋅ ( vgs – VTO ) ⋅ ------------------------ ⋅ GDSNOI⎞
8kt 1+a+a
⎝ 3 a ⎠
The following equations calculate values for the preceding equation:
vds
α = 1 – --------------------------
vgs – VTO
Linear region
α = 0
Saturation region
The flicker noise is calculated as:
AF 1 ⁄ 2
KF ⋅ ids
flickernoise = ⎛ -----------------------------⎞
⎝ f ⎠
Parameter Description
Function Parameter
DC M, TCV, XTI
grading M
mobility BEX
EG eV 1.16 Energy gap for the gate to drain and gate to source
diodes at 0 ° K
■
1.17 - silicon
■
0.69 - Schottky barrier diode
■
0.67 - germanium
■ 1.52 - gallium arsenide
GAP1 eV/× 7.02e-4 First bandgap correction factor, from Sze, alpha term
■
7.02e-4 - silicon
■
4.73e-4 - silicon
■ 4.56e-4 - germanium
■
5.41e-4 - gallium arsenide
GAP2 x 1108 Second bandgap correction factor, from Sze, beta term
■
1108 - silicon
■ 636 - silicon
■
210 - germanium
■ 204 - gallium arsenide
TCV (VTOTC) 1/× 0.0 Temperature compensation coefficient for VTO (threshold
voltage)
or
t2
eg ( t ) = 1.16 – 7.02e-4 ⋅ ------------------------
t + 1108.0
TLEV = 2
tnom 2
egnom = EG – GAP1 ⋅ -----------------------------------
tnom + GAP2
t2
eg ( t ) = EG – GAP1 ⋅ -----------------------
t + GAP2
TLEV = 0 or 1
PB ( t )
CGD ( t ) = CGD ⋅ 1 + M ⋅ ⎛ 4.0e-4 ⋅ Δt – -------------- + 1⎞
⎝ PB ⎠
eg ( t )
PB ( t ) = PB ⋅ ⎛⎝ -------------⎞⎠ – vt ( t ) ⇒ 3ln ⎛⎝ -------------⎞⎠ + ---------------------- – ------------
t t egnom
tnom tnom vt ( tnom ) vt ( t )
TLEVC = 1
CGS ( t ) = CGS ⋅ ( 1 + CTS ⋅ Δt )
CGD ( t ) = CGD ⋅ ( 1 + CTD ⋅ Δt )
PB M
CGD ( t ) = CGD ⋅ ⎛ --------------⎞
⎝ PB ( t )⎠
Δt
CGD ( t ) = CGD ⋅ ⎛ 1 – 0.5 ⇒dpbdt ⇒-------⎞
⎝ PB⎠
The next equation calculates values for the preceding equations:
PB ( t ) = PB + dpbdt ⋅ Δt
TLEV = 0 or 1
– egnom + 3 ⋅ vt ( tnom ) + ( 1.16 – egnom ) ⋅ ⎛⎝ 2 – -------------------------------⎞⎠ – PB
tnom
tnom + 1108
dpbdt = ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
tnom
TLEV = 2
– egnom + 3 ⋅ vt ( tnom ) + ( EG – egnom ) ⋅ ⎛ 2 – -----------------------------------⎞ – PB
tnom
⎝ tnom + GAP2⎠
dpbdt = -----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
tnom
t BEX
BETA ( t ) = BETA ⋅ ⎛⎝ -------------⎞⎠
tnom
If BETATCE=0; otherwise (TriQuint model):
BETATCE ( t – tnom )
BETA ( T ) = BETA ⋅ 1.01
Several alias TOM parameters are defined for existing LEVEL=3 parameters to
make the conversion easier. An alias allows the original name or the alias name
to be used in the .MODEL statement. However, the model parameter printout is
in the original name. Please note that in two cases, a sign reversal is needed,
even when using the TOM parameter name.
Q VGEXP
TRG1 TRG
TRD1 TRD
TRS1 TRS
Name
(Alias) Description
DELTA Ids feedback parameter of the TOM model. This parameter is not used
if its value is zero. DELTA may be negative or positive.
i ds
i ds ⇒--------------------------------------------------------------------------------------------------
max [ ( – 1 + v ntol ) ,( DELTA + v ds ⋅ i ds ) ]
CAPDS Weff
Drain-to-source capacitance: CAPDSeff = CAPDS ⋅ ------------ ⋅ M
Leff
to reproduce the TOM model. Using non-zero values for these parameters with
nonzero BETATCE, DELTA, or CAPDS, results in a hybrid model.
Model Description
This section describes the DC and Capacitance equations for the HSPICE
LEVEL=7 JFET model.
DC Equations
The DC equations for the HSPICE LEVEL=7 JFET model are:
Drain to Source Current (IDS)
I DS = I O ⋅ ( 1 + LAMBDA ⋅ V DS )
Q
IO = β ⋅ VG ⋅ fK
α ⋅ V DS
f K = ----------------------------------------------------
K 1,1/K
[ 1 + ( α ⋅ V DS ) ]
V G = Q ⋅ V ST ⋅ log [ exp ( u ) + 1 ]
V GS – V TO + γ ⋅ V DS
u = ----------------------------------------------------
Q ⋅ V ST
V ST = V STO ⋅ ( 1 + M STO ⋅ V DS )
Transconductance
Q–1
⎛ Q ⋅ β ⋅ fK ⋅ VG ⎞
GM = ⎜ ------------------------------------------------⎟ ⋅ ( 1 + LAMBDA ⋅ V DS )
⎝ 1 + exp ( – u ) ⎠
Output Conductance
( V GS – V TO + γ ⋅ V DS ) ⋅ M ST0
G DS = LAMBDA ⋅ I 0 + G M ⋅ γ – ⎛ -----------------------------------------------------------------------------⎞
⎝ 1 + M ST0 ⋅ V DS ⎠
Q
Q ⋅ I 0 ⋅ M ST0 ⎞ ⎛
⎛ -------------------------------------- α ⋅ β ⋅ VG ⎞
- + ⎜ -----------------------------------------------------
-⎟ ⋅ ( 1 + LAMBDA ⋅ V DS )
⎝ 1 + M ST0 ⋅ V DS⎠ ⎝ K 1+1/K⎠
[1 + (α ⋅ V ) ] DS
– V GS
I LS = ILK ⋅ ⎛⎝ 1 – exp ------------⎞⎠
PLK
– V GS
G LS = ⎛ -----------⎞ ⋅ ⎛ exp ------------⎞
ILK
⎝ PLK⎠ ⎝ PLK ⎠
– V GD
I LD = ILK ⋅ ⎛ 1 – exp -------------⎞
⎝ PLK ⎠
– V GD
G LD = ⎛ -----------⎞ ⋅ ⎛ exp -------------⎞
ILK
⎝ PLK⎠ ⎝ PLK ⎠
Temperature and Geometry Dependence
BETATCE ⋅ ( T – T NOM )
β = AREA ⋅ BETA ⋅ 1.01
ALPHATCE ⋅ ( T – T NOM )
α = ALPHA ⋅ 1.01
V TO = VTO + VTOC ⋅ ( T – T NOM )
Capacitance Equations
The capacitance equations for the HSPICE LEVEL=7 JFET model are:
Combined Gate Charge
Q GG = Q GL ⋅ f T + Q GH ⋅ ( 1 – f T ) + QGG0 ⋅ ( V GSI + V GDI )
∂f T
C GS = C GSL ⋅ f T + C GSH ⋅ ( 1 – f T ) + ( Q GL – Q GH ) ⋅ -------------- + QGG0
∂V GSI
∂f T
C GD = C GDL ⋅ f T + C GDH ⋅ ( 1 – f T ) + ( Q GL – Q GH ) ⋅ --------------- + QGG0
∂V GDI
f T = exp ( – Q GGB ⋅ I DS ⋅ V DS )
∂f T
--------------
- = – QGGB ⋅ [ I DS + ( g m + g ds ) ⋅ V DS ] ⋅ f T
∂V GDI
∂f T
------------
- = – QGGB ⋅ [ I DS + g ds ⋅ V DS ] ⋅ f T
∂V GD
Lower Power Gate Charge
Q GL = qgI + QGCL ⋅ ( V GS + V GD )
Name (Alias)
Units Default Description
Name (Alias)
Units Default Description
K - 2 Knee-function parameter
Name (Alias)
Units Default Description
DC Model
V GS 2 α1*V DS
I D = I DSS ⎛ 1 – ---------⎞ tanh ⎛ ----------------------⎞
⎝ VP ⎠ ⎝ V GS – V P⎠
∂I D V GS α1 ⋅ V DS
g m = ------------ = I DSS – ------ ⎛⎝ 1 – ---------⎞⎠ ⋅ tanh ⎛⎝ -----------------------⎞⎠
2
∂V GS VP VP V GS – V P
V GS 2⎞
⎛ 1 – -------- α1 ⋅ V DS⎞ – α1 ⋅ V DS
- ⋅ sec h 2 ⎛ ----------------------
- --------------------------
⎝ VP ⎠ ⎝ V GS – V P ⎠ ⋅ 2
V GS – V P
∂I D V GS 2γ V GS α1 ⋅ V DS
g DS = ------------ = I DSS ⋅ ⎛ 1 – ---------⎞ ⋅ – ---------------- ⋅ tanh ⎛ -----------------------⎞
∂V DS ⎝ VP ⎠ 2 ⎝ V GS – V P ⎠
VP
V GS⎞
⎛ 1 – -------- α1 ⋅ V DS⎞ α1 ⋅ ( V GS – V PO )
- ⋅ sec h 2 ⎛ ----------------------
- ⋅ --------------------------------------------
⎝ VP ⎠ ⎝ V GS – V P ⎠ 2
( V GS – V P )
V P = V TO + γ V DS
Table 41 DC Model Parameters
Name
(Alias) Units Default Description
Name
(Alias) Units Default Description
⎛ ⎞
CGD ⎜ vds ⎟
------------- ⎜ 1 – ----------------------------------------------------
-⎟
2 ⎜ 2⎟
vds + ⎛ -----------------------⎞ ⎠
2 1
⎝ ⎝ ALPHA1⎠
⎛ ⎞
CGD ⎜ vds ⎟
------------- ⎜ 1 + -----------------------------------------------------⎟
2 ⎜ 2⎟
vds + ⎛ -----------------------⎞ ⎠
2 1
⎝ ⎝ ALPHA1⎠
1 2 2
vnew = --- [ veff + vte + ( veff – vte ) + ( 0.2 ) ]
2
Table 42 Gate Capacitance Model Parameters
Name
(Alias) Units Default Description
Noise Model
Two current sources model the thermal noise generation in the drain and
source regions (RD and RS resistances):
■
inrd
■
inrs
inrd and inrs are modeled by:
1⁄ 2
inrs = ⎛⎝ --------⎞⎠
4kt
rs
4kt 1 ⁄ 2
inrd = ⎛ --------⎞
⎝ rd ⎠
Channel thermal and flicker noise are modeled by the ind current source, and
defined by the equation:
ind = channel thermal noise+ flicker noise
8kt ⋅ g m 1 ⁄ 2
channel thermal noise = ⎛ ---------------------⎞
⎝ 3 ⎠
AF 1 ⁄ 2
KF ⋅ ids
flicker noise = ⎛ -----------------------------⎞
⎝ f ⎠
Table 43 Noise Model Parameters
Name
(Alias) Units Default Description
Example
.MODEL NCH NJF LEVEL =8
+ IDSS = 69.8e-3 VTO = -2 GAMMA = 0
+ ALPHA1 = 1 RS = 0 RD = 0
+ CGS = 1e-15 CGD = 2e-16 PB = 0.8
+ IS = 5e-16 AF = 1 KF = 0
+ FC = 0.5
.END
References
[1] GaAs FET Device and Circuit Simulation in SPICE, IEEE Transactions on
Electron Devices, Vol. ED-34.
[2] A MESFET Model for Use in the Design of GaAs Integrated Circuits, IEEE
Transactions on Microwave Theory, Vol. MTT-28 No. 5.
[3] H. Statz, P.Newman, I.W.Smith, R.A. Pucel, and H.A. Haus, GaAs FET
Device and Circuit Simulation in Spice.
[4] Tsivids, Yanis P., Operation and Modeling of the MOS Transistor, McGraw-
Hill, 1987, p. 340.
[5] A.J. McCamant, G.D. Mc Cormack, and D.H.Smith, An Improved GaAs
MESFET Model for SPICE, IEEE.
[6] W.Curtice, “A MESFET Model for Use in the Design of GaAs Integrated
Circuits,” IEEE Tran, Microwave, and H.Statz, P.Newman, I.W.Smith, R.A.
Pucel, and H.A. Haus, “GaAs FET Device And Circuit Simulation in SPICE”.
■
LEVEL=6 Philips Bipolar Model (MEXTRAM Level 503)
■
LEVEL=8 HICUM Model
■
HICUM Model Advantages
■
LEVEL=9 VBIC99 Model
■
LEVEL=9 Model Parameters
■
LEVEL=10 Philips MODELLA Bipolar Model
■
LEVEL=11 UCSD HBT Model
■
LEVEL=13 HICUM0 Model
Note: To modify high-injection effects for BiCMOS devices, use the IKF
and IKR high-current Beta degradation parameters.
The SUBS model parameter facilitates the modeling of both vertical and lateral
geometrics.
Selecting Models
To select a BJT device, use a BJT element and model statement. The element
statement uses the name of the simulation device model to reference the
model statement. The following example uses the reference name MOD1. This
example uses an NPN model type to describe an NPN transistor.
Q3 3 2 5 MOD1 <parameters>
.MODEL MOD1 NPN <parameters>
You can specify parameters in both element and model statements. If you
specify the same parameter in both an element and a model, then the element
parameter (local to the specific instance of the model) always overrides the
model parameter (global default for all instances of the model, if you did not
define the parameter locally). The model statement specifies the type of
device—for example, for a BJT, the device type might be NPN or PNP.
You can override global depletion capacitance equation selection that uses
the .OPTION DCAP=<val> statement in a BJT model by including
DCAP=<val> in the BJTs .MODEL statement.
Convergence
Adding a base, collector, and emitter resistance to the BJT model improves its
convergence. The resistors limit the current in the device so that the forward-
biased pn junctions are not overdriven.
Parameter Description
Example
.MODEL t2n2222a NPN
+ ISS=0. XTF= 1. NS = 1.00000
+ CJS=0. VJS= 0.50000 PTF= 0.
+ MJS=0. EG = 1.10000 AF = 1.
+ ITF=0.50000 VTF= 1.00000
+ BR =40.00000 IS = 1.6339e-14 VAF=103.40529
+ VAR=17.77498 IKF= 1.00000
+ NE =1.31919 IKR= 1.00000 ISC= 3.6856e-13
+ NC =1.10024 IRB= 4.3646e-05 NF = 1.00531
+ NR =1.00688 RBM= 1.0000e-02 RB =71.82988
+ RC =0.42753 RE = 3.0503e-03 MJE= 0.32339
+ MJC=0.34700 VJE= 0.67373 VJC= 0.47372
+ TF =9.693e-10 TR =380.00e-9 CJE= 2.6734e-11
+ CJC=1.4040e-11 FC = 0.95000 XCJC=0.94518
Parameters
Transient model parameters for BJTs are composed of two groups: junction
capacitor parameters and transit time parameters.
■
CJE, VJE, and MJE model the base-emitter junction.
■
CJC, VJC, and MJC model base-collector junction capacitance.
■
CJS, VJS, and MJS model the collector-substrate junction capacitance.
TF is the forward transit time for base charge storage. TF can be modified to
account for bias, current, and phase, by XTF, VTF, ITF, and PTF. The base
charge storage reverse transit time is set by TR. To select from several sets of
temperature equations for the BJT model parameters, set TLEV and TLEVC.
Parameter Description
DC BF, BR, IBC, IBE, IS, ISS, NF, NR, NS, VAF, VAR
Parameter Description
noise KF, AF
Name (Alias)
Unit Default Description
BULK (NSUB) 0.0 Sets the bulk node to a global node name. A
substrate terminal node name (ns) in the
element statement overrides BULK.
Name (Alias)
Unit Default Description
otherwise:
ISE (C2, JLE) amp 0.0 Base-emitter leakage saturation current. If ISE
is greater than 1e-4, then:
otherwise:
IKF (IK, JBF) amp 0.0 Corner for forward Beta high-current roll-off.
Use zero to indicate an infinite value.
IKR (JBR) amp 0.0 Corner for reverse Beta high-current roll-off.
Use zero to indicate an infinite value.
IRB (JRB, amp 0.0 Base current, where base resistance falls half-
IOB) way to RBM. Use zero to indicate an infinite
value.
RCeff =RC/(AREA ⋅ M)
Function Parameter
emission coefficients TNC1, TNC2, TNE1, TNE2, TNF1, TNF2, TNR1, TNR2,
TNS1, TNS2
Function Parameter
GAP1 eV/× 7.02e-4 First bandgap correction factor (from Sze, alpha
term):
■
7.02e-4 - silicon
■ 4.73e-4 - silicon
■
4.56e-4 - germanium
■ 5.41e-4 - gallium arsenide
Scaling
Scaling is controlled by the element parameters AREA, AREAB, AREAC, and M.
The AREA parameter, the normalized emitter area, divides all resistors and
multiplies all currents and capacitors. AREAB and AREAC scale the size of the
base area and collector area. Either AREAB or AREAC is used for scaling,
depending on whether vertical or lateral geometry is selected (using the SUBS
model parameter). For vertical geometry, AREAB is the scaling factor for IBC,
ISC, and CJC. For lateral geometry, AREAC is the scaling factor. The scaling
factor is AREA for all other parameters.
The following formula scales the DC model parameters (IBE, IS, ISE, IKF,
IKR, and IRB) for both vertical and lateral BJT transistors:
Ieff = AREA ⋅ M ⋅ I
In the preceding equation, I can be IBE, IS, ISE, IKF, IKR, or IRB.
For both the vertical and lateral, the resistor model parameters, RB, RBM, RE,
and RC are scaled by the following equation.
R
Reff = ----------------------------
AREA ⋅ M
In the preceding equation, R can be RB, RBM, RE, or RC.
nc
(collector node)
I1 (Q1)
nb
ns
(base node)
(substrate node)
I2 (Q1)
I4 (Q1)
ne
(emitter node)
I3 (Q1)
Collector Conductance
∂ic ∂ic
g o = ----------- = – -----------
∂vce vbe = const.
∂vbc vbe = const.
Transconductance
∂ic
gm = -----------
∂vbe vce = const.
∂ic ∂ic
gm = ----------- + -----------
∂vbe ∂vbc
∂ic
gm = ----------- – g o
∂vbe
collector
CBCP rc
rb
Base ice
cbe ibe
ibs cbs
substrate re
CCSP
CBEP
emitter
collector
CBCP rc CCSP
cbe ibe
re
CBEP
emitter
collector
CBCP rc
rb
Base gm ice go
vbe
cbe gbe
gbs cbs
substrate re
CCSP
CBEP
emitter
collector
CBCP rc CCSP
cbe gbe
re
CBEP
emitter
collector
CBCP rc inrc
rb
Base gm ice go inc
vbe
inrb
cbe gbe inb
gbs cbs
substrate re inre
CCSP CBEP
emitter
collector
CBCP CCSP
rc inrc
inrb
cbe gbe inb
re inre
CBEP
emitter
Variable Definitions
f Frequency
Variable Definitions
gm Transconductance
go Collector conductance
rb Base resistance
Variable Definitions
Quantities Definitions
Δt t - tnom
vt(t) k ⋅ t/q
vt(tmon) k ⋅ tnom/q
Quantities Definitions
ib base current
ic collector current
is substrate current
Quantities Definitions
vs substrate voltage
power power
gm transconductance
rx base resistance
ro collector resistance
emitter
base
collector
buried collector
substrate
substrate
In Figure 31, the views from the top demonstrate how IBE is multiplied by
either base area, AREAB or collector area, AREAC.
substrate
C vertical B lateral transistor
E E C
Area Area AreaC
AreaB
AreaB
DC Model Equations
DC model equations are for the DC component of the collector current (ic) and
the base current (ib).
Current Equations: IS Only
If you specify only IS, without IBE and IBC:
vbe vbc vbc vbc
-------------------- -------------------- -------------------- --------------------
ic = ------------ ⋅ ⎛ e NF ⋅ vt – e NR ⋅ vt⎞ – ------------ ⇒⎛ e NR ⋅ vt – 1⎞ – ISCeff ⇒⎛ e NC ⋅ vt – 1⎞
ISeff ISeff
qb ⎝ ⎠ BR ⎝ ⎠ ⎝ ⎠
vbe vbc
-------------------- --------------------
ib = ------------ ⋅ ⎛⎝ e NF ⋅ vt – 1⎞⎠ + ------------ ⇒⎛⎝ e NR ⋅ vt – 1⎞⎠
ISeff ISeff
BF BR
vbe vbc
-------------------- --------------------
+ ISEeff ⇒⎛ e NE ⋅ vt – 1⎞ + ISCeff ⋅ ⎛ e NC ⋅ vt – 1⎞
⎝ ⎠ ⎝ ⎠
Current Equations: IBE and IBC
If you specify IBE and IBC, instead of IS:
vbe vbc
-------------------- --------------------
ic = ----------------- ⋅ ⎛ e NF ⋅ vt – 1⎞ – ----------------- ⇒⎛ e NR ⋅ vt – 1⎞
IBEeff IBCeff
qb ⎝ ⎠ qb ⎝ ⎠
vbec vbe
-------------------- --------------------
– ----------------- ⇒⎛⎝ e NR ⋅ vt – 1⎞⎠ – ISCeff ⇒⎛⎝ e NC ⋅ vt – 1⎞⎠
IBCeff
BR
vbe vbc
-------------------- --------------------
ib = ----------------- ⋅ ⎛⎝ e NF ⋅ vt – 1⎞⎠ + ----------------- ⇒⎛⎝ e NR ⋅ vt – 1⎞⎠
IBEeff IBCeff
BF BR
vbe vbc
-------------------- --------------------
+ISEeff ⇒⎛ e NE ⋅ vt – 1⎞ + ISCeff ⋅ ⎛ e NC ⋅ vt – 1⎞
⎝ ⎠ ⎝ ⎠
vsc > – 10 ⋅ NS ⋅ vt
isc = – ISSeff
vsc ≤– 10 ⋅ NS ⋅ vt
Lateral Transistors
vbs
-------------------
ibs = ISSeff ⋅ ⎛ e NS ⋅ vt – 1⎞
⎝ ⎠
vbs > – 10 ⋅ NS ⋅ vt
ibs = – ISSeff
vbs ≤– 10 ⋅ NS ⋅ vt
If you do not specify either IBE or IBC:
ISSeff = ISS ⋅ AREA ⋅ M
If you specify both IBE and IBC:
ISSeff = ISS ⋅ AREAC ⋅ M
vertical
ISSeff = ISS ⋅ AREAB ⋅ M
lateral
1
q1 = --------------------------------------------
⎛ 1 – ----------
vbc- ----------
–
vbe-⎞
⎝ VAF VAR⎠
Otherwise, if UPDATE=1 and
vbe-
vbc- ----------
---------- + ≥0
VAF VAR
then
vbc vbe
q1 = 1 + ----------- + -----------
VAF VAR
vbe vbc
-------------------- --------------------
q2 = ----------------- ⋅ ⎛ e NF ⋅ vt – 1⎞ + ----------------- ⋅ ⎛ e NR ⋅ vt – 1⎞
ISeff ISeff
IKFeff ⎝ ⎠ IKReff ⎝ ⎠
With IBE and IBC, the preceding equation is:
vbe vbc
-------------------- --------------------
q2 = ----------------- ⋅ ⎛ e NF ⋅ vt – 1⎞ + ----------------- ⋅ ⎛ e NR ⋅ vt – 1⎞
IBEeff IBCeff
IKFeff ⎝ ⎠ IKReff ⎝ ⎠
In the preceding equation:
■ IBE=IS if IBE=0
■
IBC=IS if IBC=0
q1
qb = ------ ⋅ [ 1 + ( 1 + 4 ⋅ q2 ) NKF ]
2
tan ( z ) – z
rbb = RBMeff + 3 ⋅ ( RBeff – RBMeff ) ⋅ -----------------------------------------------
z ⋅ tan ( z ) ⋅ tan ( z )
– 1 + [ 1 + 144 ⋅ ib / ( π2 ⋅ IRBeff ) ] 1 / 2-
z = ----------------------------------------------------------------------------------------------
24 ⎛ ib ⎞ 1 / 2
------ ⋅ -----------------
π2 ⎝ IRBeff⎠
The following equation calculates the argtf value for the preceding equation:
vbc
2 -----------------------------
argtf = XTF ⋅ ⎛ ------------------------⎞ ⋅ e 1.44 ⋅ VTF
ibe
⎝ ibe + ITF⎠
vbe – MJE
cbedep = CJEeff ⋅ ⎛ 1 – ----------⎞
⎝ VJE⎠
vbe ≥ FC ⋅ VJE
vbe
1 – FC Þ ( 1 + MJE ) + MJE ⋅ ----------
VJE
cbedep = CJEeff ⋅ ------------------------------------------------------------------------------------
-
( 1 – FC ) ( 1 + MJE )
DCAP=2
The base-emitter depletion capacitance is determined as follows:
vbe < 0
vbe – MJE
cbedep = CJEeff ⋅ ⎛ 1 – ----------⎞
⎝ VJE⎠
vbe ≥ 0
DCAP=3
Limits peak depletion capacitance to FC ⋅ CJCeff or FC ⋅ CJEeff with proper
fall-off when forward bias exceeds PB (FC ≥ 1).
vbc – MJC
cbcdep = XCJC ⋅ CJCeff ⋅ ⎛⎝ 1 – -----------⎞⎠
VJC
vbc ≥ FC ⋅ VJC
vbc
1 – FC ⇒1 ( + MJC ) + MJC ⋅ -----------
VJC
cbcdep = XCJC ⋅ CJCeff ⋅ --------------------------------------------------------------------------------------
-
( 1 – FC ) ( 1 + MJC )
DCAP=2
Specify DCAP=2 to select one of the following equations:
vbc < 0
vbc – MJC
cbcdep = XCJC ⋅ CJCeff ⋅ ⎛ 1 – -----------⎞
⎝ VJC⎠
vbc ≥ 0
vbcx –MJC
cbcx = CJCeff ⋅ ( 1 – XCJC ) ⋅ ⎛ 1 – ------------⎞
⎝ VJC ⎠
vbcx ≥ FC ⋅ VJC
vbcx
1 – FC Þ ( 1 + MJC ) + MJC ⋅ ------------
VJC-
cbcx = CJCeff ⋅ ( 1 – XCJC ) ⋅ --------------------------------------------------------------------------------------
( )
( 1 – FC ) 1 + MJC
DCAP=2
Specify DCAP=2 to select one of the following equations:
vbcx < 0
– MJC
cbcx = CJCeff ⋅ ( 1 – XCJC ) ⋅ ⎛⎝ 1 – ------------⎞⎠
vbcx
VJC
vbcx ≥ 0
Substrate Capacitance
The function of substrate capacitance is similar to that of the substrate diode.
To switch it from the collector to the base, set the SUBS model parameter.
vbs –MJS
cbs = CJSeff ⋅ ⎛⎝ 1 – ----------⎞⎠
VJS
Forward Bias vbs ≥ 0
vsc – MJS
csc = CJSeff ⋅ ⎛ 1 – ----------⎞
⎝ VJS⎠
Forward Bias vsc ≥ 0
4 ⋅ k ⋅ t 1/2
inrb = ⎛ ---------------------⎞
⎝ rbb ⎠
The mean square short-circuit collector resistance noise current equation is:
4 ⋅ k ⋅ t 1/2
inrc = ⎛ ---------------------⎞
⎝ RCeff ⎠
The mean square short-circuit emitter resistance noise current equation is:
4 ⋅ k ⋅ t 1/2
inre = ⎛ ---------------------⎞
⎝ REeff ⎠
The noise associated with the base current is composed of two parts: shot
noise and flicker noise. Typical values for the flicker noise coefficient, KF, are
1e-17 to 1e-12. They are calculated as:
2 ⋅ q ⋅ fknee
In the preceding equation, fknee is the noise knee frequency (typically 100 Hz
to 10 MHz), and q is electron charge.
KF ⋅ ib AF
inb = ( 2 ⋅ q ⋅ ib ) + ⎛ --------------------------⎞
2
⎝ f ⎠
2 2 2
inb = shot noise + flicker noise
shot noise = ( 2 ⋅ q ⋅ ib ) 1 / 2
KF ⋅ ib AF 1 / 2
flicker noise = ⎛⎝ --------------------------⎞⎠
f
The noise associated with the collector current is modeled as shot noise only.
inc = ( 2 ⋅ q ⋅ ic ) 1 / 2
Parameter Description
t2
eg ( t ) = 1.16 – 7.02e-4 ⋅ ------------------------
t + 1108.0
TLEV=2
tnom 2
egnom = EG – GAP1 ⋅ -----------------------------------
tnom + GAP2
t2
eg ( t ) = EG – GAP1 ⇒-----------------------
t + GAP2
t XTB
BF ( t ) = BF ⋅ ⎛ -------------⎞
⎝ tnom⎠
t XTB
BR ( t ) = BR ⋅ ⎛ -------------⎞
⎝ tnom⎠
facln
ISE - -------------
ISE ( t ) = --------------------------
XTB
⋅ e NE
⎛ ------------
t -⎞
⎝ tnom⎠
facln
ISC - -------------
ISC ( t ) = --------------------------
XTB
⋅ e NC
⎛ -------------⎞
t
⎝ tnom⎠
facln
ISS - -------------
ISS ( t ) = --------------------------
XTB
⋅ e NS
⎛ -------------⎞
t
⎝ tnom⎠
TLEV=0, 1 or 3
IS ( t ) = IS ⋅ e facln
facln
-------------
IBE ( t ) = IBE ⋅ e NF
facln
-------------
IBC ( t ) = IBC ⋅ e NR
The following equation calculates the facIn value for the preceding equations:
BR ( t ) = BR ⋅ ( 1 + TBR1 ⋅ Δt + TBR2 ⋅ Δt 2 )
TF ( t ) = TF ⋅ ( 1 + TTF1 ⋅ Δt + TTF2 ⋅ Δt 2 )
TR ( t ) = TR ⋅ ( 1 + TTR1 ⋅ Δt + TTR2 ⋅ Δt 2 )
NF ( t ) = NF ⋅ ( 1 + TNF1 ⋅ Δt + TNF2 ⋅ Δt 2 )
NR ( t ) = NR ⋅ ( 1 + TNR1 ⋅ Δt + TNR2 ⋅ Δt 2 )
NE ( t ) = NE ⋅ ( 1 + TNE1 ⋅ Δt + TNE2 ⋅ Δt 2 )
NC ( t ) = NC ⋅ ( 1 + TNC1 ⋅ Δt + TNC2 ⋅ Δt 2 )
NS ( t ) = NS ⋅ ( 1 + TNS1 ⋅ Δt + TNS2 ⋅ Δt 2 )
VJC ( t )
CJC ( t ) = CJC ⋅ 1 + MJC ⋅ ⎛ 4.0e-4 ⋅ Δt – ----------------- + 1⎞
⎝ VJC ⎠
VJS ( t )
CJS ( t ) = CJS ⋅ 1 + MJS ⋅ ⎛ 4.0e-4 ⋅ Δt – ---------------- + 1⎞
⎝ VJS ⎠
eg ( t )
VJE ( t ) = VJE ⋅ ------------- – vt ( t ) ⇒ 3 ⋅ ln ⎛ -------------⎞ + ---------------------- – ------------
t t egnom
tnom ⎝ tnom ⎠ vt ( tnom ) vt ( t )
eg ( t )
VJC ( t ) = VJC ⋅ ------------- – vt ( t ) ⇒ 3 ⋅ ln ⎛ -------------⎞ + ---------------------- – ------------
t t egnom
tnom ⎝ tnom⎠ vt ( tnom ) vt ( t )
egnom eg ( t )
VJS ( t ) = VJS ⋅ ------------- – vt ( t ) ⇒ 3 ⋅ ln ⎛⎝ -------------⎞⎠ + ---------------------- – ------------
t t
tnom tnom vt ( tnom ) vt ( t )
TLEVC=1
CJE ( t ) = CJE ⋅ ( 1 + CTE ⋅ Δt )
VJC MJC
CJC ( t ) = CJC ⋅ ⎛ -----------------⎞
⎝ VJC ( t )⎠
VJS MJS
CJS ( t ) = CJS ⋅ ⎛ ----------------⎞
⎝ VJS ( t )⎠
Δt
CJC ( t ) = CJC ⋅ ⎛ 1 – 0.5 ⇒dvjcdt ⇒-----------⎞
⎝ VJC⎠
Δt
CJS ( t ) = CJS ⋅ ⎛⎝ 1 – 0.5 ⇒dvjsdt ⇒----------⎞⎠
VJS
If TLEV=0, 1, or 3, then:
RB ( t ) = RB ⋅ ( 1 + TRB1 ⋅ Δt + TRB2 ⋅ Δt 2 )
RC ( t ) = RC ⋅ ( 1 + TRC1 ⋅ Δt + TRC2 ⋅ Δt 2 )
t BEX
RC ( t ) = RC ⋅ ⎛ -------------⎞
⎝ tnom⎠
t BEXV
VO ( t ) = VO ⋅ ⎛⎝ -------------⎞⎠
tnom
S C
Cx
+ -
Vbcx
lepl
Ci
csc D + -
Ibs=BRS(ibc-isc) Vbc
isc
RB
B npn
RE
Cx
+ -
Vbcx
lepi
Ci
+ -
Vbc
RB
B npn
cbs ibs
RE
S E
ki – kx – ln ⎛ ---------------⎞ + ---------------------------
1 + ki vbc – vbcx
⎝ 1 + kx⎠ NEPI ⋅ vt
Iepi = ------------------------------------------------------------------------------------
⎛ -------------------------
RCeff -⎞ ⎛ vbc – vbcx-⎞
-----------------------------
⎝ NEPI ⋅ vt⎠ ⋅ ⎝ 1 + VO ⎠
If you set the GAMMA model parameter to zero, then the ki and kx values both
become one, and:
vbc – vbcx
Iepi = -------------------------------------------------------------------
RCeff ⋅ ⎛⎝ 1 + ------------------------------⎞⎠
vbc – vbcx
VO
qi = QCOeff ⋅ ⎛ ki – 1 – -----------------------⎞
GAMMA
⎝ 2 ⎠
qx = QCOeff ⋅ ⎛ kx – 1 – -----------------------⎞
GAMMA
⎝ 2 ⎠
∂ GAMMA ⋅ QCOeff
Ci = ( qi ) = ⎛⎝ --------------------------------------------------⎞⎠ ⋅ e vbc / ( NEPI ⋅ vt )
∂vbc 2 ⋅ NEPI ⋅ vt ⋅ kx
∂ ( qx ) = ⎛ GAMMA ⋅ QCOeff⎞
Cx = -------------------------------------------------- ⋅ e vbcx / ( NEPI ⋅ vt )
∂vbcx ⎝ 2 ⋅ NEPI ⋅ vt ⋅ kx ⎠
If GAMMA=0, then the Ci and Cx values become zero.
Example
This example is based on demonstration netlist quasisat.sp, which is available
in directory $<installdir>/demo/hspice/bjt:
Note: All parameter values in the following model must come from
either a .PARAM statement or the subcircuit call.
Example
The following is a subcircuit definition that converts the National Semiconductor
scaled BJT model to a form usable in HSPICE.
History of VBIC
VBIC was developed by engineers at several companies. The detailed
equations for all elements are given in the referenced publication. Recent
information and source code can be found on the web site:
http://www-sm.rz.fht-esslingen.de/institute/iafgp/neu/VBIC/index.html
The HSPICE implementation complies with standard VBIC. Starting in release
2001.4 of the VBIC model, self-heating and excess phase have been
implemented or enabled.
The large signal equivalent circuit for VBIC is shown in Figure 35. Capacitors
CBCO, CBEO and resistors RCX, RBX, RE, and RS are linear elements, all other
elements of the equivalent circuit are nonlinear.
VBIC Parameters
Figure 35 lists the parameters that you can set for the model, and shows the
default values for each parameter. The same parameter names are used in the
table and the previous referenced publication.
Starting in Version 2003.03, the BJT LEVEL=4 model prints FT in the .OP
output.
CBCO Collector
Substrate
RS RCX
I Qbcp
t bcp RBIP/qBP
Iccp t
Qbep RCi
Ibep
CBEO Re
If values of parameters given by the user are beyond their ranges, those
parameters will be reset to new values and warnings will be printed unless you
set .OPTION NOWARN.
TREF (TNOM) o
C 27 Nominal measurement temperature of
parameters (do not use TNOM alias)
Noise Analysis
The following sources of noise are taken into account:
■
The thermal noise of resistors RBX, RCX, RE, RS, RBP, RCI, RBI
■
Shot noise of currents IBE, IBEP, ICC, ICCP
■ Flicker noise due to currents IBE, IBEP
Noise due to IBEX and IGC is not included.
■
If you do not specify either the DTEMP or the DTA parameter, then
DTEMP=0.0.
■ If you specify DTA but not DTEMP, then DTEMP uses the DTA value.
■
If you specify DTEMP, then simulation uses the DTEMP value, and ignores the
DTA value.
For a description of the MEXTRAM model, refer to:
http://www-us.semiconductors.com/Philips_Models/
Parameter Description
Parameter Description
OUTFLAG
RBV ohm 100. Variable part of the base resistance at zero bias
Example
This example is located in the following directory:
$installdir/demo/hspice/bjt/mextram.sp
■
Split base-collector and base-emitter depletion capacitance
■
Current crowding and conductivity modulation of the base resistance
■
First order approximation of distributed high frequency effects in the intrinsic
base (high frequency current crowding and excess phase shift)
■
Ohmic resistance of epilayer
■
Velocity saturation effects on the resistance of the epilayer
■
Recombination in the base (meant for SiGe transistors)
■
Early effects in the case of a graded bandgap (SiGe)
■
Thermal noise, shot noise, and 1/f-noise
■
Self-heating
You can use either of two parameters to specify the difference between the
circuit temperature and the ambient temperatures in the MEXTRAM model:
■
DTEMP instance parameter as specified in the element statement
■
DTA (difference between the device temperature and the ambient analysis
temperature) global model parameter.
DTA and DTEMP both default to zero. DTEMP overrides DTA locally, if you specify
both. Simulation uses the value of DTEMP to de-rate the temperature in model
equations and other parameters.
■
If you do not specify either the DTEMP or the DTA parameter, then
DTEMP=0.0.
■
If you specify DTA but not DTEMP, then DTEMP uses the DTA value.
■
If you specify DTEMP, then simulation uses the DTEMP value, and ignores the
DTA value.
This model is described at:
http://www.semiconductors.philips.com/Philips_Models/newsflashmextram504
The following topics are covered:
■
Notes on Using MEXTRAM 503 or 504 Devices
■
LEVEL=6 Model Parameters (504)
■
Mextram 504 Update Releases
TREF o
C 25.0 Reference temperature
XEXT - 0.63 Part of Iex, Qex, Qtex, and Isub that depends on
the base-collector voltage Vbc1
RTH o
C/W 0 Thermal (self-heating) resistance
■
Accurate description of the high-current operating region (including quasi-
saturation and saturation).
■ Distributed modelling of external base-collector region.
■
Proper handling of emitter periphery injection and charge storage.
■
Internal base resistance as a function of operating point (conductivity
modulation and emitter current crowding), and emitter geometry.
■
Sufficiently physical model equations allowing predictions of temperature
and process variations, as well as scalability, even at high current densities.
■
Parasitic capacitances, independent on operating point, are available in the
equivalent circuit, representing base-emitter and base-collector oxide
overlaps, that become significant for small-size transistors.
■
Weak avalanche breakdown is available.
■
Self-heating effects are included. Non-quasi-static effects, resulting in a
delay of collector current AND stored minority charge, are modelled as
function of bias.
■
Collector current spreading is included in minority charge and collector
current formulation.
■
Extensions for graded-base SiGe HBTs have been derived using the
Generalized Integral Charge-Control Relation (GICCR); the GICCR also
permits modelling of HBTs with (graded) bandgap differences within the
junctions.
■
Base-emitter tunneling model is available (for example, for simulation of
varactor leakage).
■
Simple parasitic substrate transistor is included in the equivalent circuit.
■
Simple parallel RC network taking into account the frequency dependent
coupling between buried layer and substrate terminal.
■
Parameter extraction is closely related to the process enabling parametric
yield simulation; parameter extraction procedure and list of test structures
are available; HICUM parameters can be determined using standard
measurement equipment and mostly simple, decoupled extraction
procedures.
■
Simple equivalent circuit and numerical formulation of model equations
result in easy implementation and relatively fast execution time.
■
Vertical NQS effects have been implemented in HICUM through “additional
delay times” for both minority charge and forward transfer current.
■ The correlation between base and collector noise is included in the HICUM
compact model.
If you use these features with easily-measurable basic variables (such as
junction capacitances and transit time), the results are more accurate than if
you use SGPM. This improved accuracy applies to digital circuit, small-signal
high-frequency, and especially high-speed large-signal transient simulation.
Also, you can laterally scale HICUM over a wide range of emitter widths and
lengths, up to high collector current densities. The scaling algorithm is generic,
and has been applied to the SGPM (within its validity limits).
In summary, HICUM’s major advantages over other bipolar compact models
are:
■
Scalability
■
Process-based and relatively simple parameter extraction
■
Predictive capability in terms of process and layout variations
■
Fairly simple numerical formulation facilitating easy implementation and
resulting in still reasonable simulation time compared to the (too) simple
SGPM at high current densities
Parameter Description
nt Temperature node
DTEMP Difference between the element temperature and the circuit temperature in
degrees Celsius. The default is 0.0.
Example
The following is an example of a BJT Q1 model:
Q1 1 2 0 4 QM area=1*0.5*5 dtemp=0.002
■
Emitter is connected to node 0.
■
Substrate is connected to node 4.
■
QM references the name of the BJT model.
rsu
² ² ² ² S
CjS Ø
iTS
Csu rCx
ijSC C’ C
² ² ² ² ² ² ² ²
ijBCx CdS ijBCi iAVL
C’BCx C”BCx CjCi CdC
CrBi Ø
B B* B’ Ø
² ² ² ² ² ² ² ² iT
rBx CjEp *
Ø rBi CjEi CdE
CEox
ijBEp iBEt ijBEi
² ² ² ² ²
E’
Tj
rE
²
Ø
P Rth Cth
²
E ²
(a) (b)
Notes:
(a) The external BC capacitance consists of a depletion and a bias-independent capacitance (for example,
oxide) capacitance with the ratio C’BCx / C”BCx being adjusted with respect to proper modeling of the HF
behavior.
(b) Thermal network used for self-heating calculation.
Input Netlist
.DATA test_data vbe vce vsub
0.0 0.0 0.0
0.1 0.0 0.0
0.2 0.0 0.0
0.3 0.0 0.0
0.4 0.0 0.0
0.5 0.0 0.0
0.6 0.0 0.0
0.7 0.0 0.0
0.8 0.0 0.0
0.9 0.0 0.0
1.0 0.0 0.0
.ENDDATA
.OPTION
.TEMP 26.85
VIN 2 0 vbe
VC 1 0 vce
VS 4 0 vsub
VE 3 0 0
Q1 1 2 3 4 hicum
.DC data=test_data
.PRINT DC I(VIN) i2(q1) I(VC) i1(q1) I(VCS) i4(q1)
.MODEL hicum NPN Level=8
+ tref =26.85
+ c10=.3760000E-31 qp0=.2780000E-13 ich=.2090000E-01
+ hfc=.1000000E+01
+ hfe=1.0000000E+00 hjei=.000000E+00
+ hjci=.100000E+01 tr=1.00000000E-9
+ cjei0=.81100E-14 vdei=.950000E+00 zei=.5000000E+00
+ aljei=.18000E+01
+ cjci0=.11600E-14 vdci=.800000E+00 zci=.3330000E+00
+ vptci=.41600E+03
+ rci0=.127800E+03 vlim=.700000E+00 vpt=.5000000E+01
+ vces=.100000E+00
+ t0=.47500000E-11 dt0h=.210000E-11 tbvl=.400000E-11
+ tef0=.180000E-11 gtfe=.140000E+01 thcs=.300000E-10
+ alhc=.750000E+00
+ fthc=.600000E+00
+ latb=.376500E+01 latl=.342000E+00 fqi=.9055000E+00
+ alit=.450000E+00 alqf=.225000E+00
+ favl=.118600E+01 qavl=.111000E-13 alfav=.82500E-04
+ alqav=.19600E-03
+ ibeis=.11600E-19 mbei=.101500E+01 ibeps=.10000E-29
+ mbep=.200000E+01
+ ireis=.11600E-15 mrei=.200000E+01 ireps=.10000E-29
+ mrep=.200000E+01
Internal Transistors
This section lists the HICUM LEVEL=8 internal transistor parameters. This
includes parameter names, descriptions, units, default values, factors, and
notes.
LATB - 0 Scaling factor for Qfc in l_E (“I” is the letter L—not
the number 1)
Peripheral Elements
This section lists the HICUM LEVEL=8 model peripheral element parameters.
This includes parameter names, descriptions, units, default values, factors, and
notes.
Table 88 BJT LEVEL=8 BE Depletion Capacitance
External Elements
This section lists the HICUM LEVEL=8 model external element parameters.
This includes parameter names, descriptions, units, default values, factors, and
notes.
To use the self-heating HICUM feature (in BJT LEVEL=8), set VERS=2.1 and
set an RTH parameter value other than 0. If you use vers=2.0 or RTH=0, then
self-heating is OFF.
The self-heating effect also applies to the circuit temperature as an increased
self-heating temperature. T=Tckt(circuit temperature) + Tsh(self heating
temperature) + dtemp (difference between circuit temperature and ambient
temperature).
The default parameter values for HICUM version 2.2 are follow:
****
.DATA test_data vbe vce vsub
0.0 0.0 0.0
0.1 0.0 0.0
0.2 0.0 0.0
0.3 0.0 0.0
0.4 0.0 0.0
0.5 0.0 0.0
0.6 0.0 0.0
0.7 0.0 0.0
0.8 0.0 0.0
0.9 0.0 0.0
1.0 0.0 0.0
.ENDDATA
.OPTION post
.TEMP 26.85
VIN 2 0 vbe
VC 1 0 vce
VS 4 0 vsub
VE 3 0 0
Q1 1 2 3 4 hicum
.DC data= test_data
.PRINT DC I(VIN) i2(q1) I(VC) i1(q1) I(VCS) i4(q1)
.MODEL hicum NPN Level=8
+ tref = 26.85
+ c10=.3760000E-31 qp0=.2780000E-13 ich=.2090000E-01
+ hfc=.1000000E+01
+ hfe=1.0000000E+00 hjei=.000000E+00
+ hjci=.100000E+01 tr=1.00000000E-9
+ cjei0=.81100E-14 vdei=.950000E+00 zei=.5000000E+00
+ aljei=.18000E+01
+ cjci0=.11600E-14 vdci=.800000E+00 zci=.3330000E+00
+ vptci=.41600E+03
+ rci0=.127800E+03 vlim=.700000E+00 vpt=.5000000E+01
+ vces=.100000E+00
+ t0=.47500000E-11 dt0h=.210000E-11 tbvl=.400000E-11
+ tef0=.180000E-11 gtfe=.140000E+01 thcs=.300000E-10
+ alhc=.750000E+00
+ fthc=.600000E+00
+ latb=.376500E+01 latl=.342000E+00 fqi=.9055000E+00
+ alit=.450000E+00 alqf=.225000E+00
+ favl=.118600E+01 qavl=.111000E-13 alfav=.82500E-04
+ alqav=.19600E-03
+ ibeis=.11600E-19 mbei=.101500E+01 ibeps=.10000E-29
+ mbep=.200000E+01
+ ireis=.11600E-15 mrei=.200000E+01 ireps=.10000E-29
+ mrep=.200000E+01
+ rbi0=.000000E+00 fdqr0=.00000E+00 fgeo=.730000E+00
+ fcrbi=.00000E+00
+ cjep0=.00000E+00 vdep=.105000E+01 zep=.4000000E+00
+ aljep=.24000E+01
+ ceox=.000000E+00
+ cjcx0=.00000E+00 vdcx=.700000E+00 zcx=.3330000E+00
+ vptcx=.10000E+03
+ ccox=.000000E+00 fbc=.1526000E+00
+ ibcxs=.10000E-29 mbcx=.200000E+01 ibcis=.11600E-19
+ mbci=.101500E+01
+ cjs0=.000000E+00 vds=.6000000E+00 zs=.44700000E+00
+ vpts=.100000E+04
+ rcx=.0000000E+00 rbx=.0000000E+00 re=.00000000E+00
+ kf=.00000000E+00 af=.00000000E+00
+ vgb=.1170000E+01 alb=.6300000E-02 alt0=.000000E+00
+ kt0=.0000000E+00
+ zetaci=.1600E+01 alvs=.100000E-02 alces=.40000E-03
+ zetarbi=0.5880E+00 zetarcx=0.2230E+00
+ zetarbx=0.2060E+00 zetare=0.0000E+00
+ rth=0.0 cth=0.0
+ ibets=.00000E+00 abet=.000000E+00
+ itss=.000000E+00 msf=.0000000E+00 tsf=0.000000E+00
+ iscs=.000000E+00
+ msc=.0000000E+00
+ rsu=.0000000E+00 csu=.0000000E+00
.END
Usage Notes
The following information applies to the HSPICE device model for the VBIC99
device:
■
Set LEVEL to 9 to identify the model as a VBIC99 bipolar junction transistor
model.
■
The reference temperature, TREF, equals 27 degrees.
■
The VBIC99 model supports AREA and M factor scaling.
■
This model supports self-heating. Model parameters are RTH and CTH.
Parameter Description
AREA The normalized emitter area. VBIC99 LEVEL=9 model has no area
effect. Default value=1. Area is used only as an alias of the
multiplication factor (M).
OFF Sets the initial condition to OFF for this element in DC analysis. You
cannot use OFF with VBE or VCE.
Parameter Description
Effects of VBIC99
The VBIC99 model includes several effects that are improved compared to the
VBIC95 model:
■
Addition of temperature dependency for several parameters.
■
Base-emitter breakdown model.
■
Reach-through model for base-collector depletion capacitance.
■
High-current beta rolloff effect.
■
Fixed collector-substrate capacitance,
■
Reverse transport saturation current.
Model Implementation
The following parameters were added to the VBIC99 model and are not in the
VBIC95 model.
Equivalent Circuits
E C
IRLAT IFLAT
CFLAT CRLAT
IRVER
IRE IRC
ILE ILC
RBE RBC
ISE ISC
CET CCT
CFVER CRVER
CFN CRN
CπL
CμL
dlπL dlμL
dlRVER
dlFVER
CSB
dlSE dlSC
GSB
Model Equations
Early Factors
The Early factors for the components of the main current I p are derived from
the variation of the depletion widths in the base relative to the base width itself.
Early factor of the lateral current components
⎧ V1 ⎞ 2 ⎫
⎪ ⎛ ⎛ 1 – ------- V-⎞
+δ ⎛ 1 – ------- - + δ ⎞ ⎪
⎪ ⎜ ⎝ VD ⎠ ⎝ VD ⎠ ⎟ ⎪
FLAT = hyp ⎨ 1 – ⎜ ------------------------------------- – --------------------------------------- ⎟ ⋅ δE ⎬
1⎪ ⎜ 1 + EARL EAFL
--------------- 1 + --------------- ⎟ ⎪
⎪ ⎝ 2VD 2VD ⎠ ⎪
⎩ ⎭
Early factor of the forward vertical current component
⎧ ⎛ ⎛ V E2B1 ⎞ 2 V CLB ⎞ 2 ⎞ ⎫
⎪ -------------
- δ ⎛ ------------ + δ⎟ ⎪
⎜4 ⎝ 1 –
⎠
+ 4
⎝
1 –
⎠
⎪ VD VD ⎪
FFVR = hyp ⎨ 1 – ⎜⎜ ---------------------------------------------- – -------------------------------------------- ⎟⎟ ⋅ δE ⎬
T T
1⎪ EARV EAFV ⎪
⎜ 1 + --------------- 1 + --------------- ⎟
⎪ ⎝ 2VD T 2VD T ⎠ ⎪
⎩ ⎭
Early factor of the reverse vertical current component
⎧ ⎛ ⎛ V E1B ⎞ 2 V C2B2 ⎞ 2 ⎞ ⎫
⎪ ⎜ 4 ⎝ 1 – ----------- - + δ 4 ⎛ 1 – -------------- + δ⎟ ⎪
⎪ VD ⎠ ⎝ VD ⎠ ⎪
FRVER = hyp ⎨ 1 – ⎜ ------------------------------------------- – ---------------------------------------------- ⎟ ⋅ δE ⎬
T T
1⎪ ⎜ EARV EAFV ⎟
⎜ 1 + --------------- 1 + --------------- ⎟ ⎪
⎪ ⎝ 2VD T 2VD T ⎠ ⎪
⎩ ⎭
Model Parameters
■
EAFL
■
EAFV
■
EARL
■
EARV
Currents
The ideal diode equations are as follows.
Ve1b/Vt
If 1 = I S ( e – 1)
Ve2b1/Vt
If 2 = I S ( e – 1)
Vc1b/Vt
Ir 1 = I S ( e – 1)
Vc2b2/Vt
Ir 2 = I S ( e – 1)
model parameter: Is
The Ip main current is as follows.
Ip = I flat + I fver – I rlat – I rver
⎛ ⎞
⎜ 4 × Xifv × If 2 ⎟
I fver = ⎜ ---------------------------------------------⎟ ÷ Ffver
⎜ If 2⎟
⎝ 3 + 1 + 16 × ------- Ik ⎠
-
Model parameters:
■
Xifv
■ Ik
Base Current
Forward components
The total forward base current is composed of an ideal and a non-ideal
component. Both components depend on the bottom part of the emitter-base
junction.
Ideal component:
If 2
Ire = --------
Bf
Non-ideal component:
Ve2b1/Vt
Ibf × ( e – 1 )-
Ile = ------------------------------------------------
Ve2b1/Vt Vlf/2Vt
e +e
Model parameters:
■
Bf
■
Ibf
■
Vif
Reverse components
The total reverse base current is composed of an ideal and a non-ideal
component. Both components depend on the bottom part of the collector-base
junction.
Ideal component:
Ir 2
Irc = --------
Br
Non-ideal component:
Vc2b2/Vt
Ibr × ( e – 1)
Ilc = -------------------------------------------------
Vc2b2/2Vt Vlt/2Vt
e +e
Model parameters:
■ Br
■ Ibr
■
Vlr
Substrate current
Forward components
The forward substrate component depends on the bottom part of the emitter-
base junction. It consists of an ideal component, and a component subject to
high injection effects. The XHES parameter determines the fraction that is
subject to high injection.
4 × Xhes × Xes × If 2
Ise = ( 1 – Xhes ) × Xes × If 2 + -------------------------------------------------------
If 2
3 + 1 + 16 × --------
Ik
Model parameters:
■
Xes
■ Xhes
Reverse components
The reverse substrate component depends on the bottom part of the collector-
base junction. It consists of an ideal component, and a component subject to
high injection effects. The XHCS parameter determines the fraction that is
subject to high injection.
4 × Xhcs × Xcs × Ir 2
Isc = ( 1 – Xhcs ) × Xcs × Ir 2 + -------------------------------------------------------
Ir 2
3 + 1 + 16 × --------
Ik
Model parameters:
■
Xcs
■
Xhcs
Charges
Depletion Charges
The Poon-Gummel formulation models the depletion charges.
Emitter-base depletion charge
⎧ ⎫
⎪ ⎪
– Cje ⎪ Vde – Ve2b1 ⎪
Qte = --------------- × ⎨ -----------------------------------------------------
Pe ⎬
1 – Pe ⎪ ------ ⎪
2 2
⎪ ⎛ 1 – Ve2b1 ----------------⎞ + δ ⎪
⎩ ⎝ Vde ⎠ ⎭
Model parameters:
■
Cje
■
Vde
■
Pe
Collector-base depletion charge
⎧ ⎫
⎪ ⎪
– Cjc ⎪ Vdc – Vc2b2 ⎪
Qtc = --------------- × ⎨ -----------------------------------------------------
Pc ⎬
1 – Pc ⎪ ------ ⎪
2 2
⎪ ⎛ 1 – Vc2b2 ----------------⎞ + δ ⎪
⎩ ⎝ Vdc ⎠ ⎭
Model parameters:
■ Cjc
■
Vdc
■
Pc
Substrate-base depletion charge
⎧ ⎫
⎪ ⎪
– Cjs ⎪ Vds – Vsb ⎪
Qts = --------------- × ⎨ ---------------------------------------------
Ps
-⎬
1 – Ps ⎪ ------ ⎪
2 2
⎪ ⎛ 1 – ---------⎞ + δ ⎪
Vsb
⎩ ⎝ Vds⎠ ⎭
Model parameters:
■
Cjs
■
Vds
■ Ps
Series Resistances
The emitter includes the following series resistance:
■
Reex—constant
■
Rein—constant
The collector includes the following series resistance:
■
Rcex—constant
■
Rcin—constant
The conductivity modulation of the base resistances is derived from the fact
that the voltage drop across the epitaxial layer, is inversely proportional to the
electron concentration under the emitter and collector.
Base resistance under the emitter:
2 × Rbev
Rbe = Rbec + -----------------------------------
If 2
1 + 16 × --------
Ik
Base resistance under the collector:
2 × Rbcv
Rbc = Rbcc + -----------------------------------
Ir 2
1 + 16 × --------
Ik
The Rb resistance models the ohmic leakage, across the substrate-base
junction.
Noise Equations
For noise analysis current sources are added to the small signal equivalent
circuit. In these equations:
■ f represents the operation frequency of the transistor.
■
Df is the bandwidth.
When measured at 1 Hz, a noise density is obtained.
Thermal Noise
2 ⋅ k ⋅ Tk-
4-------------------------
iN REEX= ⋅ Δf
REEX
2 ⋅ k ⋅ Tk-
4-------------------------
iN REIN= ⋅ Δf
REIN
2 ⋅ k ⋅ Tk-
4-------------------------
iN RCIN= ⋅ Δf
RCIN
2 ⋅ k ⋅ Tk-
4-------------------------
iN RCEX= ⋅ Δf
RCEX
2 ⋅ k ⋅ Tk-
4-------------------------
iN RBE= ⋅ Δf
RBE
2 ⋅ k ⋅ Tk-
4-------------------------
iN RBC= ⋅ Δf
RBC
2 ⋅ k ⋅ Tk-
4-------------------------
iN RSB= ⋅ Δf
RSB
Lateral Collector Current Shot Noise
2
iN CLAT= 2 ⋅ q ⋅ I FLAT – I RLAT ⋅ Δf
Vertical Collector Current Shot Noise
2
iN CVER= 2 ⋅ q ⋅ I FVER – I RVER ⋅ Δf
Forward-base Current Shot Noise and 1/f Noise
1 – AF AF
2 KF ⋅ MULTI ⋅ I RE ⋅ I LE
iN B= 2 ⋅ q ⋅ I RE – I LE ⋅ Δf + -------------------------------------------------------------------------------------- ⋅ Δf
f
Series Resistance
SPC
RCIN T = RCIN × T N
SNBN
RBCC T = RBCC × T N
SNB
RBCV T = RBCV × T N
SNBN
RBEC T = RBEC × T N
SNB
RBEV T = RBEV × T N
The BJT LEVEL=10 model assumes that REEX and RCEX are temperature
independent.
Depletion Capacitances
TEMP
VDxt = – 3k ---------------- ⋅ 1n ( T N ) + VDx ⋅ T N + ( 1 – T N ) ⋅ Vgap
q
VDx PX
CJ XT = CJx ⋅ ⎛ --------------⎞
⎝ VD XT⎠
Emitter-base Junction
Vgap=VGEB, x=E
Collector-base Junction
Vgap=VGCB, x=C
Substrate-base Junction
Vgap=VGSB, x=S
VD T
EAFL T = EAFL ⋅ -----------
VD
VD T
EARL T = EARL ⋅ -----------
VD
VD T
EAFV T = EAFV ⋅ -----------
VD
VD T
EARV T = EARV ⋅ -----------
VD
TI
BF T = BF ⋅ T N ( AE – SPB ) ⋅ exp {q ⋅ ( VGB – VGE ) ⋅ ⁄ k}
2 TI
IBF T = IBF ⋅ T N ⋅ exp {q ⋅ ( VGJE ⁄ 2 ) ⋅ ⁄ k}
( 1 – SPB )
IK T = IK ⋅ T N
BF T
BR T = BR ⋅ ----------
BF
IBF T
IBR T = IBR ⋅ ------------
IBF
2 TI
ISST = ISS ⋅ T N ⋅ exp {q ⋅ VGSB ⋅ ⁄ k}
( SPB – 1.0 )
TLAT T = TLAT ⋅ T N
TLAT T
TFVR T = TFVR ⋅ -----------------
TLAT
( SX – 1.0 )
TFN T = TFN ⋅ T N
TLAT T
TRVR T = TRVR ⋅ -----------------
TLAT
TFN T
TRN T = TRN ⋅ --------------
TFN
All other model parameters are temperature-independent.
Usage Notes
The following information applies to the HSPICE device model for the UCSD
HBT device:
1. Set BJT LEVEL=11.
2. The default room temperature is 25o C in the HSPICE, but is 27o C in most
other simulators. When comparing to other simulators, do one of the
following:
• set the simulation temperature to 27, or
• set TEMP 27, or
• set .OPTION TNOM=27
3. The set model parameter should always include the model reference
temperature, TREF. The default value for TREF is 27.
4. You can use DTEMP with this model to increase the temperature of individual
elements, relative to the circuit temperature. Set its value on the element
line.
5. The HBT (BJT LEVEL=11) model includes self-heating effects. If you turn
on self-heating, then set RTH to more than zero and SELFT to 1 in the model
card.
Parameter Description
Qxxx BJT element name. Must begin with Q, which can be followed by up
to 1023 alphanumeric characters.
OFF Sets initial condition to OFF for this element in DC analysis. You
cannot use OFF with VBE or VCE.
Model Equations
This section describes the model equations for the HSPICE BJT LEVEL=11
model.
Current Flow
There are seven different current flow calculations for the BJT LEVEL=11
device model.
■
Intrinsic collector current contributions. This model computes the electron
flow between Ei and Ci nodes by using equations similar to the Gummel-
Poon model with modifications to take into account the potential spike that
can appear at the base-emitter or base-collector junctions of HBTs. This
model separates the electron current into forward and reverse components,
Icf and Icr.
Icf=IS * [exp(qVbei/NF/KT) - 1] / D
Icr=IS * [exp(qVbci/NR/KT) -1] / D
In these equations:
D= qb + IS* exp(qVbei/NA/KT) / ISA + IS*exp(qVbci/NB/KT)/ISB
ISA, ISB, NA and NB are new parameters. ISA and ISB approximate the
transition currents, from base-transport controlled to potential-barrier
controlled, current flow.
qb partially retains the standard BJT model form (a fractional increase in the
base charge associated with the bias changes).
qb= q1/2 * [1+(1+4*q2)0.5]
q1=1/ [1 - Vbci/VAF -Vbei/VAR]
q2=IS/IK*[exp(qVbei/NF/KT)-1]
In the preceding equations, qb omits the reverse knee current contribution.
As noted in the following, qb is not used to define the ac model in the fashion
of the Gummel-Poon model.
The total collector current Icc is:
Icc=Icf - Icr
This formulation uses the IS, NF, VAF, VAR, and IK parameters, established
in the SPICE BJT model in addition to the ISA, ISB, NA, and NB parameters
described above.
■
Intrinsic Base-Emitter Diode. Ideal and non-ideal components are included:
■
Extrinsic Base-Collector Diode. This diode has customary I-V
characteristics with its own saturation current and ideality factor.
Ibcx= ISCX * [exp(q Vbcx/ NCX /KTop) - 1]
■
Substrate-Extrinsic Collector Diode. This diode allows for conducting
substrates. Use it primarily for SiGe HBTs.
Ics= ICS * [exp(-q Vcs/ NCS /KTop) - 1]
In accordance with the model topology, the external currents through the
E,B, and C nodes are:
Ib= Ibei + Ibex - Ibk + Ibci + Ibcx
Ic= Icc + Ibk -Ibci-Ibcx - Ics
Charge Storage
This section describes the following different charge storage calculations for the
HSPICE BJT LEVEL=11 device model.
■
Base-Emitter Charge. The overall charge stored at the base-emitter junction
has components associated with the base-emitter depletion layer:
• Qbej, which is current-independent.
• Qbediff, a collector current-dependent charge. Qbediff corresponds to a
portion of the base charge, and the (collector current-dependent) base-
collector charge.
Qbe= Qbej + Qbediff
■
Base-Emitter Depletion Charge, Qbej. The depletion charge, Qbej, follows
equations standard for SPICE, modified to allow specification of a minimum
capacitance CEMIN (corresponding to reach-through to an n+ layer).
As studied by Chris Grossman, there is often an extra component of charge
storage at the base-emitter heterojunction of HBTs, associated with a
minimum in the conduction band energy profile.
Qbej is computed using DepletionCapMod.
Define:
Vmin= VJE*[1-(CJE/CEMIN)(1/MJE)]
(the critical voltage for attaining the minimum capacitance value)
If Vbei<FCE*VJE and Vbei<Vmin:
Qbej=CEMIN*(Vbei-VJE)+CEMIN*VJE*MJE/(MJE-1) *
(CJE/CEMIN)(1/MJE)
Cbej=dQbej/dVbei=CEMIN
If Vbei<FCE*VJE and Vbei>Vmin:
Qbej= -CJE*VJE*(1-Vbei/VJE)(1-MJE) / (1-MJE)
Cbej= CJE*(1-Vbei/VJE)(-MJE)
If Vbei>FCE*VJE, and CJE>CEMIN*(1-FCE)MJE:
Qbej=-CJE*VJE/(1-FCE)MJE*[(1-FCE)/(1-MJE)+FCE-Vbei/VJE
-MJE*(FCE-Vbei/VJE)2/2/(1-FCE)]
Cbej=CJE/(1-FCE)MJE*[1+MJE*(Vbei/VJE-FCE)/(1-FCE)]
If Vbei>FCE*VJE, and CJE<CEMIN*(1-FCE)MJE,
Qbej=CEMIN*(Vbei-VJE)+CEMIN*VJE*MJE/(MJE-1)
*(CJE/CEMIN)(1/MJE)+CJE*VJE*(Vbei/VJE-FCE)2*MJE/2/(1-
FCE)(MJE+1)
Cbej=CEMIN + CJE*VJE*MJE*(Vbei/VJE-FC)/(1-FCE)(MJE+1)
■
Base-Emitter Diffusion Charge, Qbediff. The diffusion charge in HBTs is
associated with contributions from minority carriers in the base, and from
mobile charge in the collector depletion region. In homojunction transistors,
diffusion charge storage in the emitter is also present. The LEVEL=11
model evaluates the base and collector-depletion region contributions
separately (if necessary, the emitter charge storage can be associated with
the base contribution).
• Specify the base charge through the base transit time, TFB. This transit
time varies with bias through several mechanisms:
• The Early effect causes a change in transit time with junction voltage.
• In heterojunction transistors, there is frequently a minimum in the
conduction band, on the base side of the base-emitter (and potentially
base-collector) heterojunction. Minority carriers tend to accumulate in
these potential wells.
The stored charge adds to the base charge (to a good approximation).
In the lowest order, the charge stored is directly proportional to the
collector current, and thus contributes to TFB. For a greater degree of
accuracy, the depth of the potential well on the emitter side varies with
Vbe. Similarly, the amount of charge stored at the base-collector side
varies with Vbc.
The equations used to describe the effects are:
TFBt=TFB*(1+Vbei/VAR+Vbci/VAF) +
TBEXS*exp(-q(Vbei-VJE)/NA/KTop) +
TBCXS*exp(q(Vbci-VJC)/NB/KTop)
Note: Different signs are associated with the BE and BC
junction effects. The value of the T temperature to
describe these effects is assumed to be Top.
You can use any of these methods to specify collector
charge:
■
Extrinsic Base-Collector Charge, Qbcx. The Qbcx stored charge consists of
a depletion charge and a diffusion charge.
Standard SPICE does not use the diffusion charge component. However,
this component can be an important contribution to saturation stored-charge
in many HBTs (in addition to the contribution associated with the intrinsic
base-collector junction).
The corresponding charge storage time, TRX, might be different from the
intrinsic time, TRI. This difference occurs because of implant-induced
recombination, surfaces, or other structural changes.
The depletion charge corresponds to a standard depletion region
expression (without considering charge density modulation due to current),
modified to allow for a minimum value of capacitance under a reach-through
condition.
Furthermore, as indicated in the following, if you assign a value other than
unity to the XCJC variable, then the depletion charge is partitioned between
the Bx-Cx capacitance and the B-Cx capacitance.
Qbcx= TRX*Ibcx + XCJC*Qbcxo
In the preceding equation, Qbcxo is the depletion charge.
As a result the dependences of Ibcx on Vbcx, a diffusion capacitance results
from the formulation:
Cbcxdiff= TRX* dIbcx/dVbcx
Base-Extrinsic Collector Charge (Qbcxx), and Treatment of XCJC
In standard SPICE, XCJC indicates the fraction of overall Cbc depletion
capacitance that should be associated with the intrinsic base node. The
remaining fraction (1-XCJC) is attached to the base terminal. HBT Spice
uses a similar assignment: the depletion charge associated with the
extrinsic base-collector junction is partitioned between the Bx node and the
B node:
Qbcx= TRX*Ibcx + XCJC* Qbcxo
has been defined above, between the Bx and Cx nodes, and charge
Qbcxx= (1-XCJC)*Qbcxxo
Noise
The LEVEL=11 model includes noise current generators, similar to those in
standard Spice. The noise current generators have magnitudes in units of A2/
Hz, and are computed based on 1Hz bandwidth. The noise sources are placed
in parallel with corresponding linearized elements in the small signal model.
Sources of 1/f noise have magnitudes that vary with the frequency (f); you can
use a BFN exponent, if you do not see the exact f-1 behavior.
inc2=2*q*Icc
inb2= 2 * q *Ibe + KFN * IbeAFN / f BFN
inre2= 4 * K* Td / RE
inrbx2= 4 * K * Td / RBX
inrbi2= 4 * K * Td / RBI
inrcx2= 4 * K * Td / RCX
inrci2= 4 * K * Td / RCI
inrex2= 4 * K * Td / REX
Equivalent Circuit
Rcz
Qos
Cx
S
Rci
-
- Qbcz los
lbcz
+ T
+ Qbczz Ci
-
Qbci ibci ibk +
B Bx Rbi + Qth
-
loc Ith Rth
Rbz Rez +
Qbe ibei
Ex -
Ibez
Ei
Re
Rcz
Ccs
Cz
S
ibczq
Rci
Cbcz gcs
gbcz T
Cbzz ibcq gbci ibk ibca Ci
Cbci
+
B Bz Rbi Qth
-
icc Ith Rth
Rbz Rez
gez Ez
Cbe
ibeq gbei ibea
Ei
Re
■
makes use of the advanced capabilities of HICUM LEVEL=2 and the related
knowledge base for generating geometry scalable and statistical models
■ offers an easy migration path from a conventional, single-transistor-based,
to a process-based geometry scalable parameter extraction and model
usage to meet today’s requirements for advanced integrated circuit design.
The model parameters have a clear (physical) meaning and many of them are
similar to HICUM LEVEL=2 parameters.
Parameter Description
Parameter Description
References
[1] C. McAndrew, J. Seitchik, D. Bowers, M. Dunn, M. Foisy, I. Getreu, M.
McSwain, S. Moinian, J. Parker, D. Roulston, M. Schroter, P. van Wijnen,
and L. Wagner, “VBIC95: The vertical bipolar intercompany model,” IEEE
Journal of Solid State Circuits, vol.31, p.1476-1483, 1996.
HSPICE ships hundreds of examples for your use; see Listing of Input
Demonstration Files for paths to demo files.
For libraries with multiple models of a specific active or passive device element,
you can use the automatic model selector in HSPICE to automatically find the
proper model for each transistor size.
This chapter lists device libraries that you can use. It includes the following
topics:
■
Overview of Library Listings
■
Analog Device Models
■ Behavioral Device Models
■
Bipolar Transistor Models
■
Diode Models
■ JFET and MESFET Models
ssm2210 ssm2220
rcfilt tline
Diode Models
The search path for diode models is:
$installdir/parts/dio
The required element syntax is:
Xyyyyy anode cathode modelname
Optional parameters are:
isat=value, tt=value
Table 114 Diode Model Names
dsk4a3
331
Index
C
332
Index
D
element 32 ferromagnetic 41
equation selector option 51 model parameters 39
equations 114 magnetic 37
models, SPICE 112 element output 39
parameters model parameters 37
junction 168 model parameters 39
metal and poly 59 A 42
temperature, equations 15 ALPHA 42
CAPOP model parameter 111 MS 42
charge current
base-collector 303, 305 BJT Level 10 281
base-emitter 300 breakdown 299
substrates 306 collector 298
VBIC99 269 convention
charge stoarge, HBT model 300 BJTs 177
charges, BJT Level 10 284 JFETs 114
circuit schematics diodes 63
BJT Level 10 277 epitaxial 208
HBT model 307 HBT model 298
circuits Curtice model 111, 127
BiCMOS 160
BJT 176, 178
ECL 160
D
DC
TTL 160
BJT 163
CJBR model parameter 81 equations
CJGR model parameter 81 JFETs and MESFETs 135–137
CJSR model parameter 81 equations, BJTs 189
CMC R2 resistor model 20 parameters
collector current 298 JFETs LEVEL=1 126
collector diode, substarte 300 JFETs LEVEL=2 126
JFETs LEVEL=3 128
collector-substrate junction 163
conductance DC equations, TOM3 148
BJTs 178, 179 DC model, JFET Level 8 153
diodes 64 DC operating point, BJT Level 10 278
GMIN 52 DCAP 51, 161
GMINDC 52 equation selector 161
JFETs 113 option 51
JFETs and MESFETs 113, 115 JFETs and MESFETs capacitance 113
control options 113 JFETs and MESFETs models 113
convergence 52 overriding in BJTs models 161
setting 50 DCCAP 51, 161
convergence option 51, 193
BJT model 161 JFETs and MESFETs capacitance 113
JFETs and MESFETs 113 depletion
problems capacitance
diodes 52 DCAP equation selector 52, 113, 161
core equations 73
Jiles-Atherton 37 depletion capacitance, BJT Level 10 289
333
Index
E
334
Index
F
335
Index
H
336
Index
K
337
Index
M
338
Index
N
339
Index
O
340
Index
S
341
Index
U
342
Index
W
W distribution 9
W capacitor parameter 53 model RC 8
warnings resistance 12
capacitance too large 34, 40 WM capacitor parameter 53
invalid value for CRATIO 9 WP capacitor parameter 53
IS parameter too small 57
resistance smaller than RESMIN 13 X
wire
capacitance 13 X statement 4
model XCJC 305
effective length and width 12
parameters 8, 12 Z
resistance calculation 12
model capacitance Zener diodes 50
calculation 13
343
Index
Z
344