M B - P V C S G (CSG) Mosfet: Odeling of Uilt IN Otential Ariations OF Ylindrical Urrounding ATE S
M B - P V C S G (CSG) Mosfet: Odeling of Uilt IN Otential Ariations OF Ylindrical Urrounding ATE S
M B - P V C S G (CSG) Mosfet: Odeling of Uilt IN Otential Ariations OF Ylindrical Urrounding ATE S
5, October 2012
DOI : 10.5121/vlsic.2012.3506 67
MODELING OF BUILT-IN POTENTIAL VARIATIONS
OF CYLINDRICAL SURROUNDING GATE (CSG)
MOSFETS
Santosh Kumar Gupta
1
and S. Baishya
2
Department of Electronics & Communication Engineering, National Institute of
Technology Silchar, Assam 788 010, INDIA
1
santoshty@gmail.com,
2
sb@nits.ac.in
ABSTRACT
Due to aggressive scaling of MOSFETs the parasitic fringing field plays a major role in deciding its
characteristics. These fringing fields are now not negligible and should be taken into account for deriving
the MOSFET models. Due to this fringing field effect there are some charges induced in the source/drain
extension regions which will change the potential barrier at the source-channel from its theoretical
nominal values. In this paper an attempt has been made to model variation of built-in potential variation
for a cylindrical surrounding gate MOSFET. The model has been verified to be working in good agreement
with the variations of gate length and channel radius.
KEYWORDS
Barrier lowering, cylindrical surround gate (CSG) MOSFET, fringing field, short channel effects (SCEs)
1. INTRODUCTION
Silicon technology is spreading widely and supports the IT revolution which is now reshaping
society. The technology keeps improving year on year as the chip sizes are going on decreasing
with the simultaneous increase in the transistor speeds. However as we reduce the dimensions,
short channel effects (SCEs) cause several problems such as threshold voltage lowering,
increased substrate bias effect etc. [1]-[3]. To continue the scaling of Si MOSFETs in the
nanometre regime, innovative device structures have been proposed. One such structure is the
surrounding cylindrical gate MOSFET where the gate has control over the channel from all the
directions reducing the short channel effects and improving the sub-threshold characteristics.
However this structure is also not free from the SCEs.
But due to scaling of the MOSFETs the fringing field effects are becoming comparable to the
gate capacitance and cant be neglected anymore [4], [5]. Due to these fringing fields originating
from the sidewalls of the gates are capable of changing the source-channel and drain-channel
barriers as well. This change in the barrier heights may contribute to significant change in the
electrical characteristics of a highly scaled MOS transistor [6]-[15].
In this paper, an attempt has been made to account for these fringing fields and see the effect on
the barrier lowering. The results obtained are encouraging but not as accurate. Even though, this
gives an insight about the changes in the barrier potentials due to these fringing fields. In the
section II of the paper the preliminary concept of the variation of the built in potential variation
has been discussed and a general equation has been derived. In section III, we have derived
International Journal of VLSI design & Communication Systems (VLSICS) Vol.3, No.5, October 2012
68
expressions accounting for the fringing field effect of cylindrical surrounding gate MOSFETs. In
section IV, the model has been verified to be working almost in agreement with the 3D-TCAD
simulations of cylindrical surrounding gate MOSFETs. Variations in the potential barrier due to
the change in the radius and length of the channel have been taken into account for the
verification of the model developed.
2. MODEL DEVELOPMENT
In this section first we discuss about the general discussion about the built in potential variation.
In the latter part of the section the mathematical formulations for the case of cylindrical
surrounding gate MOSFET has been derived.
2.1. Variation of Built-in potential at Source-Channel junction
In weak inversion we assume that potential at the source/drain-channel boundary (junction) is
fixed at built in potential ( )
b
and is given by (for a p-type body, with doping concentration
equal to
A
N ).
| |
= +
|
\
log
2
g
A
b
i
E
N kT
q n
(1)
Here it is assumed that Fermi energy level in the source is aligned with the source conduction
band.
g
E is the silicon energy band gap,
i
n is the intrinsic carrier concentration, q is the unit
charge, T is the absolute temperature in degree Kelvin, and k is the Boltzmann constant. For the
doping levels of the order
20 -3
10 cm , the value of the
b
may be different than predicted by the
Fermi-Dirac (FD) statistics. The above equation (1) also assumes that there is no spatial variation
of carrier density inside the source/drain. But, the actual case may be very different than the
above and
b
may depend upon many other parameters like gate length, film thickness and
biasing voltages etc.
For accurate prediction of the source/drain - channel boundary potential, we need to solve the 1-D
Poissons equation inside the source,
( )
=
2
2
s D
Si
d q
n N
dy
(2)
Where,
s
n is the electron concentration in the source,
Si
is the permittivity of silicon and
D
N is
the source doping density.
In equilibrium for an undoped body (channel), the electron concentration in the source assuming
the Maxwell-Boltzmann (MB) statistics can be given by
=
T
V
s i
n n e (3)
From (2) and (3) we get
( )
| | | |
| |
= = =
| |
\ \
2
2
1
T T
V V D i
s D i D
Si Si Si D
qN n d q q
n N n e N e
N
dy
(4)
When =
bi
V , =
s D
n N and we have =
bi
T
V
V
D i
N n e . Substituting this into (4), we get
International Journal of VLSI design & Communication Systems (VLSICS) Vol.3, No.5, October 2012
69
( )
| |
| |
|
|
| = =
|
|
|
|
\
\
2
2
1 1
bi
T T
bi
T
V
V V D i D
V
Si Si
V
i
qN n qN d
e e
dy
n e
(5)
Integrating the above equation with the boundary conditions
= 0
d
dy
at =
s
y y and =
bi
V
Multiplying both sides by 2d dy we get-
( ) ( )
| | (
| |
| (
= = +
|
| (
\
|
(
\
2
2
0
2
2
2 2 1
bi bi
T T
V V
V V D D
T
Si Si
qN qN d d d d
e V e C
dy dy dy
dy
(6)
Where
0
C is a constant of integration and may be found by using the boundary conditions as
follows
Using
= 0
d
dy
at we =
s
y y get-
[ ]
=
0
2
D bi T
Si
qN V V
C (7)
Substituting above value of
0
C in to (6) we get
( )
[ ]
(
| |
(
= +
|
(
\
(
2
2
2
bi
T
V
D bi T V D
T
Si Si
qN V V
qN d
V e
dy
(8)
Above equation can be simplified to
( )
( )
| |
(
=
|
(
\
(
2
2
1
bi
T
V
bi V D T
Si T
V
qN V d
e
dy V
(9)
Above equation do not has an analytical solution and needs to be numerically evaluated to get an
accurate value of the source-channel boundary potential. Simpler, yet reasonably accurate value
of the boundary potential can be obtained using the depletion approximation. Using the boundary
potential in terms of the electric field at the boundary may be given by [2]
=
2
2
Si b
b bi
D
E
V
qN
(10)
Where
b
E is the electric field at the source-channel boundary.
2.2. Evaluation of E
b
for cylindrical surrounding gate (CSG) MOSFETs
A cylindrical surrounding gate (CSG) MOSFET is symmetrical about its axis so we can use only
2D Poissons equation for the calculation of the potential variations. For a CSG MOSFET, in the
weak inversion the channel potential can be calculated using the following 2-D Poissons
equation-
International Journal of VLSI design & Communication Systems (VLSICS) Vol.3, No.5, October 2012
70
( ) ( ) ( ) ( )
| |
+ = |
|
\
2
2
, ,
1
A
Si
r z r z
qN
r
r r r
z
(11)
Where r is the distance in the radial direction of the cylindrical channel region, ( ) , r z is the
potential variations in the channel region, N
A
is the doping concentration in the channel region.
Figure 1 Cylindrical surrounding gate MOSFET (here only channel region with gate oxide is
shown for simplicity)
Assume parabolic potential profile in the vertical direction of the channel.
( ) ( ) ( ) ( ) = + +
2
1 2 3
, r z c z c z r c z r (12)
Where ( )
1
c z , ( )
2
c z and ( )
3
c z are the constants and can be determined the following boundary
conditions-
1. Surface potential is a function of z only at the Si/SiO
2
interface.
( ) ( ) = ,
S
R z z (13)
2. The electric field in the center of the silicon pillar is zero due to symmetry along z-axis.
( )
=
0
,
0
r
r z
r
(14)
3. The electric field at = r R (Si/SiO
2
interface) is continuous.
( )
( )
= (
,
f
GS s FB
Si
r R
C r z
V z V
r
(15)
Where
ox
t - gate oxide thickness, ( )
S
z - surface potential,
GS
V - gate to source voltage,
FB
V - flat
band voltage,
| |
= +
|
\
ln 1
ox
f ox
t
C R
R
is the gate capacitance per unit area,
4. The potential at the source end is
( ) ( ) = = 0, 0 0
s b
(16)
International Journal of VLSI design & Communication Systems (VLSICS) Vol.3, No.5, October 2012
71
5. The potential at the drain end is
( ) ( ) = = + , 0
s b DS
L L V (17)
Where
DS
V is the drain to source voltage.
Solving (13) (15) constants ( )
1
c z , ( )
2
c z and ( )
3
c z are found as:-
( ) ( ) [ ]
( )
( ) ( ) [ ]
(
= +
(
=
(
= +
(
1
2
3
1
2 2
0
1
2 2
f f
S GS FB
Si Si
f f
S GS FB
Si Si
C R C R
c z z V V
c z
C R C R
c z z V V
(18)
Substituting above values in (12) the 2-D potential in the channel is given by-
( ) ( ) [ ] ( )
(
= + + (
(
2
, 1
2 2 2
f f f
S GS FB GS s FB
Si Si Si
C R C R C
r z z V V V z V r
R
(19)
Substituting this value of ( ) , r z in to (1), we get-
( )
( )
=
2
2
2
S
S
d z
z
dz
(20)
Where =
2
2
f Si
C R and ( ) =
2
A Si GS FB
qN V V .
Solution of the above equation is of the form as given by
( )
= +
2
z z
S
z Ae Be (21)
Where A and B are constants, found using (16) and (17), as given by
( )
( )
( )
( )
| |
+ +
|
\
=
| |
+
|
\
=
2
2
1
2sinh
1
2sinh
L
b DS
L
b DS
e V
A
L
e V
B
L
(22)
Now, the electric field at the barrier (source-channel boundary), i.e. at = 0 z , can be given by the
potential of the channel at the boundary, i.e.,
International Journal of VLSI design & Communication Systems (VLSICS) Vol.3, No.5, October 2012
72
( ) ( )
=
=
(
| |
( = = + +
` |
(
\
)
2
0
0
,
1
2 2
f f S
b
Si Si
z
z
C R C r d r z d z
E
dz R dz
(23)
Or
( )
( )
=
=
(
| |
( = = + +
` |
(
\
)
2
0
0
,
1
2 2
f f z z
b
Si Si
z
z
C R C r d r z
E Ae Be
dz R
( )
| |
| =
|
\
2
1
2 2
f f
b
Si Si
C r C R
E A B
R
(24)
And at the surface, i.e., = r R , the electric field will be
( ) ( ) = =
b
E A B B A (25)
Now, substituting the above value of the electric field at the boundary in to (10), we get
( )
+
=
2 2 2
2
2
Si
b bi
D
A B AB
V
qN
(26)
Substituting the values of A and B from (22) into (26), we get
( ) ( )
{ }
| |
+ + + + =
|
\
2
2 2 2 2 2 2 2
2
4 2 4 2 0
Si Si DS D Si DS D bi
u m n u m n V qN d V qN d V (27)
where
| |
+ = =
|
\
2 2 b b
u u ,
( )
= 1
L
e m,
( )
= 1
L
e n , ( ) =2sinh d L
Solution of the above quadratic equation is given by
=
2
4
2
u (28)
Where
( ) =
2
2
Si
m n , ( ) = +
2 2
4 2
Si DS D
m n V qN d , and
| |
= +
|
\
2 2 2
2
4 2
Si DS D bi
V qN d V .
The above equation can be re-written as
= + =
2 2
2
4 4
2 2
b
u (29)
From which the value of
b
can be given by
International Journal of VLSI design & Communication Systems (VLSICS) Vol.3, No.5, October 2012
73
=
2
2
4
2
b
(30)
3. RESULTS AND DISCUSSIONS
For the validation of the above model a typical cylindrical surrounding gate (CSG) MOSFET was
simulated using the Sentaurus TCAD tool from the Synopsys as shown in Fig. 2. The CSG
MOSFET parameters used for the simulation and validation of the model developed are shown in
the Table 1. The structures have been generated using the structure editor and the device
simulations have been carried using the Sentaurus device tool of Sentaurus TCAD [16]. The
potential has been extracted using the Tecplot SV and the current voltage characteristics using the
Inspect tool.
The CSG MOSFETs modelled are of n-type with lightly doped channel (Boron,
16 3
1 10 cm
) and
heavily doped source/drain regions (Arsenic,
20 3
1 10 cm