Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Книги+1 / 2013 [Chandan_Kumar_Sarkar]_Technology_CAD

.pdf
Скачиваний:
176
Добавлен:
11.03.2016
Размер:
34.78 Mб
Скачать

MOSFET Characterization for VLSI Circuit Simulation

283

7.3.4.1 Simulation Setup

For all the simulation results provided in this chapter, conventional bulk NMOS transistor has been selected. A 45-nm technology node has been selected, and the drawn channel length of the transistor is taken to be 65 nm, if not mentioned otherwise. The supply voltage is taken to be 1 V. The physical oxide thickness is 1.1 nm and the electrical oxide thickness is 1.75 nm, considering poly-depletion effect and inversion layer thickness. The substrate is uniformly doped with concentration equal to 3.24E18/cm3. The source/drain concentration is 2E20/cm3. The source/drain junction depth is 14 nm. The HSPICE simulation tool has been used to obtain all the simulation results, with BSIM4 as the compact model. The model parameters of the corresponding predictive technology model [4] have been taken for simulation purposes.

7.3.4.2  Threshold Voltage Characterization with Substrate Bias Effect

The variation of threshold voltage of a large geometry MOS transistor (W = L = 10 m) is shown in Figure 7.6(a). It is observed that the threshold voltage increases from its zero substrate bias value as the substrate voltage is increased. The value of zero substrate bias long geometry threshold voltage VT 0 as observed from Figure 7.6(a) is 0.466 V. The threshold voltage variation with substrate bias follows (7.5). Noting the values of the necessary model parameters from the model file, the theoretical curve is drawn and compared with the simulation results. It is observed that the theoretical curve closely follows the simulation results.

The variation of the substrate sensitivity of threshold voltage with the substrate bias is shown in Figure 7.6(b). The sensitivity obtained from theoretical formulation is also shown in Figure 7.6(b).

7.3.4.3  Threshold Voltage Characterization for Short Channel Transistors

The variation of threshold voltage with channel lengths of a MOS transistor is shown in Figure 7.7(a). It is observed that as the channel length is reduced from 100 nm onwards, the threshold voltage reduces. This is referred to as the threshold voltage roll-off. The effect is more pronounced when the applied drain bias is high. This is demonstrated in Figure 7.7(b). The amount of threshold voltage roll-off as observed from the simulation results for different substrate bias and drain bias are summarized in Table 7.1. It is observed that with the increase in substrate bias, the short channel effect increases. This is easy to understand considering the fact that with the increase in substrate bias, the depletion width increases and consequently the short channel effect increases as suggested in Brew’s relation.

The DIBL effect is shown in Figure 7.8. The DIBL coefficient η is defined as

the slope of the curve (η = VT (VDS =1V)VT (VDS =0.1V)). From simulation results, its

0.9

284

Technology Computer Aided Design: Simulation for VLSI MOSFET

 

0.66

VDS = 50 mV

eoretical

 

0.64

 

VGS = 1 V

Simulation

 

0.62

 

 

 

0.60

 

 

VT (V)

/dVBS

0.58

 

 

 

 

 

0.56

 

 

 

 

 

0.54

 

 

 

 

 

0.52

 

 

 

 

 

0.50

 

 

 

 

 

0.48

 

 

 

 

 

0.46

 

 

 

 

 

0.0

–0.2

–0.4

–0.6

–0.8

–1.0

 

 

 

VBS (V)

 

 

 

 

 

(a)

 

 

–0.14

VDS = 50 mV

 

 

 

 

 

 

eoretical

 

 

–0.15

VGS = 1 V

 

Simulation

 

 

 

 

 

 

 

–0.16

 

 

 

 

 

–0.17

 

 

 

 

 

dVT

–0.18

–0.19

–0.20

0.0

–0.2

–0.4

–0.6

–0.8

–1.0

 

 

 

VBS (V)

 

 

(b)

FIGURE 7.6

(a) Variation of threshold voltage with substrate bias for n-channel MOS transistor. (b) Substrate sensitivity of threshold voltage.

MOSFET Characterization for VLSI Circuit Simulation

285

VT (mV)

 

 

VDS = 50 mV

VBS = 0 V

 

 

VBS = –0.5 V

660

 

 

 

 

 

 

VBS = –1 V

640

 

 

 

 

 

 

 

 

620

 

 

 

 

 

