US5880838A  System and method for optically measuring a structure  Google Patents
System and method for optically measuring a structure Download PDFInfo
 Publication number
 US5880838A US5880838A US08659168 US65916896A US5880838A US 5880838 A US5880838 A US 5880838A US 08659168 US08659168 US 08659168 US 65916896 A US65916896 A US 65916896A US 5880838 A US5880838 A US 5880838A
 Authority
 US
 Grant status
 Grant
 Patent type
 Prior art keywords
 structure
 light
 beam
 amplitude
 incident
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Expired  Lifetime
Links
Images
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
 G01B11/00—Measuring arrangements characterised by the use of optical means
 G01B11/24—Measuring arrangements characterised by the use of optical means for measuring contours or curvatures
 G01B11/2441—Measuring arrangements characterised by the use of optical means for measuring contours or curvatures using interferometry

 G—PHYSICS
 G01—MEASURING; TESTING
 G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
 G01B11/00—Measuring arrangements characterised by the use of optical means
 G01B11/02—Measuring arrangements characterised by the use of optical means for measuring length, width or thickness

 G—PHYSICS
 G01—MEASURING; TESTING
 G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
 G01B9/00—Instruments as specified in the subgroups and characterised by the use of optical measuring means
 G01B9/02—Interferometers for determining dimensional properties of, or relations between, measurement objects
 G01B9/02015—Interferometers for determining dimensional properties of, or relations between, measurement objects characterised by a particular beam path configuration
 G01B9/02017—Interferometers for determining dimensional properties of, or relations between, measurement objects characterised by a particular beam path configuration contacting one object several times
 G01B9/02019—Interferometers for determining dimensional properties of, or relations between, measurement objects characterised by a particular beam path configuration contacting one object several times contacting different points on same face of object

 G—PHYSICS
 G01—MEASURING; TESTING
 G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
 G01B9/00—Instruments as specified in the subgroups and characterised by the use of optical measuring means
 G01B9/02—Interferometers for determining dimensional properties of, or relations between, measurement objects
 G01B9/02015—Interferometers for determining dimensional properties of, or relations between, measurement objects characterised by a particular beam path configuration
 G01B9/02024—Measuring in transmission, i.e. light traverses the object

 G—PHYSICS
 G01—MEASURING; TESTING
 G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
 G01B9/00—Instruments as specified in the subgroups and characterised by the use of optical measuring means
 G01B9/02—Interferometers for determining dimensional properties of, or relations between, measurement objects
 G01B9/02055—Interferometers for determining dimensional properties of, or relations between, measurement objects characterised by error reduction techniques
 G01B9/02075—Interferometers for determining dimensional properties of, or relations between, measurement objects characterised by error reduction techniques of particular errors
 G01B9/02078—Caused by ambiguity
 G01B9/02079—Quadrature detection, i.e. detecting relatively phaseshifted signals
 G01B9/02081—Quadrature detection, i.e. detecting relatively phaseshifted signals simultaneous quadrature detection, e.g. by spatial phase shifting

 G—PHYSICS
 G01—MEASURING; TESTING
 G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
 G01B2210/00—Aspects not specifically covered by any group under G01B, e.g. of wheel alignment, caliperlike sensors
 G01B2210/56—Measuring geometric parameters of semiconductor structures, e.g. profile, critical dimensions or trench depth

 G—PHYSICS
 G01—MEASURING; TESTING
 G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
 G01B2290/00—Aspects of interferometers not specifically covered by any group under G01B9/02
 G01B2290/70—Using polarization in the interferometer
