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

Книги2 / 1988 Kit Man Cham, Soo-Young Oh, John L. Moll, Keunmyung Lee, Paul Vande Voorde, Daeje Chin (auth.) Computer-Aided Design and VLSI Device Development 1988

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

350

Computer-Aided Design

311.

3D

10

(fF)

 

 

....

 

 

 

 

 

......

 

 

 

 

 

....

....

"--

 

 

 

2

-.- - _-

 

 

 

 

---------

 

 

 

 

 

...

0.'

0.71

1.00

1.•

1••

 

 

VERTI CAL. SPAC I NG

 

 

 

 

(a)

 

 

311

 

 

3D

 

10

1. (fF)

.....................................................

 

-..-....

----------------

 

 

O.

0.'

 

------------

0.71

1.00

1.•

1

 

 

 

 

VERTICAL. SPACING

 

 

 

 

(b)

 

 

 

311

 

3D

 

10

1. (fF)

 

 

 

 

 

........................

 

 

 

 

_--

.............................

 

..........

 

 

 

 

 

-------------------------

0._--:;:......_--..;__-:-';;;;--_--;-'.;-_---,.....---1

 

0.'

0.71

1.00

1..

1.•

 

 

VERTICAL. SPACING

 

 

 

 

 

(c)

 

 

Fig. 15.17. Capacitance calculation for

the

three cases in Fig.

 

15.16 by FCAP3.

 

 

Examples of Parasitics Simulation

351

complete calculation is done three-dimensionally by FCAP3.

 

Fig. 15.17(a) shows the interline capacitance for the

parallel lines in

Fig. 15.16(a). Since the lines are paralle~ the two-dimensional calculation is an accurate simulation and is the same as the three-dimensional calculation. As the vertical spacing becomes larger, the capacitance decreases for both the two-dimensional and one-dimensional calculations. However, the relative error of one-dimensional calculation to the two-dimensional calculation becomes larger since the effect of the fringing field increases. For the 45° crossing lines in Fig. 15.16(b), we have to do the three-dimensional simulation to obtain a correct capacitance. As shown in Fig. 15.17(b), the twodimensional calculation underestimates the capacitance given by the threedimensional calculation, and the one-dimensional calculation underestimates it even more. Since the 9<Y' crossing has more three-dimensional characteristics than the 45° crossing, two-dimensional and one-dimensional calculations become almost useless. For the spacing of 1.5 ~m, the result in Fig. 15.17(c) shows that the two-dimensional calculation estimates the capacitance to be only half of the correct capacitance calculated by the three-dimensional calculation.

The above examples emphasize the importance of the three-dimensional simulation which most current circuit extractors do not consider. When they extract a circuit from its layout including interconnect capacitances, the extractors usually calculate them using a two-dimensional approximation, namely, areal and peripheral coefficients. But the smaller the size of interconnects and the greater the number of interconnect levels, the more inaccurate this approach is.

Fig. 15.18 shows another example of capacitance calculation of two-level interconnect lines. This time, three different dielectric materials are used to be more realistic. Due to the symmetry, half of the whole problem is simulated by cutting with a :xz-plane along the center of the second level metal (note the triad in the figure for the directions). Table 15.2 lists the amounts of charge on each conductor with a bias condition. Also, the potential plots in three ways - bird's-eye-view plot, contour plot, and one-dimensional plot - are shown in Fig. 15.19(a), (b), and (c). The bird's-eye-view and contour plots are the potential distribution on the plane, z = 2.4~, and the one-dimensional plot is

352

Computer-Aided Design

5.& 7.0 &.~

l.~

Fig. 15.18. Two-level on-chip metal lines simulated by FCAP3.

Conductor Potential

Charge (fC)

ground

OV

-1.469

Cll

OV

-0.338

C12

IV

3.023

C13

OV

-0.448

C21

OV

-0.768

Table 15.2. Capacitance values of the structure in Fig. 15.18 calcu-

lated by FCAP3.

 

 

the potential plot along the line, x

= 4.1 JJm and z = 2.4 JJm. The plots show

the high potential along the first-level center conductor (M12) and the suppressed potential due to the second-level conductor (M21) near the xzplane.

