practical method for modeling conductor surface roughness ......practical method for modeling...
TRANSCRIPT
Practical Method for Modeling Conductor
Surface Roughness Using Close Packing
of Equal Spheres
Bert Simonovich Lamsim Enterprises Inc. [email protected]
2
Acknowledgements
David Dunham, Molex Inc.
Scott McMorrow, Teraspeed Consulting. -A Division of Samtec
Dr. Yuriy Shlepnev, Simberian Software Corp.
Dr. Alexandre Guterman
3
Outline
Overview
Conductor loss
Copper foil fabrication
Modeling roughness
Hexagonal Close-packing of Equal Spheres Model
Practical method to determine rough conductor loss
Case study
With just -3.5dB delta @12.5 GHz => ½ the eye height with rough copper
4
Failure To Model Roughness Can Ruin You Day
Rough Copper Smooth Copper
Simulated with Keysight ADS
25Gb/s
Current distribution at DC is uniform through cross-sectional area of conductor
5
DC Current Distribution Through a Rectangular Conductor
t
w
DC resistance is proportional resistivity and inversely proportional to the cross sectional area
6
DC Resistance
t
w
_ / mDC condRt w
ρ = Bulk resistivity of the material in Ω-m
X-sectional Area
Above ~10MHz current flows mainly along “skin” of the conductor
7
Skin Effect
δ
Skin depth (δ) is effective thickness where AC current flows
8
Skin Depth
9
Skin Depth vs Frequency
μ0 = Permeability of free space in H/m σ = Conductivity in S/m.
0
1
f
δ
Skin depth inversely proportional to 𝑓
10
AC Resistance
δ
w
_ / m2
AC condRw
Effective X-sectional Area
Reduced cross-sectional area causes AC resistance to increase proportional to 𝑓
High frequency currents concentrated mostly along surface facing reference plane due to proximity effect
11
Current Distribution Microstrip
Reference Plane
12
Current Distribution Stripline
Reference Plane
Symmetrical
Reference Plane
Reference Plane
Asymmetrical
Reference Plane
Return current on respective reference plane ≈ +/-3H from signal conductor center
13
Return Current Distribution
H
3H 3H
H2
H1
3H1
3H2 3H2
3H1
RAC_microstrip = RAC_trace + RAC_ref
14
AC Resistance Microstrip
H
w
t
2
0_ 2
1 1 40.94 0.132 0.0062 ln /m ; when 0.5 10
5.8 0.03AC microstrip
wf w w w wHR f
w Hw H H t H
H w
Reference [1]
1. Determine RAC_microstrip1
2. Determine RAC_microstrip2
3. Combine both in parallel
4. Determine Insertion Loss:
15
Stripline Conductor Loss
_ 1 _ 2
_ 1 _ 2
_ /mAC microstrip AC microstrip
AC microstrip AC
AC str
microstrip
ipline
R f RR
f
R Rf
f f
_
10( ) 20log dB/m( )
AC stripline
smooth
R fIL f e
Zo f
H2
H1
w2
w1
t
16
Conductor Loss Model Validation
4 in
Excellent correlation!
_
10( ) 20log dB/m( )
AC stripline
smooth
R fIL f e
Zo f
VS
Keysight ADS Momentum
No such thing as a perfectly smooth PCB conductor surface
Roughness is always applied to promote adhesion to the dielectric material
17
Conductor Roughness
18
Copper Foil Manufacturing Processes
VS
Rolled Electro-deposited (ED)
Lower Cost Smoother
Copper bar fed through a series of progressively smaller rollers to achieve final thickness
Roller smoothness determines final smoothness of foil
19
Rolled Copper Foil Fabrication Process
Copper Foil
Reducing Rollers
Copper Bar
Drum speed controls foil thickness
Matte side always rougher than drum side
20
Electrodeposited Copper Foil Fabrication Process
Cathode
Drum Side
Matte Side
Anode
Drum
Untreated Foil
Cu Sulfate
21
Common Roughness Profiles
IPC Standard Profile IPC Very Low Profile(VLP) Ultra Low Profile (ULP)Class
SEM Photos Courtesy [6]
No min/max spec < 5.2 μm max -Other names: HVLP, VSP -No IPC spec -Typically < 2 μm max
22
Electro-deposited Copper Foil Nodulation Treatment
Drum Side Untreated
Matte Side Untreated
Drum Side Treated
Matte Side Treated
Untreated Foil
Nodulation Treatment
Treated Foil
Drum Side
Matte Side Matte Side
Drum Side
SEM Photos Courtesy [2]
Profilometers are often used to measure surface roughness
23
Measuring Surface Roughness
Stylus Tip
Surface Profile
Stylus Sample
Motor Drive and Processing Unit
Measured Profile
Stylus Pickup Unit
3D Scan Profile
Faster
More reliable
More accurate
24
Optical Profilometer
Sample
Measured Profile
Optics
CCD
Processing Unit
Average roughness (Ra) typically specified for drum side on data sheet
Ra = Arithmetic average of the absolute values of deviations Yi
25
Average Roughness Parameter
N
a i
i
R YN
1
1
Ten-point mean roughness (Rz) typically specified for matte side on data sheet
Rz= Sum of the average of the five highest peaks and the five lowest valleys over the sample length
26
Ten-point Mean Roughness Parameter
z Pi Vi
i i
R Y Y
5 5
1 1
1 1
5 5
27
Modeling Copper Roughness
Assumes 2D corrugated surface
Based on mathematical fit to S.P. Morgan Power Loss Data (1948)
Lose accuracy above 5GHz for rough copper
28
Hamerstad & Jenson Model
22
1 arctan 1.4rough
HJ
flat
PK
P
RMS tooth height in meters
SF = scaling factor representing the ratio of the length of the rough surface (Lrough )to the spatial length (Lspatial) –Ref [2]
Impractical from first principles perspective – Lrough not published in data sheets
29
Modified Hamerstad & Jenson Model
2
21 arctan 1.4 1
rough
HJM
flat
PK SF
P
Lspatial
Lrough
𝐴𝑚𝑎𝑡𝑡𝑒
𝐴𝑓𝑙𝑎𝑡 = relative area of the matte base compared to a flat surface
ai = radius of the copper sphere (snowball) of the ith size, in meters
𝑁𝑖
𝐴𝑓𝑙𝑎𝑡 = number of copper spheres of the ith size per unit flat area in
sq. meters
δ (f) = skin-depth, as a function of frequency, in meters
30
Huray “snowball” Model
SEM Photo Courtesy [3]
1
2 2
21
43 ( ) ( )1
2 2
jrough matte i i
SRH
iflat flat flat i i
P A N a f fK f
P A A a a
Reference [2]
31
Huray Model Prior Art
11 spheres min; 38 spheres max of radius 1μm to fit within hex tile area and height of 5.8μm
Fit equation parameters to measured data
Assumes stacked “snowballs” arranged in hexagonal lattice
SEM Photo Courtesy [3] Plot Courtesy [2]
VNA Measurement
Model
Benefits:
Practical
Accurate
Issues:
Expertise Required
Time
Money
32
Design Feedback Method*
*Reference [2]
Design Product Channel Simulation
Dk
Df
Roughness
Extract Parameters
Fit
Model Test Design Coupon
X-section Data
Fab
Hexagonal Close-packing of Equal Spheres (HCPES) Model
33
34
Why Bother?
Helps make informed decision sooner - “Sometimes an OK answer NOW! is more important than a good answer late” – Eric Bogatin
Fast simulation time - Practical for what-if spreadsheet analysis
Minimal expertise required
Useful to sanitize CAD tools
Useful to gain intuition on what to expect with measurements and help determine root cause of differences
Similar to Huray Model
Based on close-packing of 11 equal sized spheres
Does not require SEM analysis to determine stack height (HRMS ) or hexagonal tile area
35
HCPES Model
Reference [5]
Hexagonal Tile Area
Assumes nodule treatment applied
to perfect flat surface
Sphere radius and hex tile area
determined solely on published
roughness parameters from
manufacturer’s data sheet
36
HCPES Correction Factor
2
2
2
1 66( ) ( )
12
rough Hex
HCPES
flat
r
P AK f
P f f
r r
HCPES lattice structure loosely resembles the actual SEM photo
37
HCPES Model Lattice Structure
SEM Photo Courtesy [3]
Lattice structure scales inversely to the square of the height
38
HCPES Model Scalability
H H/2
Total of 11 equal sized spheres
HRMS = height of 2 tetrahedrons plus 2 sphere radii
Hexagonal tile perimeter surrounds 7 base spheres exactly
39
HCPES Model Anatomy
Hexagonal Tile
Determine Height of Single Tetrahedron
40
H
41
1. Determine DE
Given: Each side of the tetrahedron = 2r Using Pythagorean theorem:
2
3DE DF
2 2
2 2
2
3
22
3
23
3
DE DB BF
r r
r
A
D
B
F
C E
42
2. Determine Height (AE)
Therefore:
2 2
22 2
2 33
26
3
AE AD DE
r r
r
A
D
B
F
C E
Since HRMS = height of 2 tetrahedrons + sphere dia.
Therefore sphere radius is:
43
Determine HCPES Sphere Radius
2 2
22 6 1
3
RMSH AE r
r
22 6 1
3
RMSHr
r
r
A
E
HRMS
AHex = 6 x area of equilateral triangle ADG
44
Determine Hexagonal Tile Area
2
2
6
1 16 1 3 1
3 3
16 3 1
3
HexA DF AF
r r
r
Hexagonal Tile Area
45
Method to Determine Rough Conductor Loss
Typically:
Must consider roughness of each side when determining AC resistance
Matte sides bonded to core
Drum sides bonded to prepreg
Drum sides roughened with oxide or etch treatment prior to lamination
46
Stripline Geometry with Surface Roughness Example
47
Dual Triangular Sawtooth Profile (DTSP) Model
Used to approximate RMS height of matte and drum side
2 3
aRMSa
RH
2 3
zRMSz
RH
Matte Side Rz
Drum Side Ra
PCB Conductor DTSP
Mean
Mean
Not to Scale
Ra
Rz
48
1. Determine RMS Tooth Height of Matte and Drum Sides
_2 3
zRMS matte
RH
_2 3
aRMS drum
RH
*Use Micro-etch Roughness for Drum Side
*
49
2. Determine HCPES Matte & Drum Correction Factors
_
22 6 1
3
RMS matte
matte
Hr
2
2
_
16 3 1
3Hex matte matteA r
_
22 6 1
3
RMS drum
drum
Hr
2
2
_
16 3 1
3Hex drum drumA r
2
_
_2
2
1 66( ) ( )
12
matte
Hex matte
HCPES matte
matte matte
r
AK f
f f
r r
2
_
_2
2
1 66( ) ( )
12
drum
Hex drum
HCPES drum
drum drum
r
AK f
f f
r r
50
3. Determine AC Resistances of Each Surface
Reference [1]
1
2
0 1 1 1_ 2
1 1 1
41 1ln 0.94 0.132 0.0062 /mAC w
f w w wR f
w t H H
2
2
0 2 2 2_ 2
2 2 2
41 1ln 0.94 0.132 0.0062 /mAC w
f w w wR f
w t H H
1
0 1_ _1
1 1 1
1 1
/m
5.8 0.03
AC ret
w
f HR f
w w H
H w
2
0 2_ _ 2
2 2 2
2 2
/m
5.8 0.03
AC ret
w
f HR f
w w H
H w
* when 0.5 10n
n
w
H
*
*
51
4. Determine Stripline Rough Conductor Loss
2 2_ _ _ _ /mAC drum HCPES drum AC w AC retR f K f R f R f
1 1_ _ _ _ /mAC matte HCPES matte AC w AC retR f K f R f R f
_ _
_ _
_ _
/mAC matte AC drum
AC stripline rough
AC matte AC drum
R f R fR f
R f R f
_ _
10( ) 20log dB/m( )
AC stripline rough
rough
o
R fIL f e
Z f
Case Study
52
53
Test Platform
Case 1 Megtron-6
HVLP Cu
Case 2 N4000-13EP
VLP Cu
Photos courtesy [2]
12 Layer test boards designed, built and tested by Molex Inc., courtesy of David Dunham
Generalized Modal S-parameters (GMS) data courtesy Scott McMorrow, Teraspeed Consulting Group
Generalized Modal S-parameters (GMS) data courtesy Yuriy Shlepnev, Simberian Software Corp
Parameter Case 1 Megtron-6 Case 2 N4000-13EP
Dk 3.62 @50GHz 3.6-3.7 @10GHz[i]
Df 0.006 @ 50GHz 0.008-0.009 @ 10GHz[ii]
Rz HVLP 1.50 μm -
Rz VLP - 2.50 μm
Ra w/Micro-etch [iii] 1.44 μm 1.44 μm
Trace Thickness, t 15.23 μm 15.23 μm
Trace Widths w1, w2 251 μm, 236 μm 251 μm, 236 μm
Dielectric Heights, H1,H2 249 μm, 231 μm 249 μm, 231 μm
GMS trace length 10.15 cm (4.00 in) 10.15 cm (4.00 in)
Zo(fo) ohms [iv] 52.29 @ 50GHz 52.07 @ 10GHz
[i] Dk = 3.65 used [ii] Df = 0.0085 used [iii] CO-BRA BOND® SM is an example of a hydrogen peroxide/sulfuric acid micro-etch treatment often used by PCB
fabricators to improve the adhesion of copper surface to dielectric materials. [iv] Zo(fo) = Characteristic impedance determined by 2D field solver at frequency fo
54
Board Parameters From Data Sheets and Design
Stripline Geometry Reference
Determining Total Insertion Loss
55
_ Total diel cond roughIL f IL f IL f
Keysight ADS
• Svensson/Djordjevic wideband Debye model used to ensure causality for dielectric loss
• Conductivity parameter set to a value much-much greater than the normal conductivity of copper ensures the conductor is lossless for the simulation
_ _
10( ) 20log dB/m( )
AC stripline rough
rough
o
R fIL f e
Z f
Simulation Correlation Results
56
Excellent correlation!
VLP HVLP
Model Comparisons
57
2
2
2
1 66( ) ( )
12
rough Hex
HCPES
flat
r
P AK f
P f f
r r
HCPES HJ
22
1 arctan 1.4rough
HJ
flat
PK
P
HJM
2
21 arctan 1.4 1
rough
HJM
flat
PK SF
P
Tuned SF=1.65
Correction Factor Comparisons (Ra = 1.44μm; Rz =2.5 μm)
58 KHJM SF = 1.65 KHJM SF = 1.65
59
Summary and Conclusions
1. Using the concept of hexagonal close-packing of equal spheres, a novel method to accurately calculate sphere size and hexagonal tile area was devised for use in the Huray model.
2. By using published roughness parameters and dielectric properties from manufacturers’ data sheets, we show the need for further SEM analysis or experimental curve fitting, may no longer be required for preliminary design and analysis.
3. HCPES model looks promising as a practical alternative to previous modeling methods.
60
Ongoing Research
Test the HCPES model to see how well this method applies to other material and copper roughness
61
References [1] S. Hall, H. Heck, “Advanced Signal Integrity for High-Speed Digital Design”, John Wiley & Sons, Inc., Hoboken, NJ, USA., 2009.
[2] E. Bogatin, D. DeGroot, P. G. Huray, Y. Shlepnev, “Which one is better? Comparing Options to Describe Frequency Dependent Losses”, DesignCon2013 Proceedings, Santa Clara, CA, 2013.
[3] Huray, P.G.; Hall, S.; Pytel, S.; Oluwafemi, F.; Mellitz, R.; Hua, D.; Peng Ye, "Fundamentals of a 3-D “snowball” model for surface roughness power losses," Signal Propagation on Interconnects, 2007. SPI 2007. IEEE Workshop on , vol., no., pp.121,124, 13-16 May 2007 doi: 10.1109/SPI.2007.4512227.
[4] Huray, P. G. (2009) “The Foundations of Signal Integrity”, Chapter 6, John Wiley & Sons, Inc., Hoboken, NJ, USA., 2009
[5] L. Simonovich, “Method for Modeling Conductor Surface Roughness”, Provisional Patent Application # US 61/977,138
[6] ISOLA-Group, “Copper Foil 102 Presentation”, 2012