Fast measurement of X-ray diffraction from tilted layers

Abstract

A method for analysis includes directing a converging beam of X-rays toward a surface of a sample having multiple single-crystal layers, including at least a first layer and a second layer that is formed over and tilted relative to the first layer. The X-rays that are diffracted from each of the first and second layers are sensed simultaneously while resolving the sensed X-rays as a function of angle so as to generate a diffraction spectrum including at least a first diffraction peak due to the first layer and a second diffraction peak due to the second layer. The diffraction spectrum is analyzed so as to identify a characteristic of at least the second layer.

Description

FIELD OF THE INVENTION

The present invention relates generally to X-ray analysis, and specifically to X-ray measurement of thin film properties.

BACKGROUND OF THE INVENTION

X-ray diffractometry (XRD) is a well-known technique for studying the crystalline structure of matter. In XRD, a sample is irradiated by a monochromatic X-ray beam, and the locations and intensities of the diffraction peaks are measured. The characteristic diffraction angles and the intensity of the diffracted radiation depend on the lattice planes of the sample under study and the atoms that make up the crystalline material. For a given wavelength λ and lattice plane spacing d, diffraction peaks will be observed when the X-ray beam is incident on a lattice plane at angles θ that satisfy the Bragg condition: nλ=2d sin θ_B, wherein n is the scattering order. The angle θ_B, that satisfies the Bragg condition is known as the Bragg angle. Distortions in the lattice planes due to stress, solid solution, or other effects lead to observable changes in the XRD spectrum.

XRD has been used, inter alia, for measuring characteristics of single-crystal layers produced on semiconductor wafers. In some of these measurements, relative angular tilt between the layers is taken into account. For example, Pesek et al. describe XRD techniques of this sort in “Lattice Misfit and Relative Tilt of Lattice Planes in Semiconductor Heterostructures,” Semiconductor Science and Technology 6 (1991), pages 705-708, which is incorporated here by reference. The authors used XRD to investigate the tilt in various compound semiconductor layers grown on miscut GaAs substrates.

Cohen et al. describe similar techniques in “High-Resolution X-Ray Diffraction for Characterization and Monitoring of Silicon-on-Insulator Fabrication Processes,” Journal of Applied Physics 93 (2003), pages 245-250, which is also incorporated herein by reference. The authors found XRD to be applicable to multilayered silicon-on-insulator (SOI) structures fabricated by wafer bonding. The tilt and rotation of each crystalline layer with respect to the substrate allowed them to obtain independent measurements of each SOI film.

SUMMARY

Embodiments of the present invention that are described hereinbelow provide improved methods and devices for fast, high-resolution XRD (HRXRD).

There is therefore provided, in accordance with an embodiment of the present invention, a method for analysis, that includes directing a converging beam of X-rays toward a surface of a sample having multiple single-crystal layers, including at least a first layer and a second layer that is formed over and tilted relative to the first layer. The X-rays that are diffracted from each of the first and second layers are sensed simultaneously, while resolving the sensed X-rays as a function of angle so as to generate a diffraction spectrum including at least a first diffraction peak due to the first layer and a second diffraction peak due to the second layer. The diffraction spectrum is analyzed so as to identify a characteristic of at least the second layer.

Typically, sensing the X-rays includes deploying a detector array having elements configured to capture and resolve the X-rays over a range of elevation angles simultaneously, wherein the range is at least 2 degrees.

In some embodiments, sensing the X-rays includes detecting the X-rays over a range of elevation angles at two or more different azimuthal angles so as to generate at least first and second diffraction spectra at first and second azimuthal angles, respectively. Analyzing the diffraction spectrum may then include comparing the first and second diffraction spectra so as to quantify a tilt of the second layer relative to the first layer. In one embodiment, comparing the diffraction spectra includes finding angular locations of the diffraction peaks as a function of the azimuthal angles, and fitting the angular locations to a parametric formula in order to find at least one of a magnitude and a direction of the tilt. Alternatively or additionally, directing the converging beam includes adjusting an azimuthal orientation between the sample and the beam so as to nullify an effect of the tilt on the diffraction spectrum. Typically, the tilt causes a separation between angular locations of the first and second diffraction peaks to vary as a function of the azimuthal angles, and comparing the diffraction spectra may further alternatively or additionally include correcting a measured value of the separation responsively to the quantified tilt.