It is difficult to benchmark the above calculations against real measurements, since the values are too small to be measured directly. In fact, this is one of the main reasons that we have to resort to simulations to estimate the values. To verify the accuracy of three-dimensional calculations by FCAP3, we can build a scaled up structure and measure the capacitance of it. Fig. 15.20 shows a structure with metal bars and Table 15.3 lists the measured capacitance values for different configurations. The capacitance values simulated by FCAP3 and percentage difference between the measured and simulated values

Examples of Parasitics Simulation

353

x

y

o

L-__-W____ __

(b)

u

o'--------------------------

'y

(c)

Fig. 15.19. (a) Bird's-eye-view plot and (b) contour plot of potential at the plane z=2.4J,'m of the structure in Fig. 15.18 calculated by FCAP3. (c) Potential plot along the line x=4.1J,'m z=2.4J,'m.

354

I ground

Computer-Aided Design

ground

T2 = 0.5 inch

Fig. 15.20. Side views of the structure used to benchmark the FCAP3 calculations.

HI

Sl

H~

S~

 

Co

Cl•

Cll

Cl~

Co-Cl•

 

(inches)

 

Measured

 

(pF/l1inches)

 

-Cll-C12

0.30

1.00

1.00

1.00

19.94

13.13

2.50

4.06

0.25

 

 

 

 

FCAP3

19.40

12.97

2.35

4.08

0.00

 

 

 

 

DiJJ(%)

-2.71

-1.22

-6.00

0.49

 

1.03

1.00

0.25

1.00

Measured

18.81

4.21

2.88

11.38

0.34

 

 

 

 

FCAP3

18.88

4.25

2.44

12.19

0.00

 

 

 

 

DiJJ(%)

0.37

0.95

-15.28

7.12

 

0.37

1.00

1.00

4.50

Measured

17.77

11.66

3.17

2.29

0.65

 

 

 

 

FCAP3

16.84

11.02

3.25

2.57

0.00

 

 

 

 

DitI(%)

-5.23

-5.49

2.52

12.23

 

1.27

5.00

0.60

4.50

Measured

10.74

5.78

0.17

4.39

0.40

 

 

 

 

FCAP3

10.51

5.28

0.22

5.01

0.00

 

 

 

 

DitI(%)

-2.14

-8.65

29.41

14.12

 

1.27

5.00

0.60

10.00

Measured

9.58

6.38

0.44

2.45

0.31

 

 

 

 

FCAP3

9.35

5.87

0.60

2.88

0.00

 

 

 

 

DiJJ(%)

-2.40

-7.99

36.36

17.55

 

1.27

10.00

0.60

10.00

Measured

9.44

6.64

0.09

2.50

0.21

 

 

 

 

FCAP3

9.35

6.18

0.17

3.00

0.00

 

 

 

 

DitI(%)

-0.95

-6.93

88.89

20.00

 

Table 15.3. Comparison of measured and simulated capacitances of the structure in Fig. 15.20.

are also listed in Table 15.3. The conservation of charge requires C11 +C12 +C19-C0 to be zero, where C11, C12, C19, and Co are the capacitances between the first-level center metal and the two adjacent first-level metals, the capacitance between the first-level center metal and the three second-level metals, the capacitance between the first-level center metal and the ground, and the total capacitance of the first-level center metal, respectively. As shown in Table 15.3, this is consistently satisfied by all the values from FCAP3. On

Examples of Parasitics Simulation

355

the contrary, the values from the experiment have some offsets. We can interpret the offsets are due to a measurement error. This measurement error causes large percentage difference between the simulation and measurement for the small capacitance values. However, if the measurement error is considered, Table 15.3 shows an excellent agreement between the measurement and the simulation by FCAP3.

Via Hole Capacitance Calculation by FCAP3

The capacitance between a conductor plane and a line which goes through a hole in the plane is important to calculate accurately in modeling of the integrated circuit packages and printed circuit boards. We can easily simulate the structure using the second level input scheme of FCAP3. We need only simulate one-eighth of the structure, taking advantage of the symmetry of the structure as shown in Fig. 15.21(a). Fig. 15.21(b) shows the calculated capacitance values with "respect to the line diameter. The two-dimensional approximation is obtained by multiplying the thickness of the plane and the

