ISPD 2001 ISPD 2001
- K. Banerjee, M. Pedram and A. H. Ajami
Analysis and Optimization of Analysis and Optimization of Thermal - - PowerPoint PPT Presentation
Analysis and Optimization of Analysis and Optimization of Thermal Issues in High Thermal Issues in High Performance VLSI Performance VLSI Kaustav Banerjee Stanford University Massoud Pedram and Amir H. Ajami University of Southern
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
n O Die
and Cooling Technologies (Same Rn) TDie at Other Technology Nodes Calculated
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
∆ ∆ ∆ ∆T increases with increasing tox Self Heating under DC Stress (IRPS 96)
Thermal impedance θ θ θ θj, defined by ∆ ∆ ∆ ∆T = P x θ θ θ θj
100 200 300 400 500 600 1 2 3 4
Power [W] ∆ ∆ ∆ ∆T [0C]
M1 M3 M2 M4 M1: Metal1 M2: Metal2 M3: Metal3 M4: Metal4
increasing tox
Thickness of AlCu in M4 is doubled
L W= = = =1000 µ µ µ µ µ µ µ µ m 3 m TiN Silicon
θ θ θ θj
tox
AlCu
Oxide
Cross section
ISPD 2001 ISPD 2001
Impact of Scaling Using Low-k (IEDM 96)
Metal
SiO2 SiO2 Low-K (Gap Fill)
Metal
DC Conditions
As W decreases SH increases. Low-k increases SH by 10-15%
100 200 300 400 500 1 2 3
Power [W] ∆ ∆ ∆ ∆T [
0C]
Metal 1 Standard Dielectric 0.75 µm 1.5 µm 3 µm
Standard Dielectric
100 200 300 400 500 1 2 3
Power [W] ∆ ∆ ∆ ∆T [
0C]
Metal 1 Low-k Dielectric 0.75 µm 1.5 µm 3 µm
Low-K Dielectric
ISPD 2001 ISPD 2001
Chip Power and Area increases Negligible Change in Power Density Current Density in Metal Lines Increases Number of Metal Levels Increases Electromigration (EM) Time to Failure Decreases Increased ρ ρ ρ ρ (T) Wire Delay Increases
ISPD 2001 ISPD 2001
1.0 1.5 2.0 2.5 3.0 3.5 4.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2
1.0 1.5 2.0 2.5 3.0 3.5 4.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2
( ITRS ’99 ) Air HSQ SiO2 Polyimide
180 nm 100 nm 100-130 nm 130 nm 70 nm <50 nm
ISPD 2001 ISPD 2001
Conduction
Generation (q’’’), Constant Properties (k)
2
( Interconnect )
2
( Others )
– Uniform jrms for all Metal Lines ( ITRS ’99 )
Si Substrate Si Substrate ILD ILD Cu IMD ILD ILD Interconnect Interconnect
TDie
ISPD 2001 ISPD 2001
Change in Power Density
( ITRS ’99 )
Due to Joule Heating of Interconnects ( FEM Simulation )
50 100 150 200 250 180 130 100 70 50 35
Technology Node [nm] Temperature [°C]
0.0 0.2 0.4 0.6 0.8 1.0
Power Density [W/mm2]
Tmax TDie
50 100 150 200 250 50 100 150 200 250 180 130 100 70 50 35 180 130 100 70 50 35
Technology Node [nm] Temperature [°C]
0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0
Power Density [W/mm2]
Tmax TDie
(IEDM 2000)
ISPD 2001 ISPD 2001
Temperature [°C]
100 nm 130 nm
180 nm
35 nm 70 nm
1 2 3 4 5 6 7 8 9 10 1 2 3 4 5 6 7 8 9 10 1 2 3 4 5 6 7 8 9 10 120 140 160 180 200 220 120 140 160 180 200 220
Distance from Substrate [µ µ µ µm]
50 nm
209 °C 126 °C
Global Wires Global Wires
(IEDM 2000)
ISPD 2001 ISPD 2001
20 40 60 80 100 70 180 130 100 50 35
% Decrease in TTF Technology Node [nm]
20 40 60 80 100 20 40 60 80 100 70 180 130 100 50 35 70 180 130 100 50 35
% Decrease in TTF Technology Node [nm]
100 ) ( ) ( 1 %
max
× × × ×
− − − = = = =
Die
T TTF T TTF TTF in D
10 20 30 40 50 180 130 100 70 50 35
% Increase in RC Delay Technology Node [nm]
10 20 30 40 50 10 20 30 40 50 180 130 100 70 50 35 180 130 100 70 50 35
% Increase in RC Delay Technology Node [nm]
100 1 ) ( ) ( %
max
× × × ×
− − − = = = =
Die
T T Delay RC in I ρ ρ ρ ρ ρ ρ ρ ρ
(IEDM 2000)
ISPD 2001 ISPD 2001
Scaling trends that cause increasing thermal
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
Transport of mass in metal interconnects under an applied current density EM lifetime reliability modeled using Black’s equation given by,
− − − − m B n
TTF is the time-to-fail A is a constant that depends on line geometry and microstructure j is the DC or average current density Q is the activation energy for EM ( ~ 0.7 eV for AlCu) Tm is the metal temperature
ISPD 2001 ISPD 2001
Accelerated EM stress data yields A, Q, and n in Black’s equation, and a value of log-normal σ σ σ σLN Typical goal: achieve a 10 year lifetime EM stress data + Black’s equation gives a technology limit to the maximum allowed current density ( javg ) for the required failure rate and a desired lifetime at a reference temperature Tref ( ~ 100 0C) The javg limit does not comprehend self heating
ISPD 2001 ISPD 2001
Peak, Average, and RMS current densities: For an unipolar waveform: EM is determined by javg, and self-heating by jrms
A I j
peak peak =
= = =
( ( ( ( ) ) ) )dt
t j T j
T avg
= = = 1
( ( ( ( ) ) ) )
= = =
T rms
dt t j T j
2
1
A is the cross sectional area of interconnect, T is the time period of the current waveform.
T ton Ipeak Irms Iavg = I
DC
r =
ton/T
peak rms
j r j = = = =
peak avg
j r j = = = =
r is the duty factor
ISPD 2001 ISPD 2001
= = =
− − − − m B n
T k Q exp j A TTF
θ θ θ θ
R R I T T T
rms ref m heating self 2
= = = = − − − − = = = = ∆ ∆ ∆ ∆
− − − −
Rθ
θ θ θ is the effective thermal impedance
given by,
eff ins ins
W L K t R = = = =
θ θ θ θ
Weff is the effective metal width to account for quasi-2D heat conduction.
L
Silicon tins
Metal
insulator Wm Kins Tm Tref q tm
ISPD 2001 ISPD 2001
Typically, design rules specify javg from EM and jrms from self-heating separately. Self-consistent approach: comprehends EM and self-heating simultaneously. The lifetime at any javg and metal temperature Tm, should be equal to or greater than the lifetime value (e.g., 10 year) under the design rule current density (j0).
2 ref B 2 avg m B
j T k Q exp j T k Q exp
≥ ≥ ≥
ISPD 2001 ISPD 2001
r j j
2 rms 2 avg =
= = =
( ( ( ( ) ) ) ) ( ( ( ( ) ) ) )
eff ins ref m m m m m ins ref B m B 2
W K T T T W t t T k Q exp T k Q exp j r − − − −
= = = ρ ρ ρ ρ
This is a single equation in the single unknown temperature Tm
ISPD 2001 ISPD 2001
1.0E+05 1.0E+06 1.0E+07 1.0E+08 1.0E+09 0.0001 0.001 0.01 0.1 1
Duty Cycle r Maximum jrms [A/cm2]
1.0E+05 1.0E+06 1.0E+07 1.0E+08 1.0E+09
Maximum jpeak [A/cm2]
Kins X 10-6 W/m0K 1.15 (SiO2) 0.60 (HSQ) 0.25 (Polyimide k) 0.024 (Air)
materials are introduced.
ISPD 2001 ISPD 2001
ineffective in increasing jpeak.
100 200 300 400 500 600 700 800 0.0001 0.001 0.01 0.1 1
Duty Cycle r Self-Consistent Metal Temperature Tm [0C]
1.0E+04 1.0E+05 1.0E+06 1.0E+07 1.0E+08 1.0E+09
Maximum jpeak [A/cm2] Air j0 X 105 A/cm2
6.0 9.0 13.5 20.3
SiO2 SiO2
ISPD 2001 ISPD 2001
Comparison with AlCu
1 2 3 4 5 6 7 Oxide HSQ Polyimide Air
Dielectric Maximum Allowed j peak (MA/cm
2)
Cu: javg = 1.8 MA/cm2 Cu: javg = 0.6 MA/cm2 AlCu: javg = 0.6 MA/cm2
0.1-µm/Metal 8 r = 0.1
materials are introduced
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
Thermal effects predominant in semi-global and global interconnects which are: Away from the Si substrate Long Typically split into buffered segments Long interconnects can be optimally buffered.
sopt s lopt
ISPD 2001 ISPD 2001
distributed nature of the interconnect.
0.25 µ µ µ µm and 0.1 µ µ µ µm technology Full 3-D Interconnect capacitance extracted Accurate lopt and sopt values determined by SPICE simulations
jrms (max)
sopt sopt lopt
ISPD 2001 ISPD 2001
Effect of Thermal coupling Included (NTRS Based)
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 Oxide HSQ Polyimide Air
Dielectric javg [x105 A/cm 2]
Performance Reliability Reliability-Multilevel
0.1-µm/Metal 8
even after considering thermal coupling
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
Electrostatic Discharge (ESD) A short duration (< 200 ns), high current (> 1 A) event Can cause open circuit failure of metals and latent damage that impact EM reliability
ISPD 2001 ISPD 2001
Failure current densities are much higher than under normal circuit conditions
1 2 3 4 2 4 6 8
Current Density, J [x107 A/cm2] Resistance Rise, ∆Rmax/R0
200 400 600 800 1000 1200
Temperature Rise, ∆Tmax [0C]
Metal 1 Metal 2 Metal 3 Metal 4
∆t = 200 ns
γcrit
Melt Threshold
Self-heating characteristics of AlCu lines under short- pulse stress conditions (Electron Device Letters 97)
identical SH
due to smaller surface area to volume ratio
temperature is ~ 1000 0C
ISPD 2001 ISPD 2001
Metal 4 ~ 12 µ µ µ µm Metal 1 ~ 12 µ µ µ µm
ISPD 2001 ISPD 2001
Significant Electromigration Performance Degradation
1 µ µ µ µm 0.15 µ µ µ µm Unstressed AlCu Stressed AlCu
ISPD 2001 ISPD 2001
Thermal Analysis including Interconnect Joule Heating based on ITRS ’99
Peak Temperatures in ICs Increase with Technology Scaling in Spite of Constant Power Density Significant Implications for Performance and Reliability Advanced Chip Cooling Techniques may be Necessary
Thermal Effects and Reliability
Thermal Effects Strongly Impacts EM Self-Consistent Analysis: Thermal + EM Point-to-Point Interconnects Optimized for Performance Meets Reliability Based Current Density Limits High-Current Design Rules Must be Followed for I/O and ESD Protection Circuit Interconnects
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
L tox tm
Interconnect Oxide
Via Via
Substrate
(Tref )
(Tline )
Heat equation in Interconnect (DAC 2001)
2 2 line m
2 2 2 2
line line ref
λ λ λ λ and θ θ θ θ are constants
ISPD 2001 ISPD 2001
200 400 600 800 1000 1200 1400 1600 1800 2000
108 109 110 111 112 113 114 115 116 117
Position x (micron)
Thermal profile for different technology nodes
T(x) (c)
d
d
Tref = 100 °C
T(x=0) = 100 °C T(x=2000) = 100 °C
L=2000 µ µ µ µm
ISPD 2001 ISPD 2001
DPM, Functional block clock gating Thermal time constant is much higher that signal propagation constant
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
Introduces error in wire-planning and
ISPD 2001 ISPD 2001
Unit length capacitance is not affected
ISPD 2001 ISPD 2001
∆ ∆ ∆ ∆x rd CL L
L L L d L L x
L L L
° ° °C
ISPD 2001 ISPD 2001
µ µ µm): 5-6% increase for each 20-degree increase in long global lines
5 10 15 20 25 30 35 40 45 50 30 50 70 90 110 130 150 Temperature ( C ) Delay Increase (%) L=100 L=400 L=700 L=1000 L=2000
ISPD 2001 ISPD 2001
200 400 600 800 1000120014001600 18002000 20 40 60 80 100 120 140 160 180 200
Position (X)
T(x) Exponential Thermal Profile TH TL
8 10 12 14 16 18 20 22 24 30 50 70 90 110 130
P e a k Te m pe ra ture ( C )
L=2000,T1 L=1000,T1 L=2000,T2 L=1000,T2
TL=30 °C
Direction of Thermal Gradient is Important
ISPD 2001 ISPD 2001
T(x) T(x) T(x) T(x) x x x x
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
Have long global segments in the highest metal layers delay variations affect skew
ISPD 2001 ISPD 2001
1 2 2 3 3 3 3
ISPD 2001 ISPD 2001
l L-l
p
q x 1 2 2
ISPD 2001 ISPD 2001
*
*
l
ISPD 2001 ISPD 2001
T(x) T(x) T(x) T(x) x x x x L/2 L/2 L/2 L/2
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
9.57 911
µ µ µ µ=300, σ σ σ σ=700
0.0 1000
µ µ µ µ=1000, σ σ σ σ=400
2.40 979.5
TH=170, TL=130
3.63 968.66
TH=170, TL=110
2.65 1021
TH=170, TL=130
3.98 1032
TH=170, TL=110
5.42 1042 TH=170, TL=90 7.78 1210
µ µ µ µ=2000, σ σ σ σ=1000
5.24 957.5
TH=170, TL=90
( ) T x ax b = +
H L
T T a L − =
L
b T = ( )
bx
T x a e− = ⋅
H
a T =
1 ln( )
H L
T b L T =
2 2
( ) 2 m ax
( )
x
T x T e
µ σ − − −
= ⋅
l=L/2 skew%
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
ISPD 2001 ISPD 2001
Delay model for non-uniform line temperature presented Delay based on uniform worst case line temperature is not sufficient Direction of thermal gradients is important
Non-uniform substrate temperature introduces skew in the clock tree Bottom-up merging techniques must consider non- uniform interconnect thermal profile Skew can be minimized by suitable merging