Abstract
Description
The invention described herein was made in the performance of work under a grant from the Department of the Army, and is subject to the provisions of Public Law 96517 (35 U.S.C. §202), in which the Contractor has elected to retain title. The U.S. Government has certain rights in this invention pursuant to Grant No. DAAH0493G0192 awarded by the Army.
1. Field of the Invention
The present invention is a system and method for measuring the dimensions of a small structure, such as a microelectronic device. More particularly, the present invention is an optical measuring device that uses a polarized light beam, reflected off a very small structure, to measure the dimensions of the structure.
2. Background and Summary
Current microelectronics fabrication processes can produce structures as small as 0.25 microns. Indeed, structures even smaller than 0.25 microns will soon be fabricated commercially. (For the remainder of this description, the term "structure" will be used to mean any type of feature to be measured by the present invention, including grooves, pits, lines, bumps, and other three dimensional objects.)
The quality control inspection required at each step of the fabrication process becomes a problem with structures having such small dimensions. The small size allows many structures to be fabricated on a wafer, meaning that more layers of interconnections must also be fabricated on the wafer. Thus, the method for inspecting such devices must be fast and efficient in order to inspect every wafer on a fabrication line. In addition, the method must not damage the wafer.
Several methods currently exist for measuring such small structures, but these methods suffer from disadvantages that prevent their use or at least make them impractical for inspecting microelectronic wafers online. One such method employs a scanning electron microscope to measure the structures. This method is expensive and slow, however, because each wafer must be placed in a vacuum chamber, which must then be evacuated, before the measurement can be performed. Moreover, most electron microscopy techniques require that the wafer be cut so that the wafer's profile can be measured. These techniques, therefore, cannot be used to inspect wafers online.
Conventional optical microscope methods also cannot be used for such small structures. Such methods can barely resolve structures of this scale and thus do not provide accurate dimensional data.
Another method for measuring such small structures, which employs an atomic force microscope, is even slower than electron microscopy. This method, moreover, is impractical for structures with high aspect ratios, i.e., structures in which the height or depth is much larger than the smallest dimension. An example of a high aspect ratio structure is a contact hole etched through oxide layers for contacting a circuit device with metallization. Another example of a high aspect ratio structure is a capacitor trench etched into a silicon substrate for DRAM cell fabrication.
Optical diffraction methods have also been used to measure small microelectronic structures. These optical methods measure the scatter of plane waves and assume that the structures form gratings (i.e., periodic structures). This assumption, however, generally does not hold. Further, these methods cannot be used to measure individual or aperiodic structures.
Accordingly, the inventor recognized a need for an optical system and method for measuring small individual or aperiodic structures that is fast, efficient, and inexpensive and that overcomes the disadvantages and drawbacks of the prior art. The present invention provides such a system and method.
An object of the present invention is to provide an optical system and method for measuring a lateral dimension or a vertical dimension of a structure by focusing a polarized beam of light onto the structure and measuring the phase and amplitude of the reflected polarization components, where the lateral dimension is smaller than the wavelength of the light beam.
Another object of the present invention is to provide methods for determining the amplitude and phase of the reflected polarization components of the light beam.
Another object of the present invention is to provide methods for testing the stability and accuracy of the method of calculating the amplitude and phase of the polarization components.
Yet another object of the present invention is to provide an apparatus for detecting the reflected polarized light beam.
The system of the present invention employs a light source to generate a beam of light that is polarized and focused onto a structure to be measured. Preferably, the light beam is approximately monochromatic. The structure is illuminated with both a transverse electric (TE) polarized field and a transverse magnetic (TM) polarized field. The structure has a plurality of lateral dimensions (e.g., length and width) and at least one vertical dimension (e.g., depth or height). These lateral and vertical dimensions will be generally referred to throughout this description as "structural parameters." The term "structural parameter" should also be understood to include the type of material making up the structure, e.g., gold or silver. One of these structural parameters is preferably substantially larger than the wavelength of the incident light beam. Another one of the structural parameters is small enough so that, when the incident light is diffracted off the structure, the amplitude and phase of the diffracted TM and TE are affected differently by the structure. Preferably, the light beam has polarization components substantially parallel (TE) and perpendicular (TM) to the larger structural parameter.
The focused light beam is diffracted off the structure and is detected by a detector. The TM and TE fields are affected differently by the diffraction off the structure. The TE field can be used as a reference to analyze the phase and amplitude changes in the TM field. The relationship between phases and amplitudes of the TE and TM fields is dependent on the structural parameters of the structure. Thus, the parameters of the structure can be determined by examining the amplitude and phase relationship between the reflected TE and TM fields.
The system uses beam splitters, wave plates, and detector elements to measure the reflected light beam. Preferably, the light beam propagates through a polarizer and a beam splitter and is focused onto the structure. The light beam then preferably is reflected off the structure and back through the beam splitter. Alternatively, the light beam can be transmitted through the structure. The inphase and quadrature components of the TM field, using the TE field as a reference, are then measured. The system may have an inphase and quadrature leg, so that both components can be measured in parallel. Alternatively, the system may have only one leg for measuring the inphase and quadrature components, with a time delay being employed to measure the components in series. Preferably, the system includes four separate detector elements, two each for the inphase and quadrature legs. The TE and TM polarized waves can be "interfered" in both the inphase and quadrature legs.
Another aspect of the invention involves determining the amplitude and phase of the diffracted electromagnetic far fields for each polarization component. One method for determining the amplitude and phase is by a numerical method, which is preferably embodied in a computer program and includes several steps. First, several variables are determined and entered into the computer. The variables include: the wavelength of the focused light beam, the index of refraction of the incident medium (e.g., air), the index of refraction of the diffracting structure, the numerical aperture of the focusing lens, a description of the shape of the structure, and a choice of the type of polarization (e.g., transverse electric or magnetic field polarization). In addition, a maximum extent is selected for the surface of the structure which confines the numerical calculation to a truncated region. Generally, the truncated region is substantially larger than the wavelength divided by the numerical aperture of the focusing lens. An arbitrary spacing of grid of points along the surface of the structure is selected, and a set of basis functions is selected.
Second, lookup tables may be created. This step increases the speed of the computer program, but is not necessary, as the calculations in the third step can be performed (albeit more slowly) without the lookup tables. One such lookup table may include entries representing the geometry of the structure's surface contour. Another lookup table may include the distances along the surface contour from a reference point to each grid point. A third lookup table may include the abscissas and weights used repeatedly in the numerical integrations performed in the next step.
Third, numerical integrations are performed for surface contours of the structure. Each integration provides the tangential electric and magnetic fields along the structure's surface contour at some point. The integrand includes basis functions (preferably, sinc functions, where sinc(x)=sinπx/πx)) and a Green's function. The Green's function is a function of several of the variables input in the first step, such as the wavelength of the light beam and the indexes of refraction of the incident medium and the refracting medium. The number of sinc functions is preferably equal to the number of grid points.
Fourth, a matrix is constructed. Each matrix element is a numerical line integral along the whole surface contour of the structure. Thus, each matrix element includes a numerical integration of a 2dimensional Green's function and a sinc function. The rows of the matrix represent the grid point on which the Green's function is centered, and the columns represent the grid point on which the sinc function is centered. The matrix accounts for electromagnetic fields and for the surface contour of the structure.
Fifth, the incident electromagnetic field at each grid point along the structure's surface contour is calculated, in known fashion. These incident field values are due only to the incident light beam at the surface contour and thus do not account for diffracted light. These incident field values become the elements of an incident field vector.
Sixth, the matrix is inverted and multiplied by the incident field vector to obtain the elements comprising a previously unknown vector. The unknown vector (a tangential electromagnetic field vector) accounts for diffraction of the light beam. The elements of the electromagnetic field vector represent the actual tangential electric and magnetic fields that are present along the surface contour.
Finally, the tangential electromagnetic field vector is used to calculate the diffracted electromagnetic far field. Both the amplitude and phase of the electromagnetic far field are a function of angle. It is not necessary to calculate both the electric and magnetic far fields. If either the electric or magnetic far field is calculated, as those skilled in the art will appreciate, the other far field can be derived using Maxwell's equations.
This numerical method provides the scattered electromagnetic far field for one polarization component (e.g., the TE component). The method can be repeated to obtain the scattered far field for the other polarization component. The final step is to calculate the phase difference in the Fourier Transform responses of the TE and TM fields and the corresponding inphase and quadrature response.
An alternative method for determining the phase and amplitude of the electromagnetic far field involves the use of a sample with structures of known dimensions. Here, a set of reference structures is taken, for which the structural parameters are (preferably) precisely known. This reference structure is then measured using the present invention, and the electromagnetic far field polarization responses for the reference structures are determined. Other structures having at least one unknown structural parameter can then be measured in accordance with this invention, and the electromagnetic far field polarization response for each such structure can be compared to that for the reference structures. The unknown structural parameters can then be determined by comparing the measured response of the unknown parameters with the measured responses of the known reference structures. The unknown structural parameters will be most similar to the dimensions of that known reference structure for which the polarization quadrature measurements are most similar.
The present invention also includes methods for testing the stability and accuracy of the numerical method described above. The preferred method for testing the stability is known as a convergence test. The convergence test involves repeating the method for the same structure with an increasing and/or decreasing number (and thus spacing) of grid points. The truncation size of the structure can also be varied while repeating the method. Preferably, however, the truncated length of the structure is kept relatively short to reduce the total number of grid points and thus the time it takes to perform all the calculations. The numerical program is stable if the numerical solution converges (i.e., if the result does not change significantly when varying the number of grid points). The numerical accuracy can be estimated from the convergence data.
Another method for testing the accuracy of the numerical method uses the Fast Fourier Transform (FFT) of surface currents along the contour. The use of sinc functions for the basis functions of the surface current makes the implicit assumption that the surface current is strictly bandlimited in space. The Fourier Transforms should show a cutoff frequency if the surface currents are in fact bandlimited. The highest frequency bin will correspond to the cutoff frequency of the bandlimited functions, if a calculation is made of the FFT of the surface current coefficients, which estimate the tangential electromagnetic field along the contour. The FFTs of the surface current coefficients are then examined to determine whether they have significant amplitudes in the high frequency bins (e.g., by comparing those amplitudes to some expected amplitude). If so, it can be assumed that there is aliasing (i.e., that grid sampling is not being performed fast enough) and that the grid spacing must be decreased.
These and other aspects of the invention will now be described in detail with reference to the accompanying drawings, wherein:
FIG. 1 shows a simple coordinate system for a plane parallel waveguide.
FIG. 2 is a graph plotting the pit phase depth for different pit widths of a structure, with the pit depth equal to a quarterwave.
FIGS. 3A3B show two embodiments of a detector for polarization quadrature measurement in accordance with the present invention, in which light is reflected off a structure.
FIG. 4 shows another embodiment of a detector for polarization quadrature measurement in accordance with the present invention, in which lenses are added to allow imaging as in a conventional microscope.
FIG. 5 shows another detector for polarization quadrature measurement in accordance with the present invention, in which light is transmitted through the structure.
FIG. 6 is a graph showing the TE/TM quadrature response for a groove, in which the depth varies from zero to four wavelengths.
FIG. 7 shows a diffraction model.
FIG. 8 is a flow diagram for the preferred numerical integration method for determining the amplitude and phase of the diffracted electromagnetic far field for each polarization component.
FIG. 9 shows a generalized geometry for diffraction used in the numerical integration method of the present invention.
FIG. 10 shows a diffracting surface of a structure with the incident field generated by distant sources.
FIG. 11 is a graph showing the coordinate systems for a twodimensional diffracting structure.
FIG. 12 is a graph showing the coordinate system for the far field.
FIG. 13 is a graph showing a groove contour of a supergaussian profile of the form z(x)=a exp(x/σ).sup.υ, where a=1, σ=0.5, and υ=14.
FIG. 14 is a plot of the successive change to the far field calculation as the solution for the numerical integration method converges.
FIG. 15 is a plot showing the amplitude of the Fast Fourier Transform of J_{s},m for two different grid spacings.
FIG. 16 shows an apparatus used in experiments for measurement of the diffracted far field for focused incident fields.
FIG. 17 is a CCD image of the reflected Fourier plane for an aluminum grating with TM incidence, showing the Wood's anomalies.
FIG. 18 shows plots of the measured far field and the far field calculated by the numerical integration method for a sinusoidal aluminum grating, with the Wood's anomalies evident as a null in the reflected field.
FIGS. 19A and 19B are scanning electron microscope photographs of a groove etched in GalliumArsenide and used in a diffraction measurement, where FIG. 19A is a top view of the groove showing its width and FIG. 19B is a side view showing the groove's shape.
FIG. 20 is a plot of a contour used to represent the groove shown in FIGS. 19A and 19B.
FIG. 21 shows plots of the measured far field and calculated far field for the groove of FIGS. 19A and 19B for the illuminating light centered on the groove and for the illuminating light away from the groove.
FIG. 22A and 22B show threedimensional far field intensity profiles as the incident light source is shifted away from the groove of FIGS. 19A and 19B, where FIG. 21A is calculated using the numerical integration method and FIG. 22B is measured data from the groove.
Throughout this description, the preferred embodiments and the examples shown should be considered as exemplars, rather than as limitations on the present invention.
I. Principles of the Invention
The system and method of the present invention involve an optical technique for measuring small structures, such as grooves and bumps on an integrated circuit or compact disc. The structure is illuminated with a beam of polarized, incident light that is diffracted off the structure. The structures have a plurality of lateral dimensions and at least one vertical dimension. As stated above, both the lateral and vertical dimensions of the structure will be generally referred to as "structural parameters." The term "structural parameter" may also refer to material (e.g., gold or silver) comprising the structure. Those dimensions may also be specifically referred to as lateral or vertical dimensions. The lateral dimensions may be the length and width of the structure, and the vertical dimension can be the height or depth of the structure.
At least one of the structural parameters is substantially greater than the wavelength (λ) of the light. For example, the length of the structure may be greater than the or equal to 2λ. In addition, at least one other structural parameter is small enough, in comparison to λ, such that the amplitude and phase of the TE and TM fields of the diffracted incident light depend on the polarization of the incident light. This dependence may result, for example, when one of the structural parameters is approximately the same size as, or less than, the wavelength of the incident light. It should be understood that this latter structural parameter could be the height, width, depth, etc., of the structure.
The present invention can be used to measure structures like the deep capacitor trenches fabricated in a silicon substrate of a dynamic random access memory (DRAM) process. Integrated circuit contact holes are another example of structures that can be measured using the present invention. Contact holes are etched through oxide layers for the purpose of contacting a circuit device with metallization. Both DRAM capacitor trenches and contact holes are structures that have at least one very small lateral dimension that can be accurately and efficiently measured using the present invention.
Contact holes have a high "aspect ratio." In general, aspect ratio refers to the ratio of a longer structural parameter to a shorter structural parameter. Thus, the aspect ratio increases as a long structural parameter increases and/or a shorter structural parameter decreases. Specifically, the width of a contact hole can be as small as 0.25 microns, while the depth is several microns, meaning that this exemplary contact hole has a high aspect ratio.
DRAM capacitor trenches also have a high aspect ratio. A DRAM cell consists of a capacitor and a transistor. Some DRAM processes form the capacitor by etching a deep and narrow trench, depositing a thin oxide along the walls of the trench, and then filling the trench with a conducting material. The capacitor is formed along the walls of the trench. Thus, the depth is significantly larger than the width, giving the DRAM capacitor a high aspect ratio.
Quality control is important for both contact holes and DRAM capacitor trenches to ensure that they are of the desired dimensions. The depth of the contact hole is important, because the hole must make contact with the proper metal layer. The depth and shape of the DRAM trench capacitor are important, and the uniformity of the oxide layer is also important. These structures can currently be measured by electron beam microscopy. Electron microscopy is slow and expensive, however, making it impractical for inspecting every wafer during a fabrication process. Moreover, many electron beam microscope techniques require cutting the wafer to measure the structure's profile. Current optical microscopes can barely resolve these structures, and they do not provide accurate depth information. Interference methods are ineffective at measuring structures with high aspect ratios and do not work well at measuring the depth of structures in which the resolution limit is larger than the short structural parameter of the structure. As a result, a defective wafer cannot be detected until electrical tests are performed later in the fabrication process.
The optical technique of the present invention can be used to measure small structures, such as DRAM capacitor trenches and contact holes, with high aspect ratios. The technique employs a beam of incident light that is focused onto a structure to be measured. The incident light is preferably approximately monochromatic and preferably has polarization components approximately parallel and perpendicular to the long structural parameter of the structure. The electromagnetic field component, where the electric field is approximately parallel to the long structural parameter of the structure, will be referred to throughout this description as the transverse electric (TE) polarization. The approximately perpendicular polarization will be called the transverse magnetic (TM) polarization. Those skilled in the art will recognize, however, that the TE and TM polarizations can be reversed, and that the polarizations need not be parallel and perpendicular to the long structural parameter. The effect of diffraction of the illumination with the structure modulates the phase, amplitude, and polarization state of the reflected light. This technique can be used to measure individual and aperiodic structures.
A basic physical principle of this measurement technique is as follows: The diffraction properties for the polarization components of the incident light beam can be predicted if it is known that at least one of the structural parameters is small enough, with respect to the wavelength of the incident light, that the amplitude and phase of the diffracted TE and TM fields depend on the polarization of the incident light. The amplitude and phase behavior of each polarization component can be predicted numerically or measured from a known sample, as will be described in detail below. Therefore, the parameters of the structure can be determined by measuring the amplitude and phase relationship between the diffracted polarizations.
This technique can be used during fabrication of integrated circuits to measure small, high aspect ratio structures. For example, the technique can be used in situ for etching processes. The system of the present invention can also be designed to perform coherent or incoherent imaging of the structure. Incoherent imaging can be accommodated through the present technique because the phase measurement is selfreferencing, allowing several structures to be measured in parallel.
II. Details of the System and Method
A simple example involving a plane parallel waveguide will first be explained to gain an intuitive understanding of the source of polarization differences. Then, this basic understanding is related to the measurement of aluminum grooves using the results of rigorous numerical diffraction computation. (The numerical computation method is described in detail below, in section III.) The system for detecting the diffracted light is then described, which detects the exact diffracted polarization amplitudes and phase differences. Finally, a description is provided for the results of a rigorous simulation of diffraction of focused deep ultraviolet (DUV) light focused onto a 0.25 micron silicon groove and the corresponding measured output.
A. Propagation of Polarizations in Waveguides
Two sets of allowable modes exist for electromagnetic propagation inside a plane parallel waveguide with perfectly conducting walls. The transverse electric (TE) modes have the electric field parallel to the walls. The transverse magnetic (TM) modes have the magnetic field parallel to the walls. FIG. 1 shows a plane parallel waveguide 102 with a simple coordinate system that is chosen for this analysis. The direction of propagation of the light beam 104 is in the zdirection, and the xdirection is perpendicular to the walls of a waveguide 102. The fields in this waveguide 102 consist of discrete modes, as restricted by the boundary conditions of the tangential electric field at the walls. The solutions for the electric field inside the waveguide for the TE case (electric field parallel to the wells) are: ##EQU1## where d is the width of the waveguide 102, and p is an integer. Alternately, the solutions when the magnetic field is parallel to the walls (TM case), are: ##EQU2##
Each mode, corresponding to an integer value of p, has a cutoff width, d_{c} ^{p} =pπ/k. The value k_{z} is imaginary for widths smaller than d_{c} ^{p}, and the field amplitude weakens exponentially in the zdirection. The TE solution for the mode p=0 is E_{y} =0; there can be no TE fields. There is, however, a TM solution for p=0, namely: H=e^{jkz}. The cutoff width for the p=1 mode is d_{c} ^{1} =λ/2, and both TE and TM have propagating solutions for this mode when d>d_{c} ^{1}. Therefore, when d<λ/2, only evanescent TE waves can exist, while the TM zero mode can still propagate. Here lies the essential difference between the TE fields and the TM fields when considering the interaction of light with small trench structures. The TE fields can no longer penetrate the trench when the trench width falls below λ/2, while the TM fields can still propagate down the trench and back. The basic principle of the present invention is that, by measuring the amplitude and phase of one polarization with respect to another, the parameters of the trench can be deduced.
B. Aluminum Grooves with Varying Width
The above consideration of the waveguide explains why it might be expected that light of different polarizations would interact differently with an etched trench. Yet, it is clearly an abstraction and does not indicate what the diffraction of an illuminating spot focused on the trench will be. The present invention includes a rigorous, numerical method for calculating the diffracted electromagnetic fields for light focused onto arbitrary two dimensional structures. (This numerical method, which is described in detail below in section III, is preferably emodied as a computer program.)
Grooves of varying width in aluminum will be analyzed to provide an example of optical diffraction. An actual aluminum groove could be a contact hole on which a thin film of aluminum is sputtered. The groove phase depth in this example was calculated as a function of groove width for a groove having a depth of λ/4 to study the effect of changing the groove width.
FIG. 2 is a plot of the pit phase depth for different pit widths, where the pit depth is equal to a λ/4. The phase depth is defined as the phase of the electric field at the bottom of the groove. It is expected, from scalar diffraction considerations, that the phase depth of a groove would simply be the groove depth compared to the wavelength of the illuminating light. The results depicted in FIG. 2, however, show that the polarization state has a significant effect on the diffraction. This is apparent in the region where the pit width varies between λ/2 and λ, where λ is adjusted for the index of refraction of the incident medium. The phase depth of the TM wave in that region remains fairly constant, while the phase depth of the TE wave increases close to the quarterwave depth of the pit. Therefore, the pit width appears to be crucial to the penetration of the TE wave into the pit.
C. Polarization Quadrature Measurement (PQM)
The present invention utilizes a method referred to in this description as "polarization quadrature measurement" or "PQM." The fundamental idea behind PQM is to measure the amplitude and phase of the TM diffracted far field with respect to the TE field. The previous simulation showed that TE polarized light and the TM polarized light each see a structure (e.g., a groove) with a different effective depth when illuminating the groove. Therefore, respective far fields of the TM and TE polarized lights will have different amplitudes and phases. The amplitude and phase of the diffracted far field depends on the phase depth of the groove and the proportion of incident light that is diffracted off the bottom of the groove and the surface around the groove. The numerical method described below in section III can also be used to calculate the diffracted far field for each polarization.
Normally, the phase of an optical field cannot be detected without interfering it with a reference optical field. The diffracting structure is illuminated in PQM measurement with both TE polarized and TM polarized fields. The TE field can be used as a reference to analyze the phase and amplitude change in the TM field, because the TM and TE fields are differently affected by the diffraction. Those skilled in the art will recognize that, conversely, the TM field can be used as a reference to analyze the TE field. Measuring the relationship between the two orthogonally polarized fields is similar to the principle of conoscopic holography, which is known in the art and described in Sirat & Psaltis, "Conoscopic Holography," Optics Letters, Vol. 10, p.4 (1985).
The advantage of using the TE field as a reference, rather than constructing a traditional interferometer by splitting off part of the TM field before the diffraction and interfering the two parts, is that coherence length, stability, and vibrations are not an issue. Another advantage is the ease of measuring both the inphase and quadrature components. To measure both inphase and quadrature components of an electromagnetic field using traditional interferometry would require a separate reference leg with an optical path length onequarter wave different from the first reference leg.
Consider the electric field at the center of the far field propagating in the zdirection (with the zdirection oriented as shown in FIG. 1). It will be assumed that a single plane wave is present, because, in reality, an aperture can always be placed on the detector in the Fourier plane to restrict spatial components. The electric fields and corresponding detector voltages used in the following formulation, if desired, can be considered as functions of far field angle. In addition, the convention will be used that the TE reflected field is polarized in the xdirection, so that the TM electric field points in the ydirection. Consequently, the total electric far field is:
E=α.sub.ε ε.sup.jφ.sbsp.ε x+α.sub.m e.sup.jφ.sbsp.m y Eq. 3
where the factor of e^{j}(kzωt) has been ignored. The convention will also be used that the long structural parameter of the groove is in the xdirection. Thus, the TE electric field is α.sub.ε ε^{j}φ.sbsp.ε, and the TM electric field is α_{m} ε^{j}φ.sbsp.m
FIGS. 3A3B show two embodiments of a PQM detector system 300 of the present invention, in which the incident light beam 302 is reflected off the structure 304 to be measured. The reflecting system 300 includes a polarized light source 306, a lens 308 for focusing the incident light beam 302 onto the structure 304, and a nonpolarizing beam splitter 310. The system also includes a quadrature polarization detector 312.
FIG. 3A shows a relatively simple embodiment of the PQM system 300 for measuring the reflected polarization states of the incident light beam 302. The system 300 includes a polarized light source 306 that generates a beam of polarized light that is directed onto a collimating lens 314. The polarized light is transmitted through the collimating lens 314 onto a nonpolarizing beam splitter 310, which redirects the light beam onto the objective lens 308. The objective lens 308, in turn, focusses the light beam 302 onto the structure 304 being measured. The light beam then reflects off the structure 304, back through the objective lens 308, and through the beam splitter 310. The light beam then passes through an iris 316, through a quarterwave plate 318, and through a halfwave plate 320. The polarizer/analyzer is after the halfwave plate 320 and before a detector element 322. The light then enters the detector element 322. The inphase and quadrature components are derived from four detector measurements. Each detector measurement corresponds to a different setting of the λ/2 plate 320 and the λ/4 plate 318. These four settings lead to Equations 4, 5, 6 and 7(a) and 7(b).
Another embodiment of the PQM detector system 300 is shown in FIG. 3B, in which both the inphase and quadrature components of the TM far field are detected using the TE field as a reference. The beam of polarized light in this system is generated by the polarized light source 306 and is transmitted through a beam splitter 310. The light then passes through an objective lens 308, which focusses the light beam 302 onto the structure 304. The light beam then reflects off the structure 304, back through the lens 308, and is diffracted through the beam splitter 310. The light then enters another nonpolarizing beam splitter 324, where the light is split into inphase and quadrature phase components. The inphase leg 326 of the system includes a polarizing beam splitter 328 rotated 45° with respect to the zaxis, as shown in FIG. 3B. This rotated beam splitter 328 "interferes" the TE and TM polarized waves, which are then detected by two detector elements 330a, 330b.
The quadrature leg 332 also includes a polarizing beam splitter 334 rotated 45° with respect to the zaxis as well as a quarterwave plate 336. The quarterwave plate 336 retards the TE component by 1/4 wave with respect to the TM component. As a result, the electric field in the quadrature leg 332 is:
E.sub.2 =α.sub.ε ε.sup.j(φ.sbsp.ε.sup.+π/2) x+α.sub.m e.sup.jφ.sbsp.m y Eq. 4
The π/2 phase shift in the reference results from the quarterwave retardation and leads to the quadrature detection. The sign bias due to measuring the intensity in both legs 326, 332 is removed by subtracting the complementary outputs of the polarizing beam splitters 328, 334. The voltages at detectors 330a and 338a are:
V.sub.A,C αE.sub.1,2 ·(x+y)/√2.sup.2 Eq. 5
and the voltages at detectors 330b and 338b are:
V.sub.B,D αE.sub.1,2 ·(x+y)/√2.sup.2 Eq. 6
The differential amplifiers 340, 342 in each leg 326, 332 subtract the bias terms inherent in measuring intensity. The resulting inphase and quadrature voltages appearing at the outputs of the differential amplifiers 340, 342 are:
V.sub.I =V.sub.A V.sub.B =α.sub.e α.sub.m cos (φ.sub.e φ.sub.m) Eq. 7a
V.sub.Q =V.sub.C V.sub.D =α.sub.e α.sub.m sin (φ.sub.e φ.sub.m) Eq. 7b
As long as a_{e}, a_{m}, and φ_{e} φ_{m} remain stable, this PQM detector 300 will measure both the inphase and quadrature components of the TM field with respect to the TE field. The terms a_{e} or a_{m} can be normalized by adding to the detector a third leg with an untilted polarizing beam splitter, if those terms are somehow noisy. This causes a_{e} and a_{m} to be measured independently, and V_{I} and V_{Q} can be normalized. There is no simple way to also account for instability in the phases of the far fields. No reason exists, however, for the TE and TM fields to become incoherent, because those fields travel identical paths.
FIG. 4 shows another embodiment of a PQM detector 400, illustrating only one leg 402. The detector 400 has a light source 306, a collimating lens 314, a nonpolarizing beam splitter 310, and an objective lens 308, as in the detectors of FIGS. 3A and 3B. The light beam 302 reflects off a structure 304, back through the beam splitter 310, into an exit pupil (or iris) 316, and into a detector assembly 404. The detector assembly 404 has a tilted beam splitter 406 and two detector elements 408, 410. Lenses 412, 414 are added between the beam splitter 406 and the detector elements 408, 410 to allow the PQM detector 400 to image like a conventional microscope. As a result, the PQM response can be imaged onto a detector in an array of detectors, allowing measurement for as many structures in parallel as can be handled by the optical imaging system and detector array.
Many alternative embodiments exist for the PQM detector 300. For example, the nonpolarizing beam splitter 324 and one polarizing beam splitter (328 or 334) can be eliminated if the quarterwave plate is cut in half and only retards the xdirection polarization of half the beam. Then, the detector planes would also have to be cut in half; one half would be for the inphase signal and the other for the quadrature signal.
Another embodiment of the PQM detector is a transmission system 500, as shown in FIG. 5. There, the detector assembly 502 is placed behind a sample with a diffracting structure 504. A polarized light source 506, a collimating lens 508, and an objective lens 510 are placed in front of the structure 504. The light beam 512 is focused onto the structure 504, through which the light beam is transmitted to an objective lens 514 on the other side of the structure 504. The beam then passes through an iris 516, and through a quarterwave plate 518 and a halfwave plate 520. A detector 522 detects the polarized light. The transmitting detector 500 allows the first nonpolarizing beam splitter 324 of FIGS. 3A and 3B to be eliminated.
D. Considerations and Simulation for Silicon Structures
The previous example in section II.B of diffraction from a groove was for an aluminum structure with visible red light. Silicon is only about 30% reflective at the wavelength of visible red light. As a result, the walls of the silicon structure absorb a great deal of the light. Thus, most of the incident light of a red spot that is focused onto a silicon structure having a depth greater than about one wavelength is never reflected back out of the structure. There can be no PQM signal if the structure absorbs all of the incident light.
A solution is to illuminate the groove with deep ultraviolet (DUV) light. As an example, a KrF excimer laser emits light at a wavelength of 248 nm, and a high pressure mercury arc lamp can emit light at 257 nm. Both of these sources are commercially available and currently being used in microelectronic lithography. Silicon is very conductive at these wavelengths. In addition, its reflectivity increases to about 70%. Another benefit of using the DUV light is that the numerical aperture requirements of the illuminating system can be relaxed, because the illuminating spot size is approximately the wavelength of the DUV light divided by the numerical aperture.
A simulation was performed for the TE/TM quadrature detector response to a silicon groove having a width of 0.25 micron, where the groove depth is increasing. The simulation was performed with light having a wavelength of 248 nm, corresponding to a KrF excimer laser, because silicon is very reflective at DUV wavelengths. Illuminating the groove with light at a wavelength where silicon is most reflective gives the best measurement signal, because the reflected light is being measured. The index of refraction of silicon at a wavelength of 248 nm is about 1.6+j3.6.
FIG. 6 shows the TE/TM quadrature response for the simulation for the 0.25 micron wide groove, in which depth is varied from 0λ to 4λ (0.992 micron). The depth was incremented in steps of λ/4. The numerical aperture for the incident light was 0.6. The numerical aperture for collecting the reflected light was 0.1. The horizontal axis of the graph shows the inphase voltage, V_{I}, from the PQM detector, and the vertical axis shows the quadrature voltage output, of V_{Q}. Each of the plotted points, represented by a triangle, is for a different trench depth in the calculation. Thus, the (V_{I}, V_{Q}) coordinate corresponds to the detector response for that trench depth.
The graph of FIG. 6 shows how the phase difference between the diffracted TE and TM far fields rotates as the trench deepens, causing a spiral curve. The starting point is for the groove depth equal to 0λ, and the end point is for the groove depth equal to 4λ. Both polarizations have identical reflection when no structure is present, resulting in maximum inphase voltage and zero quadrature voltage. (This point is labeled "depth=0".) The TE electric field does not penetrate as far into the groove as the TM component as the groove deepens. Therefore, the reflected TE field does not suffer as large a phase retardation as the reflected TM field. In addition, the amount of light of each polarization that is absorbed by the structure is different, resulting in different reflected amplitudes for the two polarizations. Each point in the graph represents the voltage pair (V_{I} and V_{Q}) for a different groove depth.
The voltage pairs spiral counterclockwise as the groove deepens, because the TM field is phase retarded more than the TE field. The spiral is inward because deeper grooves absorb more light, reducing the reflected amplitudes. These differences of absorption and phase are measured by the PQM detector 300, resulting in the voltage pattern shown.
The operator of an integrated circuit assembly line could precisely stop the etching process at a particular groove depth by watching the detector output if the measurement of the trench could be made in situ. An instrument would measure the PQM voltage pair for each trench if the measurement is not made in situ. The trench passes if the measured voltage pair equals the expected voltage pair. Otherwise, the etching process has been defective.
There is no requirement for a coherent source for the PQM measurement. The TE and TM fields travel exactly the same optical path. Therefore, coherence length is not an issue.
In fact, there may be an advantage to using an incoherent source. The reflection of each trench can be imaged in parallel with coherent imaging.
As described above, FIG. 4 shows a PQM detector system 400 for measuring the parameters of a trench. The added lenses 412, 414 before the detectors 408, 410 allow imaging as in a traditional microscope. The PQM response for each trench can be imaged onto a detector (such as the one shown in FIG. 4) in an array of detectors by placing the lenses 412, 414 in the PQM system. The PQM response can be measured in this way for as many trenches in parallel as the optical imaging system and detector array can handle.
Optical disc memories consist of marks, or pits, embossed on a plastic substrate and then covered with a thin aluminum film. The marks are arranged in circular tracks so that the rotation of the disc causes the marks to move under a tightly focused illumination spot. Currently, the marks are designed so that when one is illuminated, a small amount of reflected light is directed back toward the objective lens. All the incident light is reflected back when no mark is present. Therefore, each mark on a current optical disc memory has two possible stateshigh reflectivity and low reflectivity. Therefore, each mark represents one bit of information. The overall capacity of an optical disc memory, in terms of bits of information, is approximately the total number of marks that can fit on one disc.
If marks could be designed to have different effects on the two different polarizations (polarizations in this case being defined by the direction of the electric field being either parallel or perpendicular to the track direction), then polarization quadrature could be used to determine the presence of a mark. For example, instead of a mark having two possible statesmark or no marka mark could be one of four possible states. The four different states would each correspond to a different (V_{i}, V_{q}) pair in the polarization quadrature detector. In this case, each mark would now represent two bits of information, and the overall storage capacity of the disc would double.
In fact, marks with different effects on polarization already exist. The width of the pits used in the new digital video disc (DVD) is less than half the wavelength of the illuminating light. As described in section II.A above, when the width of a reflecting groove falls below half the wavelength, the phase of the reflected light for each polarization is differently affected by the pit's depth. Therefore, three different pit depths (plus the possibility of no pit) could be chosen which have the following polarization properties:
______________________________________ Parallel Perpendicular Polarization Polarization Reflected Power Reflected Power______________________________________no pit large large1st pit depth small large2nd pit depth small small3rd pit depth large small______________________________________
Since each mark can be one of four possibilities, each represents two bits of information. This basic principle can be enlarged so that each mark can be one of many possibilities, representing many bits of information. Each possible state for a mark will correspond to a different constellation point for the polarization quadrature detector. This scheme is very analogous to quadrature amplitude modulation schemes used in telecommunications.
III. Numerical Method for Determining the Amplitude and Phase of the Diffracted Far Fields
An accurate tool for modeling diffraction would be a useful tool for the present invention. An accurate diffraction model enables the design of structures with desirable far field properties. For example, patterns can be located that have far fields that are maximally separated. Such patterns might be useful in the design of an optical memory. Another use of a model is to explain various diffraction phenomena. Other uses include the inspection of submicron features fabricated in the microelectronics industry and for the analysis of optical memory formats.
FIG. 7 shows a diffraction model 700. The model includes substrate surface 702 with a known contour and an incident field 704 resulting from a beam 706 focused through an objective lens 708 with a known numerical aperture. Most numerical diffraction approaches, such as the coupledwave and coupledmode methods, require a periodic diffracting structure and a plane wave incident field. These approaches could be used for a focused incident spot on a single structure by making the grating period large and separately accounting for each plane wave that makes up the focused spot. The number of orders, however, required in the diffraction calculation, because of the large period, will cause a burden on time and memory of the computation. This computation time must be multiplied by the number of plane waves that are to be considered in composing the focused spot. Therefore, a more direct integration technique was pursued in the present invention. This technique will allow the diffraction to be calculated directly from the focused incident field and a single groove.
Other alternatives to the integral method are finite differences, finite difference time domain, and finite element analysis. The advantages of these methods are that numerical integrations are not required, and, for finite differences, the formulation is somewhat simpler. The disadvantages of these methods are that grid points are placed throughout the whole volume of the structure, rather than just on the surface, as in the integral method. The matrices become very large, although sparse, with such a large number of grid points. Other problems also arise with these methods, such as questions concerning the convergence properties and concerning how to restrict the structural boundaries, which are at infinity.
Thus, a numerical diffraction method was developed in the present invention based on coupled integral equations. The mathematical formulation begins with Maxwell's equations and derives expressions for the electric and magnetic fields inside a volume. Each expression involves a convolution of the tangential fields on the volume's surface with the electromagnetic Green's function. These same expressions are written a second time to use the tangential surface fields to calculate the fields outside the volume. Coupled integral equations are formed for the tangential fields themselves after matching boundary conditions by equating the tangential components of the fields on either side of the surface (inside and outside the volume).
Several methods of testing the results of this numerical integration method have been developed and will be described below. In addition, the results of this numerical method will be compared below to experimental diffraction measurements.
A. The Basic Steps of the Numerical Integration Method
The numerical method can be used to determine the amplitude and phase of the diffracted electromagnetic far fields for each polarization component. The numerical method is preferably embodied in a computer program and includes several steps, as shown in FIG. 8.
The first step 802 of the method shown in FIG. 8 is to determine several variables, which are preferably entered into the computer. The variables include: the wavelength of the focused light beam, the index of refraction of the incident medium (e.g., air), the index of refraction of the diffracting structure, the numerical aperture of the focusing lens, a description of the shape of the structure, and a choice of the type of polarization. An arbitrary spacing of grid of points along the surface of the structure is also selected, and a set of basis functions is selected. The maximum extent of the grid truncates the structure's surface.
A description of the shape of the structure is also determined in the first step 802. This is done using mathematical functions, for example, cosine or superGaussian functions. FIG. 10 shows that the light beam is traveling in the zdirection toward the structure. A good approximation of the surface contour can be obtained by selecting different values of x, selecting a function that best approximates the structure, and setting z=function(x). For example, z=cos(x) may be selected if that function best approximates the crosssectional shape of the structure. A superGaussian function works well for approximating a groove. Alternatively, some actual measurement of the structure could be taken.
An optional second step 804 involves creating lookup tables. This step 804 increases the speed of the computer program, but is not necessary, because the calculations in the integration step 806 can be performed (albeit more slowly) without lookup tables. One such lookup table may include entries representing the geometry of the structure's surface contour. For example, this table could include the geometry of the structure's surface contour in the x and zdirections. Another lookup table may include the distances along the surface contour from a reference point to each grid point. FIG. 11 shows these distances represented by "s", which is the distance along the structure from the zaxis to a point on the structure, where s is a function of x. A third lookup table may include the abscissas and weights used repeatedly in the numerical integrations.
Numerical integrations are performed in the third step 806. The numerical integrations are for surface contours of the structure. Each integration provides the tangential electric and magnetic fields at some point along the structure's surface contour. The integrand includes basis functions (preferably, sinc functions) and Green's functions. The Green's functions include various parameters, such as the wavelength of the light beam and the indices of refraction of the incident medium and the refracting medium. The integrand may also include the derivatives of the basis and Green's functions. These derivatives implicitly involve the geometry of the structure. The number of sinc functions is preferably equal to the number of grid points, where sinc(x) is centered on each grid point.
A matrix is constructed in the fourth step 808. Each matrix element is a numerical line integral along the whole surface contour of the structure. Thus, each matrix element includes a numerical integrations of 2dimensional Green's functions and the sinc functions. The rows of the matrix represent the grid point on which the Green's functions are centered. The columns represent the grid point on which the sinc functions are centered.
The incident electromagnetic field at each grid point along the structure's surface contour is calculated in the fifth step 810, in known fashion. These incident field values are due only to the incident light beam at the surface contour and thus do not account for diffracted light. These incident field values become the elements of an incident field vector.
The sixth step 812 involves solving for a previously unknown tangential electromagnetic field vector by inverting the matrix and multiplying it by the incident field vector. The tangential electromagnetic field vector accounts for diffraction of the light beam. The elements of the vector represent the actual tangential electric and magnetic fields that are present along the surface contour.
The seventh (and final) step 814 uses the tangential electromagnetic field vector to calculate the diffracted electromagnetic far field. The far fields must be measured because the electric and magnetic fields cannot be measured at the surface contour. The amplitude and phase of the electromagnetic far field are both a function of angle. Both the electric and magnetic far fields need not be calculated, however, because, if one such field is calculated, the other field can be derived using Maxwell's equations.
Performing steps 802814 provides the scattered electromagnetic far field for one polarization component. The method can be repeated (816) if the scattered far field for the other polarization component has not yet been calculated. The method is complete (818) when the far field for both components has been calculated.
B. Details of the Numerical Integration Method
The formulation for the integral method of numerical diffraction is described in detail by A. T. de Hoop, "Modern Topics in Electromagnetics and Antennas," PPL Conference Publication 13, peter Peregrinus, Ltd., Stevenage, England, 1977, Ch. 6, and J. Pasman, "Principles of Optical Disc Systems," Adam Hilger Ltd., 1985, Ch. 6. Specific operator matrices for the two polarizations in two dimensional diffraction problems are given in the appendix of J. G. Dil and B. A. Jacobs, "Apparent Size of Reflecting Polygonal Obstacles of the Order of One Wavelength," J.Opt.Soc.Am., Vol. 69, No. 7, 1979, at 950.
The numerical integration method of the present invention starts with Maxwell's equations for the fields with an arbitrary structural geometry, as shown in FIG. 9. The fields in each region are solved in terms of the tangential fields at the boundary between the two regions. The tangential fields can also be thought of as surface electric and magnetic currents, which act as sources for the scattered fields in each region. By integrating Maxwell's equations, the fields in each region are expressed with integral operators on the surface currents. The operators include convolution integrals of the surface currents with the free space Green's function. The boundary conditions are continuity of the tangential fields from one region to the other, and this condition is easily enforced with the surface currents. Finally, two coupled integral equations remain, because the unknown fields are set and expressed by the operators and surface currents on the surface.
The following identity is needed for the Fourier transform of the curl of a vector field. The Fourier transform of a vector field, A(r), is ##EQU3## It is implied that the inverse Fourier transform of A(k) returns a valid vector field only over V(r), because the Fourier transform of the vector field is defined only over a volume, V(r). That is, ##EQU4## An identity is now formed for the Fourier transform of the curl of a vector by using the vector identity for the curl of a scalar times a vector along with Gauss's theorem. First,
e.sup.jk·r (∇×A)=∇×(e.sup.jk·r A)∇e.sup.jk·r ×A Eq. 10
and then integrating and applying Gauss's theorem, the result is: ##EQU5## where n is the unit normal vector on S pointing away from V.
The method starts with Maxwell's equations for timeharmonic (e^{j}ωt dependence) fields in a linear, isotropic, homogeneous, and sourcefree medium. The method then moves to the Fourier domain and solves for the fields in all space in terms of the tangential fields at the surface of the structure. The inverse Fourier transform is taken at the end, using the convolution theorem. The result will be two coupled equations for the electric and magnetic fields in terms of integrals of the Green's function with the surface currents. The surface currents will act as sources, so that once these currents are known, the fields at any point in space can be calculated by integrating over the sources with the appropriate Green's function. Diffraction would not only occur on the surface, but also inside the volume, if the medium inside the volume were inhomogeneous. Volume source currents as well as the surface currents are required in such a case.
The interface between two homogeneous media is then considered. For example, an optical disc memory consists of pits in an aluminum coated disc with a polycarbonate covering. Aluminum is highly conductive even though the aluminum covering is thin, so that the fields will not penetrate very deeply. Thus, a volume filled with aluminum can be used in the calculation instead of having to consider the more difficult problem of a thin aluminum coating on polycarbonate. The incident medium in this example is the polycarbonate. There is no need to consider the volume source currents, because both the aluminum and polycarbonate are homogeneous. The formulation of the present invention will only use the surface currents.
A complex ε is used, as a convention, to account for induced current in the medium. If ε=ε'+jε", then jωε'E+J=(jωε'+σ)E=jωεE, with σ=ωε". Taking the Fourier transform of Maxwell's equations,
∇×E=jwμH Eq. 12a
∇×H=jwεE Eq. 12b
and using Eq. 11, we obtain: ##EQU6## The term r_{s} has been used to indicate the independent variable along the surface of the volume. The vector fields n×E(r_{s}) and n×H(r_{s}) are equal to the tangential components of E(r_{s}) and H(r_{s}), respectively, at the surface, except that they point tangentially to the surface and perpendicularly to E(r_{s}) and H(r_{s}). These tangential components can be thought of as surface currents, and the boundary conditions require continuity of these surface currents. These surface currents are defined, as a convenience, as:
J.sub.s,m (r.sub.s)=n×E(r.sub.s) Eq. 14a
J.sub.s,e (r.sub.s)=n×H(r.sub.s) Eq. 14b
The terms with surface integrals in Eqs. 13a and 13b can now be written as Fourier transforms of the surface currents. Substituting the definitions into Eqs. 13a and 13b, the following coupled equations are obtained:
J.sub.s,m (k)=jwμH(k)jk×E(k) Eq. 15a
J.sub.s,e (k)=jwεE(k)jk×H(k) Eq. 15b
The next step is to solve Eqs. 15a and 15b for E and H in terms J_{s},m and J_{s},e. Consider J_{s},m and J_{s},e as known. Then, Eqs. 15a and 15b form two equations with two unknowns. An equation can be written for E(H) by substituting Eq. 15a into Eq. 15b to eliminate H (or, alternatively, Eq. 15b into Eq. 15a to eliminate E) . Using the vector identity jk×jk×A=k^{2} A+jk(jk·A) and the two additional relations,
jk·H=(jwμ).sup.1 jk·J.sub.s,m Eq. 16a
jk·E=(jwε).sup. jk·J.sub.s,e Eq. 16b
formed by taking jk·Eqs. 15a and 15b, expressions are obtained for E(E). The resulting equations are: ##EQU7## In the above expressions, the notation k^{2} =ω^{2} με, which is a constant inside the volume, is substituted. This constant is not to be confused with k^{2}, which is the squared magnitude of the independent variable in Fourier space.
The final step is to take the inverse Fourier transforms of Eqs. 17a and 17b, making use of the convolution theorem. The inverse Fourier transform of the first term in each of Eqs. 17 is the Green's function for free space propagation: ##EQU8## where, as a convention, Im{k}>0 is chosen. The convolution operator is notated, for convenience, as: ##EQU9## The inverse Fourier transforms can now be performed on Eqs. 17 to obtain our desired result for the fields in terms of the surface currents. Using the Fourier transform identities,
ℑ.sup.1 {jk(jk·A(k))}=∇(∇·A(r))Eq. 20
and
ℑ.sup.1 {jk×A(k))}=∇A(r) Eq. 21
the inverse Fourier transform of Eqs. 17 are ##EQU10## The derivatives inside the convolution integrals cannot, in general, be moved in forming Eqs. 22a and 22b, because of the singularity of the Green's function. In fact, the singularity will only be of concern when rεS.
The diffraction problem is for a boundary between two optically different, homogeneous media. FIG. 10 shows a geometry consisting of two volumes, V_{0} and V_{1}. Note that the sources of the incident field are placed in V_{0}. A source free region was assumed in the initial statement of Maxwell's equations. Therefore, expressions can only be written for fields that are sourcefree in a particular region. Following de Hoop (cited above), the fields summarized in the following table will be used in this description:
TABLE 1______________________________________The Fields and Their SourceFree RegionsField description source free in______________________________________F.sub.1 total field in V.sub.1 V.sub.1F.sub.0 total field in V.sub.0F.sub.i incident field V.sub.1F.sub.s scattered field V.sub.0______________________________________
An expression is written for each field in the region where the field is sourcefree. The boundary conditions require equality of the tangential components of the total fields in each region, F_{0} and F_{1}, at the surface. Note, however, that F_{0} is not source free in any region. F_{0} can be decomposed into its incident, F_{i}, and scattered, F_{s}, parts. Each of these components has a source free region. The final result, after equating the tangential components of F_{0} and F_{1} at the surface, is a set of coupled integral equations for the unknown surface currents.
Eqs. 22 express the fields inside a volume in terms of operators on the tangential fields at a surface. The notation can be simplified significantly by using operators. Thus, Γ is the operator matrix, J is the vector of surface currents, and F is the vector of fields. Then, ##EQU11## and Eqs. 22 become: ##EQU12## In region 0, k_{0} =ω√μ_{0} ε_{0} , G_{0} (r;r_{s})=(1/4πrr_{s} )e^{jk}.sbsp.(0).sup.rr.sbsp.s.sup., and Γ.sup.(0) is the corresponding operator. Similarly, Γ.sup.(1) is the operator for region
A convenient formula is now obtained with which to express the fields in each region in terms of the tangential fields on the surface of the diffracting structure. Solving first for the fields in V^{1}, Eq. 26 is applied, providing the following result: ##EQU13## The fields in region 0 can be separated into an incident field and a scattered field, F_{0} (r)=F_{s} (r)+F_{i} (r). The incident field, by definition, has its sources somewhere in V_{0}, and it is the field that would be present if the medium in region 1 were optically identical to the medium in region 0. The sources for the scattered field are the surface currents on S. The form of Maxwell's equations for a sourcefree region was chosen in the formulation for the solution of the fields in a region (resulting in Eq. 26). Therefore, care must be taken as to how that formula is applied to the incident and scattered fields. The formula can be solved in the incident and scattered fields with the following equation, because the scattered field is sourcefree in V_{0} : ##EQU14## The minus sign for the operator in the above formula indicates retention of the convention that the surface normal points into V_{0}. The incident field, however, is not sourcefree in V_{0}. Thus, an expression cannot be directly written for the incident field in V_{0}. Instead, the incident field is expressed in V_{1}, ##EQU15## The operator Γ.sup.(0) can be used for the incident field in V_{1}, because the incident field was defined as the field that would be present in V_{1} if this medium had no optical contrast with the medium in V_{0}. Subtracting Eq. 29 from Eq. 28, and substituting F_{0} =F_{i} +F_{s}, the following relation is obtained: ##EQU16##
Eqs. 22 and 30 constitute a number of equations equal to four times the number of unknowns. Wirgin, "Optics Communications," Vol. 27, No. 2, 1978, at 189, and J. Pasman, (cited above), when solving for diffraction from periodic gratings for planewave incidence, use an approach that avoids both explicit integral equations and integrals across singularities. In the third part of each of Eqs. 27 and 30, the right hand sides are known. The fields and surface currents can be expanded into their Fourier components. In addition, a simple matrix equation results for the solution of the coefficients, because the grating of structures is periodic and only one planewave is incident. The last step is to use the coefficients in the forward matrix equation that results from the first parts of Eqs. 27 and 30 to solve for the diffracted planewave coefficients. In addition, this whole process can be performed with either the electric fields or the magnetic fields. Three sampling grids would have to be generated to use this method to solve the problem, where the diffracting surface is not periodic (i.e., has either a plurality of aperiodic structures or has a sinble structure) and the incident field is not a planewave. The three sampling grids are: (1) the surface S, (2) a surface somewhere in V_{1}, and (3) a surface somewhere in V_{0}. Then, basis functions for each grid would have to be chosen.
The approach followed by de A. T. de Hoop, (cited above), and J. G. Dil and B. A. Jacobs, (cited above), and the one that will be follow here, is to concentrate on the fields at the surface for the right hand sides of Eqs. 27 and 30. Selecting the tangential components of those Equations (i.e., multiplying those equations by nhatx), Fredholm integral equations of the first kind appear. The four integral equations, along with the boundary condition, need to be combined to form two equations with two unknowns, ##EQU17## and known right hand sides ##EQU18## The boundary condition, equality of the tangential components on the surface, are written as n×F_{0} =n×F_{1} =J.
One solution is to select the tangential components of Eqs. 27 and 30 and then add the two equations together (using Γ=n×Γ:
(Γ.sup.(1) Γ.sup.(0))J=1/2J.sub.1 +1/2J.sub.0 J.sup.(i)Eq. 31
Setting J_{1} =J_{0} =J (the boundary condition), the final coupled integral equations are obtained:
(Γ.sup.(1) Γ.sup.(0) 1)J=J.sup.(i) Eq. 32
for the unknown surface currents, J.
The integrals of Eq. 32 are numerically solved by first approximating the surface currents with basis function expansions. As a result, each of the unknown surface currents in Eq. 32 becomes a vector of coefficients. The vector of basis functions is placed with the operators, so that each of the four operators in Eq. 32 becomes a row vector of operators, where each element includes a convolution of the Green's function and a basis function. The expansion itself introduces no error if a basis function is chosen so that the expansion converges uniformly to the surface current. The expansion must be truncated, so that the numerical program can work with finite vectors. This truncation, however, introduces error.
The integrals in the operator matrix cannot be numerically evaluated until values for r are chosen. Choosing a value of rvec corresponds to choosing a point on the surface where fields should be calculated with the operator matrix. Values of r are chosen by representing the surface as a set of grid points. Weighting functions on the surface could be chosen and the fields calculated according to those weighting functions, in the most general approach, rather than calculating the fields at specific grid points on the surface. Weighting functions made up of a set of evenly spaced delta functions is equivalent to choosing a set of grid points on the surface and calculating the fields at these grid points. Solving the fields only at a fixed set of grid points does not introduce error, except that the incident field on the right hand side of the equation must also be sampled and truncated according to the weighting functions or grid points. The four operators in Eq. 32 become four matrices, because each operator for each basis function must be numerically calculated for each grid point. Each column of each matrix is for a particular basis function, and each row is for a specific grid point.
The diffracted fields can now be computed. Integral equations were formulated above in a manner independent of geometry and without any assumptions about the fields. A series of choices will now be made, however, that will tailor the numerical calculations to the specific needs of the present invention. The preferred embodiment of the present invention uses tightly focused incident fields and single (or aperiodic) deep diffracting structures with at least one structural parameter on the order of the wavelength of the light source. These structures are assumed to lie in an infinite half plane filled with a medium of complex refractive index n_{1}. The operator matrices will be rewritten after each simplification. Each operator matrix will be handled separately, with the labels ##EQU19##
The integrals can be immediately performed along the yaxis when the surface and all the fields are independent of a particular direction, say the y direction, because only the Green's function has any dependence on y. The Green's function becomes ##EQU20## which is the Hankel function of the first kind and order zero, and R=rr_{s} . The remaining integration required to calculate each matrix element is a contour integral in the xz plane.
The independent variable is better described as a scalar, s, because all the integrations are now contour integrations. The scalar, s, is the length along the contour starting from the point on the contour where x=0, which is shown in FIG. 11. The value of s at any point along the contour is calculated in the standard manner, namely, ##EQU21## The unit vectors τ×n=y and n, which point tangentially and normally, respectively, to the contour, will ease the notation for the formulation in two dimensions. The direction of τ has been chosen, so that τ×n=y.
Another consequence of the independence in the ydirection is that the problem splits into two independent polarizations: (1) E=Ey(TE) and (2) H=Hy(TM). These polarizations are independent because an incident field with one of these polarizations will scatter only to fields of the same polarization. Assuming, for example, the TE polarization, the magnetic field has the direction H=H.sub.τ τ+H_{n} n, and the surface currents have the following directions:
J.sub.s,m (s)=n(s)×E(s)y=E(s)τ(s)=J.sub.s,m (s)τ(s)Eq. 34a
J.sub.s,e (s)=n(s)×(H.sub.τ (s)τ(s)+H.sub.n (s)n(s))=H.sub.τ (s)y=J.sub.s,e (s)y Eq. 34b
The solution fields in Eqs. 22 have directions consistent with the TE field, with the surface currents fixed to these directions and independent of y. Therefore, the integral equations are completely selfconsistent with only TE fields present. The principle is the same for the TM polarization.
For the reduction of the integral equations, TE polarization will be assumed. The matrix of operators can be simplified by knowing the specific directions of the surface currents. Specifically, it can be seen that in Γ_{12} there is a term, ##EQU22## The term y must be zero, because the assumption was made that all functions are independent in y. The following equations are obtained by substituting these surface currents into the matrix elements, ##EQU23## where it has been assumed that μ=μ_{1} =μ_{0}.
The surface current will now be expanded on a set of basis functions. In general, a function can be represented by the expansion, ##EQU24## where e_{l} (x) are the basis functions, and a_{l} are coefficients. The surface currents can be assumed to be approximately constant over small segments, making the basis functions rect functions (rect(x)=1 when x<1/2, and 0 otherwise) . This approach works well when the incident field has a small angular spectrum, such as a planewave. The number of grid points required to satisfactorily describe the incident field becomes very large, however, when the incident field is a tightly focused spot. This is the case even when the diffracting surface is a flat plane.
The incident field is known to have a large but finite angular spectrum. That the angular spectrum of the incident field is finite is guaranteed by the fact that the incident field comes from far away so that only propagating planewaves can reach the diffracting surface. A good choice of basis functions might then be sinc functions sinc(x)=sin(x)/x. The Nyquist sampling theorem states that a function of finite bandwidth (angular spectrum) is identically equal to a sum of evenly spaced samples weighted by sinc functions, if the sample rate is at least twice the highest frequency (planewave angle) of the function. That is, f(x)=Σ_{l} f(lπ/W)sinc(πWxπl), if W is greater than or equal to the bandwidth of f(x). As a result, the incident field and surface currents can be perfectly represented by a set of sinc basis functions if the diffracting surface is a flat plane. The incident field is no longer a simple Fourier transform of its angular spectrum, however, for diffracting surfaces that are not a flat plane.
Moreover, a sampling rate greater than the Nyquist frequency is normally required. This higher sampling rate corresponds to allowing for evanescent waves in other numerical approaches, such as coupledwave. Intuitively, the sampled rate required for good results would be expected to increase as the surface becomes less like a flat surface. The numerical results from the calculations of the present invention support this rule.
The disadvantage of sinc basis functions is that the integral in each of the matrix elements now has infinite limits. In contrast, the use of rect basis functions places the integral only over the region where the rect is nonzero. Thus, sinc basis functions are superior for flat and nearly flat surfaces. It is not entirely clear, however, which basis function (sinc or rect) is better for deep and narrow grooves in the diffracting surface.
The expansion of the surface currents is of the following form: ##EQU25## The unaliased bandwidth of the surface currents is determined by W. The coefficients, a_{l} and b_{l}, are samples of the surface currents at the grid points s=πl/k_{0} W. The expansion has been truncated to 2N+1 terms. The following Equation is obtained by substituting the basis functions into the matrix elements and placing the coefficients in the vector of unknowns: ##EQU26## The vector of unknowns is ##EQU27##
The derivatives in the matrix elements will now be resolved. The derivatives are with respect to the field point, r_{s}, not the independents variable of integration, s'. As described above, the order of differentiation and integration cannot simply be reversed, because the integration includes the point, s'=s, where the Hankel function has a logarithmic singularity of the nature, ##EQU28## where γ is Euler's constant. The singularity cancels because each integrand contains the subtraction, Γ_{1} (R)Γ_{0} (R), and the integrand has the finite limit: ##EQU29## This cancellation does not occur, however, for the first term of Γ_{21}.sbsb.l, where ##EQU30## which has an integrable singularity at R=0. Fortunately, this term does not involve a derivative. The term with the second derivative will have the integrable singularity in the TM case. Thus, the form of the singularity must be subtracted from the integrand before differentiating. (Subtracting singularities will be discussed in the next section.)
To calculate the derivatives, the chain rule is used, which is expressed as follows:
∇H.sub.0.sup.(1) (kR)=kH.sub.1.sup.(1) (kR)(∇R)Eq. 43
∇R=((xx')x+(zz')z)/R Eq. 44
The second term of Γ_{21}.sbsb.1 has two derivatives. Thus, the nature of the singularity of kH_{1}.sup.(1) (kR) must be examined. Fortunately, the subtraction in the integrand again causes the singularity to cancel: ##EQU31## Taking the second derivative of this term will produce a singularity similar to the singularity in the first term of Γ_{21}.sbsb.l. (Integrating across these singularities will be discussed in the next section.) The equations for the matrix elements are simplified by introducing the functions,
H.sub.0 (R)=k.sub.1.sup.2 H.sub.0.sup.(1) (k.sub.1 R)k.sub.0.sup.2 H.sub.0.sup.(1) (k.sub.0 R) Eq. 46a
H.sub.1 (R)=k.sub.1 H.sub.1.sup.(1) (k.sub.1 R)k.sub.0 H.sub.1.sup.(1) (k.sub.0 R) Eq. 46b
H.sub.2 (R)=k.sub.1.sup.2 H.sub.2.sup.(1) (k.sub.1 R)k.sub.0.sup.2 H.sub.2.sup.(1) (k.sub.0 R) Eq. 46c
Eqs. 39 result in the following by completing the derivatives in those equations: ##EQU32## All the integrands above can be numerically calculated with the aid of a few additional relations. The following relations are needed for the terms with n(s) and τ(s): ##EQU33## Also,
τ(s)·∇(τ(s')·∇R)=1/R τ(s).multidot.τ(s')(τ(s')·∇R)(τ(s)·.gradient.R)! Eq. 50
which is easily calculated. If the contour has a region where dz/dx becomes infinite (or large enough that numerical roundoff errors will cause problems), then Eqs. 48 and 49 must be replaced in those regions with forms using dx/dz and ds/dz.
Finally, the units of the electric and magnetic fields will be set by choosing √μ/ε_{0} =1. This choice fixes ωμ=k_{0}. Eqs. 44, 46, 47, 48, 49, and 50 now give the completed recipe for calculating the integrands for any pair of s, s' on the contour.
Each matrix element for each m,l pair can be calculated by choosing as one of the grid points, s_{m} =πm/k_{0} W, and by numerically integrating the appropriate term from Eq. 47 along the contour. Any of the various numerical quadrature routines would work well in principle. The simplest approach could be to use a standard Gaussian quadrature routine to integrate from grid point to grid point, because the grid points are evenly spaced in s. With this idea, all the grid points, quadrature points, and their x, z, dz, ds, and Hankel function values can be tabulated at the start of the calculation, for example, in lookup tables. The numerical integration cannot be performed to the infinite limits. The integration, however, need not go too far before the error becomes small, because H_{0}.sup.(1) (kx)sinc(kx)α(kx)^{3/2} sin(kx)e^{jkx} as x→0.
The singularity in the matrix elements Γ_{} {21_{} {lm} must now be integrated. These singularities occur only when m=1. Otherwise, s_{m} falls on a null of the sinc function, and the sinc function goes to zero sufficiently fast to eliminate the singularity. One method to integrate those diagonal terms with the singularity is to subtract the form of the singularity and integrate it separately, as is known. The term becomes: ##EQU34## The subtraction is not necessary outside the interval containing the singularity, a<s'<b. In the formula, R_{l} indicates that R is calculated from the grid point s_{l} to s', and s_{l} =s_{m} is the grid point for both the basis function and the source point. The TM case is treated slightly differently because a singularity also appears in the term with the double derivative. The correct approach is to subtract the singularity before differentiating, and then differentiating the extra asymptotic terms along with the rest.
The right hand side of the general matrix Eq. 32 is the tangential component of the incident field calculated at the grid points s_{m}. Many different models for the incident fields can be used. Restrictions must be placed, however, on the incident field if the numerical problem is to remain practical. As noted above, one source of error is the truncation of the grid points. The truncation is only reasonable if the power flow across the surface at the missing grid points is negligible. Therefore, only incident fields with finite power concentrated over a relatively small region can be allowed. For example, plane waves cannot be handled with this method, but focused spots are permitted.
FIG. 7 shows a lens 708 focusing light onto the diffracting surface 702. A perfect lens transforms the light distribution at its front focal plane into a distribution of plane waves. If B(u) represents the amplitude and phase of the TE electric field at the front focal plane of the lens, where u is the transverse coordinate in this plane, then the incident field behind the lens 708 is: ##EQU35## The incident magnetic field follows from Maxwell's equations, ##EQU36## The integrals for the incident field are easily calculated numerically, as long as x and z are sufficiently small. The following relationship was used in most of the simulations: ##EQU37## where NA corresponds to the numerical aperture of the lens. The incident surface currents for the right hand side of the matrix equation are then
J.sub.s,e.sup.(i) (s)=n×E.sup.(i) =τ(s)E.sub.y.sup.(i) (x.sub.s,z.sub.s) Eq. 55a
J.sub.s,m.sup.(i) (s)=n×H.sup.(i) =y(n.sub.x (s)H.sub.z.sup.(i) (x.sub.s,z.sub.s)n.sub.z (s)H.sub.x.sup.(i) (x.sub.s,z.sub.s))Eq. 55b
The incident fields for the TM case are completely analogous: H_{y} (x,z) is calculated from the front focal plane distribution, and the electric field components are derived from the curl of H.
The matrix equation is solved using a standard numerical routine, once the matrix elements and the incident field vector are calculated. The solution vector contains the coefficients a_{1}, b_{1}, which can be substituted into Eq. 38 to find the actual surface currents anywhere on the contour. Unfortunately, the tangential fields at the contour cannot be directly measured. Useful results of this diffraction calculation would be the amplitude and phase in the far field, because these quantities can be experimentally measured. The first part of Eq. 30 is used to calculate the scattered electric field (magnetic field in the TM case). A far field approximation for the Hankel function is used for this calculation. Also, the integrand in Eq. 30 is simplified by ignoring terms that have a dependence on r stronger than 1√r. The electric far field can then be calculated by a simple numerical integration.
The Hankel function for large arguments will first be examined. That function is as follows: ##EQU38##
FIG. 12 shows a coordinate system for the far field calculation. The far field in FIG. 12 is r>>r_{s}, so that using cylindrical coordinates for r, the following is obtained: ##EQU39## The Green's function for the far field is obtained by substituting this value of R into the asymptotic form of the Hankel function: ##EQU40## Equation 4.23 for the TE case is: ##EQU41## The integration is not across a singularity, because the field point is now in the far field, and the curl may be resolved directly: ##EQU42## where terms with r^{3/2} dependency have been ignored. The electric far field is determined by making the appropriate substitutions and ignoring the e^{jk} _{0} ^{r} √r factor: ##EQU43##
The integration in Eq. 61 is easily performed numerically because all the grid points and surface currents have already been tabulated. Values for θ can be chosen arbitrarily, although values that lie within the objective lens are generally used. The far field intensity is proportional to the magnitude squared of the electric field.
C. Methods for Testing the Results of the Numerical Integration Method
Several methods are known in the art for testing whether the numerical integration method is both numerically stable and an accurate representation of the mathematical formulation. These testing methods include comparing numerical results against analytic solutions for the small class of known solutions, checking that the numerical solution satisfies the appropriate power conservation laws, testing reciprocity relations with the numerical method, testing the numerical results with the extinction crosssection theorem, and checking the convergence of the numerical results for decreasing grid spacing and increasing matrix rank.
Different checks are useful at different stages of forming the numerical integration method. For the problem of a planar interface, a problem for which analytic solutions are known, only the matrix elements of Γ_{12} and Γ_{21} are nonzero. The planar interface problem can then be used to scrutinize the calculation of the matrix elements Γ_{21}. These elements are the most complicated calculations and, therefore, the most likely to contain errors. Power conservation relations are useful when assembling the various components of the overall calculation. For example, power conservation would not follow if a sign error were present in one of the incident fields. These two tests (i.e., the planar interface problem and power conservation relations) ensure that the numerical integration method is an errorfree representation of the formulated equations. These two tests do not, however, in general, test the numerical stability and accuracy of the computer program.
Convergence testing is the best method to evaluate the numerical stability of the numerical integration method of the present invention. The convergence test involves simply performing the method repeatedly for the same problem with an increasing number of grid points. (As described above, the numerical integration method is preferably embodied in a computer program, which can be run repeatedly with an increasing number of grid points for convergence testing.) The numerical program is stable if the numerical solution converges in a reasonable manner, at least for the particular problem that was computed. The numerical accuracy of the solution can be estimated from the convergence data.
Another testing method that will be discussed is a new check specifically for the parameter W, the scale of the sinc functions, and, thus, the inverse of the grid spacing. This method uses the Fast Fourier Transform (FFT) of the surface currents to check the validity of the assumption that they are bandlimited to W.
The power conservation rule is simply that the time averaged power flow transmitted across the contour boundary plus the power reflected from the boundary is equal to the power in the incident field. This method can only use incident fields with finite power, unlike formulations that use a plane wave as the incident field. First, it will be shown that the power flow across a boundary can be calculated with the surface currents. Second, the formulas for calculating the transmitted, reflected, and incident power from the surface current coefficients will be derived. This is done by applying the finite power incident field formulation and the orthogonality of the sinc functions.
The time averaged electromagnetic power flow across a surface is found by integrating the normal component of the real part of the Poynting vector, as is known: ##EQU44## The integration in the y direction is ignored, because the problem has been formulated in two dimensions. The normal component of the Poynting vector, however, can be represented using surface currents. The following relation is obtained from the definition of the surface currents:
J.sub.s,m ×J*.sub.s,e =(n) n·(E×H*)! Eq. 63
The directions of the surface current vectors are already known from the two dimensional treatment. For TE:
J.sub.s,m τ×J*.sub.s,e y=J.sub.s,m J*.sub.s,e n Eq. 64
and the time averaged power flow across the surface is ##EQU45## Note that the sign of the integrand has been changed so that the power flow is from V_{0} into V_{1}, opposite the direction of n.
The power flow in terms of the surface coefficients is obtained by substituting the basis function expansions of the surface currents, Eqs. 38, into the expression for the average power of Eq. 65: ##EQU46## This last formula results from the orthogonality relation, ##EQU47## where l and l' are integers and δ_{ll'} =1⃡l=l'.
The power calculated in Eq. 66 gives the total power flow across the boundary, which is the transmitted power. To test the power conservation relation,
<P>=<P.sup.(i) ><P.sup.(r) > Eq. 68
The incident power, <P.sup.(i) >, and reflected power, <P.sup.(r) >, must be calculated. The incident power is simply the power flow across the boundary for the case where no optical contrast exists between the two regions. In that case, a_{l} =a_{l}.sup.(i) and b^{l} =b_{l}.sup.(i). Substitution of these coefficients, already calculated as the right hand side of the matrix equation, into Eq. 66 gives <P.sup.(i) >.
The tangential fields must first be separated into incident and scattered fields in order to calculate the reflected power. The total fields in V_{0}, as stated in the formulation of the integral equations, can be separated into the incident fields and the reflected scattered fields. The boundary condition ensures that the surface currents are representations at the boundary of the tangential fields of V_{0} as well as V_{1}. Therefore, the surface currents can also be separated into incident and reflected components:
J.sub.s,m (s)=J.sub.s,m.sup.(i) (s)+J.sub.s,m.sup.(r) (s) J.sub.s,e (s)=J.sub.s,e.sup.(i) (s)+J.sub.s,e.sup.(r) (s) Eq. 69
The orthogonality of the basis functions gives the following equation when the surface currents for the reflected fields are expanded on the same sinc basis set:
a.sub.l.sup.(r) =a.sub.l a.sub.l.sup.(i) b.sub.l.sup.(r) =b.sub.l b.sub.l.sup.(i) Eq. 70
The reflected power can now be calculated using the coefficients a_{l}.sup.(r) and b_{l}.sup.(r) in Equation 4.59.
Convergence of the numerical results means that, if a calculation is repeated several times with smaller grid spacing each time, then the numerical results will become more accurate with each repetition. Here, the convergence should be tested for both decreasing grid spacing (increasing W) and increasing number of grid points with a fixed grid spacing. The latter test is employed because this invention is applied to individual or aperiodic structures, rather than periodic structures (e.g., gratings), and the integral is approximated to infinity by integrating only as far as the grid points extend.
Convergence testing is important and is preferably performed before accurate results can be obtained. There is a tradeoff between accuracy and computation time in every numerical calculation. A grid spacing can be chosen that gives satisfactory accuracy for a reasonable computation time, by considering the convergence properties. Naturally, the convergence estimation itself will require relatively long computation time because the test should extend to grid spacing smaller than required. Also, the test should be carried to grid spacings small enough to determine if the results are oscillating or converging uniformly. Finally, the numerical results for the surface currents and the far field will converge at different rates for different contour points or far field angles, depending on the particular contour.
The choice of sinc functions for the basis functions of the surface current allows for another type of convergence check. As described above, expanding a function on a basis set of sinc functions makes the implicit assumption that the function is strictly bandlimited. Also, the coefficients, a_{l}, b_{l}, are samples of the surface currents at the points s=πl/k_{0} W.
The Fourier transforms of the surface currents can be calculated based on exactly knowing the functions everywhere. The Fourier transforms should show a cutoff frequency, if the surface currents are in fact bandlimited. The highest frequency bin will correspond to the cutoff frequency of the bandlimited functions, if the FFTs of the surface current coefficients are calculated. Therefore, one can reasonably assume that aliasing is occurring and that the grid spacing should be shortened (i.e., that the bandwidth, W, should be increased), if the FFTs of the surface current coefficients have significant amplitudes in the high frequency bins.
This example simulates measurement of the depth of a groove etched in silicon. Semiconductor fabrication can involve etching single pits and grooves into a substrate, usually silicon, where the lateral dimension is as small as a quarter micron and the depth is one micron or more. Conceivably, a useful quality control measurement would be to measure the diffraction of a spot of laser light focused onto the silicon groove or pit. The measured diffraction could then be compared to a "perfect" example and a quality determination made. The sensitivity of this type of measurement could be numerically estimated by calculating the changes in diffraction patterns for grooves with varying parameters.
A single groove 1 micron deep and 1 micron wide was selected for this example to illustrate the convergence properties of the numerical integration method. An objective lens with a numerical aperture of 0.6 focuses light onto the groove and captures the reflected far field with the same aperture. The wavelength of the light is 0.6328 micron. The index of refraction of silicon at this wavelength is approximately 3.88+j0.02. The numerical integration method of the present invention does have a requirement that the boundary have a continuous first derivative. Therefore, a rectangular groove with sharp corners cannot be used. In consequence, a supergaussian profile is used. The superGaussian profile is of the form, z(x)=aexp((x/σ).sup.υ). FIG. 13 shows the groove shape 1302, in which the contour was truncated at s=±4 micron for these calculations. The groove contour 1302 is for a=1, σ=0.5, and υ=14.
TABLE 2 below shows the convergence for the reflected, transmitted, and total power for an increasing number of grid points.
TABLE 2______________________________________POWER CONSERVATION TEST FOR THE INTEGRAL METHODCONVERGES FOR A SINGLE GROOVE IN SILICON# of GridPoints R T R + T______________________________________83 3.171 4.1515 0.980299 0.2984 0.7019 1.0002115 0.3006 0.6996 1.0002131 0.2998 0.6993 0.9991147 0.3000 0.6994 0.9994165 0.3005 0.6991 0.9997______________________________________
The reflected and transmitted power values for the calculation with 83 grid points are both nonphysical, but the sum is actually close to one. Also, R and T oscillate as they converge when more grid points are used.
FIG. 14 shows a plot of the successive change to the far field calculation as the numerical solution converges, as an absolute limitation to the value of using the numerical integration method to predict the far field intensity. Plot 1402, which is very different from the other plots 1404, 1406, 1408, 1410, 1412 is for the trial with 83 grid points (i.e., the curve N=83 in FIG. 14). The other far field plots (i.e., N=99, 115, 131, 147, 165) converge, but not uniformly for all the points. Notice the strong oscillation of the far field as sin(θ) approaches the numerical aperture of 0.6. These oscillations are similar to the Gibbs phenomenon for truncated Fourier transforms. They appear here because of the similarity between the far field calculation (Equation 4.54) and a Fourier transform. Also, the integration for the far field calculation is truncated as the contour is truncated. The Gibbs phenomenon oscillations can be reduced by weighting the contour integration for the far field calculation with an appropriate window function.
The final check is the FFT of the surface currents. FIG. 15 shows the amplitude of the FFT of J_{s},m for two different values of the grid spacing, W. A sample rate of W=5 corresponds to 83 grid points, and W=10 corresponds to 165 grid points. The highest contour spatial frequency for the first case is 5k_{0}, and the FFT shows large amplitudes there. Therefore, the surface currents are highly aliased with this sample rate. The amplitude of the FFT goes to zero, however, for the greater sample rate of 10k_{0} at the higher contour spatial frequencies. The result of aliasing brings energy that should be at higher frequencies into the lower contour spatial frequencies, distorting the calculated surface currents. This principle explains why the far field calculation for 83 grid points was significantly different in shape from the others.
D. Experimental Verification of the Numerical Integration Method
The diffracted far field for each sample was measured in a similar manner to a Fourier plane measurements, as is known. FIG. 16 shows an detector apparatus 1600 used in a laboratory setup to measure the diffracted far field for focused incident fields. The detector apparatus 1600 includes an HeNe laser 1606 with polarized output, a spatial filter 1608, a first collimating lens 1610, an iris 1612, a first focusing lens 1614, a second collimating lens 1616, a beam splitter 1618, and an objective lens 1620. These elements are used to direct the light beam 1602 onto the test sample structure 1604. The detector apparatus 1600 also includes a reference mirror 1622, a second focusing lens 1624, a third collimating lens 1626, and a CCD camera 1628 to detect the reflected light beam 1602.
The iris 1612 is imaged onto the entrance pupil of the objective lens 1620, reducing the numerical aperture of the incident field. The full numerical aperture of the objective lens 1620 is used to capture the reflected field. The back focal plane, or Fourier plane, of the objective lens 1620 is imaged by lens 1624 and lens 1626 onto the CCD camera 1628, where an intensity profile of the reflected light is measured. The reference mirror 1622 is used to focus the illumination spot, and a shutter is closed when diffraction measurements are recorded.
The measurement of the far field is for a three dimensional spot, but the numerical calculation is only for two dimensions. A true comparison would be for the incident illumination being formed by a cylindrical lens rather than the spherical objective lens 1620 that is used. A crosssection of the measured far field is used, and the effect of the spot profile in the other direction is ignored. It can be argued that, if the test pattern is truly independent of one direction, no plane wave with an angle in that dimension can be diffracted to its angle in that dimension. Therefore, a cross section in the direction of interest should contain no diffraction of plane waves that would not be present if the objective lens 1620 were a cylinder.
The first test sample consisted of sinusoidal gratings in Aluminum. The gratings were fabricated by the Rochester Photonics Corporation (Rochester, N.Y.) using a rastered laser spot to expose a photoresist. A coating of approximately 500 Angstroms of aluminum was evaporated onto the grating after the photoresist was developed to form the grating. Finally, the actual depth and profile of the grating was measured with an atomic force microscope at the Jet Propulsion Laboratory, Pasadena, Calif. The surface measurement showed that the grating was quite sinusoidal with a peak to valley height of 0.185 micron. The grating period was 1.0 micron.
The Wood's anomaly for this grating should be quite evident. The angular location can be estimated where one expects to see strong absorption by calculating which plane wave angle is coupled into the plane wave traveling parallel to the surface by the grating. The Floquet condition, sinθ_{r} =sinθ_{i} +iλ/d, was used, where θ_{r} is the reflected angle, θ_{i} is the incident angle, and i is an integer. TABLE 3 below demonstrates that, by setting the reflected angle to sinθ_{r} =1, anomalous absorption is expected at the incident angles,
TABLE 3______________________________________ i sinθ.sub.i = 1  iλ/d______________________________________ +1 0.367 +2 0.266______________________________________
for this aluminum grating with λ=0.6328 μm. Consequently, notches in the reflection are expected to be measured at these angles.
FIG. 17 shows the CCD image of the reflected Fourier plane for an aluminum grating with TM incidence, clearly showing the Wood's anomalies. This image contains several interesting features. First, the numerical aperture for the incident spot was about 0.6, which is equal to the maximum sin(θ) of the incident plane waves. In addition, the +1 and 1 diffracted orders are shifted sin(θ)=λ/d =0.6328/1.0=0.6328, thus forming the two partial circles that do not quite meet in the center. The faint complete circle that fills in the background is the zero order reflected spot. Small bands at the top and bottom of the zero order spot are diffracted light that was captured because the numerical aperture for the collection of the reflected light was about 0.7. This numerical aperture is larger than the numerical aperture for the incident light. Finally, the two dark arcs within each diffracted order result from the absorption of the Wood's anomalies. The darker band is for the +1 angle in the Table 3, and the other is for the +2 angle.
The numerical simulation for this problem was calculated for a sinusoidal surface in a medium with index of refraction 1.5+j7.3. The optical field is only about 2.7% of the incident amplitude at a depth of 500 Angstroms, because the imaginary part of the index of refraction is so large. Therefore, it is reasonable to approximate the thin film of aluminum by filling the whole volume with aluminum.
FIG. 18 compares the numerical, calculated result to a vertical crosssectional measurement taken from the image for the aluminum grating. There is good agreement between the calculated results and the measured results at the smaller angles. The nulls in the reflected far field power due to the Wood's anomalies are quite evident and occur very close to the predicted angles. The small difference in the null angle between the numerical result and the measured cross section may be the result of a difference in the actual numerical aperture and the one used in the numerical calculation. This is because the horizontal scale of the measured data was scaled to the calculated data, assuming the numerical aperture used for each was identical. The vertical scale for each plot was derived by matching the measured and calculated reflection for a planar aluminum interface.
Experimental verification was also performed for a set of grooves etched into galliumarsenide. FIGS. 19A and 19B are photographs of one of the etched grooves and were taken by a scanning electron microscope. FIG. 19A is a top view of the groove and shows the groove width, which is about 0.35 micron. FIG. 19B shows a side view of the groove. The side view of FIG. 19B was obtained by forming a mask across the grooves and then etching the substrate away, thereby opening the groove to the side. The photograph of FIG. 19B was taken at a 60° angle. Direct measurement of the groove's depth from the photograph of FIG. 19B is difficult. The groove height was measured by comparing the groove to the large etch height, which could be measured with a surface profiler. The groove depth was measured to be about 0.65 micron.
FIG. 20 is a graphical representation of a surface contour that was used to estimate the shape of the galliumarsenide groove. The contour of FIG. 20 was constructed from a raised cosine pulse.
The measurement of the diffracted far field was performed by sliding the groove underneath the illuminating spot. The intensity profile at the CCD 1628 was recorded at each step of 0.1 micron. FIG. 21 compares the measured and calculated far fields for the illuminating spot centered on the groove and for the illuminating spot away from the groove. The overall change in intensity level between on and off the groove seems to agree between the calculated and measured data. The peak in the measured intensity in the middle of the far field is an artifact of light reflected from a lens surface that comes to a focus in the CCD plane; these peaks can be ignored. The reflected field calculated for the absence of a groove is simply the Fresnel reflection coefficients, and this reflection increases with far field angle for TE illumination.
FIGS. 22A and 22B are graphical representations comparing the measurement with the complete calculation for the groove shifting underneath the illuminating spot. The graph of FIG. 22A is calculated using the numerical integration method, while the graph of FIG. 22B is the measured data. As before, the peak running down the middle of the measured data of FIG. 22B can be ignored. Features in common with the calculated and measured data are the waves on the left side of the plots. One difference between the two plots is that the slope on the front left part is steeper for the measured data than for calculated data.
The numerical integration method described in section III is only one way to determine the amplitude and phase of the diffracted far field. Another method uses a sample having a set of reference structures, for which the structural parameters are already precisely known. This method will be referred to as the reference structure method.
The structural parameters of the reference structures are first measured using any known method. For example, a scanning electron microscope or atomic force microscope can be used to measure the reference structures.
The reference structures are then measured in accordance with the present invention, i.e., by the PQM technique explained in this description. Thus, the electromagnetic far field polarization responses for the reference structures are determined using PQM. The results of the PQM measurement for the structural parameters of the reference structures are then compared to the measurements taken for those parameters. The PQM method is thus calibrated, so as to account for differences between the measured and PQM results.
Next, a structure of unknown parameters is measured by the PQM method. The PQM results are then adjusted using the calibration factor derived from the differences between the measured and PQM results for the reference structures. An accurate measurement of unknown structures can thus be obtained without using the numerical integration method described above.
Although only a few embodiments have been described in detail above, those having ordinary skill in the art will readily realize that many modifications are possible without departing from the advantageous teaching therein. For example, several methods for testing the results of the numerical integration method have been described. It should be understood by those having ordinary skill in the art, however, that other methods could be used to test the results. Moreover, several embodiments of the apparatus for carrying out the PQM technique have been described. Those skilled in the art should understand, however, that other embodiments of this apparatus also exist. A hardware embodiment has been described above for carrying out various functions. Many advantages of the present invention are obtained from the hardware operation. However it should be understood by those having ordinary skill in the art that certain of these advantages and functions can also be carried out in software. For example, a special purpose digital signal processor could be appropriately programmed with firmware to carry out the functions of this assembly. Other similar modifications would also be apparent to those having ordinary skill in the art, and are intended to be encompassed within the following claims.
Claims (10)
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

US08659168 US5880838A (en)  19960605  19960605  System and method for optically measuring a structure 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

US08659168 US5880838A (en)  19960605  19960605  System and method for optically measuring a structure 
Publications (1)
Publication Number  Publication Date 

US5880838A true US5880838A (en)  19990309 
Family
ID=24644329
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

US08659168 Expired  Lifetime US5880838A (en)  19960605  19960605  System and method for optically measuring a structure 
Country Status (1)
Country  Link 

US (1)  US5880838A (en) 
Cited By (209)
Publication number  Priority date  Publication date  Assignee  Title 

US6160251A (en) *  19980605  20001212  Sandia Corporation  Method and apparatus for distinguishing actual sparse events from sparse event false alarms 
US6185323B1 (en) *  19971023  20010206  International Business Machines Corporation  Method characterizing a feature using measurement imaging tool 
US6304330B1 (en)  19991006  20011016  Metrolaser, Inc.  Methods and apparatus for splitting, imaging, and measuring wavefronts in interferometry 
US20020015146A1 (en) *  19970922  20020207  Meeks Steven W.  Combined high speed optical profilometer and ellipsometer 
US20020039184A1 (en) *  20001003  20020404  Sandusky John V.  Differential numerical aperture methods and device 
US6383824B1 (en)  20010425  20020507  Advanced Micro Devices, Inc.  Method of using scatterometry measurements to control deposition processes 
US6433871B1 (en)  20010525  20020813  Advanced Micron Devices, Inc.  Method of using scatterometry measurements to determine and control gate electrode profiles 
US6458605B1 (en)  20010628  20021001  Advanced Micro Devices, Inc.  Method and apparatus for controlling photolithography overlay registration 
US20020145740A1 (en) *  20010326  20021010  Meeks Steven W.  Material independent optical profilometer 
US20020163634A1 (en) *  20010326  20021107  Meeks Steven W.  Multiple spot size optical profilometer, ellipsometer, reflectometer and scatterometer 
US6479309B1 (en)  20010525  20021112  Advanced Micro Devices, Inc.  Method and apparatus for determining process layer conformality 
US6479200B1 (en)  20010419  20021112  Advanced Micro Devices, Inc.  Method of controlling stepper process parameters based upon scatterometric measurements of DICD features 
US6504612B2 (en)  20001114  20030107  Georgia Tech Research Corporation  Electromagnetic wave analyzer 
US20030025905A1 (en) *  20010326  20030206  Meeks Steven W.  Method of detecting and classifying scratches, particles and pits on thin film disks or wafers 
US6529278B2 (en) *  19991227  20030304  Canon Kabushiki Kaisha  Optical interference apparatus and position detection apparatus 
US6529282B1 (en)  20010611  20030304  Advanced Micro Devices, Inc.  Method of controlling photolithography processes based upon scatterometric measurements of photoresist thickness, and system for accomplishing same 
US20030043372A1 (en) *  20010830  20030306  Bernd Schulz  Apparatus for determining an overlay error and critical dimensions in a semiconductor structure by means of scatterometry 
US20030042579A1 (en) *  20010830  20030306  Bernd Schulz  Semiconductor structure and method for determining critical dimensions and overlay error 
US20030044702A1 (en) *  20010830  20030306  Bernd Schulz  Semiconductor structure and method for determining critical dimensions and overlay error 
US6535288B1 (en)  20010710  20030318  Advanced Micro Devices, Inc.  Machine readable code to trigger data collection 
US6537833B1 (en)  20010619  20030325  Advanced Micro Devices, Inc.  Method and apparatus for characterizing an interconnect structure profile using scatterometry measurements 
US20030058443A1 (en) *  19980306  20030327  Yiping Xu  Spectroscopic scatterometer system 
US6556291B2 (en)  20000531  20030429  Nidek Co., Ltd.  Defect inspection method and defect inspection apparatus 
US6556303B1 (en)  20010710  20030429  Advanced Micro Devices, Inc.  Scattered signal collection using strobed technique 
US6579733B1 (en)  20020116  20030617  Advanced Micro Devices, Inc.  Using scatterometry to measure resist thickness and control implant 
US6582863B1 (en)  20010611  20030624  Advanced Micro Devices, Inc.  Method of controlling photolithography processes based upon scatterometric measurements of subnominal grating structures 
US6597712B2 (en) *  20010226  20030722  Hitachi, Ltd.  Laser diode module 
WO2003060422A2 (en) *  20020121  20030724  Scinex Ag  Interferometric optical system 
US6602723B1 (en)  20010402  20030805  Advanced Micro Devices, Inc.  Method of integrating scatterometry metrology structures directly into die design 
US6602727B1 (en)  20010427  20030805  Advanced Micro Devices, Inc.  Scatterometry based active control of exposure conditions 
US20030147086A1 (en) *  20020205  20030807  Allan Rosencwaig  Analysis of isolated and aperiodic structures with simultaneous multiple angle of incidence measurements 
US6617087B1 (en)  20010627  20030909  Advanced Micro Devices, Inc.  Use of scatterometry to measure pattern accuracy 
US6618149B1 (en)  20010406  20030909  Advanced Micro Devices, Inc.  Method of identifying film stacks based upon optical properties 
US6623994B1 (en) *  20020926  20030923  Advanced Micro Devices, Inc.  Method for calibrating opticalbased metrology tools 
US6630362B1 (en) *  20010613  20031007  Advanced Micro Devices, Inc.  Method and apparatus for performing trench depth analysis 
US6639663B1 (en)  20010523  20031028  Advanced Micro Devices, Inc.  Method and apparatus for detecting processing faults using scatterometry measurements 
US6643008B1 (en)  20020226  20031104  Advanced Micro Devices, Inc.  Method of detecting degradation in photolithography processes based upon scatterometric measurements of grating structures, and a device comprising such structures 
US20030210760A1 (en) *  20020513  20031113  Nelson Alan C.  Method and apparatus for emission computed tomography using temporal signatures 
US6650423B1 (en)  20010702  20031118  Advanced Micro Devices Inc.  Method and apparatus for determining column dimensions using scatterometry 
US20030224261A1 (en) *  20020531  20031204  Bernd Schulz  Twodimensional structure for determining an overlay accuracy by means of scatterometry 
US20030223072A1 (en) *  20020531  20031204  Bernd Schulz  Scatterometer including an internal calibration system 
US6662125B2 (en) *  20010626  20031209  Fujitsu Limited  Electromagnetic wave analyzer and program for same 
US6660542B1 (en)  20010406  20031209  Advanced Micro Devices, Inc.  Method of controlling stepper process parameters based upon optical properties of incoming process layers, and system for accomplishing same 
US6665078B1 (en)  19970922  20031216  Candela Instruments  System and method for simultaneously measuring thin film layer thickness, reflectivity, roughness, surface profile and magnetic pattern in thin film magnetic disks and silicon wafers 
US20040017561A1 (en) *  19970922  20040129  Meeks Steven W.  Method of detecting and classifying scratches and particles on thin film disks or wafers 
US6704478B2 (en) *  20010823  20040309  Hitachi Cable, Ltd.  Wavelength separation optical device and multiple wavelength light transmission module 
US20040046968A1 (en) *  20000329  20040311  Jon Opsal  Critical dimension analysis with simultaneous multiple angle of incidence measurements 
US6708127B2 (en) *  20010315  20040316  Lucent Technologies Inc.  Beam propagation method for stepindex waveguides 
US6707562B1 (en)  20010702  20040316  Advanced Micro Devices, Inc.  Method of using scatterometry measurements to control photoresist etch process 
US6713753B1 (en)  20010703  20040330  Nanometrics Incorporated  Combination of normal and oblique incidence polarimetry for the characterization of gratings 
US6717671B1 (en)  19970922  20040406  Candela Instruments  System for simultaneously measuring thin film layer thickness, reflectivity, roughness, surface profile and magnetic pattern 
US20040070772A1 (en) *  20011219  20040415  Shchegrov Andrei V.  Parametric profiling using optical spectroscopic systems 
US20040076319A1 (en) *  20020419  20040422  Fauver Mark E.  Method and apparatus of shadowgram formation for optical tomography 
US6735333B1 (en) *  19980730  20040511  Nidek Co., Ltd.  Pattern inspection apparatus 
US6741730B2 (en) *  20010810  20040525  Visiongate, Inc.  Method and apparatus for threedimensional imaging in the fourier domain 
US6742168B1 (en)  20020319  20040525  Advanced Micro Devices, Inc.  Method and structure for calibrating scatterometrybased metrology tool used to measure dimensions of features on a semiconductor device 
US6757056B1 (en) *  20010326  20040629  Candela Instruments  Combined high speed optical profilometer and ellipsometer 
US20040133362A1 (en) *  20021021  20040708  Emit Systems Co.  Method for significantly increasing computational speed, accuracy and versatility of electromagnetic reconstruction of shapes and composition of complex targets containing lossy materials 
US6773939B1 (en)  20010702  20040810  Advanced Micro Devices, Inc.  Method and apparatus for determining critical dimension variation in a line structure 
US6785009B1 (en)  20020228  20040831  Advanced Micro Devices, Inc.  Method of using high yielding spectra scatterometry measurements to control semiconductor manufacturing processes, and systems for accomplishing same 
US20040169850A1 (en) *  19990702  20040902  Meeks Steven W.  System and method for double sided optical inspection of thin film disks or wafers 
US6790570B1 (en)  20011108  20040914  Advanced Micro Devices, Inc.  Method of using scatterometry measurements to control stepper process parameters 
US6791697B1 (en)  20020321  20040914  Advanced Micro Devices, Inc.  Scatterometry structure with embedded ring oscillator, and methods of using same 
US20040196469A1 (en) *  20030401  20041007  Seagate Technology Llc  Optical profiler for ultrasmooth surface with normal incident beam deflection method 
US6804014B1 (en)  20010702  20041012  Advanced Micro Devices, Inc.  Method and apparatus for determining contact opening dimensions using scatterometry 
US20040212812A1 (en) *  20020618  20041028  Timbre Technologies, Inc.  Optical metrology of single features 
US20040223137A1 (en) *  20011010  20041111  Littau Michael E.  Determination of center of focus by crosssection analysis 
US6836691B1 (en)  20030501  20041228  Advanced Micro Devices, Inc.  Method and apparatus for filtering metrology data based on collection purpose 
US20050010108A1 (en) *  20020419  20050113  Rahn John Richard  Method for correction of relative objectdetector motion between successive views 
US20050014043A1 (en) *  20030718  20050120  Matsushita Electric Industrial Co., Ltd.  Power supply unit 
US6859746B1 (en)  20030501  20050222  Advanced Micro Devices, Inc.  Methods of using adaptive sampling techniques based upon categorization of process variations, and system for performing same 
US20050041258A1 (en) *  20020313  20050224  Jon Opsal  Optical scatterometry of asymmetric lines and structures 
US6867866B1 (en) *  20010810  20050315  ThermaWave, Inc.  CD metrology analysis using green's function 
US6882437B2 (en)  20020419  20050419  KlaTencor Technologies  Method of detecting the thickness of thin film disks or wafers 
US6881594B1 (en)  20021028  20050419  Advanced Micro Devices, Inc.  Method of using scatterometry for analysis of electromigration, and structures for performing same 
US20050085721A1 (en) *  20020419  20050421  University Of Washington  System and method for processing specimens and images for optical tomography 
US20050085708A1 (en) *  20020419  20050421  University Of Washington  System and method for preparation of cells for 3D image acquisition 
US20050105089A1 (en) *  20010828  20050519  Jon Opsal  Detector configurations for optical metrology 
US6898537B1 (en)  20010427  20050524  Nanometrics Incorporated  Measurement of diffracting structures using onehalf of the nonzero diffracted orders 
US20050134847A1 (en) *  20031219  20050623  Ray Hebert  Apparatus and method for optical characterization of a sample over a broadband of wavelengths with a small spot size 
US20050134848A1 (en) *  20031219  20050623  Ray Hebert  Apparatus and method for optical characterization of a sample over a broadband of wavelengths while minimizing polarization changes 
US6933158B1 (en)  20021031  20050823  Advanced Micro Devices, Inc.  Method of monitoring anneal processes using scatterometry, and system for performing same 
US6944322B2 (en)  20010328  20050913  Visiongate, Inc.  Optical tomography of small objects using parallel ray illumination and postspecimen optical magnification 
US6949462B1 (en)  20020404  20050927  Nanometrics Incorporated  Measuring an alignment target with multiple polarization states 
US6972853B1 (en)  20020927  20051206  Advanced Micro Devices, Inc.  Methods of calibrating and controlling stepper exposure processes and tools, and system for accomplishing same 
US6980300B1 (en)  20010411  20051227  Advanced Micro Devices, Inc.  Method and apparatus for generating a polishing process endpoint signal using scatterometry 
US6982793B1 (en)  20020404  20060103  Nanometrics Incorporated  Method and apparatus for using an alignment target with designed in offset 
US20060012582A1 (en) *  20040715  20060119  De Lega Xavier C  Transparent film measurements 
US6989900B1 (en) *  20010402  20060124  Advanced Micro Devices, Inc.  Method of measuring implant profiles using scatterometric techniques 
US6992764B1 (en)  20020930  20060131  Nanometrics Incorporated  Measuring an alignment target with a single polarization state 
US6991738B1 (en)  20041013  20060131  University Of Washington  Flowthrough drum centrifuge 
US6992773B1 (en)  19990830  20060131  Advanced Micro Devices, Inc.  Dualdifferential interferometry for silicon device damage detection 
US20060023219A1 (en) *  20010328  20060202  Meyer Michael G  Optical tomography of small objects using parallel ray illumination and postspecimen optical magnification 
US20060028962A1 (en) *  20040805  20060209  Oron Zachar  Method and system for encoding and detecting optical information 
US20060033921A1 (en) *  20040816  20060216  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20060066878A1 (en) *  20040614  20060330  Samsung Electronics Co., Ltd.  System and method for measuring dimension of patterns formed on photomask 
US20060096358A1 (en) *  20041028  20060511  University Of Washington  Optical projection tomography microscope 
US7061615B1 (en)  20010920  20060613  Nanometrics Incorporated  Spectroscopically measured overlay target 
US20060126074A1 (en) *  20041214  20060615  Asml Netherlands B.V.  Inspection apparatus, sample and inspection method 
US20060139592A1 (en) *  20041122  20060629  Asml Netherlands B.V.  Latent overlay metrology 
US7075741B1 (en)  20040614  20060711  Kla Tencor Technologues Corporation  System and method for automatically determining magnetic eccentricity of a disk 
US20060158659A1 (en) *  20050120  20060720  Xavier Colonna De Lega  Interferometer for determining characteristics of an object surface 
US20060183220A1 (en) *  20041109  20060817  Nelson Alan C  Automated cell sample enrichment preparation method 
US20060192936A1 (en) *  20050225  20060831  Asml Netherlands B.V.  Lithography measurements using scatterometry 
US7113284B1 (en)  20010326  20060926  KlaTencor Technologies Corporation  Material independent optical profilometer 
US20060215289A1 (en) *  20040614  20060928  KlaTencor Technologies Corp.  Servo pattern characterization on magnetic disks 
US7115858B1 (en)  20000925  20061003  Nanometrics Incorporated  Apparatus and method for the measurement of diffracting structures 
US20060250609A1 (en) *  20050506  20061109  KlaTencor Technologies Corp.  Wafer edge inspection 
US20060250611A1 (en) *  20050506  20061109  KlaTencor Technologies Corp.  Wafer edge inspection 
US20060250612A1 (en) *  19970922  20061109  Meeks Steven W  Detecting and classifying surface features or defects by controlling the angle of the illumination plane of incidence with respect to the feature or defect 
US20060256325A1 (en) *  20050511  20061116  KlaTencor Technologies Corp.  Illumination energy management in surface inspection 
US20060256345A1 (en) *  20050512  20061116  KlaTencor Technologies Corp.  Interferometry measurement in disturbed environments 
US20060256324A1 (en) *  20050510  20061116  Asml Netherlands B.V.  Inspection apparatus and method of inspection 
US20060262291A1 (en) *  20050517  20061123  Hess Harald F  Dynamic reference plane compensation 
US20070003840A1 (en) *  20050630  20070104  Asml Netherlands B.V.  Focus determination method, device manufacturing method, and mask 
US20070002336A1 (en) *  20050630  20070104  Asml Netherlands B.V.  Metrology apparatus, lithographic apparatus, process apparatus, metrology method and device manufacturing method 
US20070064572A1 (en) *  20050727  20070322  Polarizonics Corporation  Optical data encoding and detection 
US7197355B2 (en)  20020419  20070327  Visiongate, Inc.  Variablemotion optical tomography of small objects 
US7201799B1 (en)  20041124  20070410  KlaTencor Technologies Corporation  System and method for classifying, detecting, and counting micropipes 
US20070086013A1 (en) *  20051011  20070419  Zygo Corporation  Interferometry method and system including spectral decomposition 
US20070093044A1 (en) *  20051025  20070426  Asml Netherlands B.V.  Method of depositing a metal layer onto a substrate and a method for measuring in three dimensions the topographical features of a substrate 
US20070115483A1 (en) *  19970922  20070524  Oak Dave S  Surface finish roughness measurement 
US20070153274A1 (en) *  20051230  20070705  Asml Netherlands B.V.  Optical metrology system and metrology mark characterization device 
US20070153273A1 (en) *  19970922  20070705  Meeks Steven W  Material independent profiler 
US20070182964A1 (en) *  20060209  20070809  Asml Netherlands B.V.  Lithographic system, sensor, and method of measuring properties of a substrate 
US7262864B1 (en) *  20010702  20070828  Advanced Micro Devices, Inc.  Method and apparatus for determining grid dimensions using scatterometry 
US20070206480A1 (en) *  20060303  20070906  Seagate Technology Llc  Polarization state detector systems and methods for readout of multilevel oriented nanostructure optical media 
US20070222960A1 (en) *  20060321  20070927  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method with reduced scribe lane usage for substrate measurement 
US20070229785A1 (en) *  20060328  20071004  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method using overlay measurement 
US20070229837A1 (en) *  20060329  20071004  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method using overlay measurement quality indication 
US7286229B1 (en)  20050906  20071023  KlaTencor Technologies Corporation  Detecting multidomain states in perpendicular magnetic media 
US7295300B1 (en)  20050928  20071113  KlaTencor Technologies Corporation  Detecting surface pits 
US20070262771A1 (en) *  20060511  20071115  KlaTencor Technologies Corporation  Method and System for Perpendicular Magnetic Media Metrology 
US20070279742A1 (en) *  20060531  20071206  Asml Netherlands B.V.  Metrology tool, system comprising a lithographic apparatus and a metrology tool, and a method for determining a parameter of a substrate 
US20070291269A1 (en) *  20060620  20071220  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20070296960A1 (en) *  20060613  20071227  Asml Netherlands B.V.  Method of characterising the transmission losses of an optical system 
US20070296973A1 (en) *  20060622  20071227  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20080002207A1 (en) *  20060630  20080103  Asml Netherlands B.V.  Apparatus for angularresolved spectroscopic lithography characterization and device manufacturing method 
US20080018901A1 (en) *  20060721  20080124  Zygo Corporation  Compensation of systematic effects in low coherence interferometry 
US20080030701A1 (en) *  20060803  20080207  Asml Netherlands B.V.  Individual wafer history storage for overlay corrections 
WO2008015230A1 (en) *  20060802  20080207  Commissariat A L'energie Atomique  Device for characterizing unique objects 
US20080037134A1 (en) *  20060810  20080214  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device method 
US20080036984A1 (en) *  20060808  20080214  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20080043239A1 (en) *  20060815  20080221  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20080049233A1 (en) *  20020909  20080228  Zygo Corporation  MultipleAngle MultipleWavelength Interferometer Using HighNA Imaging and Spectral Analysis 
US20080065350A1 (en) *  20030306  20080313  De Groot Peter J  Profiling complex surface structures using scanning interferometry 
US20080068609A1 (en) *  20060914  20080320  Asml Netherlands B.V.  Inspection apparatus, an apparatus for projecting an image and a method of measuring a property of a substrate 
US20080069430A1 (en) *  20060808  20080320  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20080068614A1 (en) *  20030915  20080320  De Groot Peter J  Methods and systems for interferometric analysis of surfaces and related applications 
US20080074666A1 (en) *  20060925  20080327  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20080088849A1 (en) *  20050120  20080417  Zygo Corporation  Interferometry for determining characteristics of an object surface, with spatially coherent illumination 
US20080088854A1 (en) *  20061013  20080417  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080117434A1 (en) *  20061122  20080522  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080128644A1 (en) *  20061130  20080605  Asml Netherlands  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080144036A1 (en) *  20061219  20080619  Asml Netherlands B.V.  Method of measurement, an inspection apparatus and a lithographic apparatus 
US20080148875A1 (en) *  20061220  20080626  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080151228A1 (en) *  20061208  20080626  Asml Netherlands B.V.  Scatterometer, a lithographic apparatus and a focus analysis method 
US7396022B1 (en)  20040928  20080708  KlaTencor Technologies Corp.  System and method for optimizing wafer flatness at high rotational speeds 
US7397553B1 (en)  20051024  20080708  KlaTencor Technologies Corporation  Surface scanning 
US20080170780A1 (en) *  20070116  20080717  Asml Netherlands B.V.  Inspection method and apparatus,lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080174753A1 (en) *  20070122  20080724  Asml Netherlands B.V.  Method of measurement, an inspection apparatus and a lithographic apparatus 
US20080174784A1 (en) *  20061222  20080724  Zygo Corporation  Apparatus and method for measuring characteristics of surface features 
US20080180656A1 (en) *  20070126  20080731  KlaTencor Technologies Corporation  Surface characteristic analysis 
US20080180685A1 (en) *  20070131  20080731  Zygo Corporation  Interferometry for lateral metrology 
US20080198380A1 (en) *  20070221  20080821  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080212097A1 (en) *  20070301  20080904  Asml Netherlands B.V.  Method of inspection, a method of manufacturing, an inspection apparatus, a substrate, a mask, a lithography apparatus and a lithographic cell 
US20080239265A1 (en) *  20070330  20081002  Asml Netherlands B.V.  Angularly resolved scatterometer, inspection method, lithographic apparatus, lithographic processing cell device manufacturing method and alignment sensor 
US20080239277A1 (en) *  20070327  20081002  Asml Netherlands B.V.  Method of forming a substrate for use in calibrating a metrology tool, calibration substrate and metrology tool calibration method 
US20080239318A1 (en) *  20070330  20081002  Asml Netherlands B.V.  Method of measuring asymmetry in a scatterometer, a method of measuring an overlay error in a substrate and a metrology apparatus 
US20080266574A1 (en) *  20030306  20081030  Groot Peter De  Interferometer and method for measuring characteristics of optically unresolved surface features 
US20080266560A1 (en) *  20060630  20081030  Asml Netherlands B.V.  Apparatus for angularresolved spectroscopic lithography characterization and device manufacturing method 
US20080279442A1 (en) *  20070509  20081113  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7460237B1 (en)  20070802  20081202  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20090015844A1 (en) *  20020909  20090115  De Groot Peter J  Interferometry Method for Ellipsometry, Reflectometry, and Scatterometry Measurements, Including Characterization of Thin Film Structures 
US20090021723A1 (en) *  20070719  20090122  Zygo Corporation  Generating model signals for interferometry 
US20090027691A1 (en) *  20060321  20090129  Asml Netherlands B.V.  Lithographic Apparatus and Device Manufacturing Method with Reduced Scribe Lane Usage for Substrate Measurement 
EP1610089A3 (en) *  20040621  20090218  Trioptics GmbH  Method and device for measuring the angle of optical surfaces 
US20090059236A1 (en) *  20070830  20090305  KlaTencor Technologies Corporation  Wafer Edge Inspection 
EP2053349A2 (en)  20071025  20090429  ASML Netherlands B.V.  Method and apparatus for determining properties of a lithographic pattern on a substrate 
US7532331B2 (en)  20060914  20090512  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20090128827A1 (en) *  20071113  20090521  De Groot Peter  Interferometer utilizing polarization scanning 
US20090147268A1 (en) *  20071012  20090611  Zygo Corporation  Interferometric analysis of underresolved features 
US7557934B2 (en)  20061207  20090707  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20090208072A1 (en) *  20080218  20090820  Visiongate, Inc.  3d imaging of live cells with ultraviolet radiation 
US7599064B2 (en)  20070307  20091006  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method, substrate for use in the methods 
US20090262349A1 (en) *  20060616  20091022  Institut National Polytechnique De Grenoble  Electrooptic probe for vector measurement of an electromagnetic field 
US7684032B1 (en)  20050106  20100323  KlaTencor Corporation  Multiwavelength system and method for detecting epitaxial layer defects 
EP2172934A1 (en)  20081006  20100407  Thomson Licensing  Optical storage medium comprising marks with different orientation, and respective apparatus for reading of data. 
US7700247B2 (en)  20031219  20100420  International Business Machines Corporation  Differential critical dimension and overlay metrology apparatus and measurement method 
US20100128280A1 (en) *  20081126  20100527  Zygo Corporation  Scan error correction in low coherence scanning interferometry 
US20100214639A1 (en) *  20090223  20100826  Visiongate, Inc.  Optical tomography system with highspeed scanner 
US7787112B2 (en)  20071022  20100831  Visiongate, Inc.  Depth of field extension for optical tomography 
US20100232664A1 (en) *  20090312  20100916  Visiongate, Inc.  Pattern noise correction for pseudo projections 
US7835561B2 (en)  20070518  20101116  Visiongate, Inc.  Method for image processing and reconstruction of images for optical tomography 
US7852459B2 (en)  20070202  20101214  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20110032500A1 (en) *  20080320  20110210  Asml Netherlands B.V.  Inspection apparatus for lithography 
US20110032534A1 (en) *  20090519  20110210  Camtek Ltd  System and a method for broadband interferometry 
US7907765B2 (en)  20010328  20110315  University Of Washington  Focal plane tracking for optical microtomography 
US7911612B2 (en)  20070613  20110322  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7916284B2 (en)  20060718  20110329  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20110137625A1 (en) *  20091208  20110609  Asml Netherlands B.V.  Methods and Apparatus for Determining Electromagnetic Scattering Properties and Structural Parameters of Periodic Structures 
US20110218789A1 (en) *  20090924  20110908  Asml Netherlands B.V.  Methods and Apparatus for Modeling Electromagnetic Scattering Properties of Microscopic Structures and Methods and Apparatus for Reconstruction of Microscopic Structures 
US8126677B2 (en)  20071214  20120228  Zygo Corporation  Analyzing surface structure using scanning interferometry 
US20120250026A1 (en) *  20110330  20121004  Boer Guido Do  Interferometer module 
US8294907B2 (en)  20061013  20121023  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US8468471B2 (en)  20110923  20130618  KlaTencor Corp.  Process aware metrology 
US8792096B2 (en)  20080324  20140729  Asml Netherlands B.V.  Inspection apparatus for lithography 
WO2014128710A1 (en) *  20130221  20140828  Nova Measuring Instruments Ltd.  Optical phase measurement method and system 
US20150241200A1 (en) *  20120927  20150827  Mapper Lithography Ip B.V.  Multiaxis differential interferometer 
Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

US4236823A (en) *  19780918  19801202  Rca Corporation  Diffractometer for measuring signal depth and width 
US5432607A (en) *  19930222  19950711  International Business Machines Corporation  Method and apparatus for inspecting patterned thin films using diffracted beam ellipsometry 
US5486919A (en) *  19920427  19960123  Canon Kabushiki Kaisha  Inspection method and apparatus for inspecting a particle, if any, on a substrate having a pattern 
Patent Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

US4236823A (en) *  19780918  19801202  Rca Corporation  Diffractometer for measuring signal depth and width 
US5486919A (en) *  19920427  19960123  Canon Kabushiki Kaisha  Inspection method and apparatus for inspecting a particle, if any, on a substrate having a pattern 
US5432607A (en) *  19930222  19950711  International Business Machines Corporation  Method and apparatus for inspecting patterned thin films using diffracted beam ellipsometry 
NonPatent Citations (18)
Title 

Bach, et al., "Modern Topics in Electromagnetics and Antennas," PPL Conference Publication 13, Based on Lectures delivered at the 1976 Summer Inst. at the Technical University, Eindhoven, pp. 6.16.59, 1977. 
Bach, et al., Modern Topics in Electromagnetics and Antennas, PPL Conference Publication 13, Based on Lectures delivered at the 1976 Summer Inst. at the Technical University, Eindhoven, pp. 6.1 6.59, 1977. * 
Bouwhuis, et al., "Principles of Optical Disc Systems," Philips Research Laboratories, Eindhoven, Adam Hilger Ltd, Bristol and Boston, pp. 88124, 1985. 
Bouwhuis, et al., Principles of Optical Disc Systems, Philips Research Laboratories, Eindhoven, Adam Hilger Ltd, Bristol and Boston, pp. 88 124, 1985. * 
Dil and Jacobs, "Apparent size of reflecting polygonal obstacles of the orde of one wavelength," J. Opt. Soc. Am., 69(7):950, Jul. 1979. 
Dil and Jacobs, Apparent size of reflecting polygonal obstacles of the orde of one wavelength, J. Opt. Soc. Am., 69(7):950, Jul. 1979. * 
Giovannini, et al., "Angleresolved polarimetric phase measurement for the characterization of grating", Optic Letters, 21(20):1619, Oct. 15, 1996. 
Giovannini, et al., Angle resolved polarimetric phase measurement for the characterization of grating , Optic Letters, 21(20):1619, Oct. 15, 1996. * 
Prins, et al., "Scatterometric sensor for PEB process control," 1996 SPIE Microlithography Conference in Santa Clara, CA, Proceedings SPIE 2725, pp. 110. 
Prins, et al., Scatterometric sensor for PEB process control, 1996 SPIE Microlithography Conference in Santa Clara, CA, Proceedings SPIE 2725, pp. 1 10. * 
Raymond, et al., "Multiparameter CD measurements using scatterometry," 1996 SPIE Microlithography Conference in Santa Clara, CA, Proceedings SPIE 2725, pp. 112. 
Raymond, et al., Multi parameter CD measurements using scatterometry, 1996 SPIE Microlithography Conference in Santa Clara, CA, Proceedings SPIE 2725, pp. 1 12. * 
Roger F. Harrington, "TimeHarmonic Electromagnetic Fields", Introduction to Waves, McGrawHill Book Co. pp. 6691, 1961. 
Roger F. Harrington, Time Harmonic Electromagnetic Fields , Introduction to Waves, McGraw Hill Book Co. pp. 66 91, 1961. * 
Sirat and Psaltis, "Conoscopic holography", Optics Letters 10(1):4, Jan. 1985. 
Sirat and Psaltis, Conoscopic holography , Optics Letters 10(1):4, Jan. 1985. * 
Wirgin, A., "A New Theoretical Approach To Scattering From a Periodic Interface," Optics Communications, 27(2):189, Nov. 1978. 
Wirgin, A., A New Theoretical Approach To Scattering From a Periodic Interface, Optics Communications, 27(2):189, Nov. 1978. * 
Cited By (415)
Publication number  Priority date  Publication date  Assignee  Title 

US20060250612A1 (en) *  19970922  20061109  Meeks Steven W  Detecting and classifying surface features or defects by controlling the angle of the illumination plane of incidence with respect to the feature or defect 
US6665078B1 (en)  19970922  20031216  Candela Instruments  System and method for simultaneously measuring thin film layer thickness, reflectivity, roughness, surface profile and magnetic pattern in thin film magnetic disks and silicon wafers 
US20040017561A1 (en) *  19970922  20040129  Meeks Steven W.  Method of detecting and classifying scratches and particles on thin film disks or wafers 
US20020015146A1 (en) *  19970922  20020207  Meeks Steven W.  Combined high speed optical profilometer and ellipsometer 
US6956660B2 (en)  19970922  20051018  KlaTencor Technologies  System and method for measuring properties of an object using a phase difference between two reflected light signals 
US6956658B2 (en)  19970922  20051018  KlaTencor Technologies Corporation  System and method for measuring object characteristics using phase differences in polarized light reflections 
US7688435B2 (en)  19970922  20100330  KlaTencor Corporation  Detecting and classifying surface features or defects by controlling the angle of the illumination plane of incidence with respect to the feature or defect 
US20040046959A1 (en) *  19970922  20040311  Meeks Steven W.  System and method for simultaneously measuring thin film layer thickness, reflectivity, roughness, surface profile and magnetic pattern on thin film magnetic disks and silicon wafers 
US6717671B1 (en)  19970922  20040406  Candela Instruments  System for simultaneously measuring thin film layer thickness, reflectivity, roughness, surface profile and magnetic pattern 
US7630086B2 (en)  19970922  20091208  KlaTencor Corporation  Surface finish roughness measurement 
US20070153273A1 (en) *  19970922  20070705  Meeks Steven W  Material independent profiler 
US7714995B2 (en)  19970922  20100511  KlaTencor Corporation  Material independent profiler 
US7123357B2 (en)  19970922  20061017  Candela Instruments  Method of detecting and classifying scratches and particles on thin film disks or wafers 
US20070115483A1 (en) *  19970922  20070524  Oak Dave S  Surface finish roughness measurement 
US6185323B1 (en) *  19971023  20010206  International Business Machines Corporation  Method characterizing a feature using measurement imaging tool 
US20070091327A1 (en) *  19980306  20070426  Yiping Xu  Spectroscopic Scatterometer System 
US7859659B2 (en)  19980306  20101228  KlaTencor Corporation  Spectroscopic scatterometer system 
US7173699B2 (en) *  19980306  20070206  KlaTencor Technologies Corporation  Spectroscopic scatterometer system 
US7898661B2 (en)  19980306  20110301  KlaTencor Corporation  Spectroscopic scatterometer system 
US20030058443A1 (en) *  19980306  20030327  Yiping Xu  Spectroscopic scatterometer system 
US20110125458A1 (en) *  19980306  20110526  KlaTencor Corporation  Spectroscopic Scatterometer System 
US6160251A (en) *  19980605  20001212  Sandia Corporation  Method and apparatus for distinguishing actual sparse events from sparse event false alarms 
US6735333B1 (en) *  19980730  20040511  Nidek Co., Ltd.  Pattern inspection apparatus 
US20040169850A1 (en) *  19990702  20040902  Meeks Steven W.  System and method for double sided optical inspection of thin film disks or wafers 
US7061601B2 (en)  19990702  20060613  KlaTencor Technologies Corporation  System and method for double sided optical inspection of thin film disks or wafers 
US6992773B1 (en)  19990830  20060131  Advanced Micro Devices, Inc.  Dualdifferential interferometry for silicon device damage detection 
US6552808B2 (en)  19991006  20030422  Metrolaser, Inc.  Methods and apparatus for splitting, imaging, and measuring wavefronts in interferometry 
US20060132795A1 (en) *  19991006  20060622  Millerd James E  Methods and apparatus for splitting, imaging, and measuring wavefronts in interferometry 
US7298497B2 (en)  19991006  20071120  Metrolaser, Inc.  Methods and apparatus for splitting, imaging, and measuring wavefronts in interferometry 
US7170611B2 (en)  19991006  20070130  Metrolaser, Inc.  Methods and apparatus for splitting, imaging, and measuring wavefronts in interferometry 
US6304330B1 (en)  19991006  20011016  Metrolaser, Inc.  Methods and apparatus for splitting, imaging, and measuring wavefronts in interferometry 
US6529278B2 (en) *  19991227  20030304  Canon Kabushiki Kaisha  Optical interference apparatus and position detection apparatus 
US6972852B2 (en)  20000329  20051206  ThermaWave, Inc.  Critical dimension analysis with simultaneous multiple angle of incidence measurements 
US7248375B2 (en)  20000329  20070724  ThermaWave, Inc.  Critical dimension analysis with simultaneous multiple angle of incidence measurements 
US20060012803A1 (en) *  20000329  20060119  Jon Opsal  Critical dimension analysis with simultaneous multiple angle of incidence measurements 
US20040046968A1 (en) *  20000329  20040311  Jon Opsal  Critical dimension analysis with simultaneous multiple angle of incidence measurements 
US6829057B2 (en)  20000329  20041207  ThermaWave, Inc.  Critical dimension analysis with simultaneous multiple angle of incidence measurements 
US20050057760A1 (en) *  20000329  20050317  Jon Opsal  Critical dimension analysis with simultaneous multiple angle of incidence measurements 
US6556291B2 (en)  20000531  20030429  Nidek Co., Ltd.  Defect inspection method and defect inspection apparatus 
US7372565B1 (en)  20000925  20080513  Nanometrics Incorporated  Spectrometer measurement of diffracting structures 
US7115858B1 (en)  20000925  20061003  Nanometrics Incorporated  Apparatus and method for the measurement of diffracting structures 
US20040246481A1 (en) *  20001003  20041209  Accent Optical Technologies, Inc.  Differential numerical aperture methods 
US6750968B2 (en) *  20001003  20040615  Accent Optical Technologies, Inc.  Differential numerical aperture methods and device 
US20020039184A1 (en) *  20001003  20020404  Sandusky John V.  Differential numerical aperture methods and device 
US7053991B2 (en) *  20001003  20060530  Accent Optical Technologies, Inc.  Differential numerical aperture methods 
US6504612B2 (en)  20001114  20030107  Georgia Tech Research Corporation  Electromagnetic wave analyzer 
US6597712B2 (en) *  20010226  20030722  Hitachi, Ltd.  Laser diode module 
US6708127B2 (en) *  20010315  20040316  Lucent Technologies Inc.  Beam propagation method for stepindex waveguides 
US20040233419A1 (en) *  20010326  20041125  Meeks Steven W.  Combined high speed optical profilometer and ellipsometer 
US20020163634A1 (en) *  20010326  20021107  Meeks Steven W.  Multiple spot size optical profilometer, ellipsometer, reflectometer and scatterometer 
US20020145740A1 (en) *  20010326  20021010  Meeks Steven W.  Material independent optical profilometer 
US20050057747A1 (en) *  20010326  20050317  Meeks Steven W.  Method of detecting and classifying scratches, particles and pits on thin film disks or wafers 
US7075630B2 (en)  20010326  20060711  KlaTencor Technologies Corporation  Combined high speed optical profilometer and ellipsometer 
US7113284B1 (en)  20010326  20060926  KlaTencor Technologies Corporation  Material independent optical profilometer 
US6930765B2 (en)  20010326  20050816  KlaTencor Technologies  Multiple spot size optical profilometer, ellipsometer, reflectometer and scatterometer 
US6897957B2 (en)  20010326  20050524  Candela Instruments  Material independent optical profilometer 
US6909500B2 (en)  20010326  20050621  Candela Instruments  Method of detecting and classifying scratches, particles and pits on thin film disks or wafers 
US20030025905A1 (en) *  20010326  20030206  Meeks Steven W.  Method of detecting and classifying scratches, particles and pits on thin film disks or wafers 
US6757056B1 (en) *  20010326  20040629  Candela Instruments  Combined high speed optical profilometer and ellipsometer 
US7907765B2 (en)  20010328  20110315  University Of Washington  Focal plane tracking for optical microtomography 
US6944322B2 (en)  20010328  20050913  Visiongate, Inc.  Optical tomography of small objects using parallel ray illumination and postspecimen optical magnification 
US20060023219A1 (en) *  20010328  20060202  Meyer Michael G  Optical tomography of small objects using parallel ray illumination and postspecimen optical magnification 
US6602723B1 (en)  20010402  20030805  Advanced Micro Devices, Inc.  Method of integrating scatterometry metrology structures directly into die design 
US6989900B1 (en) *  20010402  20060124  Advanced Micro Devices, Inc.  Method of measuring implant profiles using scatterometric techniques 
US6660542B1 (en)  20010406  20031209  Advanced Micro Devices, Inc.  Method of controlling stepper process parameters based upon optical properties of incoming process layers, and system for accomplishing same 
US6618149B1 (en)  20010406  20030909  Advanced Micro Devices, Inc.  Method of identifying film stacks based upon optical properties 
US6980300B1 (en)  20010411  20051227  Advanced Micro Devices, Inc.  Method and apparatus for generating a polishing process endpoint signal using scatterometry 
US6479200B1 (en)  20010419  20021112  Advanced Micro Devices, Inc.  Method of controlling stepper process parameters based upon scatterometric measurements of DICD features 
US6383824B1 (en)  20010425  20020507  Advanced Micro Devices, Inc.  Method of using scatterometry measurements to control deposition processes 
US6602727B1 (en)  20010427  20030805  Advanced Micro Devices, Inc.  Scatterometry based active control of exposure conditions 
US6898537B1 (en)  20010427  20050524  Nanometrics Incorporated  Measurement of diffracting structures using onehalf of the nonzero diffracted orders 
US6639663B1 (en)  20010523  20031028  Advanced Micro Devices, Inc.  Method and apparatus for detecting processing faults using scatterometry measurements 
US6479309B1 (en)  20010525  20021112  Advanced Micro Devices, Inc.  Method and apparatus for determining process layer conformality 
US6433871B1 (en)  20010525  20020813  Advanced Micron Devices, Inc.  Method of using scatterometry measurements to determine and control gate electrode profiles 
US6582863B1 (en)  20010611  20030624  Advanced Micro Devices, Inc.  Method of controlling photolithography processes based upon scatterometric measurements of subnominal grating structures 
US6529282B1 (en)  20010611  20030304  Advanced Micro Devices, Inc.  Method of controlling photolithography processes based upon scatterometric measurements of photoresist thickness, and system for accomplishing same 
US6630362B1 (en) *  20010613  20031007  Advanced Micro Devices, Inc.  Method and apparatus for performing trench depth analysis 
US6537833B1 (en)  20010619  20030325  Advanced Micro Devices, Inc.  Method and apparatus for characterizing an interconnect structure profile using scatterometry measurements 
US6662125B2 (en) *  20010626  20031209  Fujitsu Limited  Electromagnetic wave analyzer and program for same 
US6617087B1 (en)  20010627  20030909  Advanced Micro Devices, Inc.  Use of scatterometry to measure pattern accuracy 
US6458605B1 (en)  20010628  20021001  Advanced Micro Devices, Inc.  Method and apparatus for controlling photolithography overlay registration 
US6707562B1 (en)  20010702  20040316  Advanced Micro Devices, Inc.  Method of using scatterometry measurements to control photoresist etch process 
US6773939B1 (en)  20010702  20040810  Advanced Micro Devices, Inc.  Method and apparatus for determining critical dimension variation in a line structure 
US6804014B1 (en)  20010702  20041012  Advanced Micro Devices, Inc.  Method and apparatus for determining contact opening dimensions using scatterometry 
US6650423B1 (en)  20010702  20031118  Advanced Micro Devices Inc.  Method and apparatus for determining column dimensions using scatterometry 
US7262864B1 (en) *  20010702  20070828  Advanced Micro Devices, Inc.  Method and apparatus for determining grid dimensions using scatterometry 
US6713753B1 (en)  20010703  20040330  Nanometrics Incorporated  Combination of normal and oblique incidence polarimetry for the characterization of gratings 
US6556303B1 (en)  20010710  20030429  Advanced Micro Devices, Inc.  Scattered signal collection using strobed technique 
US6535288B1 (en)  20010710  20030318  Advanced Micro Devices, Inc.  Machine readable code to trigger data collection 
US6867866B1 (en) *  20010810  20050315  ThermaWave, Inc.  CD metrology analysis using green's function 
US6741730B2 (en) *  20010810  20040525  Visiongate, Inc.  Method and apparatus for threedimensional imaging in the fourier domain 
US20050137809A1 (en) *  20010810  20050623  Chang Yia C.  CD metrology analysis using green's function 
US7038850B2 (en)  20010810  20060502  ThermWave, Inc.  CD metrology analysis using green's function 
US6704478B2 (en) *  20010823  20040309  Hitachi Cable, Ltd.  Wavelength separation optical device and multiple wavelength light transmission module 
US6711326B2 (en) *  20010823  20040323  Hitachi Cable, Ltd.  Wavelength separation optical device and multiple wavelength light transmission module 
US20060066868A1 (en) *  20010828  20060330  Jon Opsal  Detector configurations for optical metrology 
US7456964B2 (en)  20010828  20081125  KlaTencor Corporation  Detector configurations for optical metrology 
US7206071B2 (en)  20010828  20070417  ThermaWave, Inc.  Detector configurations for optical metrology 
US20080018895A1 (en) *  20010828  20080124  Jon Opsal  Detector configurations for optical metrology 
US20090066954A1 (en) *  20010828  20090312  Jon Opsal  Detector configurations for optical metrology 
US7667841B2 (en)  20010828  20100223  KlaTencor Corporation  Detector configurations for optical metrology 
US20050105089A1 (en) *  20010828  20050519  Jon Opsal  Detector configurations for optical metrology 
US6995842B2 (en)  20010828  20060207  ThermaWave, Inc.  Detector configurations for optical metrology 
US20030043372A1 (en) *  20010830  20030306  Bernd Schulz  Apparatus for determining an overlay error and critical dimensions in a semiconductor structure by means of scatterometry 
US6765282B2 (en)  20010830  20040720  Advanced Micro Devices, Inc.  Semiconductor structure and method for determining critical dimensions and overlay error 
US20030042579A1 (en) *  20010830  20030306  Bernd Schulz  Semiconductor structure and method for determining critical dimensions and overlay error 
US20030044702A1 (en) *  20010830  20030306  Bernd Schulz  Semiconductor structure and method for determining critical dimensions and overlay error 
US6767680B2 (en)  20010830  20040727  Advanced Micro Devices, Inc.  Semiconductor structure and method for determining critical dimensions and overlay error 
US6816252B2 (en)  20010830  20041109  Advanced Micro Devices, Inc.  Apparatus for determining an overlay error and critical dimensions in a semiconductor structure by means of scatterometry 
US7061615B1 (en)  20010920  20060613  Nanometrics Incorporated  Spectroscopically measured overlay target 
US20040223137A1 (en) *  20011010  20041111  Littau Michael E.  Determination of center of focus by crosssection analysis 
US7110099B2 (en)  20011010  20060919  Accent Optical Technologies, Inc.  Determination of center of focus by crosssection analysis 
US6790570B1 (en)  20011108  20040914  Advanced Micro Devices, Inc.  Method of using scatterometry measurements to control stepper process parameters 
US7826071B2 (en)  20011219  20101102  KlaTencor Corporation  Parametric profiling using optical spectroscopic systems 
US20040070772A1 (en) *  20011219  20040415  Shchegrov Andrei V.  Parametric profiling using optical spectroscopic systems 
US20090135416A1 (en) *  20011219  20090528  KlaTencor Technologies Corporation  Parametric Profiling Using Optical Spectroscopic Systems 
US7280230B2 (en)  20011219  20071009  KlaTencor Technologies Corporation  Parametric profiling using optical spectroscopic systems 
US6579733B1 (en)  20020116  20030617  Advanced Micro Devices, Inc.  Using scatterometry to measure resist thickness and control implant 
WO2003060422A2 (en) *  20020121  20030724  Scinex Ag  Interferometric optical system 
WO2003060422A3 (en) *  20020121  20040205  Scinex Ag  Interferometric optical system 
US7161683B2 (en) *  20020121  20070109  Campus Technologies Ag  Polarizationdependent grating interferometer for measuring optical profile depth and spectral properties of a sample 
US20050117169A1 (en) *  20020121  20050602  Thilo Weitzel  Interferometric optical system 
US20040212810A1 (en) *  20020205  20041028  Jon Opsal's  Analysis of isolated and aperiodic structures with simultaneous multiple angle of incidence measurements 
US20030147086A1 (en) *  20020205  20030807  Allan Rosencwaig  Analysis of isolated and aperiodic structures with simultaneous multiple angle of incidence measurements 
US6842259B2 (en)  20020205  20050111  ThermaWave, Inc.  Analysis of isolated and aperiodic structures with simultaneous multiple angle of incidence measurements 
US6813034B2 (en) *  20020205  20041102  ThermaWave, Inc.  Analysis of isolated and aperiodic structures with simultaneous multiple angle of incidence measurements 
US6643008B1 (en)  20020226  20031104  Advanced Micro Devices, Inc.  Method of detecting degradation in photolithography processes based upon scatterometric measurements of grating structures, and a device comprising such structures 
US6785009B1 (en)  20020228  20040831  Advanced Micro Devices, Inc.  Method of using high yielding spectra scatterometry measurements to control semiconductor manufacturing processes, and systems for accomplishing same 
US20050041258A1 (en) *  20020313  20050224  Jon Opsal  Optical scatterometry of asymmetric lines and structures 
US7061627B2 (en)  20020313  20060613  ThermaWave, Inc.  Optical scatterometry of asymmetric lines and structures 
US6742168B1 (en)  20020319  20040525  Advanced Micro Devices, Inc.  Method and structure for calibrating scatterometrybased metrology tool used to measure dimensions of features on a semiconductor device 
US6791697B1 (en)  20020321  20040914  Advanced Micro Devices, Inc.  Scatterometry structure with embedded ring oscillator, and methods of using same 
US7236244B1 (en)  20020404  20070626  Nanometrics Incorporated  Alignment target to be measured with multiple polarization states 
US6949462B1 (en)  20020404  20050927  Nanometrics Incorporated  Measuring an alignment target with multiple polarization states 
US7230705B1 (en)  20020404  20070612  Nanometrics Incorporated  Alignment target with designed in offset 
US6982793B1 (en)  20020404  20060103  Nanometrics Incorporated  Method and apparatus for using an alignment target with designed in offset 
US20050010108A1 (en) *  20020419  20050113  Rahn John Richard  Method for correction of relative objectdetector motion between successive views 
US7260253B2 (en)  20020419  20070821  Visiongate, Inc.  Method for correction of relative objectdetector motion between successive views 
US7542597B2 (en)  20020419  20090602  Visiongate, Inc.  Method for correction of relative objectdetector motion between successive views 
US6882437B2 (en)  20020419  20050419  KlaTencor Technologies  Method of detecting the thickness of thin film disks or wafers 
US20050085708A1 (en) *  20020419  20050421  University Of Washington  System and method for preparation of cells for 3D image acquisition 
US7738945B2 (en)  20020419  20100615  University Of Washington  Method and apparatus for pseudoprojection formation for optical tomography 
US20050085721A1 (en) *  20020419  20050421  University Of Washington  System and method for processing specimens and images for optical tomography 
US20040076319A1 (en) *  20020419  20040422  Fauver Mark E.  Method and apparatus of shadowgram formation for optical tomography 
US7197355B2 (en)  20020419  20070327  Visiongate, Inc.  Variablemotion optical tomography of small objects 
US7811825B2 (en)  20020419  20101012  University Of Washington  System and method for processing specimens and images for optical tomography 
US6770893B2 (en)  20020513  20040803  Visiongate, Inc.  Method and apparatus for emission computed tomography using temporal signatures 
US20030210760A1 (en) *  20020513  20031113  Nelson Alan C.  Method and apparatus for emission computed tomography using temporal signatures 
US6795193B2 (en)  20020531  20040921  Advanced Micro Devices, Inc.  Scatterometer including an internal calibration system 
US20030224261A1 (en) *  20020531  20031204  Bernd Schulz  Twodimensional structure for determining an overlay accuracy by means of scatterometry 
US20030223072A1 (en) *  20020531  20031204  Bernd Schulz  Scatterometer including an internal calibration system 
US7099010B2 (en)  20020531  20060829  Advanced Micro Devices, Inc.  Twodimensional structure for determining an overlay accuracy by means of scatterometry 
US7586623B2 (en)  20020618  20090908  Tokyo Electron Limited  Optical metrology of single features 
US7379192B2 (en)  20020618  20080527  Timbre Technologies, Inc.  Optical metrology of single features 
US20060187468A1 (en) *  20020618  20060824  Timbre Technologies, Inc.  Optical metrology of single features 
US20040212812A1 (en) *  20020618  20041028  Timbre Technologies, Inc.  Optical metrology of single features 
US20080259357A1 (en) *  20020618  20081023  Timbre Technologies, Inc.  Optical metrology of single features 
US7030999B2 (en) *  20020618  20060418  Timbre Technologies, Inc.  Optical metrology of single features 
US20080049233A1 (en) *  20020909  20080228  Zygo Corporation  MultipleAngle MultipleWavelength Interferometer Using HighNA Imaging and Spectral Analysis 
US7869057B2 (en)  20020909  20110111  Zygo Corporation  Multipleangle multiplewavelength interferometer using highNA imaging and spectral analysis 
US7812963B2 (en)  20020909  20101012  Zygo Corporation  Interferometry method for ellipsometry, reflectometry, and scatterometry measurements, including characterization of thin film structures 
US20090015844A1 (en) *  20020909  20090115  De Groot Peter J  Interferometry Method for Ellipsometry, Reflectometry, and Scatterometry Measurements, Including Characterization of Thin Film Structures 
US6623994B1 (en) *  20020926  20030923  Advanced Micro Devices, Inc.  Method for calibrating opticalbased metrology tools 
US6972853B1 (en)  20020927  20051206  Advanced Micro Devices, Inc.  Methods of calibrating and controlling stepper exposure processes and tools, and system for accomplishing same 
US6992764B1 (en)  20020930  20060131  Nanometrics Incorporated  Measuring an alignment target with a single polarization state 
US20040133362A1 (en) *  20021021  20040708  Emit Systems Co.  Method for significantly increasing computational speed, accuracy and versatility of electromagnetic reconstruction of shapes and composition of complex targets containing lossy materials 
US6881594B1 (en)  20021028  20050419  Advanced Micro Devices, Inc.  Method of using scatterometry for analysis of electromigration, and structures for performing same 
US6933158B1 (en)  20021031  20050823  Advanced Micro Devices, Inc.  Method of monitoring anneal processes using scatterometry, and system for performing same 
US20080266574A1 (en) *  20030306  20081030  Groot Peter De  Interferometer and method for measuring characteristics of optically unresolved surface features 
US20080065350A1 (en) *  20030306  20080313  De Groot Peter J  Profiling complex surface structures using scanning interferometry 
US7684049B2 (en)  20030306  20100323  Zygo Corporation  Interferometer and method for measuring characteristics of optically unresolved surface features 
US20100265516A1 (en) *  20030306  20101021  De Groot Peter  Interferometer and method for measuring characteristics of optically unresolved surface features 
US7466429B2 (en)  20030306  20081216  Zygo Corporation  Profiling complex surface structures using scanning interferometry 
US7948636B2 (en)  20030306  20110524  Zygo Corporation  Interferometer and method for measuring characteristics of optically unresolved surface features 
US7916308B2 (en) *  20030401  20110329  Seagate Technology Llc  Method and optical profiler 
US20040196469A1 (en) *  20030401  20041007  Seagate Technology Llc  Optical profiler for ultrasmooth surface with normal incident beam deflection method 
US6836691B1 (en)  20030501  20041228  Advanced Micro Devices, Inc.  Method and apparatus for filtering metrology data based on collection purpose 
US6859746B1 (en)  20030501  20050222  Advanced Micro Devices, Inc.  Methods of using adaptive sampling techniques based upon categorization of process variations, and system for performing same 
US20050014043A1 (en) *  20030718  20050120  Matsushita Electric Industrial Co., Ltd.  Power supply unit 
US7586620B2 (en)  20030915  20090908  Zygo Corporation  Methods and systems for interferometric analysis of surfaces and related applications 
US7456975B2 (en)  20030915  20081125  Zygo Corporation  Methods and systems for interferometric analysis of surfaces and related applications 
US20080068614A1 (en) *  20030915  20080320  De Groot Peter J  Methods and systems for interferometric analysis of surfaces and related applications 
US8107085B2 (en)  20030915  20120131  Zygo Corporation  Methods and systems for interferometric analysis of surfaces and related applications 
US7327457B2 (en)  20031219  20080205  N&K Technology, Inc.  Apparatus and method for optical characterization of a sample over a broadband of wavelengths while minimizing polarization changes 
US20050134847A1 (en) *  20031219  20050623  Ray Hebert  Apparatus and method for optical characterization of a sample over a broadband of wavelengths with a small spot size 
US20050134848A1 (en) *  20031219  20050623  Ray Hebert  Apparatus and method for optical characterization of a sample over a broadband of wavelengths while minimizing polarization changes 
US8035824B2 (en)  20031219  20111011  International Business Machines Corporation  Differential critical dimension and overlay metrology apparatus and measurement method 
US7248364B2 (en)  20031219  20070724  N&K Technology, Inc.  Apparatus and method for optical characterization of a sample over a broadband of wavelengths with a small spot size 
US7700247B2 (en)  20031219  20100420  International Business Machines Corporation  Differential critical dimension and overlay metrology apparatus and measurement method 
US20060066878A1 (en) *  20040614  20060330  Samsung Electronics Co., Ltd.  System and method for measuring dimension of patterns formed on photomask 
US7369254B2 (en) *  20040614  20080506  Samsung Electronics Co., Ltd.  System and method for measuring dimension of patterns formed on photomask 
US20060215289A1 (en) *  20040614  20060928  KlaTencor Technologies Corp.  Servo pattern characterization on magnetic disks 
US7397621B2 (en)  20040614  20080708  KlaTencor Technologies Corporation  Servo pattern characterization on magnetic disks 
US7075741B1 (en)  20040614  20060711  Kla Tencor Technologues Corporation  System and method for automatically determining magnetic eccentricity of a disk 
EP1610089A3 (en) *  20040621  20090218  Trioptics GmbH  Method and device for measuring the angle of optical surfaces 
US20060012582A1 (en) *  20040715  20060119  De Lega Xavier C  Transparent film measurements 
US20060028962A1 (en) *  20040805  20060209  Oron Zachar  Method and system for encoding and detecting optical information 
US20110007314A1 (en) *  20040816  20110113  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US8111398B2 (en)  20040816  20120207  Asml Netherlands B.V.  Method of measurement, an inspection apparatus and a lithographic apparatus 
US7791732B2 (en)  20040816  20100907  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US8760662B2 (en)  20040816  20140624  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US8054467B2 (en)  20040816  20111108  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20060033921A1 (en) *  20040816  20060216  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US7791727B2 (en) *  20040816  20100907  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US8553230B2 (en)  20040816  20131008  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20100277706A1 (en) *  20040816  20101104  Asml Netherlands B.V.  Method of Measurement, an Inspection Apparatus and a Lithographic Apparatus 
US7396022B1 (en)  20040928  20080708  KlaTencor Technologies Corp.  System and method for optimizing wafer flatness at high rotational speeds 
US6991738B1 (en)  20041013  20060131  University Of Washington  Flowthrough drum centrifuge 
US20060096358A1 (en) *  20041028  20060511  University Of Washington  Optical projection tomography microscope 
US20060183220A1 (en) *  20041109  20060817  Nelson Alan C  Automated cell sample enrichment preparation method 
US7494809B2 (en)  20041109  20090224  Visiongate, Inc.  Automated cell sample enrichment preparation method 
US20060139592A1 (en) *  20041122  20060629  Asml Netherlands B.V.  Latent overlay metrology 
US7201799B1 (en)  20041124  20070410  KlaTencor Technologies Corporation  System and method for classifying, detecting, and counting micropipes 
US20060126074A1 (en) *  20041214  20060615  Asml Netherlands B.V.  Inspection apparatus, sample and inspection method 
US7453577B2 (en)  20041214  20081118  Asml Netherlands B.V.  Apparatus and method for inspecting a patterned part of a sample 
US7684032B1 (en)  20050106  20100323  KlaTencor Corporation  Multiwavelength system and method for detecting epitaxial layer defects 
US7446882B2 (en) *  20050120  20081104  Zygo Corporation  Interferometer for determining characteristics of an object surface 
US7428057B2 (en)  20050120  20080923  Zygo Corporation  Interferometer for determining characteristics of an object surface, including processing and calibration 
US7952724B2 (en)  20050120  20110531  Zygo Corporation  Interferometer with multiple modes of operation for determining characteristics of an object surface 
US7616323B2 (en)  20050120  20091110  Zygo Corporation  Interferometer with multiple modes of operation for determining characteristics of an object surface 
US20060158658A1 (en) *  20050120  20060720  Xavier Colonna De Lega  Interferometer with multiple modes of operation for determining characteristics of an object surface 
US20080088849A1 (en) *  20050120  20080417  Zygo Corporation  Interferometry for determining characteristics of an object surface, with spatially coherent illumination 
US20060158659A1 (en) *  20050120  20060720  Xavier Colonna De Lega  Interferometer for determining characteristics of an object surface 
US20060158657A1 (en) *  20050120  20060720  De Lega Xavier C  Interferometer for determining characteristics of an object surface, including processing and calibration 
US7884947B2 (en)  20050120  20110208  Zygo Corporation  Interferometry for determining characteristics of an object surface, with spatially coherent illumination 
US7443486B2 (en)  20050225  20081028  Asml Netherlands B.V.  Method for predicting a critical dimension of a feature imaged by a lithographic apparatus 
US20060192936A1 (en) *  20050225  20060831  Asml Netherlands B.V.  Lithography measurements using scatterometry 
US7161669B2 (en)  20050506  20070109  Kla Tencor Technologies Corporation  Wafer edge inspection 
US7532318B2 (en)  20050506  20090512  KlaTencor Corporation  Wafer edge inspection 
US20060250611A1 (en) *  20050506  20061109  KlaTencor Technologies Corp.  Wafer edge inspection 
US20060250610A1 (en) *  20050506  20061109  KlaTencor Technologies Corp.  Wafer edge inspection 
US20070127016A1 (en) *  20050506  20070607  Meeks Steven W  Wafer edge inspection 
US20060250609A1 (en) *  20050506  20061109  KlaTencor Technologies Corp.  Wafer edge inspection 
US7161667B2 (en)  20050506  20070109  KlaTencor Technologies Corporation  Wafer edge inspection 
US7161668B2 (en)  20050506  20070109  KlaTencor Technologies Corporation  Wafer edge inspection 
US20060256324A1 (en) *  20050510  20061116  Asml Netherlands B.V.  Inspection apparatus and method of inspection 
US7315384B2 (en)  20050510  20080101  Asml Netherlands B.V.  Inspection apparatus and method of inspection 
US7548308B2 (en)  20050511  20090616  KlaTencor Corporation  Illumination energy management in surface inspection 
US20060256325A1 (en) *  20050511  20061116  KlaTencor Technologies Corp.  Illumination energy management in surface inspection 
US20060256345A1 (en) *  20050512  20061116  KlaTencor Technologies Corp.  Interferometry measurement in disturbed environments 
US7505143B2 (en)  20050517  20090317  KlaTencor Corporation  Dynamic reference plane compensation 
US20060262291A1 (en) *  20050517  20061123  Hess Harald F  Dynamic reference plane compensation 
US7403293B2 (en)  20050630  20080722  Asml Netherlands  Metrology apparatus, lithographic apparatus, process apparatus metrology method and device manufacturing method 
US20070003840A1 (en) *  20050630  20070104  Asml Netherlands B.V.  Focus determination method, device manufacturing method, and mask 
US20070002336A1 (en) *  20050630  20070104  Asml Netherlands B.V.  Metrology apparatus, lithographic apparatus, process apparatus, metrology method and device manufacturing method 
US20070013921A1 (en) *  20050630  20070118  Asml Netherlands B.V.  Metrology apparatus, lithographic apparatus, process apparatus metrology method and device manufacturing method 
US7532307B2 (en)  20050630  20090512  Asml Netherlands B.V.  Focus determination method, device manufacturing method, and mask 
US20070064572A1 (en) *  20050727  20070322  Polarizonics Corporation  Optical data encoding and detection 
US7286229B1 (en)  20050906  20071023  KlaTencor Technologies Corporation  Detecting multidomain states in perpendicular magnetic media 
US7295300B1 (en)  20050928  20071113  KlaTencor Technologies Corporation  Detecting surface pits 
US20070086013A1 (en) *  20051011  20070419  Zygo Corporation  Interferometry method and system including spectral decomposition 
US7636168B2 (en)  20051011  20091222  Zygo Corporation  Interferometry method and system including spectral decomposition 
US7397553B1 (en)  20051024  20080708  KlaTencor Technologies Corporation  Surface scanning 
US20070093044A1 (en) *  20051025  20070426  Asml Netherlands B.V.  Method of depositing a metal layer onto a substrate and a method for measuring in three dimensions the topographical features of a substrate 
US20070153274A1 (en) *  20051230  20070705  Asml Netherlands B.V.  Optical metrology system and metrology mark characterization device 
US7649636B2 (en)  20051230  20100119  Asml Netherlands B.V.  Optical metrology system and metrology mark characterization device 
US7480050B2 (en)  20060209  20090120  Asml Netherlands B.V.  Lithographic system, sensor, and method of measuring properties of a substrate 
US20070182964A1 (en) *  20060209  20070809  Asml Netherlands B.V.  Lithographic system, sensor, and method of measuring properties of a substrate 
US7889616B2 (en) *  20060303  20110215  Seagate Technology Llc  Polarization state detector systems and methods for readout of multilevel oriented nanostructure optical media 
US20070206480A1 (en) *  20060303  20070906  Seagate Technology Llc  Polarization state detector systems and methods for readout of multilevel oriented nanostructure optical media 
US7486408B2 (en)  20060321  20090203  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method with reduced scribe lane usage for substrate measurement 
US7821650B2 (en)  20060321  20101026  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method with reduced scribe lane usage for substrate measurement 
US20090027691A1 (en) *  20060321  20090129  Asml Netherlands B.V.  Lithographic Apparatus and Device Manufacturing Method with Reduced Scribe Lane Usage for Substrate Measurement 
US20070222960A1 (en) *  20060321  20070927  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method with reduced scribe lane usage for substrate measurement 
US20090244538A1 (en) *  20060328  20091001  Asml Netherlands B.V.  Lithographic Apparatus and Device Manufacturing Method Using Overlay Measurement 
US8264686B2 (en)  20060328  20120911  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method using overlay measurement 
US7532305B2 (en)  20060328  20090512  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method using overlay measurement 
US20070229785A1 (en) *  20060328  20071004  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method using overlay measurement 
US7391513B2 (en)  20060329  20080624  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method using overlay measurement quality indication 
US20070229837A1 (en) *  20060329  20071004  Asml Netherlands B.V.  Lithographic apparatus and device manufacturing method using overlay measurement quality indication 
US7684145B2 (en)  20060511  20100323  Microsense, Llc  Method and system for perpendicular magnetic media metrology 
US20070262771A1 (en) *  20060511  20071115  KlaTencor Technologies Corporation  Method and System for Perpendicular Magnetic Media Metrology 
US20100002326A1 (en) *  20060511  20100107  Ade Technologies, Inc.  Method and system for perpendicular magnetic media metrology 
US20070279742A1 (en) *  20060531  20071206  Asml Netherlands B.V.  Metrology tool, system comprising a lithographic apparatus and a metrology tool, and a method for determining a parameter of a substrate 
US7961309B2 (en)  20060531  20110614  Asml Netherlands B.V.  Metrology tool, system comprising a lithographic apparatus and a metrology tool, and a method for determining a parameter of a substrate 
US20090296081A1 (en) *  20060531  20091203  Asml Netherlands B. V.  Metrology Tool, System Comprising a Lithographic Apparatus and a Metrology Tool, and a Method for Determining a Parameter of a Substrate 
US7586598B2 (en)  20060531  20090908  Asml Netherlands B.V.  Metrology tool, system comprising a lithographic apparatus and a metrology tool, and a method for determining a parameter of a substrate 
US20070279644A1 (en) *  20060531  20071206  Asml Netherlands B.V.  Metrology tool, system comprising a lithographic apparatus and a metrology tool, and a method for determining a parameter of a substrate 
US7502103B2 (en)  20060531  20090310  Asml Netherlands B.V.  Metrology tool, system comprising a lithographic apparatus and a metrology tool, and a method for determining a parameter of a substrate 
US20070296960A1 (en) *  20060613  20071227  Asml Netherlands B.V.  Method of characterising the transmission losses of an optical system 
US7791724B2 (en)  20060613  20100907  Asml Netherlands B.V.  Characterization of transmission losses in an optical system 
US20090262349A1 (en) *  20060616  20091022  Institut National Polytechnique De Grenoble  Electrooptic probe for vector measurement of an electromagnetic field 
US8264685B2 (en) *  20060616  20120911  Institut National Polytechnique De Grenoble  Electrooptic probe for vector measurement of an electromagnetic field 
US8064056B2 (en)  20060620  20111122  Asml Netherlands B.V.  Substrate used in a method and apparatus for angularresolved spectroscopic lithography characterization 
US20110122496A1 (en) *  20060620  20110526  Asml Netherlands B.V.  Substrate Used in a Method and Apparatus for AngularResolved Spectroscopic Lithography Characterization 
US20070291269A1 (en) *  20060620  20071220  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US7898662B2 (en)  20060620  20110301  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20070296973A1 (en) *  20060622  20071227  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US7692792B2 (en)  20060622  20100406  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US7999940B2 (en)  20060630  20110816  Asml Netherlands B.V.  Apparatus for angularresolved spectroscopic lithography characterization 
US7999943B2 (en)  20060630  20110816  Asml Netherlands B.V.  Device manufacturing method with angularresolved spectroscopic lithography characterization 
US20080002207A1 (en) *  20060630  20080103  Asml Netherlands B.V.  Apparatus for angularresolved spectroscopic lithography characterization and device manufacturing method 
US20080266560A1 (en) *  20060630  20081030  Asml Netherlands B.V.  Apparatus for angularresolved spectroscopic lithography characterization and device manufacturing method 
US7659988B2 (en)  20060630  20100209  Asml Netherlands B.V.  Apparatus for angularresolved spectroscopic lithography characterization and device manufacturing method 
US7916284B2 (en)  20060718  20110329  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7522288B2 (en)  20060721  20090421  Zygo Corporation  Compensation of systematic effects in low coherence interferometry 
US20080018901A1 (en) *  20060721  20080124  Zygo Corporation  Compensation of systematic effects in low coherence interferometry 
US8306306B2 (en)  20060802  20121106  Commissariat A L'energie Atomique  Device for characterizing unique objects 
FR2904690A1 (en) *  20060802  20080208  Commissariat Energie Atomique  A unique object characterization 
WO2008015230A1 (en) *  20060802  20080207  Commissariat A L'energie Atomique  Device for characterizing unique objects 
US20090316979A1 (en) *  20060802  20091224  Serge Gidon  Device for characterizing unique objects 
US7808613B2 (en)  20060803  20101005  Asml Netherlands B.V.  Individual wafer history storage for overlay corrections 
US20080030701A1 (en) *  20060803  20080207  Asml Netherlands B.V.  Individual wafer history storage for overlay corrections 
US20080069430A1 (en) *  20060808  20080320  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20080036984A1 (en) *  20060808  20080214  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US7643666B2 (en)  20060808  20100105  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20090294635A1 (en) *  20060810  20091203  Asml Netherlands B.V.  Inspection Method and Apparatus, Lithographic Apparatus, Lithographic Processing Cell and Device Manufacturing Method 
US20080037134A1 (en) *  20060810  20080214  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device method 
US8120001B2 (en)  20060810  20120221  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7589832B2 (en)  20060810  20090915  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device method 
US7564555B2 (en)  20060815  20090721  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20080043239A1 (en) *  20060815  20080221  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US7715019B2 (en)  20060914  20100511  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080068609A1 (en) *  20060914  20080320  Asml Netherlands B.V.  Inspection apparatus, an apparatus for projecting an image and a method of measuring a property of a substrate 
US7969577B2 (en)  20060914  20110628  Asml Netherlands B.V.  Inspection apparatus, an apparatus for projecting an image and a method of measuring a property of a substrate 
US7532331B2 (en)  20060914  20090512  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7573584B2 (en)  20060925  20090811  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US20080074666A1 (en) *  20060925  20080327  Asml Netherlands B.V.  Method and apparatus for angularresolved spectroscopic lithography characterization 
US8233155B2 (en)  20061013  20120731  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US8294907B2 (en)  20061013  20121023  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080088854A1 (en) *  20061013  20080417  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080117434A1 (en) *  20061122  20080522  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7630087B2 (en)  20061122  20091208  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080128644A1 (en) *  20061130  20080605  Asml Netherlands  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7710572B2 (en)  20061130  20100504  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7557934B2 (en)  20061207  20090707  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080151228A1 (en) *  20061208  20080626  Asml Netherlands B.V.  Scatterometer, a lithographic apparatus and a focus analysis method 
US7630070B2 (en)  20061208  20091208  Asml Netherlands B.V.  Scatterometer, a lithographic apparatus and a focus analysis method 
US20080144036A1 (en) *  20061219  20080619  Asml Netherlands B.V.  Method of measurement, an inspection apparatus and a lithographic apparatus 
US20080148875A1 (en) *  20061220  20080626  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7924435B2 (en)  20061222  20110412  Zygo Corporation  Apparatus and method for measuring characteristics of surface features 
US20080174784A1 (en) *  20061222  20080724  Zygo Corporation  Apparatus and method for measuring characteristics of surface features 
US7916927B2 (en)  20070116  20110329  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080170780A1 (en) *  20070116  20080717  Asml Netherlands B.V.  Inspection method and apparatus,lithographic apparatus, lithographic processing cell and device manufacturing method 
US7619737B2 (en)  20070122  20091117  Asml Netherlands B.V  Method of measurement, an inspection apparatus and a lithographic apparatus 
US20080174753A1 (en) *  20070122  20080724  Asml Netherlands B.V.  Method of measurement, an inspection apparatus and a lithographic apparatus 
US20080180656A1 (en) *  20070126  20080731  KlaTencor Technologies Corporation  Surface characteristic analysis 
US7554654B2 (en)  20070126  20090630  KlaTencor Corporation  Surface characteristic analysis 
US9025162B2 (en)  20070131  20150505  Zygo Corporation  Interferometry for lateral metrology 
US20080180685A1 (en) *  20070131  20080731  Zygo Corporation  Interferometry for lateral metrology 
US7889355B2 (en)  20070131  20110215  Zygo Corporation  Interferometry for lateral metrology 
US7852459B2 (en)  20070202  20101214  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7839506B2 (en)  20070221  20101123  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20080198380A1 (en) *  20070221  20080821  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20100157299A1 (en) *  20070221  20100624  Asml Netherlands B.V.  Inspection Method and Apparatus, Lithographic Apparatus, Lithographic Processing Cell and Device Manufacturing Method 
US7701577B2 (en)  20070221  20100420  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7724370B2 (en)  20070301  20100525  Asml Netherlands B.V.  Method of inspection, a method of manufacturing, an inspection apparatus, a substrate, a mask, a lithography apparatus and a lithographic cell 
US20080212097A1 (en) *  20070301  20080904  Asml Netherlands B.V.  Method of inspection, a method of manufacturing, an inspection apparatus, a substrate, a mask, a lithography apparatus and a lithographic cell 
US7599064B2 (en)  20070307  20091006  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method, substrate for use in the methods 
US7605907B2 (en)  20070327  20091020  Asml Netherlands B.V.  Method of forming a substrate for use in calibrating a metrology tool, calibration substrate and metrology tool calibration method 
US20080239277A1 (en) *  20070327  20081002  Asml Netherlands B.V.  Method of forming a substrate for use in calibrating a metrology tool, calibration substrate and metrology tool calibration method 
US8031337B2 (en)  20070330  20111004  Asml Netherlands B.V.  Angularly resolved scatterometer 
US7570358B2 (en)  20070330  20090804  Asml Netherlands Bv  Angularly resolved scatterometer, inspection method, lithographic apparatus, lithographic processing cell device manufacturing method and alignment sensor 
US20080239265A1 (en) *  20070330  20081002  Asml Netherlands B.V.  Angularly resolved scatterometer, inspection method, lithographic apparatus, lithographic processing cell device manufacturing method and alignment sensor 
US20080239318A1 (en) *  20070330  20081002  Asml Netherlands B.V.  Method of measuring asymmetry in a scatterometer, a method of measuring an overlay error in a substrate and a metrology apparatus 
US20090262366A1 (en) *  20070330  20091022  Asml Netherlands B.V.  Angularly Resolved Scatterometer 
US7656518B2 (en)  20070330  20100202  Asml Netherlands B.V.  Method of measuring asymmetry in a scatterometer, a method of measuring an overlay error in a substrate and a metrology apparatus 
US20080279442A1 (en) *  20070509  20081113  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US8189195B2 (en)  20070509  20120529  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US7835561B2 (en)  20070518  20101116  Visiongate, Inc.  Method for image processing and reconstruction of images for optical tomography 
US7911612B2 (en)  20070613  20110322  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20090021723A1 (en) *  20070719  20090122  Zygo Corporation  Generating model signals for interferometry 
US7619746B2 (en)  20070719  20091117  Zygo Corporation  Generating model signals for interferometry 
US7460237B1 (en)  20070802  20081202  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell and device manufacturing method 
US20090059236A1 (en) *  20070830  20090305  KlaTencor Technologies Corporation  Wafer Edge Inspection 
US7656519B2 (en)  20070830  20100202  KlaTencor Corporation  Wafer edge inspection 
US20090147268A1 (en) *  20071012  20090611  Zygo Corporation  Interferometric analysis of underresolved features 
US8072611B2 (en)  20071012  20111206  Zygo Corporation  Interferometric analysis of underresolved features 
US7933010B2 (en)  20071022  20110426  Rahn J Richard  Depth of field extension for optical tomography 
US7787112B2 (en)  20071022  20100831  Visiongate, Inc.  Depth of field extension for optical tomography 
US20100321786A1 (en) *  20071022  20101223  Visiongate, Inc.  Depth of field extension for optical tomography 
EP2053349A2 (en)  20071025  20090429  ASML Netherlands B.V.  Method and apparatus for determining properties of a lithographic pattern on a substrate 
US20090168062A1 (en) *  20071025  20090702  Asml Netherlands B.V.  Inspection Method and Apparatus, Lithographic Apparatus, Lithographic Processing Cell, and Device Manufacturing Method to Measure a Property of a Substrate 
US8115926B2 (en)  20071025  20120214  Asml Netherlands B.V.  Inspection method and apparatus, lithographic apparatus, lithographic processing cell, and device manufacturing method to measure a property of a substrate 
US7978337B2 (en)  20071113  20110712  Zygo Corporation  Interferometer utilizing polarization scanning 
US20090128827A1 (en) *  20071113  20090521  De Groot Peter  Interferometer utilizing polarization scanning 
US8126677B2 (en)  20071214  20120228  Zygo Corporation  Analyzing surface structure using scanning interferometry 
US20090208072A1 (en) *  20080218  20090820  Visiongate, Inc.  3d imaging of live cells with ultraviolet radiation 
US8143600B2 (en)  20080218  20120327  Visiongate, Inc.  3D imaging of live cells with ultraviolet radiation 
US8368035B2 (en)  20080218  20130205  Visiongate Inc.  3D imaging of live cells with ultraviolet radiation 
US20110032500A1 (en) *  20080320  20110210  Asml Netherlands B.V.  Inspection apparatus for lithography 
US8681312B2 (en)  20080320  20140325  Asml Netherlands B.V.  Inspection apparatus for lithography 
US8792096B2 (en)  20080324  20140729  Asml Netherlands B.V.  Inspection apparatus for lithography 
EP2172934A1 (en)  20081006  20100407  Thomson Licensing  Optical storage medium comprising marks with different orientation, and respective apparatus for reading of data. 
US8902431B2 (en)  20081126  20141202  Zygo Corporation  Low coherence interferometry with scan error correction 
US7978338B2 (en)  20081126  20110712  Zygo Corporation  Compound reference interferometer 
US20100128276A1 (en) *  20081126  20100527  Zygo Corporation  Compound reference interferometer 
US20100128283A1 (en) *  20081126  20100527  Zygo Corporation  Interferometric systems and methods featuring spectral analysis of unevenly sampled data 
US8004688B2 (en)  20081126  20110823  Zygo Corporation  Scan error correction in low coherence scanning interferometry 
US20100128278A1 (en) *  20081126  20100527  Zygo Corporation  Fiberbased interferometer system for monitoring an imaging interferometer 
US8379218B2 (en)  20081126  20130219  Zygo Corporation  Fiberbased interferometer system for monitoring an imaging interferometer 
US20100128280A1 (en) *  20081126  20100527  Zygo Corporation  Scan error correction in low coherence scanning interferometry 
US8120781B2 (en)  20081126  20120221  Zygo Corporation  Interferometric systems and methods featuring spectral analysis of unevenly sampled data 
US20100214639A1 (en) *  20090223  20100826  Visiongate, Inc.  Optical tomography system with highspeed scanner 
US8254023B2 (en)  20090223  20120828  Visiongate, Inc.  Optical tomography system with highspeed scanner 
US8090183B2 (en)  20090312  20120103  Visiongate, Inc.  Pattern noise correction for pseudo projections 
US20100232664A1 (en) *  20090312  20100916  Visiongate, Inc.  Pattern noise correction for pseudo projections 
US20110032534A1 (en) *  20090519  20110210  Camtek Ltd  System and a method for broadband interferometry 
US8731882B2 (en) *  20090924  20140520  Asml Netherlands B.V.  Methods and apparatus for modeling electromagnetic scattering properties of microscopic structures and methods and apparatus for reconstruction of microscopic structures 
US20110218789A1 (en) *  20090924  20110908  Asml Netherlands B.V.  Methods and Apparatus for Modeling Electromagnetic Scattering Properties of Microscopic Structures and Methods and Apparatus for Reconstruction of Microscopic Structures 
US8645109B2 (en) *  20091208  20140204  Asml Netherlands B.V.  Methods and apparatus for determining electromagnetic scattering properties and structural parameters of periodic structures 
US20110137625A1 (en) *  20091208  20110609  Asml Netherlands B.V.  Methods and Apparatus for Determining Electromagnetic Scattering Properties and Structural Parameters of Periodic Structures 
US20120250026A1 (en) *  20110330  20121004  Boer Guido Do  Interferometer module 
US9690215B2 (en)  20110330  20170627  Mapper Lithography Ip B.V.  Interferometer module 
US9678443B2 (en)  20110330  20170613  Mapper Lithography Ip B.V.  Lithography system with differential interferometer module 
US9261800B2 (en)  20110330  20160216  Mapper Lithography Ip B.V.  Alignment of an interferometer module for use in an exposure tool 
US9069265B2 (en) *  20110330  20150630  Mapper Lithography Ip B.V.  Interferometer module 
US8468471B2 (en)  20110923  20130618  KlaTencor Corp.  Process aware metrology 
US8832611B2 (en)  20110923  20140909  KlaTencor Corp.  Process aware metrology 
US20150241200A1 (en) *  20120927  20150827  Mapper Lithography Ip B.V.  Multiaxis differential interferometer 
US9551563B2 (en) *  20120927  20170124  Mapper Lithography Ip B.V.  Multiaxis differential interferometer 
US9897553B2 (en)  20130221  20180220  Nova Measuring Instruments Ltd.  Optical phase measurement method and system 
EP2959262A4 (en) *  20130221  20160928  Nova Measuring Instr Ltd  Optical phase measurement method and system 
WO2014128710A1 (en) *  20130221  20140828  Nova Measuring Instruments Ltd.  Optical phase measurement method and system 
CN105143814A (en) *  20130221  20151209  诺威量测设备股份有限公司  Optical phase measurement method and system 
Similar Documents
Publication  Publication Date  Title 

Dubois et al.  Improved threedimensional imaging with a digital holography microscope with a source of partial spatial coherence  
US6556306B2 (en)  Differential time domain spectroscopy method for measuring thin film dielectric properties  
US4999014A (en)  Method and apparatus for measuring thickness of thin films  
US6667809B2 (en)  Scanning interferometric nearfield confocal microscopy with background amplitude reduction and compensation  
US5866935A (en)  Tunneling device  
US6775009B2 (en)  Differential interferometric scanning nearfield confocal microscopy  
US7315382B2 (en)  Interferometry method for ellipsometry, reflectometry, and scatterometry measurements, including characterization of thin film structures  
US20040239954A1 (en)  Resolution enhanced optical metrology  
US4674824A (en)  System for enhancement of optical features  
US7177030B2 (en)  Determination of thin film topography  
US6469788B2 (en)  Coherent gradient sensing ellipsometer  
Cuche et al.  Aperture apodization using cubic spline interpolation: application in digital holographic microscopy  
US20080049233A1 (en)  MultipleAngle MultipleWavelength Interferometer Using HighNA Imaging and Spectral Analysis  
Harasaki et al.  Fringe modulation skewing effect in whitelight vertical scanning interferometry  
Hettwer et al.  Three channel phaseshifting interferometer using polarizationoptics and a diffraction grating  
US7064829B2 (en)  Generic interface for an optical metrology system  
US6445453B1 (en)  Scanning interferometric nearfield confocal microscopy  
US20070091317A1 (en)  Method and apparatus for optically analyzing a surface  
Lawson et al.  Specification of optical components using the power spectral density function  
US5164790A (en)  Simple CD measurement of periodic structures on photomasks  
US6781701B1 (en)  Method and apparatus for measuring optical phase and amplitude  
Nyyssonen  Theory of optical edge detection and imaging of thick layers  
US20060274325A1 (en)  Method of qualifying a diffraction grating and method of manufacturing an optical element  
Kujawinska et al.  High accuracy Fourier transform fringe pattern analysis  
US6999178B2 (en)  Spatialheterodyne interferometry for reflection and transmission (SHIRT) measurements 
Legal Events
Date  Code  Title  Description 

AS  Assignment 
Owner name: CALIFORNIA INSTITUTE OF TECHNOLOGY, CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:MARX, DAVID;PSALTIS, DEMETRI;REEL/FRAME:008036/0181 Effective date: 19960601 

FPAY  Fee payment 
Year of fee payment: 4 

FPAY  Fee payment 
Year of fee payment: 8 

REMI  Maintenance fee reminder mailed  
FPAY  Fee payment 
Year of fee payment: 12 

SULP  Surcharge for late payment 
Year of fee payment: 11 