600

 

 

 

 

 

580

 

 

 

 

 

560

 

 

 

 

 

540

 

 

 

 

 

520

 

 

 

 

 

500

 

 

 

 

 

480

 

 

 

 

 

460

–0.2

–0.4

–0.6

–0.8

–1.0

0.0

L (µm)

(a)

VT (mV)

470

 

 

VBS = 0 V

 

 

VDS = 1 V

 

 

 

 

VDS = 50 mV

 

 

 

 

 

 

 

465

 

 

 

 

 

 

VBS = 0V

 

 

 

 

 

 

 

 

 

 

 

 

 

460

 

 

 

 

 

 

 

 

 

 

 

 

 

 

455

 

 

 

 

 

 

 

 

 

 

 

 

 

 

450

 

 

 

 

 

 

 

 

 

 

 

 

 

 

445

 

 

 

 

 

 

 

 

 

 

 

 

 

 

440

 

 

 

 

 

 

 

 

 

 

 

 

 

 

435

 

 

 

 

 

 

 

 

 

 

 

 

 

 

430

 

 

 

 

 

0.8

1.0

 

 

 

 

 

0.0

0.2

0.4

0.6

 

 

 

 

 

L (µm)

 

 

(b)

FIGURE 7.7

(a) Threshold voltage roll-off for low drain bias. (b) Simulation results showing that the threshold voltage roll-off increases with increased drain bias.

286

Technology Computer Aided Design: Simulation for VLSI MOSFET

TABLE 7.1

Amount of Threshold Voltage Roll-Off at Low and High Drain Bias for Different Substrate Bias

 

 

VDS = 50 mV

 

 

VDS = 1V

Roll-Off

VBS = 0 V

VBS = –0.5 V

VBS = –1 V

 

VBS = 0 V

 

 

 

 

 

 

VT(mV)

1.55

1.70

1.86

32.16

 

 

 

 

 

 

value is found to be 0.0342, and the theoretical value as calculated from the model discussed earlier is 0.0346.

7.3.4.4  Threshold Voltage Extraction

The threshold voltage extraction method through extrapolation in the linear region method is shown in Figure 7.9 and that through the second derivative method is shown in Figure 7.10. The threshold voltage as extracted from the constant current method, linear extrapolation method, and second derivative method for low drain bias and different substrate bias are summarized

VGS = 1 V

 

640

 

620

 

600

 

580

(mV)

560

540

T

520

V

 

500

 

480

 

460

440

420

0.0

–0.2

–0.4

–0.6

VDS (V)

FIGURE 7.8

Simulation results illustrating the DIBL effect.

L = 65 nm (VBS = 0 V)

L = 65 nm (VBS = –0.5 V) L = 65 nm (VBS = –1 V) L = 1 um (VBS = 0 V)

–0.8 –1.0

MOSFET Characterization for VLSI Circuit Simulation

287

ID (A)

2.0 m

VDS = 50 mV

 

VBS = –1 V

 

 

 

 

VBS = –0.5 V

 

 

1.8 m

 

 

 

 

 

 

VBS = 0 V

 

 

1.6 m

 

 

 

 

 

 

 

 

 

1.4 m

 

 

 

 

 

1.2 m

 

 

 

 

 

1.0 m

 

 

 

 

 

800.0 µ

 

 

 

 

 

600.0 µ

 

 

 

 

 

400.0 µ

 

 

 

 

 

200.0 µ

 

 

 

 

 

0.0

 

 

 

 

 

–200.0 µ

 

 

 

0.8

1.0

0.0

0.2

0.4

0.6

VGS (V)

FIGURE 7.9

Threshold voltage extraction: extrapolation in the linear region method.

dgm/dvGS (s/v)

 

 

 

 

VBS = –1 V

 

 

VDS = 50 mV

VBS = –0.5 V

25.0 m

 

 

 

VBS = 0 V

 

20.0 m

 

 

 

 

 

15.0 m

 

 

 

 

 

10.0 m

 

 

 

 

 

5.0 m

 

 

 

 

 

0.0

 

 

 

 

 

–5.0 m

 

 

 

 

 

–10.0 m

 

 

 

 

 

0.0

0.2

0.4

0.6

0.8

1.0

 

 

 

VGS (V)

 

 

FIGURE 7.10

