X-RAY DIFFRACTOMETER WITHOUT MANUAL INTERVENTION FOR DETERMINING THICKNESS OF MULTIPLE NON-METALLICA CRYSTALLINE LAYERS, FOURIER-AND AVERAGE VALUE ZERO TRANSFORMATIONS

Abstract: An x-ray diffractometry technique without manual intervention finds thickness of multiple layers of non-metallic crystalline material. A rocking curve is windowed to eliminate a main peak. The windowed curve is smoothed with a running average. The smoothed curve is subtracted from the windowed curve to yield a difference curve. The difference curve is transformed to make its average value zero and to constrain its endpoints to zero. A welch window is used. A fast Fourier transform is applied to the transformed difference curve. A thickness transform is applied to the result to yield a layer thickness. Examples: Ga As HEMT and silicon germanium thin film structures. General signal processing apparatus with the average value zero translation.

Description: X-RAY DIFFRACTOMETER WITHOUT MANUAL INTER VENΗON FOR DETERMINING THICKNESS OF MULTIPLE NON-METALLICA CRYSTALLINE LAYERS, FOURIER-AND AVERAGE VALUE ZERO TRANSFORMATIONS I. BACKGROUND OF THE INVENTION A. Field of the Invention The invention relates to the field of measuring thickness of materials. The invention relates further to a technique for processing signal data so that a Fourier transform can be performed more successfully. B. Related Art Measuring thickness of layers is of particular usefulness in the semiconductor arts. In the past, people attempted to use rocking curve outputs of x-ray diffractometers to measure thickness of layers. This was a cumbersome process. First, the user would identify bumps in a rocking curve manually and manually put cursors at adjacent maxima. From the position of the cursors, an angle of difference could be derived. From that one could derive a dominant thickness in the material under observation. What layer was dominant at a particular angle would depend on Bragg angle and selection rules. The user would have to guess based on expectations of thickness of particular layers in the material. That guess would be a starting point of simulation. The simulation parameters would then be altered until they produced a simulation curve matching the measured rocking curve. Some also tried to Fourier transform the rocking curve to yield a thickness curve, but the results were so noisy as to be useless. II. SUMMARY OF THE INVENTION The object of the invention is to use an x-ray diffractometer to determine the thickness of multiple non-metallic crystalline layers without manual intervention. Another object of the invention is to improve techniques of taking Fourier transforms. The inventors recognized that the reason that the Fourier transform did not work on the rocking curve was that the average value of the curve was far from zero. Accordingly, the inventors developed a technique for deriving a result curve from the rocking curve such that the result curve would have an average value near zero. This led to an automated method for determining thickness of multiple layers using an x-ray diffractometer. III. BRIEF DESCRIPTION OF THE DRAWING The invention will now be described by way of non-limiting example with reference to the following figures. Fig. 1 shows a rocking curve coming from an x-ray diffractometer measurement of a material composed of multiple non-metallic crystalline layers. Fig. 2 shows a smoothed version of a portion of the curve of Fig. 1 Fig. 3 shows a difference curve resulting from subtracting the curve of Fig. 2 from a portion of Fig. 1 Fig. 4 shows a Fourier transform curve of the curve from Fig. 3 Fig. 5 shows a second rocking curve. Fig. 6 shows a smoothed version of Fig. 5. Fig. 7 shows a difference curve resulting form subtracting the curve of Fig. 6 from a portion of Fig. 5. Fig. 8 shows a Fourier transform of Fig. 7. Fig 9 shows a rocking curve which is incorrectly windowed. Fig 10 shows a smoothed version of Fig. 9 within the incorrect window. Fig 11 shows the difference curve between Fig.s 9 and 10. Fig 12 shows the Fourier transform of Fig. 11. Fig 13 shows a flowchart. Fig 14 is a schematic of an x-ray diffractometer. IV. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Fig. 1 shows a rocking curve of GaAs HEMT (High Electron Mobility Transistor). The horizontal axis of the curve is scaled in units of arc seconds, from D 5275.0" to 0. At D5275" the curve reads 19.8 cps*. The vertical axis measures intensity in units of cps* (counts per second) from 0 to 250. The units on the axes of Fig.s 1-3 and 9-11 are all the same. This rocking curve is received per block 1301 of Fig. 13 from an x-ray diffractometer. The scaling of Fig. 1 has been chosen such that the main peak has been truncated, for clarity. Actual curve scaling is irrelevant. First the main peak, whose average value is VERY far from zero, is removed from the rocking curve using a window shown by line 101 in Fig. 1 per step 1302 of Fig. 13. Line 101 is at -1497.0" and the curve there has a value of 76J cps*. In the first instance, the peak must be windowed manually. However, once the window is set, the window placement in relation to the main peak can be used for all other samples of the same expected composition, so that windowing can be automated. Variations in the position of the peak can be automatically compensated for by an automated search of the points of the rocking curve to find the maximum intensity value. The windowed curve is then smoothed to yield the smoothed curve shown in Fig. 2 per block 1303 of Fig. 13. Those of ordinary skill in the art might devise any number of smoothing functions. For instance, a running average might be used. The rocking curve actually consists of a series of experimental points. In smoothing, the value of a point might be replaced by the average value of that point taken with the two adjacent points on either side. The smoothed curve is then subtracted from the windowed curve to yield the difference curve shown at Fig. 3 per block 1304 of Fig. 13. A transform or normalization is then applied to the difference curve per block 1305. This transform is such as to make the difference curve have an approximate average value near zero, with both end values forced to zero. Preferably a Welch Window is used to normalize the data. A Fast Fourier Transform of the transformed difference curve is made using angle as the transform variable per block 1306 of Fig. 13. The normalized curve data typically includes about 100 data values. This data is copied into the first elements of a large buffer, typically with 2048 points, with remaining data values set to 0. The FFT (Fast Fourier Transform) algorithm transforms this buffer into a power spectrum which is the curve shown in figure 4. The actual magnitude of the points in the power spectrum are not important, because only the position in the curve is needed to compute a thickness not a height. Therefore, for convenience, the power spectrum magnitudes are normalized so that the maximum data value in the curve is set to 100% and all the other values are scaled to their fraction of the maximum value. Thus, in the curve of Fig. 4, the horizontal axis is marked in units of thickness in Angstroms (A) and the vertical axis is marked in normalized percentage. Clear peaks are found at 93 A (A), 683 A (B), and 700 A (C). The position of a peak in the power spectrum is determined by finding a local maximum, the position of which is an index (I) into the buffer. This index is converted to a thickness value (t) , • * w -, 3600 * 180 t = ι * * 2 * N * cos(R) s * π using the following equation: Where W is the wavelength of the X-rays used to take the data in Angstroms; Ν is the number of points in the FFT B is the Bragg angle of the material S is the interval between samples in the rocking curve, in arcseconds W „ 3600 * 180 . , . . . * is the thickness constant. 2 * N * cos(R) s * π Because t is simply a constant multiple of I, the X-axis of the power spectrum is labeled by multiplying the index positions of the points by the thickness constant. Fig. 5 shows a rocking curve for a Silicon Germanium thin film structure. Fig. 5 starts at -7110.0 arc seconds with a value of approximately zero AU, where "AU" means arbitrary units. The units are not significant because the values will be normalized later. Fig. 5 has more than one main peak, i.e. more than one peak whose average value is far from zero. Accordingly, a windowing function 501 is used to remove all main peaks, i.e. all peaks with high rates of change. The window 501 is placed at -1480" where the rocking curve value is 0.05 AU, which is rounded to 0.0 for display of the label. The windowed curve is then smoothed to yield the smoothed curve of Fig. 6. The smoothed curve is then subtracted from the windowed curve to yield the difference curve of Fig. 7. The difference curve is then transformed to make its average value zero and to constrain its end values to zero. The resulting curve is then Fast Fourier transformed using angle as the transform value, then transformed as described above, to yield the transform curve of Fig. 8. Like Fig. 4, Fig. 8 is in units of A on the horizontal axis and normalized percentage on the vertical axis. The peak values are then automatically determined according to a peak finding algorithm. The thickness values established according to Figures 4 and 8 were cross- checked experimentally by the simulation parameter technique described in the background section of this application and found to be accurate to within a range of 1% to 3.5%. Figs 9-12 show the results of improper window placement at the outset. Fig. 9 shows windowing the curve of fig. 1 incorrectly at 901. Fig. 10 shows the resulting smoothed function. Fig. 11 shows the difference curve between fig. s 9 and 10. This curve is then transformed and Fast Fourier transformed as described above to yield Fig. 12. A comparison of Fig. 12 and Fig. 4 shows that the incorrect placement of the initial window has rendered Fig. 12 unusable with a single peak at an inappropriate thickness value. As stated before, the initial window function must be chosen first empirically by a user to eliminate main peaks. However, once the window position is known, all future thickness measurements on the same type of sample can be performed automatically. Fig. 13 shows a flow chart of the operation of the signal processing apparatus of invention as applied in an automated context. The reference numerals of this flow chart have been inserted in appropriate places in the text above to show the steps of the operation. Fig. 14 shows an x-ray diffractometer. The diffractometer includes a source 1401, a sample holder 1403, a detector 1404, and a processor 1405. X-rays 1406, from the source 1401, incident on the sample 1402, produce diffraction radiation 1407 which is detected at detector 1404. Motion of the source, sample, or detector, can yield the so-called rocking curve, which is then provided to the processor 1405. The processor 1405 then performs the operations of Fig. 13.