In a disclosed embodiment, analyzing the diffraction spectrum includes identifying a strain in the second layer.

In one embodiment, a third layer is formed over the second layer, so that the diffraction spectrum includes a third diffraction peak due to the third layer, and analyzing the diffraction spectrum includes finding a property of the third layer responsively to a displacement between the second and third diffraction peaks.

There is also provided, in accordance with an embodiment of the present invention, apparatus for analysis, including an X-ray source, which is configured to direct a converging beam of X-rays toward a surface of a sample having multiple single-crystal layers, including at least a first layer and a second layer that is formed over and tilted relative to the first layer. A detector assembly is configured to simultaneously sense the X-rays that are diffracted from each of the first and second layers while resolving the sensed X-rays as a function of angle so as to generate a diffraction spectrum including at least a first diffraction peak due to the first layer and a second diffraction peak due to the second layer. A processor is coupled to analyze the diffraction spectrum so as to identify a characteristic of at least the second layer.

The present invention will be more fully understood from the following detailed description of the embodiments thereof, taken together with the drawings in which:

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic side view of an X-ray metrology system, in accordance with an embodiment of the present invention;

FIGS. 2A and 2B are schematic sectional views of semiconductor wafers showing details of layers formed thereon, whose properties are analyzed in accordance with an embodiment of the present invention;

FIGS. 3A and 3B are schematic representations of HRXRD spectra obtained at different azimuthal angles from a sample having tilted layers, in accordance with an embodiment of the present invention; and

FIG. 4 is a plot that schematically shows a variation in peak position as a function of azimuthal angle in HRXRD spectra obtained from a sample having tilted layers, in accordance with an embodiment of the present invention.

DETAILED DESCRIPTION OF EMBODIMENTS

When XRD measurements are to be made on tilted layers using conventional methods of detection, painstaking adjustment and scanning of the X-ray sample and detection angles are generally required. The diffraction spectrum must typically be measured first at one rocking angle to find the substrate spectrum and, for highly tilted layers, at another, different angle to find the spectrum of the tilted layer. Finding the correct angles often entails a lengthy process of trial and error, and it is possible that spectral components due to tilted layers may be missed entirely.

Embodiments of the present invention that are described hereinbelow provide improved methods and systems for use in fast HDXRD measurement and analysis of tilted epitaxial films on a sample. “Tilt” in this context refers to non-parallel orientation of the crystal planes of the film relative to the planes of an underlying layer. In one embodiment, for example, the sample is a silicon wafer in which the silicon substrate constitutes a first layer, and a second layer comprising an epitaxial film is formed over the substrate. Thus, references to a “layer” in the description and claims that follow should be understood, where appropriate, to refer to the substrate, as well. The second layer may be tilted relative to the first layer either by design, or due to imperfections in the fabrication process. Alternatively, the methods and systems described herein may be applied to samples and crystalline layers (including, but not limited to, epitaxial films) of other types.

An X-ray source directs a converging beam of X-rays toward the surface of the sample. A detector assembly senses the X-rays that are diffracted from each of the first and second layers (including the substrate “layer,” as noted above), while resolving the sensed X-rays as a function of angle. For this purpose, the detector assembly typically comprises an array of detector elements, which detect the diffracted X-rays at different, respective elevation angles. The detector assembly thus generates a diffraction spectrum, which includes respective diffraction peaks due to at least the first and second layers. Because the system both irradiates and detects the scattered X-rays over a range of angles simultaneously, the entire spectrum, including the Bragg diffraction peaks due to two (or more) different layers, can be captured rapidly in a single position of the source and detector assemblies, without the need for mechanical scanning. A processor analyzes the diffraction spectrum in order to identify characteristics of the tilted layer or layers.