Threshold voltage extraction: second derivative method, evaluated at VDS = 50 mV.

288

Technology Computer Aided Design: Simulation for VLSI MOSFET

TABLE 7.2

Threshold Voltage Value Extracted through Different Methods

 

Extracted Value of the Threshold Voltage at

 

VDS = 50 mV. W = 10 μm L = 65 nm

 

Method

VBS = 0 V

VBS = –0.5 V

VBS = –1 V

 

 

 

 

 

Constant current

0.457 V

0.559 V

0.645

V

Linear extrapolation

0.409 V

0.502 V

0.583

V

Second derivative

0.462

0.558

0.636

 

 

 

 

 

 

in Table 7.2. For the constant current method, the constant current is taken to 1 μA, the channel width is 10 μm and the channel length is 65 nm.

7.4 I-V Characterization

Precise knowledge of the I-V characteristics of a MOS transistor is a basic requirement for a good VLSI designer. The fundamental current transport equations are introduced, followed by channel charge, mobility, and velocity saturation effects. The I-V models for long and short channel devices are derived, followed by some advanced issues.

7.4.1  Current Density Equations

The total current density is the sum of the drift current density and the diffusion current density, written as

Jn = qnµnξ + qDn

dn

(7.18a)

 

dx

 

Jp = qpµpξ − qDp

dp

(7.18b)

dx

 

 

The total conduction current density is thus J = Jn + Jp. The diffusion coefficients Dn and Dp for electrons and holes are related to the corresponding mobilities n and p through Einstein’s relationship [5]. Thus the current densities are written as follows:

Jn = qnµnξ + kTµn

dn

(7.18c)

 

dx

 

Jp = qpµpξ − kTµp

dp

(7.18d)

dx

 

 

MOSFET Characterization for VLSI Circuit Simulation

289

The electric field ξ, which is defined as the electrostatic force per unit charge, is written as ξ = −dψi/dx. It may be noted that gradual channel approximation has been assumed, according to which the variation of the electric field in the y-direction (along the channel) is much less than that in the x-direction (perpendicular to the channel). With this the conduction current densities are written as

Jn = −qnµn dφn

dx

Jp = −qpµp dφp

dx

The quasi-Fermi potentials φn and φp are defined as [5]

φn ≡ ψi

kT

 

n

q

ln

 

 

 

 

 

ni

φp ≡ ψi +

kT

 

p

q

ln

 

 

 

 

 

 

ni

The electron current density at a point (x,y) in the channel is

Jn (x, y) = −qµnn(x, y) dVCS (y)

dy

(7.19a)

(7.19b)

(7.20a)

(7.20b)

(7.21)

Here VCS(y) is the quasi-Fermi potential. The total current at any point y along the channel is

xi

dVCS

 

 

IDS (y) = qW nn(x, y)

dx

(7.22)

dy

0

 

 

 

The integration is carried out from x = 0 to x = xi, the bottom of the inversion layer where ψ = ΦF . There is a sign change as the drain current flows in the negative y direction. The inversion charge density is defined as

Qinv (y) = −qxi

n(x, y)dx

(7.23)

0

 

 

With this the drain to source current is given as

 

IDS = µn WL VDS [Qinv (V)] dVCS

(7.24)

0

 

 

290

Technology Computer Aided Design: Simulation for VLSI MOSFET

 

G

 

Inversion layer

xt

tox

 

– – – – – –

y

 

n+ S

n+ D

 

xb

x

 

Wdm

 

 

B

Depletion layer

FIGURE 7.11

Inversion layer forms one capacitor with the gate and another capacitor with the body. Surface mobility is a function of the average electric fields at the top and the bottom of the inversion charge layers.

The quasi-Fermi potential VCS (y) = 0 at y = 0(source) and VCS (y) = VDS at y = L(drain).

7.4.2  Channel Inversion Charge Density

With reference to Figure 7.11, the following assumption is made. The inversion layer in the channel of a MOS transistor is a sheet of charge and there is no potential drop or band bending across the inversion layer. This is referred to as the charge sheet approximation [5]. This inversion layer forms a capacitor with the gate, the oxide being the dielectric. Also it forms another capacitor with the body, the depletion layer being the dielectric. Thus the inversion layer is coupled with both gate and substrate of the transistor. The inversion charge density in strong inversion is given by

