The present invention relates to providing physical measurements of the thickness of a material and more particularly relates to measuring the thickness and rate of etching, and composition of such etched material while the material is being etched.
Many micro-electromechanical systems (MEMS) devices, sensors, integrated circuits, and optical and electro-optical elements require controlled removal of materials such as silicon and silicon oxides. A range of etching processes, including dry and wet etching processes, can be used to remove material to produce patterns or useful features. Moreover, many sensors and MEMS devices are required to operate in harsh chemical environments. In these cases, etching of the materials used in the sensors by the chemical environment can lead to device failure. Consequently, measurement of etching of the materials of construction to be used in the device is required. A key aspect of etching processes is monitoring the thickness of the material. The materials can be homogeneous, such as silicon metal, or heterogeneous, such as silicon oxide layers coated on silicon metal substrates. Because of the critical dimensions involved, it is advantageous to be able to accurately measure the material thickness, not only after etching, but also in situ, as the etching occurs. That is, it is desirable to be able to measure material thickness and rate of change of material thickness dynamically. It is also desirable to monitor the composition of the material. This is particularly important for heterogeneous materials where etching can lead to changes in the composition of the material.
There are inherent difficulties that complicate the measurement process in etching processes that make some conventional approaches unworkable for in situ measurement. Some current ways to characterize etching, for example of Si and SiO2, include quartz crystal microbalance techniques, profilometry, potentiometry, spectroscopic ellipsometry, and spectrophotometric methods (FTIR, UV-Vis)). Quartz crystal microbalance techniques can be used to carry out accurate etching measurements (K. T. Lee, S. Raghavan, “Etch Rate of Silicon and Silicon Dioxide in Ammonia-Peroxide Solutions Measured by Quartz Crystal Microbalance Technique” Electrochemical and Solid-State Letters, vol. 2, 172-174 (1999)). However, quartz crystal microbalance methods like this require that the material to be etched be first coated on the quartz crystal. This is not convenient and limits the application of this method to materials from which suitable coatings on the quartz crystal monitor can be made. Even for materials, which can be coated, the quartz crystal microbalance technique is limited to coatings, which can be prepared in the operating range of the quartz crystal, that is temperatures below 570° C.
Potentiometric methods can be applied to evaluation of in situ etching of materials, notably silicon such as the open circuit method such as that reported by EP Patent Application No. 0725435A2 by Schmidt et al. entitled “Electrochemical Measurements for in-situ Monitoring of Semiconductor Wafer Cleaning Processes”. However, this method provides an indirect measure of material thickness and requires calibration of the potentiometric output, such as the voltage potential, by another means such as spectroscopic ellipsometry. It is also inconvenient because in practice this method requires electrical contacts to be made on the material, such as by vacuum deposition of a metal on a surface of the material to be etched. Furthermore, this method requires use of a reference electrode, which can limit the usefulness of this method. For example, in the etching environment, the electrode can be degraded by the etching solution.
Spectroscopic ellipsometry and methods that use spectrophotometric techniques such as FTIR, UV-VIS and optical emission can be used to perform in situ measurements. U.S. Pat. No. 7,049,156 entitled “System and Method for In-Situ Monitor and Control of Film Thickness and Trench Depth” by A. Kueny describes a method for measuring thickness of a layer using spectral reflectometry and comparing to known models and requires development of complex algorithms for each new material investigated. U.S. Pat. No. 6,888,639 by Goebal et al. entitled “In-Situ Film Thickness Measurement Using Spectral Interference at Grazing Incidence” also describes a reflectance spectroscopic technique but requires measuring at large grazing angles. Both of these techniques are limited to front side applications. In U.S. Pat. No. 6,413,867 by Sarfaty et al., entitled “Film Thickness Control Using Spectral Interferometry”, a further variation of use of spectral reflectance interferometry is described. No. In this method, the observed spectral interference fringes as a function of wavelength are compared to a reference data set using pattern recognition techniques in order to determine when the appropriate etching end point has been reached. This technique does not measure etching rates and the requirement for a reference sample limits its scope of applicability.
U.S. Pat. No. 5,694,207 entitled “Etch Rate Monitoring by Optical Emission Spectroscopy” by Hung et al describes an indirect method of measuring the rate of plasma etching on a silicon wafer by measuring optical emission from the plasma. This method infers etch rates based on concentration of gases in the plasma and is only a front sided measurement.
Another spectroscopic approach is to study changes in the etching environment to infer etching rates of the material. For example, D. Chopra et al. in a paper entitled “In-situ Measurements of Ultrathin Silicon Oxide Dissolution Rates” in Thin Solid Films, Vol. 323, pp 170-173, 1998 uses a chemical probe dissolved in the etchant to enable spectroscopic evaluation of etching rates of silicon oxide. However this requires adding a tracking agent to the etchant and the method is inherently indirect.
The use of laser reflectance reflectometry is described by E. Steinsland et al., “In Situ Measurement of Etch Rate of Single Crystal Silicon”, paper 2D3.12P in Transducers 97, IEEE pages 707-710. This method measures the intensity of light reflected off of a silicon wafer while being etched as a function of time and measures the change in thickness by counting the build up of interference fringes. This technique can provide only relative rate information but not total thickness. Surface roughness on the sample will greatly affect the results of this type of measurement and it is limited in resolution to about 0.1 μm.
All of the above techniques require special view ports in order to protect their optical components from the etching environment. In some cases this is undesirable since it can complicate the etching set-up and adds cost. Moreover, these spectroscopic techniques are limited in their application to optically transparent etching environments. In practice, the signal to noise can be significantly reduced in the etching environment by the presence of optically dense materials such as dyes. Profilometry is not suitable for in situ measurements because it requires removal and manipulation of the sample.
An alternative solution to these limited in situ methods is to use a surrogate “witness plate” that can be subjected to the etching process and removed after a period in order to allow accurate measurement of etching outside the etching environment. For example, the witness plate can be measured outside of the etching environment by spectroscopic ellipsometry. However, such a solution requires space in the etching environment, requires an interface for its removal and reinsertion, introduces additional surface area and waste, and necessitates time delay so that the ability to obtain dynamic measurement data is compromised.
Although the methods described in the above listing may provide some measure of accuracy in determining etching, there is a significant need for improvement. In situ measurement would provide the most highly accurate data for determining the rate of etching, useful in maintaining precision control of the etching process and characterizing the chemical compatibility of the material. There exists a need for an improved method for measuring etching of materials including coated materials.
In accordance with the present invention, there is provided a method of measuring the thickness and the rate of change of thickness of a material having a surface while the material is being etched, comprising:
a) illuminating the material surface with low coherence light, a portion of which transmits through the material and a portion of which is reflected;
b) etching the material surface and while etching, collecting a portion of the reflected light from each optical interface of the material with a low coherence light interferometer;
c) calculating the thickness and rate of change of thickness of the material or part of the material according to the obtained interferometric data; and
d) storing or displaying the resultant thickness and rate of change of thickness of the material.
The present invention provides an effective method of measuring the rate of change of thickness of material while the material is being etched.
As a further advantage, the method of the present invention allows real-time monitoring of the etching rate, useful in a control loop that regulates the etch rate.
In another aspect of the present invention, there is provided a method of measuring the thermo optic coefficient of a material comprising:
a) illuminating the material with low coherence light, a portion of which transmits through the material and a portion of which is reflected;
b) heating or cooling the material over a defined time interval;
c) collecting a portion of the reflected light from each optical interface of the material with a low coherence light interferometer at a multiplicity of times within the defined time interval;
d) calculating the optical thickness of the material at the said multiplicity of times according to the obtained interferometric data;
e) monitoring the temperature of the material as a function of time during the defined time interval;
f) calculating the thermo optical coefficient by determining the slope of the change in optical thickness with respect to temperature during the defined time interval; and
g) storing or displaying the thermo optic coefficient of the material.
The present invention provides a unique way of calculating the thermo optic coefficient of a material. This method can be used simultaneously with etching the material so that changes to the etching rate can be made in real time.
In accordance with the present invention it has been determined that low coherence interferometry can be used to measure etching rates of materials in situ, during the etching process. The methods and apparatus of the present invention are particularly well suited to determine etching rates of homogeneous materials and materials comprised of coated substrates such as coatings on silicon wafers, on glass and on other flat substrates, and to determine the stability of these materials and coatings when subjected to an environment containing etchants.
In the present invention etching processes, which can be monitored by in situ low coherence interferometry, include any process used to remove material. Examples of etching processes include chemical polishes and wet etching processes that use acidic solutions such as hydrofluoric acid (HF) and hydrofluoric acid/nitric acid/acetic acid (HNA) mixtures, and wet etching processes that use basic solutions, such as potassium hydroxide (KOH) solutions and ethylenediamine solutions. Overviews of some wet etching processes used in micromachining applications can be found in S. Wolf and R. N. Tauber, “Silicon Processing for the VLSI Era, Vol. 1, Process Technology,” Lattice Press, Sunset Beach, pp 514-538, 1986; D. L. Kendall, R. A. Shoultz, “Wet Chemical Etching of Silicon and SiO2, and Ten Challenges for Micromachiners,” in Microlithography, Micromachining, and Microfabrication, Vol. 2, P. R. Choudhury, Ed., SPIE Optical Engineering Press, London, 1997, pp 41-97. In addition, low coherence interferometry can be used to monitor the stability of materials in fluid management systems in which the fluids can etch the surfaces of the materials in contact with the fluids.
Additional examples of methods for removing material include dry etching processes, such as reactive ion etching (RIE), deep reactive ion etching (DRIE), plasma-etching, ion milling, and sputtering. Overviews of these dry etching processes can be found in S. Wolf and R. N. Tauber, “Silicon Processing for the VLSI Era, Vol. 1, Process Technology,” Lattice Press, Sunset Beach, pp 539-585. In addition to dry- and wet-etching, low coherence interferometry can also be used to measure material removal during mechanical material removal processes, like chemical mechanical polishing (CMP), used for example in the fabrication of semiconductor integrated circuits. These examples of etching processes, such as wet and dry etching and CMP, which can be monitored by in situ low coherence interferometry, are meant to be instructive and not limiting.
All of the fibers in the apparatus shown in
Data acquisition, analysis and display of data are performed utilizing a computer, instrument control and display unit 10 containing appropriate hardware, such as National Instrument data acquisition cards. A bidirectional communication interface 20 is used to control data flow from the interferometer to the computer by sending appropriate control signals to the interferometer 70 including control of piezoelectric modulators 108 and 109, monitoring the detector signals from detectors 96 and 97 and providing data triggering signals. The periodicity of the laser light is utilized to track the optical distance that the low-coherent light interferometer modulators scan. In our examples signal processing and data analysis routines are run under a Labview program development environment (available from National Instruments) running on computer, instrument control and display unit 10 to analyze the low-coherent light interferograms resulting from reflections at optical interfaces in the sample.
The laser 101 in the interferometer is utilized to track the distance the optical path has changed during the push pull operation of piezoelectric modulators 108, 109 in the all fiber interferometer shown in
For the low-coherence broadband light source 76, constructive interference occurs when the path lengths of the two arms in the interferometer are equal within a few coherence lengths. In order for constructive interference to occur, light must be reflected back into the interferometer from the sample 5. This will occur at each optical interface in the sample 5. The distance between adjacent interference peaks represents the optical thickness (group index of refraction (n) times the physical thickness) of the materials making up the sample 5.
Since the instrument uses a stabilized laser light source for providing constant distance interval measurements, the instrument measures absolute optical path distance defined as (n) multiplied by physical thickness. The measurement configuration of the interferometer is the optical autocorrelation mode, in which light reflecting from the sample is input to both arms of the Michelson interferometer. In the autocorrelation mode, light reflecting from the sample is made to interfere with itself, and both arms of the interferometer see reflections from all of the optical interfaces in the sample. As the path lengths of the two arms of the interferometer are changed, a series of interference peaks are observed, indicating the optical path differences between adjacent optical interfaces. The self-correlation condition occurs when the two path lengths of the Michelson interferometer are equal, in which case all optical interfaces in the sample interfere constructively. The measured distance between the largest peak, at zero path length difference, and the first set of adjacent peaks is the shortest optical path difference in the sample.
An alternate configuration for an all fiber based interferometer is shown in
It is instructive to describe how the expected interferometric signals are derived and how the calculations are performed. It is assumed that there is minimal absorption and scattering in the material so that peak intensities are determined by reflection and transmission and index of refraction. Assume light intensity Io is incident on the 2 layer material structure shown in
Assuming there is no absorption and no scattering in the materials it can be assumed that the intensity on the first interface is Io the incident light intensity. The light intensity of the light transmitted into the top layer of the material L1 is given by
L
1
=I
o(1−R1) (2)
Similarly the light intensity transmitted into the second layer L2 is given by
L
2
=L
1(1−R2)=Io(1−R1)(1−R2) (3)
And the light intensity being transmitted past the third optical interface L3 is given by
L
3
=L
2(1−R3)=Io(1−R1)(1−R2)(1−R3) (4)
In an interferometer which is set up in an optical autocorrelator configuration, the light that comes back from each optical interface interferes with light from each of the other optical interfaces. The signal coming back to the interferometer from the first optical interface S1 is given by
S1=IoR1 (5),
the signal coming back to the interferometer from the second optical interface S2 is given by
S
2
=I
o
R
2(1−R1)2 (6)
and the signal coming back to the interferometer from the third optical interface S3 is given by
S
3
=I
o
R
3(1−R1)2(1−R2)2 (7).
For the interfaces in
The complete interferogram for this type of sample is given by
where λ is the central wavelength of the light source and k and the rest of the relationships are derived below.
A treatment of interference of partially-coherent light is found in Fundamentals of Photonics, 1991 by B. Saleh and M. Teich. When two partially-coherent light beams interfere, the intensity of the combined beam I(x) as a function of distance x is given by:
I(x)=Is+Ir+2√{square root over (IsIr)}|gsr(x)|cos φ(x) (9)
where Is and Ir are the intensities of the individual light beams, gsr(x) is the normalized mutual coherence function and φ(x) is the phase difference between the two light waves. For NIR SLED light sources, the coherence function is Gaussian as a function of distance. For the case where the sample and reference beams are mutually coherent at location xo, the third (interference) term in equation 9 called S(x) can be written as:
where k is a constant which is related to the source coherence length. For a Gaussian distribution, the source coherence length (LC) is given by the expression:
where Δλ is the source spectral bandwidth. The coherence length defines the full width at half maximum of the Gaussian function in Equation 2. When x−xO=LC/2 the amplitude of the normalized Gaussian function=½. The value of k, which satisfies this relationship, is
For a 1300 nm source with a 60 nm bandwidth, the coherence length is calculated to be 12.429 μm and k=1.794747×1010/m2.
Of central importance for signal processing is the development of a true peak location algorithm. The goal is to find the true envelope center of an interferogram (a Gaussian function times a cosine function) when the data are not sampled at the location of the true Gaussian maximum. This must also be performed in the presence of noise from the environment. A variety of alternatives were evaluated including use of beats from multiple wavelength sources, or choice of sampling rate, moment calculations, Gaussian peak analysis, up-conversion, envelope detection, and Hilbert Transform method and Fourier transform phase analysis. The Fourier transform phase analysis technique enables calculating the thickness of thin organic films coated on either silicon or glass substrates in the range from 10 Angstroms (1 nm) up to a few microns in thickness. The Fourier transform phase analysis technique is based on applying the Shift Theorem to a discrete Fourier transform data set. An article by B. Danielson and C. Boisrobert, entitled “Absolute Optical Ranging Using Low Coherence Interferometry”, Applied Optics, 30, 2975, 1991 describes this approach. As taken from R. Bracewell, The Fourier Transform and its Applications, Second Edition, McGraw Hill Book Company, New York, 1978, the Fourier Shift Theorem can be stated as follows:
F(s)=∫−∞∞f(x)e−2πixsdx (13)
where s is the frequency variable and x is the position coordinate. The Fourier Transform shift theorem can be written as:
∫−∞∞f(x−a)e−2πixsdx=e−2πiasF(s) (14)
where a is the shift in the x coordinate. If δx is the sampling distance interval, P the calculated phase slope per point in the FFT centered around the frequency fo of maximum magnitude in the FFTs power spectrum, and N the number of points in the FFT, then it can be shown that:
The spatial frequency fo is calculated from the expression:
In order to use the phase slope algorithm, an initial guess is made as to the x axis location of each of the interferogram peaks. This is done by choosing the location of the absolute value of the maximum amplitude of each of the peaks indicating optical interfaces in the interferogram as the location of the initial guess. A 256 point subset centered around this initial guess is taken and the first 128 points shifted to the end of the 256 point data subset are taken such that the most intense interferogram points are located at the beginning and end of this subset. To reduce noise and improve precision, data points in the middle of this array are set equal to zero (zero filling). The number of zero-filled points is dependent upon the bandwidth of the light source. For a 1550 nm laser and 1300 nm SLED with 50 nm bandwidth, we typically zero fill the central 140 points of the shifted interferogram. The complex FFT of the zero-filled data array is taken and the resulting FFT values are transformed to polar coordinates (magnitude and phase). The center spatial frequency of the FFT is determined by locating the array index value corresponding to the data point having the maximum value of the magnitude spectrum. This frequency is checked for validity based upon expected frequency values obtained from equation (16). The center spatial frequency of the FFT is verified by determining if it falls within the acceptable range, and the phase slope calculation is performed by performing a linear least squares fit on the phase around the points centered on spatial frequency fo. Phase unwrapping is required if the phase angle exceeds the range from −π to +π. The phase measured at fo is used in equation (15) to calculate the true location of the peak by determining the shift δx from the initial guess location a. The distance between each set of adjacent peaks, gives the optical thickness of the substrate-plus-coating layer at the time the peaks are measured. This process is repeated during an entire etching rate monitoring sequence. In order to determine the thickness divide the measured optical thickness by the index of refraction of the layer material.
In order to apply using low-coherence interferometry for monitoring the rate of change of thickness during in-situ etching the thickness is monitored as a function of time. The rate of change of thickness of the layer is determined by using the peak locations of adjacent maxima determined at known different times, a first time τ1 and at a second time τ2, by the interferometer to measure total optical path which corresponds to the optical thickness of the substrate and the layer and subtracting the optical thickness of the substrate plus layer at the known different times and dividing by the difference in the known times (τ2-τ1) to determine the rate of change in the optical thickness of the layer. This corresponds to taking the derivative of the change in thickness as a function of time.
The measured optical thickness of a material will also change with temperature due to thermal expansion and the thermo optical effect. The observed temperature dependence is given by the thermo optic coefficient of temperature
given by
which is the sum of two terms where α is the thermal coefficient of expansion of the material and dn/dT is the change of group index of refraction with temperature. Note that the complete thermo optic coefficient is not a constant and is proportional to the thickness of the sample as shown in equation (17). It is also dependent on the temperature range since dn/dT will also depend on temperature. In silicon, there is a slight increase in dn/dT with temperature (see G. Cocorullo et al, Appl. Phys. Lett., Vol. 74, No. 22, 31 May 1999 entitled “Temperature dependence of the thermo optic coefficient in crystalline silicon between room temperature and 550 K at the wavelength of 1523 nm”. As shown in Example 1 low coherence interferometry can be used to measure the thermo optic coefficient of a material directly without the need of independently measuring the thermal coefficient of expansion and the change of group index of refraction with temperature. In Example 5, etching data have been corrected using the thermo optic coefficient.
In the examples shown below the interferometer shown in
This example shows the effect of temperature on optical thickness measured by low coherence interferometry using the apparatus described in
This example shows how low coherence interferometry has been implemented to measure in situ etch rates for a homogeneous material, in this case silicon. A silicon coupon in the 100-orientation (Si(100)) (709.7 sq mm, 0.35 mm thick) polished on both sides was mounted in the fixture shown in
This example shows how low coherence interferometry can be used to measure in situ etch rates for a homogeneous material, in this case borosilicate glass. A borosilicate glass coupon (709.7 sq mm, 1.1 mm thick) was mounted the fixture shown in
This examples shows how low coherence interferometry can be used to measure in situ etch rates at different temperatures for a homogeneous material using an opaque etchant. A silicon coupon (Si(100)) (709.7 sq mm, 0.35 mm thick) polished on both sides was mounted the fixture shown in
This example shows how temperature corrections can be applied to in situ etching measurements made by low coherence interferometry. From optical thickness data shown in
This example shows how in situ etch rates measured by low coherence interferometry can be used to measure temperature profiles for etching of materials. The experimental data presented in
Table 1 lists etching rates for a range of materials determined by low coherence interferometry using the fixtures and recirculation systems described in the previous examples. The experimental conditions including the etchant, the pH value of the etchant, the pressure in the etching chamber, and the temperature in the etching chamber have been provided in the table. The list includes the homogeneous materials (samples A, B, C, and D), materials comprised of a substrate and one coating (samples E, F, and G), and materials comprised of a substrate and two coatings (samples H and I). The etch rates in Table 1 are given for the substrate itself in the case of homogeneous materials (A, B, C, and D). For the materials with coatings, the etch rates are for the topmost coating; for a one layer coating, the topmost coating is coating 1 (E, F, and G), for a two-layer coating the topmost coating is coating 2 (H and I). For the homogeneous materials, the optical thickness is related to the material thickness by dividing the optical thickness by the refractive index of the material, as in Examples 2-4. For the coated materials, approximate coating thicknesses prior to etching have been provided in the table. For the coated substrates, the relationship between optical thickness and material thickness, including coating thickness, is determined by the complete interferogram (equation 8), which includes the refractive indices and thicknesses of the coatings and of the substrate. The substrates in Table 1 include substrates polished on both sides (double-side polish, dsp) and substrates polished on one side only (single-side polish, ssp). Si(100) is 100-oriented silicon metal. Si(111) is 111-oriented silicon metal. The etchants are as follows; 1015 FF is Kodak Versamark 1015 FF flush fluid and FR 1014 is Kodak Versamark FR 1014 replenisher fluid, and potassium hydroxide (1 M, KOH). Both 1015 FF and FR 1014 were obtained from Kodak Versamark. The 1 M KOH was prepared from potassium hydroxide pellets and water prior to use.
The data presented in Table 1 show that low coherence interferometry can be used to measure etch rates for many materials, including different homogeneous materials such as Si(100) (samples A and B), Si(111) (sample C), and quartz (sample D), as well as for coated materials with a variety of coatings, such as silicon nitride (sample G), silicon oxynitride (sample H), and different silicon oxide glasses (samples E, F, and I). The table includes data collected at elevated pressure, such as at 60 psi and temperature, such as at 88° C. (sample G). A large range of etch rates are also demonstrated, from low rates of a few nm/h to nearly 2000 nm/h (sample B). The ability to measure double-side polished, such as sample A, and single-side polished materials, such as sample B, is also demonstrated. This expands the utility of the method. Furthermore, the ability to measure both dsp and ssp samples shows that this method can be used even when signal intensities are significantly reduced, as is observed in ssp materials relative to dsp materials.
Taken together, the data presented in Table 1, considered with the previous examples shows how different materials respond to a variety of etchants. Consequently, the etch rate can be used to characterize the composition of a material.
The invention has been described in detail with particular reference to certain preferred embodiments thereof, but it will be understood that variations and modifications can be effected within the spirit and scope of the invention.
The present application is related to U.S. Ser. No. 11/262,868, filed Oct. 31, 2005, by Michael Alan Marcus et al., entitled “Measuring Layer Thickness Or Composition Changes”.