Next: The WRspice User Interface
Up: Superconductor Devices
Previous: Josephson Junctions
Josephson Junction Model
Type Name: jj
The Josephson junction model is an extended version of the RSJ model
as used by Jewett.
Unless stated otherwise, this information also applies to the
Verilog-A Josephson junction model provided with WRspice in the
The parameters marked with an asterisk in the area column scale
with the ics parameter given in the device line, not necessarily
linearly. The present model paradigm assumes that the model
parameters apply to a ``reference'' junction, which is a typical
mid-critical current device as produced by the fouhdry.
Instantiations derive from the reference device for a desired critical
current. Appropriate scaling, not necessarily linear, will be applied
when formulating instance capacitance and conductances.
|JJ Model Parameters
||Quasiparticle branch model
||Critical current model
||Reference junction critical current
|vg or vgap
||Gap voltage spread
||Reference junction capacitance
||Reference junction capacitance per critical current
||Capacitance area/edge scaling
|rsub or r0
||Reference subgap resistance
|rnorm or rn
||Reference normal state resistance
||Conductance area/edge scaling
|icfct or icfact
||Ratio of critical to gap currents
||Voltage to specify fixed shunt resistance
||Shunt resistor series inductance, constant part
||Shunt resistor series inductance, proportional part
||Time step phase change limit
||Critical current first zero
Detailed information about these parameters is presented below.
If this flag is given a nonzero integer value, the junctions will be
modeled as a ``pi junction'' meaning that the zero-current phase is
rather than zero. Such devices have been fabricated using
ferromagnetic barrier materials. Although these devices have some
interesting behavior, they are not at this point available or used
with any frequency.
The rtype parameter determines the type of quasiparticle branch
modeling employed. Legal values are as follows:
||The junction is completely unshunted, all shunt conductances
set to zero.
||Standard piecewise-linear model.
||Analytic exponentially-derived approximation.
||Fifth order polynomial expansion model.
||``Temperature'' variation, allow modulation of the gap parameter.
Values for rtype larger than 1 are not currently supported in
the Verilog-A model supplied with WRspice in the Verilog-A
The default is rtype=1. Setting rtype=0 will disable
modeling of the quasiparticle current, effectively setting the RSJ
shunt resistance to infinity. Conditions with rtype=1 and 2 are as described by Jewett, however it is not assumed that the
normal resistance projects through the origin. The icfact
parameter can be set to a value lower than the default BCS theoretical
value to reflect the behavior of most real junctions. The
quasiparticle resistance is approximated with a fifth order polynomial
if rtype=3, which seems to give good results for the modeling of
some NbN junctions (which tend to have gently sloping quasiparticle
Rtype=4 uses a piecewise-linear quasiparticle characteristic
identical to rtype=1, however the gap voltage and critical
current are now proportional to the absolute value of the control
current set with a control=src_name entry in the device
line. This is to facilitate modeling of temperature changes or
nonequilibrium effects. For control current of 1 (Amp) or greater,
the full gap and critical current are used, otherwise they decrease
linearly to zero. If no device control source is specified, the
algorithm reverts to rtype=1. It is expected that a nonlinear
transfer function will be implemented with a controlled source, which
will in turn provide the controlling current to the junction in this
mode. For example, the controlling current can be translated from a
circuit voltage representing temperature with an external nonlinear
source. The functional dependence is in general a complicated
function, but a reasonable approximation is
1 - (T/Tc)4
. See the
examples for an example input file (ex8.cir) which illustrates
It is currently not possible to use other than the piecewise linear
model with temperature variation. If rtype=4, then legal values
for the critical current parameter are cct=0 (no critical
current) and cct=1 (fixed critical current). If another value
is specified for cct, cct reverts to 0. Thus, magnetic
coupling and quasiparticle injection are not simultaneously available.
The cct variable can take on the following values:
||No critical current.
||Fixed critical current.
||Sin(x)/x modulated supercurrent.
||Symmetric linear reduction modulation.
||Asymmetric linear reduction modulation.
Values for cct larger than 1 are not currently supported in the
Verilog-A model supplied with WRspice in the Verilog-A examples.
The control instance parameter should be used with devices using
cct 2,3, or 4. With cct=2, the first zero is equal to the
value of the model parameter icon. For cct=3, the maximum
critical current is at control current zero, and it reduces linearly
to zero at control current = icon
. Junctions with cct=4
have maximum critical current at control current = -
linear reduction to zero at control current = +
icon. If cct is specified as 2, 3, or 4, the area parameter, if given, is set
to unity. Otherwise, the model parameters are scaled appropriately by
the area before use.
range: 1nA - 0.1A
This is the critical current of the reference junction, which defaults
to 1.0mA if not given. This parameter is not used if cct is 0.
the superconducting current through a Josepjson junctions is
I = Icsin()
is the critical current. and the junction ``phase'' is
The V(t) is the junction voltage, and
is the magnetic
The icrit parameter should not be confused with the
ics instance parameter. The latter is actually a scale
factor which specifies the instantiated device critical current as
well as appropriately scaling conductances and capacitance, from
the model reference current which is icrit.
- vgap or vg
range: 0.1mV - 10.0mV
This parameter specifies the gap voltage, which in a hysteretic
Josephson junction is a voltage at which there is a large and abrupt
increase in conductivity. This parameter is material dependent, the
default 2.6mV is appropriate for niobium junctions.
This specifies the assumed width, in voltage, of the quasiparticle
step region, or gap. In this region, current increases sharply with
increasing voltage. The default value of 80uV is reasonable for
high-quality niobium/aluminum oxide Josephson junctions independent of
range: 0.0 - 1nF
This is the capacitance of the reference junction, in farads. This
will override the cpic parameter if given, setting a fixed value
for reference junction capacitance, invariant with icrit. If
not given, junction specific capacitance is set via the cpic
parameter, see below.
range: 0.0 - 1e-9
This supplies the default capacitance per critical current in F/A.
This defaults to the MIT Lincoln Laboratory SFQ5EE process <a
href="tolpygo">[Tolpygo]</a> value (0.7pF for 1.0 mA), and will set
the junction capacitance if <tt>cap</tt> is not given. With
<tt>cap</tt> not given, changing <tt>icrit</tt> will change the
assumed capacitance of the reference junction.
range: 0.0 - 1.0
This is a new parameter in the current model, which is intended to
account for nonlinearity in scaling of capacitance with area (or
critical current, we actually define ``area'' as the actual over the
reference critical current). It is anticipated that the actual
junction capacitance consists of two components: a physical area
dependent ``bulk'' term, and a perimeter-dependent fringing term. The
cmu is a real number between 0 and 1 where if 0 we assume no
perimeter dependence, and if 1 we assume that all variation scales
with the perimeter. The default value is 0. The capacitance of an
instantiated junctions is as follows:
C = cap(A(1 - cmu) + cmu)
is the ``area'' scaling factor, which is the ratio of the
junction critical current to the reference critical current.
range: 8mV - 100mV
This is the product of the reference subgap resistance and the
reference device critical current. This parameter is commonly
provided by foundries, and is a standard indicator for junction
quality (higher is better). Values tend to decrease with increasing
critical current density. This defaults to the value for the MIT
Lincoln Laboratory SFQ5EE process, which is 16.5mV, The
reference junction subgap resistance is obtained from the value of
this parameter and the critical current, unless given explicitly.
- rsub or r0
range: 8mV/icrit - 100mV/icrit
The reference junction subgap resistance can be given directly
with this parameter, and a given value will override the
vm value if also given.
range: 1.5mV - 1.9mV
This is the product of the reference junction ``normal state''
resistance and the critical current, where the normal state resistance
is the differential resistance measured well above the gap. The
default value is that provided for the MIT Lincoln Laboratory SFQ5EE
process which is 1.65mV. This too is a commonly given
parameter from Josephson foundries for characterizing junctions. If
not specified explicitly, this provides the reference junction normal
state resistance from the critical current.
- rnorm or rn
range: 1.5mV/icrit - 1.9mV/icrit
The reference junction normal state resistance can be given explicitly
with this parameter, which will override icrn if this is also
range: 0.0 - 1.0
This is analogous to cmu, and applies to the subgap and normal
conductances. The vm, in particular, may vary with junction
physical size, with small junctions having lower vm than larger
ones. This parameter should capture this effect. It is taken that a
significant part of the conductivity is due to defects or
imperfections around the periphery of the junction area, and the
contribution would therefor scale with the perimeter. The scaling for
conductivity is as follows:
Gx = Gx0(A(1 - gmu) + gmu)
refers to either the subgap or normal conductance,
is the same parameter for the reference junction. The A
the scaling parameter, that is, the ratio of instance to reference
critical currents. The default value is 0, meaning that scaling is
assumed purely linear, which will be the case until a number is
provided through additional data analysis. It may prove necessary to
have separate scaling parameters for subgap and above gap condutance,
at which time a new model parameter may be added.
- icfct or icfact
range: 0.5 - /4
This parameter sets the ratio of the critical current to the
quasiparticle step height. Theory provides the default value of
which is usually adequately close. Characterization of
fabricated junctions would provide an improved number.
range: 0.0 - 2.0mV
This parameter is unique in that it does not describe an as-fabricated
junction characteristic. Rather, it is for convenience in specifying
a shunt resistance to use globally in SFQ circuits, If given (in
volts) conductance will be added automatically so that the product of
the total subgap conductance and the critical current will equal vshunt. This avoids having to calculate the value of and add an
explicit resistor across each Josephson junction, as used for damping
in these circuits. The designer should choose a value consistent with
the process parameters and the amount of damping required. Higher
values will provide less damping, usually critical damping is desired.
This parameter defaults to 0, meaning that no additional demping is
supplied by default.
range: 0.0 - 2.0pH
range: 0.0 - 10.0pH/
These parameters specify series parasitic inductance in the external
shunt resistor. the vshunt parameter must be given a value such
that the added external conductance is positive, or these parameters
are ignored. The inductance consists of a constant part (lsh0)
assumed to come from resistor contacts, plus a value (lsh1)
proportional to the resistance in ohms, intended to capture the length
range: 0.001 - 1.0
This is mainly for compatibility with the Verilog-A Josephson junction
model provided with WRspice in the Verilog-A examples. It has the
same effect and range as the WRspice dphimax parameter, and
will override dphimax if given. If not given, dphimax
will be used. This is the maximum phase change allowed between
internal time points.
range: 1e-4 - 1.0
This parameter applies when the cct parameter is set to one of
the choices larger the 1, where critical current modulation is
modeled. The value of icon is the first value for (assumed)
full suppression of critical current.
The parameter is not currently recognized by the Verilog-A Josephson
junction model provided with WRspice, as that model does not
currently support values of cct larger than 1.
Next: The WRspice User Interface
Up: Superconductor Devices
Previous: Josephson Junctions
Stephen R. Whiteley