1.4Ca acitance (fF)

1.2

j

y III

1-x

10.O\Jm

r::t::.V""

6.O\Jm

~

.4

 

 

 

 

10.0\lIl

\I~o.S\.Ull 0

.2

2-DIMENSIONAL CALCULATION

 

_.....-_

 

 

 

- .......... -.....--......------.....---

 

 

 

f

. 1

1.5

2

Z5

3

 

 

 

 

Radius(micro-meter)

 

 

 

(a)

 

 

 

(b)

 

 

Flg.1S.l1. (a) A via hole and a line through the via. (b) Capacitance calculations of the via structure v.s. the radius of the line by FCAP3 and two-dimensional approximation.

356

Computer-Aided Design

capacitance per unit length of infinitely long cylinders. It underestimates the capacitance, especially when the diameter of the line is small, because it does not include the fringing field.

Via Resistance Calculation by FCAP3

2.6

2.6~

2.3

2.45

1.0

Fig. 15.22. A via connection of two-level on-chip metal lines. The shaded areas are contacts to the metal layers for the simulation

d (pm)

Resistance (0)

0.0 (plug)

0.1237

1.10.1292

1.350.1306

1.60.1320

1.850.1329

Table 15.4. Resistance of the via connection in Fig. 15.22 by FCAP3.

Examples of Parasitics Simulation

357

FCAP3 can calculate the resistance of a three-dimensional structure using the duality between capacitance and resistance as described in section 4.3. With this method, the resistance of a via connect in Fig. 15.22 is calculated varying the depth of the void at the center. The shape of the void depends on the multi-level metal process used. The result in Table 15.4 shows that the resistance does not change significantly. This is reasonable because the void is located where the current density would be low even if the void were filled [15.3].

References

[15.1] S. Seki and H. Hasegawa, "Analysis of Crosstalk in Very High-Speed LSI/VLSl's Using a Coupled Multiconductor MIS Microstrip Line Model," IEEE Trans. Electron Devices, ED-31, no. 12, pp.1948-1953, Dec. 1984.

[15.2] H. Hasegawa and S. Seki, "Analysis of Interconnection Delay on Very High-Speed LSI/VLSI Chips Using an MIS Microstrip Line Model," IEEE Trans. Electron Devices, ED-31, no. 12, pp.1954-1960, Dec. 1984.

[15.3] K. Lee, P. Vande Voorde, M. Varon, and Y. Nishi, "A Method to Reduce the Peak Current Density in a Via," Tech. Digest ofIEDM 1987.

Appendix

Source Information of 2-D programs

PROGRAM SUPREM

SUPRA

SOAP

SOURCE INFORMATION OF 2-D PROGRAMS SOURCE ADDRESS

Stanford Univ.

Office of tech. licensing, Stanford Univ.

 

105 Enclna Hall, Stanford, CA 94305

.

·

·

Public

I ,

GEMINI

CADDET

Hitachi

PISCES

Stanford Univ.

FCAP2

Hewlett Packard

TECAP2

.

HPSPICE

 

·

Tech. Administration, Hitachi

Technical

P.O. Box 2. KokubunJl, Tokyo, Japan

exchange

Office of Tech. licensing, StanfordUniv.

Public

105 Enclna Hall, Stanford, CA 94305

 

Hewlett Packrad Lab.

Technical

3500 Deer Creek Rd., Palo Alto, CA 94304

exchange

·

Commercial

 

 

Commercial

Table of Symbols

B

B

BfA

BVDS

C

C

Co

Co

Co

Clo C2

CA

CD

CD

Cj

Cn

Cn

C~

sensitivity matrix element angstrom unit (lE-lO m)

parabolic oxide growth-rate constant empirical parameter in mobility model linear oxide growth-rate constant source-drain breakdown voltage

capacitance

impurity concentration

oxidant concentration at the oxide interface

capacitance of a transmission line with dielectric materials removed output loading capacitance

doping concentrations areal capacitance coefficient depletion capacitance

drain junction capacitance

oxidant concentration at the silicon interface coefficient of auger recombination

n - pocket characteristic length gate oxide capacitance