FIG. 1 is a schematic side view of a system 20 for HDXRD of a sample 22, in accordance with an embodiment of the present invention. Sample 22 is mounted on a motion stage 24, allowing accurate adjustment of the position and orientation of the sample. An X-ray source 26 directs a converging X-ray beam 28 toward a small area 34 on sample 22. Typically, source 26 comprises an X-ray tube 30 with suitable optics 32 to focus and monochromatize beam 28. Beam 28 typically subtends at least 2°, and may subtend as much as 4° or even more, depending on optics 32, in order to irradiate sample 22 over a large range of angles simultaneously.

X-rays are diffracted from sample 22 in a generally diverging beam 38, which is received by a detector assembly 36. The detector assembly typically comprises a detector array 42, such as a CCD array, comprising multiple detector elements, configured so as to resolve beam 38 as a function of elevation angle θ. Typically, the angular span of array 42 is comparable to that of beam 28, i.e., at least 2°, and possibly 4° or greater. A knife edge 40 and/or other optical elements may be used to limit beam 28 and/or beam 38 and block undesired scattered radiation that might otherwise strike array 42 and interfere with the diffraction measurement.

Various types of X-ray sources and detector assemblies may be used in system 20. Details of such components are described, for example, in U.S. Pat. Nos. 7,076,024, 7,120,228 and 7,551,719, whose disclosures are incorporated herein by reference.

The positions of source 26 and detector assembly 36 are controlled by motion assemblies 44 and 46, respectively. In the simplified view shown in FIG. 1, the motion assemblies comprise curved tracks, which permit the source and detector assembly to be positioned at the appropriate elevations, typically in the vicinity of the Bragg angles of the layers that are to be analyzed. Other suitable motion assemblies will be apparent to those skilled in the art. For the sake of this example, it is assumed that the lattice planes creating the diffraction pattern are approximately parallel to the surface of sample 22, so that the incidence and takeoff angles defined by beams 28 and 38 relative to the surface are both equal to the Bragg angle. Alternatively, source 26 and detector assembly 38 may be positioned at different incidence and takeoff angles in order to measure diffraction from lattice planes that are not parallel to the surface of sample 22.

In addition, as noted above, stage 24 may be configured to translate the X-Y location on the sample that falls within area 34, as well as to rotate the azimuthal angle φ of the sample relative to beam 28. (As shown in FIG. 1, the X-Y plane is taken to be the sample surface, with the Z-axis perpendicular to the surface; θ is the elevation angle relative to the Z-axis; and φ is the azimuthal angle of rotation about the Z-axis.)

A signal processor 48 receives and analyzes the output of assembly 36, so as to measure a spectrum 50 of the flux of X-ray photons diffracted from sample 22 as a function of angle at a given energy or over a range of energies. Typically, sample 22 has one or more thin surface layers, such as thin films, at area 34, so that distribution 50 as a function of elevation angle exhibits a structure that is characteristic of diffraction effects due to the surface layer and underlying layers. Processor 48 analyzes the angular spectrum in order to determine characteristics of one or more of the layers of the sample, such as the composition, thickness, lattice strain and/or tilt angle of the layer, using methods of analysis described hereinbelow.

The components of system 20 and the techniques described herein may be integrated into systems that provide other types of measurement functionality, such as X-ray reflectometry and scattering measurements. Additionally or alternatively, these components and techniques may be integrated as process monitoring tools in manufacturing systems, such as systems for semiconductor wafer fabrication. Integrated systems of these types are described in greater detail in the above-mentioned patents.

