A METHOD FOR DESIGNATION/SORTING SEMICONDUCTORS WAFERS ACCORDING TO PREDICTED OXYGEN PRECIPITATION BEHAVIOUR

Abstract: A method for marking/sorting semiconductor wafers in accordance with predicted values of oxygen precipitation in the wafers that will occur after the performance of a thermal process step. In a preferred embodiment, this oxygen precipitation prediction is obtained by means of a correlation equation Y = A₀ + A₁ OX i + A₂/LT + A₃C + A₁₂ OX i /LT + A₁₃ OX i C + A₂₃ C/LT, where OX i is the value of interstitial oxygen in the wafer measured prior to the thermal process step, LT is the minority carrier lifetime for the wafer measured prior to the thermal process step, and C is the initial substitutional carbon in the wafer measured prior to the thermal process step. The constants A₀, A₁, A₂, A₃, A₁₂, A₂₃, and A₁₃ are initialized on a first sample of wafers representative of a desired range of crystal types. In a preferred embodiment, this initialization is performed using Y values equal to the change in interstitial oxygen after the performance of the thermal process step. This predictive method can be performed on every wafer prior to semiconductor device fabrication therein, thereby improving device yield after wafer processing.

Description: Background of Invention The present invention relates to a method for the physical designation of semiconductor wafers according to predicted oxygen precipitation behavior for a given set of thermal process steps. This method utilizes non-destructive tests on the wafer and thus can be performed on every wafer prior to semiconductor device fabrication therein. In the mass-production of large scale integrated circuit devices, it is highly desirable to detect material deficiencies that may produce either unsatisfactory LSI devices, or low yields of satisfactory devices, prior to the initiation of any fabrication process steps. Such process steps add cost to the individual wafers, and should be minimized for any wafers ultimately destined for scrap. It is now well known that the existence of oxygen in single crystal semiconductor wafers such as, for example, silicon wafers, has a profound effect on ultimate device yields and device characteristics in large scale integrated (LSI) electronic circuits, especially field effect and bipolar transistor digital circuits. The presence of this oxygen in the semiconductor wafers has a beneficial effect in terms of wafer hardening (solution hardening) and impurity gettering. The wafer hardening arises due to oxygen atoms located interstitially in the semiconductor lattice acting to prevent, to some degree, the semiconductor atoms from moving within their lattice. Accordingly, these oxygen atoms sitting interstitially in the lattice act to harden the lattice and prevent wafer warpage and lattice shear. The impurity gettering arises due to oxygen precipitates formed in the lattice. Essentially, these oxygen precipitates microscopically disrupt the semiconductor lattice and form a dense localized network of lattice dislocations. Metallic impurities arising in crystal growth, wafer fabrication, wafer cleaning, and/or during subsequent hot processing steps such as epitaxial deposition (800-1200°C), dopant drive-in (900-1200°C), or chemical vapor deposition (700-1100°C), tend to migrate to these oxygen precipitate sites where they are trapped. These oxygen precipitates are typically formed by oxygen atoms leaving the interstitial solution in the semiconductor lattice to form microscopic SiO₂ particles or Si x O y complexes. They are essential to removing the metallic impurities which are added in the above-listed later device processing steps and which would destroy device operation. The major phenomena affecting oxygen precipitation behavior for a semiconductor wafer of, for example, silicon, undergoing a given defined sequence of thermal steps, are the presence of undesired contaminants, the initial oxygen concentration in the lattice solution, and the grown-in or annealed microstructure of the lattice. Accordingly, during semiconductor device manufacturing via the growing of a monocrystalline boule of suitable doped silicon, much effort has been expended in the past to limit the contaminants in the as-grown boule. Likewise, much effort has been expended to control the growth process to obtain a uniform radial and axial distribution of oxygen compatible with subsequent process and device needs. However, despite improvements in crystal growing technique, there still exist variations in individual wafers sliced from a given boule that may render any given one of the wafers unsuitable for subsequent device processing. Likewise, there are clear differences between individual boules. Such differences manifest themselves due to differences in thermal processing temperatures, heating and cooling rates, furnace configurations, spacing of wafers in boats, non-oxidizing or oxidizing ambients, the intrusion of trace contaminants, etc. These differences in process details can produce profound differences in the resulting wafers. Accordingly, it is highly desirable to be able to test each wafer cut from a boule, before or after polish, to insure that it is suitable for the fabrication of a device thereon, and prior to performing any further device fabrication steps. This testing could then be used for either marking or sorting the wafers in order to minimize the number of cost-adding process steps performed prior to scrapping the wafers. However, at the present time there is no satisfactory method for predicting the oxygen precipitation content that will exist after the performance of a thermal process step. Current manufacturing techniques generally require the dedication for destructive test purposes of a certain number of wafers from a test lot. This wafer dedication not only precludes the use of these same wafers for device fabrication, but also places undue reliance on statistical sampling analysis. Moreover, it has been found that the use of sister wafer testing for crystals is, in many cases, inaccurate due to micro-fluctuations and along accros the crystal in the various atomic concentrations and lattice structure. A variety of studies have been performed in the art directed at oxygen precipitation content in semiconductor wafers. For example, see the study by Huber et al, Satellite Symposium, Electro Chemical Society, ESSDERC 82, Aggregation Phenomena of Point Defects In Silicon, Munich, September 17, 1982, wherein oxygen precipitation engineering in silicon is discussed. However, such studies do not adequately disclose methods of accurately predicting oxygen precipitation content of a wafer after a thermal process step. The invention as claimed is intended to remedy the above-defined drawbacks. The invention solves the problem of how to predict oxygen precipitation content in a semiconductor wafer which exists after a given set of device processing thermal cycles have been performed. This invention further solves the problem of predicting oxygen precipitation content without damaging the wafer under test. The advantages offered by the invention are that oxygen precipitation content in a semiconductor wafer can be predicted prior to the performance of thermal processing steps, thereby preventing the addition of cost-adding process steps to wafers which are to be scrapped. A further advantage to the present invention is that the oxygen precipitation content of the wafer is predicted by means of non-destructive tests performed on the wafer, thereby permitting the testing of every wafer. Thus, it is no longer necessary to dedicate certain wafers from a lot for test purposes, or to rely on statistical sampling analysis for wafer batches. A still further advantage of the present invention is that through the use of this oxygen precipitation prediction method, it is now possible to select wafers in accordance with predicted oxygen precipitation values for a given product/process application, to provide rapid feedback for the optimization of the device process and/or fabrication conditions (including crystal growth and subsequent thermal treatments), and to characterize wafers on an arbitrary standardized "figure of merit" oxygen precipitation scale, to facilitate sorting. These options can be exercised independently of a particular product process. Summary of the Invention Briefly, the present invention comprises a method for the physical designation of semiconductor wafers according to predicted oxygen precipitation behavior, as defined in claim 1. Preferred embodiments are defined in the dependent claims. Brief Description of the Drawings Fig. 1 is a schematic diagram of an apparatus which may be utilized to measure the minority carrier lifetime. Fig. 2 is a graph of the predicted oxygen precipitation vs. measured values of the change in interstitial oxygen ΔOx = Ox i - Ox f . Detailed Description of the Preferred Embodiment One method of carrying out the invention is described in detail below. This method relates to the physical designation of semiconductor wafers according to predicted oxygen precipitation behavior, for a given set of thermal process steps. This method for predicting the oxygen precipitation behavior of a wafer after a thermal process step is based on the proportionality of the oxygen precipitating amount to both the initial interstitial oxygen content realized before the performance of the thermal process steps in question, and to the microscopic structure and content of other impurities in the wafer before the performance of the thermal process steps. The following independent parameters are utilized in the present prediction method: the initial interstitial oxygen value, OX i , for the wafer under test; the minority carrier lifetime, LT; and the initial substitutional carbon value, C, in the wafer. The initial interstitial oxygen value, OX i , provides the reference or starting point for determining the change in interstitial oxygen after the performance of a thermal process step. The reciprocal of the minority carrier lifetime is utilized because it measures the degree of nucleation of the wafer lattice for oxygen precipitation, i.e., it is proportional to the density of microdefects which can act as oxygen precipitation sites. The substitutional carbon, C, value is used because carbon tends to form silicon carbide particles which constitute nucleation sites. Note however, that carbon, cannot be outdiffused from the silicon surface after wafer fabrication (unlike the oxygen), and thus acts to attract contaminants which can destroy the functionality of any devices formed in the wafer surface layer. The method comprises the first step of measuring the initial value OX i of interstitial oxygen in the wafers in a non-destructive manner prior to performing the given set of thermal processing steps. There are a variety of techniques available for taking an initial interstitial oxygen measurement. In this regard, see ASTM F121-83. However, commercially available fast fourier transform infrared (FTIR) instruments are preferred for taking interstitial oxygen measurements on thin wafer samples suitable for LSI product manufacture. The next step is to measure the minority carrier lifetime LT in the wafers in a non-destructive manner which is free of interference from surface recombination due to uncontrolled fast surface states, with the measurement being performed prior to performing the given set of thermal processing steps for device fabrication, but after any anneals to enhance the nucleation for precipitation. In a preferred embodiment, the surface recombination may be turned off via a strong d.c. accumulation bias applied by means of a wet-cell containing an easily cleaned electrolyte connected to low voltage electrodes. The LT measurement may be obtained by means of a photoconductive decay tester which operates by measuring microwave reflectance after pulsed illumination of the wafer by a diode laser. The diode laser in this test configuration may be set behind a quartz plate bearing a slotline antenna and a transparent bias electrode. A detailed description of this apparatus for testing the minority carrier lifetime, as well as a description of the FTIR device, is provided infra. The next step in the present method is to correlate, in a single correlation, the OX i values and the LT values for a calibration sample of wafers to Y values which are proportional to the predicted oxygen precipitation value for the series of wafers. This correlation is then used to predict precipitation on wafers of like nature. There are a variety of methods available in the art for correlating the OX i values and the LT values for a first set of wafers in order to obtain a correlation equation for Y. In one embodiment, a first wafer set of initial oxygen values OX i and reciprocal of minority carrier lifetime values LT are correlated to Y values which, in this example, may be proportional to the change in interstitial oxygen values, Δ0X, measured after the given set of thermal process steps have been performed on this first wafer set. This correlation yields an equation for Y with an OX i term and a LT term. It is preferred that some form of multi-variate regression analysis be utilized to perform this correlation. In general, multivariate regression analysis, of whatever functional form for F(X1, X2, ...Sn), seeks to approximate the fit of Y-data to X-data, by minimizing the errors of the fit. The most common method of optimization is that of "least squares fit", or "mean-square regression", in which the individual error terms, Yi - Fi, for n-sets of data are squared (to remove sign dependence), then summed, with the result minimized through optimization of the coefficients in the function F. This is also called the minimized "root mean square" error, RMS, if the square-root of the sum is taken. The function F(X) can take many forms, including linear, power, exponential, Fourier, transformation of variable, etc., with the form sometimes clearly suggested by the phenomena providing the data. Among the various procedures for implementing multivariate regression analysis (one dependent variable, several independent variables to exponents of one or greater, with optional cross-terms) are those of "the general linear model" (in which multiple equations are expressable in matrix form), "the forward solution of the abbreviated Doolittle method", "the backward elimination procedure", etc. In a preferred embodiment, multivariate regression linear analysis is utilized, both because it is the simplest in form, and also because it is suggested by the wafer physics. For purposes of illustration only, and not by way of limitation, a discussion of the the least squares regression analysis method is provided below. For this method, with a first wafer set of N data pairs, X i and Y i , where Y is the dependent variable and X is the independent variable, a function F(x) can be obtained to approximate Y with a linear fit to the data. In the simplest form, with only a single independent variable OX i , the function F(x) may be equal to F(X) = A₀+A₁X i , with the function F calculated such that the sum of the squares of the error terms [Y(X n ) - F(X n )]² , is minimized. In this equation for F(X), the coefficient A₀ is the system offset, while the slope coefficient, A₁, is the coefficient for the independent variable, X. In order to determine A₀, the following equations are calculated for the series of n data pairs: X =(1/n) Σ X i , Y =(1/n) Σ Y i , where the sum sign Σ denotes summation over the n terms. Accordingly, A₀ = Y - A₁ X , and A₁=(ΣX i Y i -n XY )/(ΣX i ²- X ²) . For this linear regression of Y on X, the scatter between Y and F can be used to calculate a correlation coefficient CC, which has a magnitude approaching unity for a good fit, and approaching 0 for a minimal relationship. For one independent variable, CC²=(n Σ X i Y i -ΣX i ΣY i )²/(n ΣX i ²-(Σ X i )²)(nΣY i ²-(ΣY i )²) . Y is exactly equal to X if the intercept or system offset term A₀=0, the slope or first order coefficient A₁=1, and the magnitude of CC=1. The above-described linear regression fit to the data can be extended to other relationships for F(x), such as an Nth order regression of Y on X. In the preferred embodiment, this linear regression analysis is extended to M-dimensional space, such that Y is a function of a plurality of independent variables such as the initial interstitial oxygen value OX i and the reciprocal of minority carrier lifetime LT. This type of analysis is termed multi-variate regression analysis for Y. In order to obtain the coefficients for the equation for Y using a series of N data sets, there are a variety of software programs commercially available for both main-frame and desk-top computers. By way of example, the multi-variate regression analysis can be performed by the program "REGGY", by Dr. Jack Prins, STAT 7, in the APL Statistical Libraries. This program is one of a number of programs based on a standard text analysis of multi-variate regression of the type described in Introduction To Multivariate Statistical Analysis, by T.N. Anderson, Wiley, 1966; Statistical Theory and Methodology, by K. Brownlee, Wiley, 1965; Introduction to Statistical Analysis , by W. Dixon and F. Massey, McGraw-Hill, 1969; Introduction to Theory of Statistics, by A. Mood and F. Graybill, McGraw-Hill, 1963; and Applied Regression Analysis, by N. Draper and H. Smith, Wiley, 1966. The software program can be used on both a main-frame computer (IBM installed user program SH20-1841-1) and on an IBM Personal computer. (See "The Directory of Personally Developed Software for IBM Personal computers direct from IBM"). The program will provide multi-variate regression analysis with the option of cross-terms among the various independent variables. By way of example, the typed inputs to be entered into a system using the REGGY program is set forth below. where each of OX i , OX₂, LT, and CARB is a "vector" name for a series of N values. It should be noted that in line 1 the right argument of REGGY is α = 0.05 in order to provide a confidence level of at least 95% (or 1 - α).OX2 (PPM) in line 2 represents the interstitial oxygen measurement after the thermal process steps have been performed. The variable OX2 is the dependent variable. Note that similar analysis can be performed by defining the dependent variable as ΔOX1 = OX1-OX2 . The independent variables in this program are listed in line 4 as OX1 (PPM) (independent variable 1), 1/LT (independent variable 2) (the minority carrier lifetime in microseconds), and CARB (PPM) (independent variable 3). Note that CARB is an additional independent variable which may be utilized to correlate the carbon content of the wafer with the oxygen precipitation. The line 5 relates to the degree of fit and has a 1 typed therein to indicate the use of first order terms only in this embodiment. If cross-terms between the independent variables are desired, then in line 6 a 1 is typed. In lines 9-11 the vectors of the variables to be crossed are listed. Accordingly, to obtain one cross-term between the variables 1 and 2, line 9 is entered as shown. Likewise, in order to obtain cross-terms between independent variables 1 and 3, line 10 is entered as shown. Finally, in order to obtain cross-terms between independent variables 2 and 3, line 11 is entered as shown. It is acceptable statistical practice to censor a small percentage of data called "sports", defined as values for Y falling more than ± 2 or 3 standard deviations from the least-squares fit, F(X1, X2, ...Xn), since the ± 2 and ± 3 std. dev. limits are; expected to contain, respectively, 95.49 and 99.74 percent of the statistical population. In a graphical example provided infra, censoring of the three worst data sets, less than 2.3 percent of the total sample size, caused the recalculated CC to increase from 0.89 to 0.92, and the standard deviation to decrease by 14.9 percent. (said "sports" standard deviation to decrease by 14.9 percent. (said "sports" and much of the data scatter are believed due to limitations of the prototype lifetime test apparatus used). Accordingly, if a sample of wafers representative of a desired range of crystal types in terms of oxygen content range, carbon content range, and microstructure type suitable for thin polished semiconductor slices for LSI product fabrication is chosen, and if the above-described measurements of the initial interstitial oxygen OX i and the minority carrier lifetime LT are made prior to the performance of the thermal process step or steps, and then the interstitial oxygen is measured again for each wafer after the thermal process steps have been performed, then a correlation equation for Y can be obtained as described above. This equation for Y will be calibrated for the particular thermal cycle under test, which can emulate product fabrication or else "benchmark" standardized conditions. After this calibration procedure, the equation for Y is ready for use in predicting the oxygen precipitation content of the crystal types used for the initialization. Accordingly, each wafer to be tested has a measurement of its initial value OX i of interstitial oxygen performed in approximately the same manner as in the initialization step and prior to performing the given set of thermal processing steps on the wafer. Additionally, the minority carrier lifetime LT for the wafer is measured in approximately the same manner as in the initialization step and prior to performing the given set of thermal processing steps on the wafer. Finally, a value Y indicative of the predicted oxygen precipitation value for the wafer is determined using the previously initialized equation for Y. Then, each wafer is physically designated to indicate a particular value of Y or a range of Y values. In a further embodiment, in order to obtain a greater correlation coefficient CC, a linear series equation for Y is utilized with an OX i term, a 1/LT term, and a cross-term therebetween. In a preferred embodiment, the step of correlating to obtain the linear series equation for Y comprises the step of correlating the OX i and the 1/LT values to Y values. However, in this embodiment the Y values are equal to the interstitial oxygen values OX f measured after the performance of the set of thermal processing steps subtracted from the initial interstitial oxygen values OX i , to thereby obtain an equation for Y which predicts the change in the interstitial oxygen from the initial values of oxygen OX i after the performance of the set of thermal process steps. In wafers where the carbon content of the wafer is not de minimus, a step is included of measuring the initial value C of substitutional carbon in the wafers in a non-destructive manner prior to performing the given set of thermal processing steps. (In this regard, see ASTM F123-83. Commercially available FTIR instruments are preferred for taking substitutional carbon measurements on thin wafer samples suitable for SLI product manufacture.) The correlating step then comprises the step of correlating the first set of initial interstitial oxygen values OX i and the reciprocal of minority carrier lifetime values LT and the initial carbon values C to values proportional to the interstitial oxygen OX f measured after the given set of thermal processing steps have been performed in order to obtain the linear series equation for Y. This linear series equation for Y thus has an OX i term, a 1/LT term, and a C term for the measured initial carbon value. In yet a further embodiment of the present method, the step of correlating to obtain a the linear series equation for Y includes an OX i term, a 1/LT term, a C term, and one or more cross-terms between these independent variables. From experiment, it appears that the cross-term between the OX i term and the 1/LT term is most preferred in terms of yielding an optimum correlation. Additional cross-terms between the OX i term and the C term, and between the 1/LT term and the C term, may be added to further improve the CC correlation value. In view of the above, a most preferred version of the correlation equation indicative of the oxygen precipitation in a given wafer is set forth below: Y = ΔOX = A₀+A₁OX i +A₂/LT+A₃C+A₁₂OX i /LT+A₁₃OX i C+A₂₃C/LT. Note that the terms A₀, A₁, A₂, A₃, A₁₂, A₁₃, and A₂₃ are all coefficient constants which are determined during the initialization of the equation for Y. It should be noted that only linear terms are utilized in the equation, although other higher order terms may be added as desired. It may also be desirable to add a triple cross term A₁₂₃. However, a correlation coefficient exceeding 0.9 was achieved with the preferred equation set forth above. With the above-described preferred equation, it is possible to predict the oxygen precipitation in a given wafer, and then either physically mark that wafer in accordance with its predicted oxygen precipitation value, or sort that wafer in accordance with its predicted oxygen precipitation value. Preferred Embodiment: The present method will be demonstrated utilizing boron-doped silicon single-crystal ingots grown by means of the Czochralski pulling method with automatic control, as described in U.S. Patent No. 3,493,770. The wafers resulting from this standard crystal-growing process are typically on the order of 300 to 700 microns in thickness. The first step in implementing the present method for the physical designation of semiconductor wafers according to predicted oxygen precipitation behavior is to initialize the equation for Y on a representative first sample set of wafers which cover a desired range of crystal types in terms of oxygen content range, carbon content range, and microstructure. In the present example, 132 wafers from 8 crystals representing 3 distinct crystal growth processes were utilized for the initialization. These wafers were obtained by means of standard Czochralski puller-model systems. A standard 650°C "oxygen-donor anneal" was performed in ingot form before ingot slicing. Then chem-thinning, and two-sided polishing steps were performed on the wafers. The wafers in this first wafer sample were then measured for initial interstitial oxygen values, initial carbon, and initial minority carrier lifetime. This first wafer sample was then separated into two groups, each having similar distributions in OX i , LT, and C, and one of these two wafer groups was subjected to a 750°C anneal in nitrogen for four hours to "pre-nucleate" these wafers. All wafers were then cleaned and treated by means of a two-sided brush cleaner fed with a solution of ammonium hydroxide and folowed by water rinse. The wafers were then retested for interstitial oxygen content and minority carrier lifetime. It was found that the interstitial oxygen readings made after this pre-nucleation agreed to within a few tenths ppm with the interstitial oxygen reading made before the pre-nucleation. However, it is known that such a pre-nucleation anneal can obscure the carbon readings. Accordingly, carbon readings must be taken prior to any pre-nucleation anneal step. Comparison of the minority carrier lifetime readings showed that pre-nucleating lowered the mean carrier lifetime reading by approximately 93.5%. Accordingly, LT measurements should be made after any pre-nucleation heating steps. The interstitial oxygen measurements listed above may be taken by a fast fourier transform infrared instrument (FTIR). Such an instrument can accurately measure oxygen and carbon concentrations on thin semiconductor wafers suitable for LSI product manufacture and in a manner suitable for subsequent device processing of those wafers. The use of FTIR devices is well known in the art. By way of example, see the application notes from BIO-RAD:Note 102, "Precise and Rapid Measurement of Interstitial Oxygen Concentration in Silicon", by K. Krishnan; and Note 101, "Method to Measure the Precipitated and Total Oxygen Concentration in Silicon", by Jastrzebski, Zanzucchi, Thebault, and Lagowski, reprinted from the Journal of the Electro Chemical Society, 129(7), July 1982. The use of the FTIR method permits a resolution of well below one part per million atomic (ppma) for both oxygen and carbon, on wafers below 0.5 mm in thickness. Note that the FTIR measurement instrument accomplishes a non-destructive measurement on finished wafers. Such FTIR measurements are most accurate for stress-free surfaces and plane-parallel two-side polish. However, correction terms can be added for one-side polished or chem-thinned surfaces. There are a variety of techniques which may be utilized in order to measure the minority carrier lifetime. It should be noted, though, that it appears to be imperative to exclude interference from surface recombination of the minority carriers during such a measurement. Such surface recombination exclusion may be accomplished, by way of example, by strongly accumulating with majority carriers or inverting with minority carriers the semiconductor surface by means of a d.c. bias voltage. Note that accumulation is preferred to inversion, so as to have no induced carrier-depleted region below the surface. In applying such a bias to the wafer surface, the wafer surface must not be damaged. A preferred method for applying this d.c. bias without surface destruction of the wafer is by means of an electrolyte wet-cell. See, for example, the electrolytic wet-cell disclosed in U.S. Patent No. 4,125,440. An electrolyte for the cell should be chosen that can be easily cleaned after the lifetime test and either (a) suppresses anodic oxidation, or (b) deliberately anodizes a few tens of angstroms of the wafer surface, after which the voltage can be reduced. By way of example, electrolytes with an anion which is difficult to oxidize and with a cation which is difficult to reduce such as tetramethylammonium tetrafluoroborate, or 1% ammonium sulfate, may be utilized. Further examples include very dilute solutions (on the order of several percent) of semiconductor grade acids, excluding those which could remove any oxidized silicon layers and/or are hazardous to equipment or operators. For example, several percent of hydrochloric acid in semiconductor grade high purity water could be used. Note that most or all suitable acqueous electrolytes and some suitable water-miscible organics, such as alcohols, glycols, etc., can easily be removed after test by standard cleaning methods which are compatible with subsequent product fabrication. Note also that any extra cleaning steps required due to the use of this electrolyte bath can generally be left to the front-end cleaning step performed in a standard LSI device line. This wet-cell bias arrangement for suppressing surface recombination may be combined with a photoconductive-decay lifetime apparatus of the type which utilizes, for example, microwave energy. A preferred minority carrier testing apparatus may be based on the measurement of the microwave reflectance from a wafer surface which has been irradiated by a laser pulse. A minority carrier lifetime apparatus for carrying out the present minority carrier lifetime measurement is shown in Fig. 1. A polished wafer 10 to be lifetime tested is disposed in a wet-cell 12, comprised of a tank 13 having cylindrical or rectangular movable and/or compressible walls 14 and filled with an electrolyte 18. The polished wafer 10 is disposed on the top of a cylindrical vacuum pedestal 20, including vacuum ports 21 and 20, and is shown in cross-section in the figure. The inner portion of vacuum pedestal 20 should be of quartz or other tough glass, of low dielectric constant and no conductivity, to be transparent to incident photons and microwave radiation. Electrical contact is made to the bottom surface of the polished wafer 10 by means of a ring of deformable or spring-loaded electrodes 22 inside the radius of an O-ring seal 16 which separates the electrodes 22 from the electrolyte 18. Generally the voltage applied to the deformable electrodes 22 is a ground voltage. Preferably, this electrical contact via the electrodes 22 is outside of the area to be measured for minority carrier lifetime and is separated from the wafer area directly contacting the electrolyte. The other electrode required for the application of the bias voltage will be discussed below. The photoconductive-decay-lifetime measuring apparatus may be comprised of an electrically-switchable source of photons 30 for generating electron-hole carrier pairs in the wafer 10, and a microwave sending/receiving circuit 31 for radiating mircowaves toward the wafer 10 and then measuring the microwave reflectance therefrom. In Fig. 1, the photon source 30 comprises a laser diode 30 disposed directly above the top surface of the polished wafer 10 to be measured. The microwave sending/receiving circuit 31 may comprise a microwave sending/receiving antenna 32 of the type disclosed in the reference by Y. Ogita, "Contactless Measurement of Carrier Lifetime in Semiconductor Wafers by Use of Strip Lines", JPN. J. APL. Physics, 19(11) 1980. More specifically, the antenna 32 may comprise a microwave slotline/stripline 32 disposed on a substrate 34. The slotline 32 is formed in the shape of a rectangular U and is disposed above the wafer 10 in such a manner that the laser diode 30 may inject its light radiation through the center space in the U. Accordingly, the substrate 34 must be transparent to the laser diode light radiation and must have a sufficiently low dielectric constant to allow leakage of the microwave radiation from the slotline. In a preferred embodiment, the transparent plate 34 may be formed by a plate of quartz or other tough glass with the slotline 32 disposed on the top surface thereof. In order to apply the d.c. bias to the electrolyte 18, a transparent electrode 36 may be disposed on the bottom surface of the transparent substrate 34. A bias voltage may be applied to this bottom transparent electrode via an electrical line (not shown). This transparent upper electrode 36 in the electrolyte wet-cell biasing arrangement may be formed simply by disposing a high resistance very thin coating of, for example, tin-oxide or gold, on the underside of the transparent substrate 34 which holds the microwave slot/line 32 on its upper surface. Such an electrode coating is essentially transparent to photons and to electromagnetic radiation. The use of this underside coating as the upper electrode in the biasing arrangement permits the precise control of the electrolyte solution thickness between the wafer 10 and the microwave slotline antenna 32. Additionally, this transparent substrate in contact with the electrolyte permits control of the meniscus effects in the microwave properties of the electrolyte solution and minimizes the microwave absorption by the solution if high frequency microwave radiation is utilized. In the alternative, this upper bias electrode may be implemented simply by disposing a circular wire electrode 37 around the inner wall 14 of the wet cell and under the surface of the electrolyte 18, as shown in Fig. 1. The slotline antenna 32 is then free to move separately. Note that it is possible to provide the transparent substrate 34 with a gradual upward curvature away from the center point, to thereby increase the clearance of the slotline antenna 32 from the edges of the wet-cell boundaries. The microwave apparatus to be utilized with the slotline antenna 32 comprises a microwave oscillator 40, a circulator 42 connected to receive a microwave signal from the oscillator 40, a microwave line 44 with a tuner/attenuator 46 connecting the circulator 42 to the slot line antenna 32, and an analyser 48 connected through a diode 50 to the circulator 42. The analyser 48 may simply be an oscilloscope, preferably digital. The microwave line 44 can be waveguide of coaxial, stripline, or slotline type, for efficient coupling to the sender/receiver 32, as per S.B. Cohn, IEEE - MTT-17, No. 10, Oct. 1969, and Cohn et al, ibid, No. 12, Dec. 1969. The apparatus shown in Fig. 1 may be coupled to the above-described FTIR apparatus by means of a computer controlled handler of the type well known in the art. The electrolyte wet-cell 13 may be formed simply by having the walls 14 and 16 rise up into the track, or the wafer pedestal 20 move downward with respect to wall 14, or, in the alternative, the wafer may be picked up and deposited in the electrolyte cell by means of a robot arm. The electrolyte solution 18 is then added to submerge the wafer 10. Note that the wafer 10 is positioned so that one of its 22 surfaces is in electrical contact with the spring-loaded electrodes, inside O-ring 16, which forms one electrode-polarity for the biasing apparatus. The microwave transmitter/receiver apparatus including the slotline antenna 32 disposed on the transparent substrate 34 and the laser diode 30 may then be moved to a position over the wafer 10. As noted previously, the rectangular U shaped slot line antenna 32 is disposed over the wafer 10 in such a manner as to permit the light from the laser diode 30 to propagate through the space between the rectangular U of the antenna 32 and down onto the top surface of the wafer 10. The transparent substrate 34/antenna 32 combination is then lowered such that the electrode formed on the bottom surface of the substrate 34 is in contact with the surface of the electrolyte 18. Accordingly, the electrodes 22 and the coated bottom surface 36 of the transparent substrate 34 (or the alternative electrode 37) form the electrodes which provide the bias voltage to suppress surface-state recombination. The bias voltage applied in order to effect the desired d.c. voltage biasing of the top surface of the wafer 10 is typically on the order of 1 to 10 volts. (See D.E. Hill, J. Appl. Phys. 51, No. 8, Aug. 1980.) It should be noted that it may be desirable to provide the capability for moving the laser diode 30-slotline antenna 32 configuration laterally to permit measurement of the minority carrier lifetime in a multiplicity of locations on the surface of the wafer 10. In operation, with the wafer 10 properly biased in the wet-cell 13, microwave radiation in a preferred frequency range of 1 to 10 GH z or higher (using commercially available components) is directed from the slotline antenna 32 to the wafer 10. However, because the wafer 10 is made from a semiconductor material, this microwave radiation propagates, with little or no reflection, through the wafer. When the laser 30 is pulsed, light photons at the appropriate frequency are directed to the surface of the wafer 10. These photons create electron-hole pairs in the surface layer of the wafer. Microwaves impinging on this wafer surface, with its upper layer filled with electron-hole pairs, penetrate only a skin depth into the wafer's surface before being reflected by this layer of electron-hole charges. These reflected microwaves are received by the slotline antenna 32. The slotline antenna 32 converts this microwave radiation into an electrical signal which is fed via the line 44, the circulator 42, and the diode 50, to the signal analyzer 48. This microwave reflection from the wafer surface layer continues until the electron-hole pairs formed by the laser photons recombine, thereby providing a measurement of minority carrier lifetime. The lifetime apparatus actually employed in the feasibility demonstration (a) used 2 to 10 percent sodium chloride in water, (b) operated with radio wave energy of 150 megahertz, (c) had a vertical-slot wet cell which was conducive to damage of the wafer surfaces, and (d) used very low d.c. bias (well below 0.1V "to suppress film growth"), so that the apparent lifetime values were extremely sensitive to the type of wafer cleaning which preceded the test. A brush-clean with ammonium hydroxide solution and water rinse was used and was found to leave surface ions which worked with the applied d.c. voltage to produce adequate accumulation on P-type silicon wafers. It was essential to remove all salt residues before a wafer dried, to prevent corrosion and strongly adhering films on the wafer. The prototype apparatus also lacked resolution at low lifetimes, because it used a slow rise and fall xenon flash instead of a fast response laser diode. Also, note that the use of microwave energy, rather than radio wave energy (10 GHz vs. 150 MHz) produces greater sensitivity on thin wafers, because of the decreased depth of penetration of the electromagnetic energy into the light injected charge carriers. More specifically, since the skin depth penetration of electromagnetic radiation into a good conductor varies approximately as the reciprocal of the square root of the frequency, for a high injected carrier concentration, an apparatus operated at 10 GHz is more than 8 times as sensitive as the actual 150 MHz frequency used. The actual apparatus utilized to obtain minority carrier lifetime measurements for the present example formed the antenna 32 by means of a marginal oscillator employing a tank coil with eddy current loading and set to operate at 150MHZ. This apparatus was similar to that described in the reference by Verkuil and Curtis, "A High Signal-to-Noise Oscillator For Contactless Measurement of Photoinduced Carrier Lifetimes", ASTM Special Technical Publication 712, 1980. Because the electrolyte utilized to obtain the data for this example was dilute sodium chloride, the wafers were rinsed and cleaned prior to the performance of the thermal process step. Since the feasibility lifetime apparatus used very low d.c. bias, it was found necessary to precede the lifetime test with double-brush and dilute ammonia. This contrasts with the preferred embodiment, wherein the pre-test clean may be optional and wherein the post-test wafer clean comprises the steps of draining the wet-cell 12, and blowing off the wafer 10 with filtered nitrogen, and then automatically transporting the wafer through a double-side brush cleaner, which is flushed with appropriate cleaning agents. In the alternative, the wafer 20 could be sent to an automated single-wafer cleaning apparatus of the type made by the FSI Corporation which utilizes a standardized mechanical interface in which semiconductor-grade cleaning solutions can be applied. It should be noted that it is preferred to have temperature control on the order of 1°C to obtain repeatable and accurate lifetime measurements. It should also be noted that this minority carrier lifetime technique, if it utilizes a fast rise-times and fall-time light pulse, permits a resolution of minority carrier lifetimes below 1 microsecond (fast precipitating wafers) and upper lifetime limits of several orders of magnitude longer. It should also be noted that the microwave slotline antenna disclosed in Fig. 1 is not the only form of transducer which may be utilized. For example, an open rectangular waveguide, or a co-ax antenna, or miniature dipoles may be utilized to implement the transducer 32. However, the slotline antenna is advantageous in that it is not critically sensitive to the distance between the transducer 32 and the sample wafer 10 under test, yet it offers tight lateral resolutions, further delimited by the laser spot size. It should also be noted that there are potentially various other non-wet-cell means for applying the d.c. bias to the wafer or for screening surface recombination without d.c. bias. At this point, a measurement of the interstitial oxygen value OX i and the minority carrier lifetime value LT, and, if desired, the initial carbon value C, have been obtained for each wafer in the first wafer set to be used in the initialization. Now these wafers in this first wafer set are processed with the given set of thermal process steps. In the present experiment, the wafers were to be utilized for a bipolar LSI process line. Accordingly, various thermal process steps were simulated by applying the actual process step temperature, but without the actual process step being performed, i.e., no actual film deposition or photolithography was performed. In the present experiment, a thermal process step simulation was obtained by placing wafers in a suitable furnace under an inert gas, such as nitrogen, at a series of temperatures from 825° - 1105°C for about seven hours total time at high temperatures. The data reported below for the present example is for end-of-process oxygen decreases, i.e., the change in oxygen, ΔOX, from its initial interstitial oxygen reading, OX i , to its final interstitial oxygen reading OX f , after performance of the thermal process steps. In order to obtain this change in interstitial oxygen, ΔOX, the initialization step for the equation for Y wherein OX i and LT and possibly C values are correlated for a first sample set of wafers, are correlated to a set of Y values equal to the interstitial oxygen values, OX f , measured after the performance of the set of thermal processing steps subtracted from their initial oxygen values OX i . It is, of course, understood that a correlation could be obtained for the final value of the interstitial oxygen, OX f , measured after the performance of the thermal process step simply by utilizing final interstitial oxygen measurement values, OX f , during the initialization procedure. It should be noted that since FET products tend to have fewer process steps than bipolar products, the present experiment and correlation demonstrates the feasibility for the most general case of LSI applications. It should also be noted that the measurements in the present experiment were taken primarily from the wafer centers, with occasional samples measured along one or more radii. Clearly, automated embodiments of the present invention could easily measure at a multiplicity of locations per wafer. In the initial analysis which led to the present invention, the relative influences of the different parameters on the Y value (in this instance ΔOX) were obtained by performing a series of simple linear regressions and comparing the resulting values of CC for each different independent parameter. It should be noted that the numerical values to follow are provided only for purposes of illustration. The exact values of the independent variable coefficients will change with changes in wafer growing apparatus, thermal cycles, or extreme changes in crystal type, as well as with use of better lifetime apparatus than the prototype. However, the present correlation method will operate if (a) interferences from surface effects are controlled in the measurement of LT, (b) the interstitial oxygen is measured before significant precipitation is allowed to occur at high temperature processes (e.g., much above 800°C) which would cause loss from solution of some significant fraction of the oxygen in the form of precipitates not detected at the 9 micrometer line, and (c) provided that the carbon, if it is to be measured, is measured before masking thermal treatments are performed. Referring now to Table I, there is shown an attempted correlation between the change in oxygen ΔOX versus 1/LT for 1PPM groupings of measurements of initial interstitial oxygen values obtained for the first wafer sample set of 132 wafers. In this table, the first column represents the initial oxygen range for the initial value OX i in PPM, the second column lists the number of samples in that OX i range, ΔOX vs. 1/L.T. FOR 1 PPM OXYGEN GROUPS OX i -RANGE PPM SAMPLE CORR. COEF. INTERCEPT AO SLOPE A2 under 25 12 .8118 1.2388 150.5098 25 to 26 11 .8956 2.5214 98.3379 26 to 27 15 .8446 5.2783 80.6267 27 to 28 15 .7631 8.2492 80.1154 28 to 29 14 .7157 11.1590 58.2848 29 to 30 19 .3394 14.2110 26.9463 (omit 2) (17) (.5503) (13.9906) (42.2445) 30 to 31 23 .047 20.5661 -2.9022 over 31 23 .0176 22.7175 -1.5510 the third column lists the correlation coefficient CC between the measured values of ΔOX and 1/LT, the fourth column lists the calculated intercepts A₀ for a projected correlation equation of ΔOX=A₀+A₂/LT . The fifth column represents the slope or coefficient A₂ for this proposed correlation equation. It should be noted that the samples listed include both the "raw" and pre-nucleated wafers therein. Examination of Table I clearly shows a decrease in CC for very low (below 25 PPM) initial oxygen values OX i and for high oxygen values of OX i . This is consistent with the concept that very high levels of oxygen drive precipitation due to the super saturation above equilibrium concentration, for the furnace temperature, whereas very low concentrations of oxygen have a decreased tendency to precipitate. At this point, an attempt was made to obtain a linear fit of the intercept coefficient A₀ and the slope coefficient A₂ over this range of initial oxygen OX i groupings. This attempted linear fit resulted in the equation A₀ = -79.3514 + 3.2176 (OX i ) and A₂ = 645.7953 -20.8750 (OX i ) . It was found that the correlation coefficient CC for the intercept coefficient A₀ was CC = 0.9873, while the correlation coefficient CC for the slope coefficient A₂ was CC = 0.9733. These high correlation coefficients for this linear fit of the intercept coefficient A₀ and the slope coefficient A₂ were unexpected. The unusually high correlation coefficients for A₀ and A₂ suggested that a multi-variate regression analysis was possible on the present wafer sample data. Accordingly, a multi-variate regression analysis was performed on the present wafer sample data beginning with an attempted correlation between the initial oxygen measurement OX i and the minority carrier lifetime measurement LT. This multi-variate regression analysis was then extended, one independent parameter at a time, to add cross-terms, and initial carbon values C. Table II shows this multi-variate correlation for the change in interstitial oxygen, ΔOX. It should be noted that in this table as well as in Table I, the initial and final interstitial oxygen readings, OX i and OX f , were in PPM, and the LT measurements were in microseconds, and the initial carbon readings were in PPM. Referring now to Table II, the first column represents the correlation coefficient CC; the second column is the offset coefficient, A₀; the third column is the coefficient, A₁, for the initial oxygen value, OX i ; the fourth column is the coefficient, A₂, for the reciprocal of minority carrier lifetime value, 1/LT, the fifth column is the coefficient, A₃, for the initial carbon measurement C; the sixth column is the coefficient, A₁₂, for the cross-term between the initial oxygen measurement OX i and the reciprocal of the minority carrier lifetime measurement, 1/LT; the seventh column represents the coefficient A₁₃ for the cross-term between the initial oxygen value, OX i ; and the initial carbon reading, C; and finally, the eighth column is the coefficient A₂₃ for the cross-term between the reciprocal of the minority carrier lifetime measurement, 1/LT, and the initial carbon reading, C. The first row in Table II represents a correlation between the change in interstitial oxygen, ΔOX, and the initial oxygen measurement, OX i , and the reciprocal of the minority carrier lifetime measurement 1/LT, i.e., ΔOX = A₀+A₁OX i +A₂/LT . It can be seen that there is a correlation coefficient of .829 for this correlation equation for ΔOX. Row No.2 represents a correlation between the change in interstitial oxygen, ΔOX, and the initial interstitial oxygen value, OX i , and the reciprocal of the minority carrier lifetime measurement 1/LT, and includes a cross-term, A₁₂, therebetween, i.e., ΔOX = A₀ + A₁OX i + A₂/LT + A₁₂OX i /LT . A correlation coefficient CC for this correlation equation is .850. Row No.3 is a correlation between the change in interstitial oxygen, ΔOX, and the initial oxygen value, OX i , and the initial carbon value, C, with no cross-terms, i.e., ΔOX = A₀ + A₁OX i + A₃C . It should be noted that carbon is typically present in the semiconductor wafers due to graphite heaters, shields, or other parts of the crystal pulling furnace, unless special measures are used to flush this out by streaming inert gas and proper flow design. Also, carbon can enter with the starting polycrystalline semiconductor charge. If the carbon content of the wafer is low, i.e., below detection limits or below about 0.1 ppm, then the use of carbon measurement values in the correlation is probably not necessary. It can be seen that for this first wafer set the correlation coefficient obtained for this equation is .803. Row No.4 for Table II is a correlation between the change in interstitial oxygen, ΔOX, and the initial oxygen value, OX i , and the initial carbon value, C, with a cross-term, A₁₃, therebetween, i.e., ΔOX = A₀ + A₁OX i + A₃C + A₁₃OX i C . The correlation coefficient CC for this correlation equation is .814 Row No.5 is a correlation of the change in the interstitial oxygen value, ΔOX, with the initial oxygen reading, OX i , the reciprocal of the minority carrier lifetime, 1/LT, and the initial value of carbon, C, with cross-terms A₁₂, A₁₃, and A₂₃. It can be seen that a correlation coefficient CC of .895 was realized. A comparison of rows No.1 and No.3 shows that row No.1 gives a higher correlation coefficient than row No.3, indicating that the reciprocal of the minority carrier lifetime,1/LT, has a stronger influence than the initial carbon, C, on the change in interstitial oxygen, ΔOX, after performance of the thermal process step. Note again that this change in interstitial oxygen, ΔOX, is essentially equal to the oxygen precipitation in the wafer. It can also be seen that row No.2 has a CC value which is greater than row No.1, indicating that the A₁₂ cross-term provides an improvement in correlation. A comparison of row No.4 to row No.3, which relate to the correlation of the initial oxygen value, OX i , and the initial carbon reading, C, also shows an improvement with the addition of a cross-term. Row No.5 shows a very high correlation on the order of 0.9 for the independent parameters of initial oxygen, OX i , the reciprocal of the minority carrier lifetime, 1/LT, the initial carbon value, C, and the cross-terms, A₁₂, A₁₃ and A₂₃. Accordingly, this correlation equation for Y is the preferred correlation equation. Even better correlation is expected to result from the use of more accurate lifetime measuring apparatus. Fig. 2 shows the actual data points for the calculated ΔOX values versus the actual measured ΔOX values, for the 132 wafers utilizing the correlation equation set forth for row No.5, namely ΔOX=A₀+A₁₀X i +A₂/LT+A₃C+A₁₂OX i /LT+A₁₃OX i C+A₂₃C/LT . It can be seen that for normal distributions, 95.44% of the sample lies within ± 2 standard deviations of the mean, i.e., ± 2 standard deviations of the fit of ΔOX to the independent parameters. Likewise, three standard deviations includes 99.74% of the sample. It should be noted that the graph of Fig. 2 of the calculated versus the actual ΔOX includes three values (2.3% of the 132 total) which are aberrations or "sports", as they lie well outside the 3 standard deviation limits. If these three values are omitted, the correlation coefficient CC increases to 0.9225, and the standard error of the estimate decreases to 2.4 (from 2.82). Accordingly, the new ± 2 standard deviation limits contain all but 5 of the 129 wafers (96.12%), with no wafers falling outside the new three standard deviation limits. Row No.6 in Table 2 shows this correlation with the omission of these three aberrant wafers. From the above, it can be seen that there are a variety of correlation equations which may be utilized, with the equations described for Rows No.1, 2 and 5 being the most preferred. Note that since carbon is less influential than minority carrier lifetime in enhancing oxygen precipitation, and since carbon is fixed in the crystal at crystal growth, it is probably desirable to reduce carbon wafer content to below detection levels. With very low content carbon wafers, the equation for row No.2 may be the most preferred, i.e., Y = A₀ + A₁OX i + A₂/LT + A₁₂OX i /LT. With the properly initialized correlation equation as described above, every wafer can now be non-destructively tested for its initial oxygen content OX i , its minority carrier lifetime value, LT, and its initial carbon value, C, if appropriate, and then a predicted Y value, proportional to the precipitated oxygen in the wafer, obtained. This Y value could then be laser-written in terms of a unique alpha-numeric identifier, on each wafer. This alpha-numeric identifier for the predicted oxygen precipitation could also be stored in a computerized data base which could subsequently be used to correlate LSI product yields and device characteristics to the initial wafer parameters of initial oxygen, OX i , initial carbon, C, and the reciprocal of the minority carrier lifetime value 1/LT. This predicted oxygen precipitation value Y can also be utilized for sorting the wafers in accordance with predetermined threshold Y values. In this regard, automated FTIR equipment such as ("Autofeed/Auto-Scan") the BIO-RAD QS-100 system could be set to sort in accordance with ranges of Y values, provided suitable lifetime apparatus were attached. It is again noted that the particular values of the independent parameter coefficients in the equation for Y depend on the particular thermal cycles being used, and the resolution-accuracy of the particular measurement apparati used for measuring OX i , OX f , LT, and C. In the present example, the change in interstitial oxygen, ΔOX, which is proportional to or equal to the precipitation amount, agreed with the actual measured values to a correlation coefficient of on the order of 0.9 or better, despite the presence of data scatter caused by the available lifetime apparatus used in the experiment. Such data scatter is expected to be minimized if the prototype minority carrier lifetime measurement apparatus used in this example is modified by utilizing a higher wafer biasing voltage in the wet-cell 13 of on the order of 1 to 10 volts to control the surface state effects, by utilizing a higher operating frequency of 1 to 10 GH₂ or higher to obtain more sensitivity on thin wafers, and by utilizing a faster rise-time and full-time light source such as a laser as compared to a flash tube. These device parameter changes should provide excellent resolution at lifetime below 1 - 10 microseconds, as well as for higher lifetime values. It should be noted that other forms of the present algorithms set forth above may be utilized, provided that the basic independent parameters OX i and 1/LT are utilized in the algorithm. For example, the algorithm or equation could incorporate the parameter of time, and/or the oxygen diffusion coefficient for the relevant process temperatures. clearly, this equation need not be limited to the linear terms set forth, e.g., the equation might also include exponentials, etc., where such terms have been established by the physics involved. With reference to the apparatus which may be utilized in order to carry out the present method, the most sensitive means for measuring interstitial oxygen and substitutional carbon is by means of the fast fourier transform infrared apparatus (FTIR). However, it is contemplated that other alternative apparatus will be available in the future. Likewise, the most sensitive apparatus for measuring the minority carrier lifetime in a non-destructive manner on thin finished wafers is presently by some variation of a microwave reflectance photoconductive decay apparatus, with the light source and microwave sender/received antenna at vertical incidence and on the same side of the wafer, preferably the device side of the wafer, so as to minimize back-surface effects. This minority carrier lifetime apparatus should control the surface-state recombination with either strong inversion or strong accumulation, the later being preferred in order to eliminate the added depletion layer beneath the silicon surface. The preferred wet-cell electrolyte biasing arrangement for effecting this accumulation should utilize a readily cleanable solution. Additionally, it may be desirable that the wet-cell 13 have an automated metered pumping apparatus to fill, empty and temperature-control the dilute aqueous dielectric solution or other electrolyte. The actual wet-cell wall 14 (Fig. 1) can be of polytetrafluoroethylyne, polycarbonate, polysteyrene, or other material which is compatible with the electrolyte solutions to be used. Flexible bellows-seals may be added to either the wet-cell wall 14 or to the wafer pedestal 20 to permit wafer handling by air-track apparatus without front-side contact. Such equipment permits the positioning of the wafer on the pedestal 20 such that the pedestal can then be surrounded by the wet-cell wall 14, and the electrolyte solution added. Note that in such an apparatus, either the pedestal 20 can be raised and/or lowered, or the cell wall 14 can rise up into the track. The apparatus may also include a simple multi-port vacuum chuck in combination with a vacuum system that is equipped with appropriate fluid diversion in case of leakage (e.g., a drain-valved fluid trap below the height of the vacuum connection). As noted previously, the microwave reflectance apparatus can be used to take minority carrier lifetime readings in the wafer center, and can also be controlled via a microcomputer, by way of example, with digitally controlled positioners, to test at various locations in the wafer. These various locations could be referenced to alignment flats or notches at the wafer perimeter. Note that it may also be desirable to utilize low intensity flooding of white light onto the wafer 10 in order to minimize trapping effects, as per ASTM-28. Such low intensity white light flooding could be provided by a ring of low power lamps. Because of the vertical incidence of the microwave source and the vertical incidence of the response signals, shielding of the test apparatus may not be necessary, provided that the separation of the microwave slotline antenna 32 and the wafer 10 is small. However, such shielding can be added in the form of a metal or a metal-coated enclosure, if desired. The above-described invention discloses a method for predicting oxygen precipitation behavior in semiconductor wafers such as silicon for a given set of device processing thermal steps. This method can be used to (a) select wafers for a given product/process application, (b) to provide rapid feedback for the optimization of the device process and/or wafer fabrication conditions (including crystal growth and subsequent thermal treatments), and (c) to characterize wafers on on an arbitrary standardized "benchmark" or "figure of merit" scale, which can be used independently of the particular product process. The present method can be automated so as to permit an oxygen precipitation value to be determined for every wafer cut from a boule, rather than relying on statistical inferences from "sister" wafers. Accordingly, this method is a means of obtaining high wafer yield without significant crystal cropping and for optimizing the crystal growth to a given product application. Obviously, there are many modifications and variations of the present invention which are possible in light of the above teachings. It is therefore to be understood that within the scope of the appended claims, the invention may be practiced otherwise than as specifically described.

Claims: A method for the physical designation of semiconductor wafers according to predicted oxygen precipitation behavior for a given set of thermal processing steps including the steps of : Determining a correlating equation from measurements on a given set of sample wafers, by : measuring for each of said sample wafers the initial value OX i of interstitial oxygen in a non destructive manner ; measuring for each of said sample wafers the minority carrier lifetime value LT in a non destructive manner which is free of interference from surface recombination due to uncontrolled fast surface states ; performing on said given set of sample wafers said given set of thermal processing steps ; measuring for each of said sample wafers the value OX f of interstitial oxygen after the performance of the set of thermal processing steps ; determining the coefficients (A₀, A₁, ...) of the correlating equation : Y = OX i - OX f = A₀ + A₁ OX i + A₂/LT by means of regression analysis ; predicting oxygen precipitation behavior for a given product wafer by : measuring the initial value OX i of interstitial oxygen of said given product wafer in a non destructive manner; measuring for said given product wafer the minority carrier lifetime value LT in a non destructive manner which is free of interference from surface recombination due to uncontrolled fast surface states ; and, predicting with said correlating equation the oxygen precipitation behavior of said product wafer by determining the value of interstitial oxygen after said set of thermal processing steps by OX f = OX i - Y The method of claim 1, wherein the said correlating equation is Y = A₀ + A₁ OX i + A₂/LT + A₁₂ OX i /LT The method of claim 2 further including the steps of : measuring the value C of initial carbon in the calibration steps, adding a term comprising C in the correlating equation, and also measuring C for said given sample wafer. The method of claim 3 wherein the term A₃C is added in said correlating equation. The method of claim 3, wherein the term A₁₃ OX i C + A₂₃ C/LT is added in said correlating equation.

Inventors: MILLER, DONALD A.
Applicants: INTERNATIONAL BUSINESS MACHINES CORPORATION
Application year: 1987
Main classification: G01N21
IPC codes: H01L21, G01R31