Humidity Testing Correlation in Peck Comprehensive Model
Humidity Testing Correlation in Peck Comprehensive Model
D. Stewart Peck
-
3646 Highland St., Allentown, PA 18104
(215) 398-1222
The standard test for the effect of humidity has EVALUATION PROGRAM AND RESULTS
been, for many years, the electrically-biased test
of components in 85c temperature and 85% relative A survey was made of all published median-life
humidity. At the initial use of this test, it was data on electrolytic corrosion failures of aluminum
considered a design-and processing~control type of metalization (with the inclusion of one batch of data
test. Specification writers and device manufacturers on nichrome fuse electrolysis) in epoxy packages, from
felt that, if they could meet the usual requirement of available sources[l-14J, avoiding that data which did
10% failures in 1000 hours for this test, the devices not provide median life measurement at 85/85. A few
would operate as reliably as others in general-use exceptions were made where the 85/85 median could be
conditions. Attempts were not made to determine a estimated from nearby conditions with a minimum of
relationship between 85C/85%RH (85/85) test results and probable error. The 85/85 life being the base, the
expectation of life at some other environmental ratio Ro of observed median at some other condition to
condition. that at 85/85 was determined. As lives both longer and
shorter than at 85/85 were involved, this ratio
Since about 1970 we have seen the results of provides a consistent set of numbers which do not con-
various laboratories' efforts to understand such a flict. Ro is the Acceleration Factor (AF) from 85/85
relationship for plastic-encapsulated semiconductor to lower stress condition, but l/Ro would be the AF
devices. These efforts are frequently based on a small from higher-stress to 85/85; to show AF for both
number of tests at only a few conditions of humidity conditions on one plot would be confusing.
and temperature, which is not surprising considering
the long testing time for any condition much less Acceleration and Deceleration
severe than 85/85. A mathematical model is then
postulated which seems to fit these points. Serious The IEEE Transactions on Reliability, in Nov 1970,
thought should tell us, of course, that there must be published definitions which include:
variation around each of the measured median lives, in- "Acceleration Factor: The ratio between the times
cluding sampling variation, changes in the encap- necessary to obtain a stated proportion of failures for
sulation properties from time to time, contamination, two different sets of stress conditions involving the
and many other factors. Hence, a model which perfectly same failure modes and/or lIlechanisms."
fits a few data points may still not be correct. So
there is a continuing effort to find a model which Note that this definition does not distinguish between
CH2256-6/86/0000-0044$OI.00 c 1986 IEEE!IRPS 44
(t./t2) and (t2/t.), but calls either one an 100 .,.-------.------,----..----"""7\
-
acceleration factor. In this author's opinion, this is I
a:
a good practice. If, in going from an observed life at
low stress to a shorter life at high stress by a
factor, e.g., of 175, then 175 is the acceleration u 10 l = - - - - - 4 - - - - - - l - - - - - . A r r - - - - - . j
factor(AF). One divides the long life by 175 to obtain
the life at accelerated stress. If one then has an
.."'
~
observed life at the high stress condition and wants to ~
g~
predict life at low stress, it is undesirable to have
1.0
to remember that the Deceleration Factor (DF) is
0.0057142; it is only required to remember that the AF OJ
>
is 175, and one multiplies by it rather than dividing ~
oJ
by the OF. Only a single AF between two conditions is OJ
a:
required if one uses it appropriately. w
u.
0.1
:::;
Z
Also defined is: c(
o'"
item ... (identifying a common failure mechanism) ... " '" I
DECELERATION FACTOR• •ACCELERATION FACTOR
This implies that the life at the higher stress is the .00 1 '-:-...L.J...J.J...U.':'"--.J.---l...J...U.l.ill_..l.-LLW-L"':----L...J-JL.llWJ
.01 0.1 1.0 10 100
accelerated life, and is normally obtained by dividing CALCULATED MEDIAN LIFE
the long life by the AF. In the present review, the RELATIVE TO THI\T 1'.1" 85 C I 85% RH
ratio of median life at some other condition to that at
85/855 is used, giving an AF when treating lower
stresses, and a DF when treating higher stresses. The Figure 1 Ratios of life in reported humidity tests to
total AF between a stress below 85/85 and that above that in 85C/85%RH are shown, comparing the
85/85 would be the AF below 85/85 divided by the DF observed ratio Ro to the calculated Rc from
from the higher stress. the ~odel. The ratios for various test con-
ditions are shown for perspective.
Modeling and Model Results
Linearity of the regression was gauged by eye from
The modeling involved in this progr~ is to pos- Figure 1; distributions of HO/Rc are shown later.
tulate a relationship between life-and-temperature and Linear regression analysis provides the intercept on
life-and-humidity so that the product of the two the Ro axis at 0.917 + or - 0.064, and the slope as
separable factors between a given condition and 85/85 0.994 + or - 0.04, the tolerances being determined from
provides a calculated ratio Rc, according to the model, the standard deviation of error. Removal of one
which best .atches the observed ratio, Ro. A plot of maverick point changes these numbers to an intercept of
Ho va Rc is used to determine regression parameters and 0.982 and a slope of 1.026, with the s~e tolerances as
the linearity and uniformity of the regression. before.
The lognormal distribution of life was assumed, as There was also some SusplClon that the earliest
it was in all the data reviewed. The standard data may not have represented the best-controlled test
deviation estimated in all tests ranged from 0.4 to conditions, and that later data might give better
0.5, showing remarkable uniformity. Ho and Hc also results. For data published in 1979 or later, the
have lognormal characteristics, as does the ratio combination of E.=0.81eV and n=-2.66 has a correlation
Ro/Rc. The use of the life ratio allows the consistent coefficient of 0.988, with intercept and slope of 1.013
comparison of data from 5-volt and 70-vo1t operation, and 1.012, and with the best standard deviations of
from test structures and from commercial devices, and distribution above and below 85/85 (0.175 and 0.281
from different time periods. respectively). Optimization was not explored, for the
sample size for Rc>l.O is then only seven and more data
Figure 1 shows the plot of 61 data points, with Ro should be obtained before this is studied further.
on the vertical axis and Rc on the horizontal, the 45° There is the potential for reducing the confidence
line representing a 1:1 correlation. The relationship intervals for extrapolation to low-humidity conditions.
shown here is:
It is apparent, however, even from existing data,
time-to-failure {RH)n exp (Eajkt) that a single relationship can be used to relate
stresses both above and below 85/85, making viable a
where the exponent n in Figure 1 is -2.66, and Ea direct extrapolation from autoclave test results to
O.7geV. The resulting acceleration calculation is life expectations at low-stress conditions. Further,
relatively insensitive in the range of n from -2.5 to publication of more data at autoclave stresses should
-3.0 and of Ea from 0.77 to 0.8leV, using the be strongly encouraged, in order to define the model
regression parameters and the distributions of RotRc more precisely for modern test facilities, and to
above and below 85/85 (_edians and estimates of determine a measure of repeatability. The lines at the
standard deviation) as criteria. Thus there are seven lower left of Figure 1 show the OF obtainable at some
measures of fit, and not all are optimized at anyone COnditions. The acceleration of 130/85 over 85/85 is
choice of n and Ea. The following combinations are almost 20:1, and 100%RH increases that. It is
suggested as nearly equivalent (and optimuu) in interesting to note that, for all the fears about
results: humidity testing at or near lOOe or 100%RH, the data
show only marginal influence of these limits; 16 of the
Ea n correlation coefficient 61 data points are within 5% of these limits, and all
are within the range of the other data, with the
0.79 -2.66 0.986 exception of one point at 100%RH, from 1972, at the
0.81 -2.50 0.985 lowest relative value of the distribution in Figure 2
0.77 -3.00 0.987 below.
45
Uniformity of Model Match resulting distribution has an improved median of 0.956
and an s of 0.281. Deletion of these points does not
-
An even division of Ro above and below Hc reflects change the correlation coefficient of the model, but
a uniform match of the data with the model through the improves the Ro intercept from 0.982 to 0.999, and the
entire range of the data; a visual review of Figure 1 slope from 1.026 to 1.020. This suggests that, while
shows this. Figures 2 and 3 show the distributions of 100%RH is useful for reducing the time of routine
Ro/Rc for Rc<l.O and Rc>l.O respectively. The medians tests, it may be desirable to avoid this condition if
of these distributions are at 0.919 and 0.963, both one wants to optimize extrapolation to low stress
being satisfactorily close to 1.0 A further breakdown levels. Further, testing in modern equipments should
to smaller segments of the range of data would show be aimed at determining if this fine distinction is
significant at the present time.
10 High-Temperature Limitations
Rol Rc FOR
R c <1.0 Several factors may be involved in determining a
temperature limitation for testing. One of them can be
n= 42
disposed of quickly, for it may not be a real problem;
med=0.92
s= 0.295
investigators should be aware of it, in case it shows
up. An anomalous result appeared in testing of
"., silicone-molded devices with Ti-Pd-Au metalization at
"
0::
~ temperatures in the range of 150-2500 and at 80%
o
0:: 1.0
,/ saturation [15).
0"
~
0::
.
./ "" 400
/' 350
300
~ 250
<t
(.)
fJ>
... 200
{ ([
~~
.10
10 30 50 10 90 98 99.8
0.2
CUMULATIVE PERCENTAGE
'<w"t 150
z
~ I"'----
::::i """..:::
u
Figure 2 The distribution of deviations (Ro/Rc) from
the model for stress greater than (life ::! 100 .......
::>
shorter than) at 85/85.
10
!;;:
'"
W
11.
::.
w
.... 50
" ""'
"
0:: ~ Figure 4 Examples of possible anomalous life be-
o / ~ havior at high temperature in 80% sat-
a: 1.0 uration (as seen in silicone encapsu-
0"
~ lations) .
0:: ,.
This is shown in three examples in Figure 4 as a
• reversal in the temperature effect, giving better life
at 2500 than at 2000. This may not be a problem in
present testing because of the difference in tempera-
ture capabilities of the materials. If, however,
improvements allow testing beyond 1500, this effect
.10
0.2 10 30 50 10 90 98 99.8
might possibly appear. It was not investigated at the
CUMULATIVE PERCENTAGE
time (1969), so there is no background for judging its
possible presence in other materials.
Figure 3 The distribution of deviations (Ho/Hc) from There are other evidences of notable behavior at
the model for stress less than (life great- 1500. Figure 5 shows an addition to figure 1, adding
er than) at 85185. data for 1500 results for several epoxy devices[15]. At
80% saturation these results, together with those at
distributions of similar uniformity. The data from 200e on silicone encapsulations, show both quite longer
stresses higher than 85/85 have a better distribution and quite shorter life than would be predicted by the
and a smaller estimate, s, of the standard deviation. present model which appears so effective up to that
The confidence intervals (90%) for the standard stress. These data were obtained in the late 1960's
deviation above 85/85 are 0.25-0.36, and are 0.26-0.46 and may simply reflect the lack of understanding at
for lower stress. These are not greatly different, but that time of proper testing procedures in autoclaves.
suggest that the distribution at higher temperature may One might guess, however, that anything done "wrong"
be somewhat better. This is not surprising, might tend to shorten life rather than lengthen it.
considering the shorter test times, with less Therefore, the regime of testing at 1500 and above, for
opportunity for test condition perturbation. epoxies, is one which should be re-examined when it can
In the data of figure 2, six results represent be done safely with respect to glass transition
lOO%RH conditions; if these results are eliminated, the temperature.
46
----------------~
r-
~
1 I '0
• 0°0 D
~
...:r The D.E.T. curves as shown by Koelmans[17) show
r l50Cl80'll. /0 0
that the drastic reduction in surface conductivity (vs
......> ~---<>-
0 1.0
- ; - ------l-
zoo,eo ., I f /f. ~ I RIl) takes place at about 5%; at this condition, at 150C
the life should be 36 times that at 85/85--an
~
., .,
" ~ ~t:XX0!'4-100C/8$% untestable time. Even at 20OC, the AF at 5% is 3.66,
W
II: 0.1
ZOO/80 U- --J..-. still a poor test time. At 20", however, the ratio to
w
u.
::::; -]~- ...~ -- llo,Cf8""'--- life at 85/85 is 0.91 at lSOC, 50 there is some
.
z
ffi .01 --
.4
i%
-~-~--
wl( If ,...- 120Cl'~%
Condition Acceleration
100
0.1 0.5 2 10 30 50 70 90 98 99.599.9
121/100 16
CUMULATIVE PERCENT DEFECTIVE
135/94 30
140/94 40
140/100 50
Figure 6 Composite life distribution from 85/85 150/100 77
and 121/100 tests.
The equivalent of 1000 hours at 85/85 could be 20
Establishing Humidity Test Conditions hours at 140/100, allowing this test to be completed
within one day. This would only be a measurement of
Figure 6 is an exampl~ of the combined test the freak level in the product, but could very well be
results of both 85/85 and 121/100 conditions (on the a specification requirement, since it will control the
basis of the equivalent time at the 85/85 conditions) early failure rates. In a relatively short time (5 to
from many samples from the same production period of a 20 days) one could describe the median life of any
semiconductor product. During the 1000 hours at 85/85, present epoxy system, or of one which might meet
only the freak level of the product was measured (but Military requirements.
the same result could have been obtained in 48 hours at
121/100, or in 28 hours at 135/100). Further, the CONCLUSIONS
important information that the product had a median
life of about 4000 hours at 85/85 was completely First, there is available here a relationship
unavailable except by testing up to about 3000 hours at between temperature, humidity and life for electrolytic
85/85, 200 hours at 121/100, or 84 hours at 135/100. corrosion of aluminum metalization, which can be used
to extrapolate autoclave testing results, accounting
Which test is the more desirable? The shortest for the possible effect of the epoxy Tg, directly to
test provides the quickest feedback for control, the expected life at low stress down to about 30%RH, below
greatest throughput of samples, and the greatest incen- which life may be longer than indicated by the model.
tive for frequent testing. Avoiding testing at 100%RH may give a slight
48
------------~_=s!
improvement in the confidence limits of the extrapola- lation to the low-humidity-atress region.
tion, from past data; the possible effect with modern
equipMents should be explored further. The present model provides the tightest distri-
bution of fit to all the data used.
Second, this relationship allows the establishment
of very-short-time tests to replace the present 1000- 10
hour, 85/85 tests, and this should be done i~ediate~
More test results should be reported at both below and
above 85/85 stress. Low-stress tests are needed to 'l
improve the statistics in that area, and the high- r- Re >1.0
~"3
stress tests should emphasize measures of repeat-
ability, particularly at 100%RH. There appears to be ~
~yl
no other reason not to use this condition for routine
V~
testing, except for the reproducibility of ratio to u
ex:
V
85/85, compared to that of lower RH. ;p 1.0
o
Third, the existence of epoxy packaging systems ~
for semiconductor devices with median lives greater a: ./ -'
than 20,000 hours has been suggested. Such devices '/
could have acceptable performance in Military applica- W
tions, and the autoclave results discussed here /'
indicate that they could be qualified in reasonable
tiDe periods. Formulation of humidity- test require-
ments for these applications may now be timely. .1
0.2 10 30 50 70 90 98 99.8
APPENDIX A CUMULATIVE PERCENT
A paper by O. Hallberg[18] also reviewed publish- Figure 8 Distributions of deviations from the model
ed data on the ratios of median life in other for the models of Figure 7, for stresses
conditions to median life in 85/85, and plotted 51 less than 85/85.
points of data according to the (RH)2 model of
Lawson(20,21] and the (T+RH) model of Reich and 10
Hakim(22]. This work was effective in reducing the
number of models, but did not use distributions of the
observed ratios about the calculated ratios as a I--- t- R c < 1.0
measure of fit, which is felt here to be the better "3
method of cOlllparison, together with the regression par'- Vl~
ooneters.
l:?"2
~~
a:u
REGRESSION PARAMETERS o ~ ./
log Ro = a + b (log Rc)
a;: 1.0
0"
~ '7
CORRELATION a b
II:
,......-........ l./V
".
COEFFICIENT (Ro INTERCEPT) (SLOPE) ./ V
(ACTUAL Ro) 1/ V./
REICH & HAKIM
1972 .974 .780 .911 .P'
PECK & ZIERDT
1973 .956 .670 1.18
1
0.2 10 30 50 70 90 98 99.8
LAWSON
1974·1984 .950 .775 1.09 CUMULATIVE PEIICENT
PECK
1986 .986 .982 1.03 Figure 9 Distributions of deviations from the model
for the models of Figure 7, for stresses
greater than 85/85.
Figure 7 Regression par~eters for four humidity
test models.
Figure 7 shows the regression parameters for the 2. J. L. Flood, "Reliability Aspects of Plastic En-
61 points of data collected. Figures 8 and 9 show the capsulated Integrated Circuits", Int'l ReI Phys
distributions of the deviations of the observed data Symp, 1972, pp 95-99
from each model, as numbered above.
3. M.C. Halleck, "The Ie Plastic Package: A Simple
Since all the models except 4 evolve from limited Method for Predicting Package PerforDIance" ,
data, it is not surprising that they show regions where Int'l ReI Phys Symp, 1972, pp 88-94
the model is relatively far from observed fact,
although none are off by extreme ratios. For example, 4. D.S. Peck and C.H. Zierdt, Jr., "Temperature-
the Reich and Haki. model evolved largely from the Humidity Acceleration of Metal-Electrolysis
Panama tests in high humidity, and the model fits well Failure in Semiconductor Devices", Int'l ReI
in that area, but does not provide as good an extrapo- Phys Symp, 1973, pp 146-152
49
5. N. Lycoudes, "The Reliability of Plastic Micro-
-
circuits in Moist Environments", Solid State
Technology, 1978 Oct