FIGS. 2A and 2B are schematic sectional views of sample 22, showing two examples of layers formed on the sample. The properties of these layers are analyzed in accordance with embodiments of the present invention that are described below. In these embodiments, the sample is a silicon wafer (which is referred to as a substrate 52, or alternatively as a carrier layer, as explained above). In FIG. 2A, a SiGe layer 58 is formed on substrate 52. In FIG. 2B, a buried oxide (BOX) insulating layer 54 is formed over substrate 52, and a silicon-on-insulator (SOI) layer 56 is deposited on the BOX layer. SiGe layer 58 is then formed on the SOI layer.

Sample 22 is analyzed in system 20 in order to find the diffraction angle of SiGe layer 58, which varies as a function of the layer composition and thus gives an indication of the germanium concentration. The diffraction angle will, however, be offset from the Bragg angle if the lattice planes in layer 58 are rotated (tilted) with respect to the corresponding planes in the underlying substrate 52. In the case of bulk substrates, tilt in layer 58 may arise when it is grown on a substrate whose surface is inclined relative to a major crystallographic plane, i.e., a miscut substrate. The Bragg angle may also be affected by strain in SOI layer 56, resulting in changes in the lattice spacing in both layer 56 and layer 58. These offsets are overcome by the ability of system 20 to capture high-resolution XRD spectra over a range of elevation angles, as illustrated in the figures that follow.

FIGS. 3A and 3B are schematic representations of HRXRD spectra obtained at different azimuthal angles from sample 22, having tilted layers as shown in FIG. 2B, in accordance with an embodiment of the present invention. The spectra are shown on a logarithmic scale in counts per second (CPS) captured by the elements of array 42, as a function of the respective elevation angles of the elements. The angular scale is adjusted, for the sake of convenience, so that a peak 60 due to diffraction from substrate layer 52 is taken as the origin (θ=0).

A second peak 62 due to SOI layer 56 is offset relative to peak 60 as the result of the tilt in layer 54 and possibly lattice strain in layer 56. This offset varies with azimuthal angle, as shown by the difference in peak positions between FIGS. 3A and 3B. There is little or no tilt (less than 1 arc minute), however, between layers 56 and 58. Therefore, the offset between peak 62 and a peak 64 due to SiGe layer 58 (which is about 0.5°) remains constant and indicates accurately the concentration of germanium in the layer. Thus, because layer 56 is available to serve as a reference, there is no need in this example to actually measure the tilt in order to extract the desired diffraction data with respect to layer 58.

Alternatively, the difference in peak positions in spectra taken at different azimuthal angles can be used to measure the layer tilt explicitly. Because of the change of the azimuth between FIGS. 3A and 3B, peaks 62 and 64 in the two spectra, due to SOI layer 56 and SiGe layer 58, are shifted with respect to substrate peak 60. This varying angular shift of the layer peak positions relative to that of the underlying substrate is expected, since the inclination angles of layers 56 and 58 are reversed with respect to the substrate layer 52 when the sample is rotated 180° in φ (about an axis normal to the surface of the substrate 52, corresponding to right/left reversal in FIG. 2B).

Measurements of the dependence of peak offset on the azimuth can be applied to calculate the relative tilt between the layers. The offset of either peak 62 or 64 can be used for this purpose, depending on the detailed structure of sample 22. In the present example peak 64, due to SiGe layer 58, would probably be chosen for tilt measurement, since this layer is thicker, and the peak is sharper and more intense; hence its position can be determined more precisely.

FIG. 4 is a plot that schematically shows a variation in relative peak positions as a function of azimuthal angle in HRXRD spectra obtained from the tilted layers of sample 22, in accordance with an embodiment of the present invention. Data points 70, 72, 74 and 76 indicate the position of peak 64, due to SiGe layer 58, relative to substrate peak 60, measured at four different azimuths. A curve 78 is fitted to the data points in accordance with a parametric model, which is described below. This fit gives a tilt angle of 1608 arcsec.