Qinv = −[Cox (VGS VCS (y) VT 0 ) + Cdm (VBS VCS (y))]

(7.25)

Simplification of (7.25) leads to the following expression for inversion charge density:

 

Qinv = −Cox (VGS mVCS (y) VT )

(7.26)

where

 

 

 

 

 

m 1 + α = 1 +

Cdm = 1 +

3tox

 

(7.27)

Wdm

 

 

 

Cox

 

In (7.26), VT = VT 0

Cdm

VBS = VT 0 − αVBS = VT 0 (m 1)VBS

has been taken

Cox

using (7.4b).

 

 

 

 

MOSFET Characterization for VLSI Circuit Simulation

291

As defined earlier, m is referred to as the body-effect coefficient.* The value of m is typically 1.2, however, it can be taken to be unity for simplified calculations. Clearly the bulk charge factor is closely related to the body-effect parameter, as observed also from (7.4b).

The inversion charge density in the weak inversion region is given by [7]

Qinv

qεSi NCH

 

ψs 2ΦF VCS

 

 

UT exp

 

 

(7.28)

4ΦF

UT

 

 

 

 

In (7.28), UT = kT/q is the thermal voltage, and NCH is the effective channel concentration.

7.4.3  Carrier Mobility Degradation Model

In the inversion layer of a MOS transistor, the current flow is determined by the surface mobility, whose value is much lower than the bulk mobility of the carriers. This is because of several mechanisms of scattering, primarily the phonon or the lattice scattering, the coulombic scattering, and the surface roughness scattering [5]. For good quality interfaces, phonon scattering is the dominant mechanism at room temperature.

The surface mobility is a function of the average of the electric fields at the top and the bottom of the inversion charge layer. These fields are shown in Figure 7.11.

From Gauss’s law using the depletion layer as the Gaussian box, it is possible to write

ξxb =

QB =

Cox (VT VFB 2ΦF )

(7.29a)

εSi

 

εSi

 

Considering the Gaussian box to be a box that encloses both the depletion and the inversion layer, we have

ξ = QB + Qinv

xt εSi

From (7.29a) and (7.26), we have, after substitutions,

ξxt = Cox (VGS VFB 2ΦF )

εSi

(7.29b)

(7.29c)

It is to be noted that the effect of the lateral field is ignored and m = 1 for simplicity.

* Some authors refer to this as the bulk-charge factor.

292

Technology Computer Aided Design: Simulation for VLSI MOSFET

The average electric field is defined as

 

 

ξeff =

1

(ξxb + ξxt )

(7.30)

 

 

 

2

 

 

Substituting from (7.29a) and (7.29c), and after some simplifications for n+ poly-gate n-channel MOS transistor

ξeff =

VGS + VT + 0.2V

(7.31)

6tox

 

 

Physically, ξeff means the average electric field experienced by the carriers in the inversion layer. The dependence of the surface mobility of the carriers on this average electric field and hence on the gate bias is given by the following empirical relationship [3]:

s =

0

 

 

1 + (

ξeff

)

υ

(7.32)

 

 

ξ0

 

 

In (7.32), µ0 is the low field surface mobility, ν is a constant whose value is ~1.85 for electrons at the surface, and ~1.0 for holes at the surface. ξ0 is the critical electric field (~0.9 MV/cm for electrons at surface and ~0.45 MV/cm for holes at the surface). The model proposed in (7.32) fits experimental data well, but because it involves a power function, it is difficult to integrate the model in a circuit simulation program. A Taylor series expansion of (7.32) and retaining only up to three terms, the following expression for the vertical field mobility degradation model is derived:

s =

0

(7.33)

1 + UA (VGStox+VT )+ UB (VGStox+VT )2

In (7.33), UA and UB are two parameters, whose values are to be extracted from the experimental I-V data. The substrate bias dependence of the mobility is incorporated by introducing another parameter UC in (7.33). With this, the model becomes [6]

s =

 

 

0

 

 

 

 

 

 

(7.34)

1 + (UA + UC .VBS )(

VGS +VT

)+ UB (

VGS +VT

)2

 

 

 

 

tox

tox

 

µs =

 

 

 

0

 

 

 

 

 

 

 

1 + UA (

VGS +VT

)+ UB (

VGS +VT

)

2

 

 

 

(7.35)

 

 

 

 

 

(1

+ UCVBS )

 

tox

tox