Claims: CLAIMS: 1. An x-ray system comprising: an x-ray source ( 1401 ); a sample holder (1403); a detector (1404); - a processor (140) coupled with the detector, the processor being arranged to calculate unambiguously, from an output signal of the detector, respective thicknesses of multiple layers of non-metallic, crystalline material, without manual intervention in the calculation. 2. The system of claim 1 wherein the processor is arranged to perform the following operations receiving (1301) data from the detector in the form of a signal curve; receiving (1302) predetermined empirical parameters; generating a thickness curve from the empirical parameters and the signal curve; and interpreting the thickness curve to yield thickness values. 3. The system of claim 2 wherein the empirical parameters comprise a window function for windowing a main peak of the signal curve. 4. The system of claim 2 wherein the signal curve gives an intensity as a function of an angle of either the source or the detector; the generating step comprises - windowing (1302) the curve based on the empirical parameters to eliminate at least one peak, which is relatively larger than other peaks in the curve, to yield a windowed curve; applying (1303) a smoothing function to the windowed curve to yield a smoothed, windowed curve; subtracting (1304) the smoothed, windowed curve from the windowed curve to yield a difference curve having an average value near zero; Fourier (1306) transforming the difference curve to yield the thickness curve; the interpreting operation comprises associating (1305) thickness values with peaks in the thickness curve. 5. The system of claim 3 wherein the processor is arranged to perform the following additional operations: after the subtracting step, norming (1305) the difference curve to yield a normed difference curve having an average value near zero and end values near zero, the Fourier transforming operation then being performed on the normed difference curve; and after the Fourier (1307) transform operation, thickness transforming a result of the Fourier transform to yield the thickness curve. 6. The system of claim 3 wherein the interpreting operation comprises windowing the thickness curve by applying thickness windows at expected thickness values and taking peak values within thickness windows. 7. Signal processing apparatus comprising: - means (140) for receiving a data signal from a signal source as a signal curve; a processor (1405) programmed to process the signal curve to yield an output curve having an average ordinate value near zero; means for interpreting and/or outputting the output curve to yield a measurement relative to the data signal. The apparatus of claim 6 wherein the measurement is a thickness measurement. 9. The apparatus of claim 6 wherein the processor is programmed to perform the following operations: - windowing (1302) the signal curve to eliminate at least one peak that is relatively larger than other peaks in the signal curve, a result of the windowing being a windowed curve; applying (1303) a smoothing function to the windowed curve to yield a smoothed, windowed curve; and subtracting (1304) the smoothed, windowed curve from the windowed curve to yield a difference curve having an approximate average value near zero; and Fourier (1306) transforming the difference curve to yield an output curve. 10. The apparatus of claim 7 wherein the processor is arranged to perform the following additional operations: after the subtracting step, norming (1305) the difference curve to yield a normed difference curve having an average value near zero and end values near zero, the Fourier transforming operation then being performed on the normed difference curve; and - after the Fourier transform operation, transforming ( 1307) a result of the Fourier transform to yield the output curve with appropriate units of measure. 11. The apparatus of claim 6 wherein the smoothing function is a running average.

Inventors: MOORE CHRISTOPHER J L, SUMMERS JAMES A, MOORE, CHRISTOPHER, J., L., SUMMERS, JAMES, A.
Applicants: KONINKL PHILIPS ELECTRONICS NV, KONINKLIJKE PHILIPS ELECTRONICS N.V.
Application year: 1999
Main classification: G01B15
IPC codes: G01B15, G01N23