The relative tilt between the planes in layer 58 and corresponding planes in substrate 52 is obtained by the following relationship:Tilt=[Δθ(φ)+Δθ(φ+180)]/2  (1)wherein Δθ(φ) is the angular separation between layer peak 64 (or 62) and substrate peak 60 measured at azimuth φ.

Additional measurements at other azimuthal angles can be used to further improve the estimation of the tilt and also to determine its direction with respect to a reference point, such as an orientation notch on sample 22. In this case, the measured angular separations Δθ(φ) between the peaks as a function of azimuth φ are fitted to the following parametric formula:Δθ(φ)=Δθ₀+TiltMagnitude×cos(φ+TiltDirection)  (2)The constant (Δθ₀), TiltMagnitude and TiltDirection parameters can be fitted, using a non-linear regression algorithm, to three or more measurements of the layer-to-substrate peak separations. A fit of this sort is illustrated in FIG. 4.

The tilt that is measured in this manner can be used in finding the peak separation with improved accuracy, which in turn gives a more accurate value of the strain and/or composition of layer 58. There are a number of ways in which tilt and peak separation can be measured and used, depending, inter alia, on the capabilities of the HDXRD measurement system.

Assuming, for example, that stage 24 is capable of moving sample 22 so that beam 28 can access all sites (x,y) on the sample at all azimuths φ, then the method described above can be used to measure the tilt at each and every site on the wafer. The measured peak separation Δθ(φ) is adjusted accordingly to give the corrected separation:Δθ′(φ)=Δθ(φ)−TiltHere Tilt is determined for the specific azimuth φ. This approach may be useful if the magnitude of the tilt varies significantly across the surface of sample 22, but it reduces the throughput of system 20 since at least two measurements are required at all sites.

An alternative method is to find the tilt properties at a single site and then apply this tilt to all subsequent measurement points. This approach can be employed if the capabilities of stage 24 are limited or in order to improve throughput. This method is useful, for example, if stage 24 does not allow all azimuth values φ to be accessed at all sites (x,y) on the wafer. In this case it is assumed that the tilt does not vary significantly across the surface of the sample and may be regarded as a constant value for a specific azimuth.

As a further alternative, the tilt magnitude and direction may be determined using three or more measurements at a single site on the wafer. This information is then applied in driving stage 24 to set the azimuthal orientation of sample 22 such that the effect of the tilt, as described in equation (2), is nullified. This condition is satisfied when the sample azimuth is oriented at a right angle (±90°) with respect to the fitted TiltDirection. This approach can be used in measuring symmetrical reflections when the diffraction planes are almost parallel to the sample surface, as long as the TiltDirection is constant across the wafer. Thus, in this case, the effects of tilt at all measured sites are eliminated by determining the tilt properties at a single site. When the sample azimuth is oriented in this way, it does not matter if the TiltMagnitude varies across the wafer, because the effective tilt at the proper azimuth will remain equal to zero.

Although the embodiments shown in the figures above relate to specific types of layer structures on a silicon wafer, the principles of the present invention are similarly applicable in making rapid XRD measurements on crystalline samples of other types in which the lattice angles may not be well controlled or precisely known in advance. It will therefore be appreciated that the embodiments described above are cited by way of example, and that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of the present invention includes both combinations and subcombinations of the various features described hereinabove, as well as variations and modifications thereof which would occur to persons skilled in the art upon reading the foregoing description and which are not disclosed in the prior art.

Claims

1. A method for analysis, comprising: directing a converging beam of X-rays toward a surface of a sample having multiple single-crystal layers, including at least a first layer and a second layer that is formed over and tilted relative to the first layer;simultaneously sensing the X-rays that are diffracted from each of the first and second layers while resolving the sensed X-rays as a function of angle so as to generate a diffraction spectrum comprising at least a first diffraction peak due to the first layer and a second diffraction peak due to the second layer,wherein sensing the X-rays comprises detecting the X-rays over a range of elevation angles at two or more different azimuthal angles so as to generate at least first and second diffraction spectra at first and second azimuthal angles, respectively; andanalyzing the diffraction spectrum so as to identify a characteristic of at least the second layer,wherein analyzing the diffraction spectrum comprises comparing the first and second diffraction spectra so as to quantify a tilt of the second layer relative to the first layer.

2. The method according to claim 1, wherein sensing the X-rays comprises deploying a detector array having elements configured to capture and resolve the X-rays over a range of elevation angles simultaneously.

3. The method according to claim 2, wherein the range is at least 2 degrees.

4. The method according to claim 1, wherein comparing the first and second diffraction spectra comprises finding angular locations of the first and second diffraction peaks as a function of the azimuthal angles, and fitting the angular locations to a parametric formula in order to find at least one of a magnitude and a direction of the tilt.

5. The method according to claim 1, wherein directing the converging beam comprises adjusting an azimuthal orientation between the sample and the beam so as to nullify an effect of the tilt on the diffraction spectrum.

6. The method according to claim 1, wherein the tilt causes a separation between angular locations of the first and second diffraction peaks to vary as a function of the azimuthal angles, and wherein comparing the first and second diffraction spectra comprises correcting a measured value of the separation responsively to the quantified tilt.

7. The method according to claim 1, wherein analyzing the diffraction spectrum comprises identifying a strain in the second layer.

8. The method according to claim 1, wherein a third layer is formed over the second layer, and wherein the diffraction spectrum comprises a third diffraction peak due to the third layer, and wherein analyzing the diffraction spectrum comprises finding a property of the third layer responsively to a displacement between the second and third diffraction peaks.

9. Apparatus for analysis, comprising: an X-ray source, which is configured to direct a converging beam of X-rays toward a surface of a sample having multiple single-crystal layers, including at least a first layer and a second layer that is formed over and tilted relative to the first layer;a detector assembly, which is configured to simultaneously sense the X-rays that are diffracted from each of the first and second layers while resolving the sensed X-rays as a function of angle so as to generate a diffraction spectrum comprising at least a first diffraction peak due to the first layer and a second diffraction peak due to the second layer,wherein the detector assembly is configured to sense the X-rays over a range of elevation angles at two or more different azimuthal angles so as to generate at least first and second diffraction spectra at first and second azimuthal angles, respectively; anda processor, which is coupled to analyze the diffraction spectrum so as to identify a characteristic of at least the second layer,wherein the processor is configured to compare the first and second diffraction spectra so as to quantify a tilt of the second layer relative to the first layer.

10. The apparatus according to claim 9, wherein the detector assembly comprises a detector array having elements configured to capture and resolve the X-rays over a range of elevation angles simultaneously.

11. The apparatus according to claim 10, wherein the range is at least 2 degrees.

12. The apparatus according to claim 9, wherein the processor is configured to find angular locations of the first and second diffraction peaks as a function of the azimuthal angles, and to fit the angular locations to a parametric formula in order to find at least one of a magnitude and a direction of the tilt.

13. The apparatus according to claim 9, and comprising a stage, which is configured to adjust an azimuthal orientation between the sample and the beam so as to nullify an effect of the tilt on the diffraction spectrum.

14. The apparatus according to claim 9, wherein the tilt causes a separation between angular locations of the first and second diffraction peaks to vary as a function of the azimuthal angles, and wherein the processor is configured to correct a measured value of the separation responsively to the quantified tilt.

15. The apparatus according to claim 9, wherein the processor is configured to identify a strain in the second layer by analyzing the diffraction spectrum.

16. The apparatus according to claim 9, wherein a third layer is formed over the second layer, and wherein the diffraction spectrum comprises a third diffraction peak due to the third layer, and wherein the processor is configured to find a property of the third layer responsively to a displacement between the second and third diffraction peaks.