US6661004B2 - Image deconvolution techniques for probe scanning apparatus - Google Patents

Image deconvolution techniques for probe scanning apparatus Download PDF

Info

Publication number
US6661004B2
US6661004B2 US09/789,992 US78999201A US6661004B2 US 6661004 B2 US6661004 B2 US 6661004B2 US 78999201 A US78999201 A US 78999201A US 6661004 B2 US6661004 B2 US 6661004B2
Authority
US
United States
Prior art keywords
probe
sample
image
estimate
shape
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 - Fee Related, expires
Application number
US09/789,992
Other versions
US20010038072A1 (en
Inventor
Bernardo D. Aumond
Kamal Youcef-Toumi
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Massachusetts Institute of Technology
Original Assignee
Massachusetts Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Massachusetts Institute of Technology filed Critical Massachusetts Institute of Technology
Priority to US09/789,992 priority Critical patent/US6661004B2/en
Assigned to MASSACHUSETTS INSTITUTE OF TECHNOLOGY reassignment MASSACHUSETTS INSTITUTE OF TECHNOLOGY ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: AUMOND, BERNARDO D., YOUCEF-TOUMI, KAMAL
Publication of US20010038072A1 publication Critical patent/US20010038072A1/en
Application granted granted Critical
Publication of US6661004B2 publication Critical patent/US6661004B2/en
Adjusted expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • G06T5/73
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01QSCANNING-PROBE TECHNIQUES OR APPARATUS; APPLICATIONS OF SCANNING-PROBE TECHNIQUES, e.g. SCANNING PROBE MICROSCOPY [SPM]
    • G01Q30/00Auxiliary means serving to assist or improve the scanning probe techniques or apparatus, e.g. display or data processing devices
    • G01Q30/04Display or data processing devices
    • G01Q30/06Display or data processing devices for error compensation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01QSCANNING-PROBE TECHNIQUES OR APPARATUS; APPLICATIONS OF SCANNING-PROBE TECHNIQUES, e.g. SCANNING PROBE MICROSCOPY [SPM]
    • G01Q40/00Calibration, e.g. of probes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration by the use of local operators
    • G06T5/30Erosion or dilatation, e.g. thinning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S977/00Nanotechnology
    • Y10S977/84Manufacture, treatment, or detection of nanostructure
    • Y10S977/849Manufacture, treatment, or detection of nanostructure with scanning probe
    • Y10S977/85Scanning probe control process
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S977/00Nanotechnology
    • Y10S977/84Manufacture, treatment, or detection of nanostructure
    • Y10S977/849Manufacture, treatment, or detection of nanostructure with scanning probe
    • Y10S977/86Scanning probe structure
    • Y10S977/863Atomic force probe
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S977/00Nanotechnology
    • Y10S977/84Manufacture, treatment, or detection of nanostructure
    • Y10S977/849Manufacture, treatment, or detection of nanostructure with scanning probe
    • Y10S977/86Scanning probe structure
    • Y10S977/875Scanning probe structure with tip detail
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S977/00Nanotechnology
    • Y10S977/84Manufacture, treatment, or detection of nanostructure
    • Y10S977/88Manufacture, treatment, or detection of nanostructure with arrangement, process, or apparatus for testing
    • Y10S977/881Microscopy or spectroscopy, e.g. sem, tem

Definitions

  • This invention relates generally to profilometry and more particularly, to image deconvolution techniques for probe scanning devices.
  • a structure having one characteristic dimension e.g. height
  • the others e.g. width and length
  • a high aspect ratio structure examples include probes of atomic force and scanning tunnelling microscopes, field emission probes, micro-indenters and Micro Electro-Mechanical systems (MEM's) structures.
  • Other high aspect ratio structures are found in quantum magnetic media for high-density data storage, compact disk stampers, crystalline structures, blades and biological systems such as virus particles.
  • Such high aspect ratio structures have typical heights on the order of tens of micrometers and tip radii in the range of tens of nanometers. Further, these structures may or may not be conductive.
  • Obtaining images of high aspect ratio structures poses certain challenges.
  • Desired image resolutions can be on the order of 1 nm in the vertical direction (i.e. a direction which is normal to a measurement surface) and 10 nm in the lateral direction (i.e. a direction which is parallel to the measurement surface).
  • Some profilometry methods are unable to provide such resolutions and thus such imaging criteria prevent the use of certain types of profilometry methods.
  • a stylus having a sharp tip and coupled to a hinged arm is mechanically dragged along the sample surface.
  • the deflection of the hinged stylus arm is measured and recorded as the surface profile.
  • the use of a hinged stylus arm allows measurement of very rough surfaces, for example those with peak-to-peak heights greater than 1 mm.
  • Probe-to-surface contact forces range from 10 ⁇ 3 N to 10 ⁇ 6 N.
  • physical rigidity limits the minimum stylus tip radius and hence the lateral resolution to about 0.1 mm.
  • optical profilometry many different optical phenomena (such as interference and internal reflection) can be utilized.
  • the most popular technique is based on phase-measuring interferometry, in which a light beam reflecting off the sample surface is interfered with a phase-varied reference beam. The surface profile is deduced from the resulting fringe patterns. With a collimated light beam and a large photodetector array, the entire surface can be profiled simultaneously.
  • This and other conventional optical profilometry methods are limited in lateral resolution by the minimum focussing spot size of about 0.5 ⁇ m (for visible light).
  • measurement values are dependent upon the surface reflectivity of the material being profiled.
  • a scanning probe microscope In the current state of the art, only scanning probe microscopes can meet a 10 nm lateral resolution requirement. In these microscopes, an atomically sharp (or nearly so) tip at a very close spacing to the sample surface is moved over the surface using a piezoactuator.
  • One type of scanning probe microscope is the atomic force microscope (AFM), which measures the topography of a surface with a probe that has a very sharp tip.
  • a probe assembly includes a cantilever beam from which the probe, or microstylus extends. The probe terminates at the probe tip having a typical tip radius of less than 0.1 ⁇ m.
  • the probe typically has a length on the order of a couple of micrometers and the cantilever beam typically has a length between 100 ⁇ m and 200 ⁇ m.
  • the probe In a contact mode atomic force microscope, the probe is moved relative to the surface of a sample and deflection of the cantilever is measured to provide a measure of the surface topography. More particularly, a laser beam is directed toward, and reflects off the back surface of the cantilever to impinge upon a sensor, such as a photodetector array. The electrical output signals of the photodetector array provide a topographical image of the sample surface and, further, provide feedback signals to a fine motion actuator, sometimes provided in the form of a piezoelectric actuator. In a constant force contact AFM, the fine motion actuator is responsive to the feedback signals for maintaining a substantially constant force between the probe tip and the sample, such as forces on the order of 10 ⁇ 8 N to 10 ⁇ 11 N.
  • Contact atomic force microscopy offers high lateral and vertical resolutions, such as less than 1 nm vertical resolution and less than 50 nm lateral resolution. Further, since the contact AFM relies on contact forces rather than on magnetic or electric surface effects, advantageously the contact AFM can be used to profile conductive and non-conductive samples. However, the maximum surface roughness that can be profiled is much less than that of conventional stylus profilometers which use a linear variable differential transducer (LVDT).
  • LVDT linear variable differential transducer
  • van der Waals forces are measured by vibrating the cantilever near its resonance frequency and detecting the change in the vibrational amplitude of a laser beam reflected off the cantilever due to a change in the force gradient caused by changes in the surface profile.
  • the non-contact atomic force microscope offers non-invasive profiling.
  • the technique has some disadvantages when compared to contact atomic force microscopy.
  • van der Waals forces are hard-to-measure weak forces, rendering the microscope more susceptible to noise.
  • the probe tip must be maintained at a fixed height above the sample, typically on the order of a few nanometers, and the feedback control necessary to maintain this spacing must operate slowly to avoid crashing the probe tip on the sample.
  • the effective tip radius is increased and hence the achievable lateral resolution is decreased.
  • AFM was primarily developed for high-resolution 3-D imaging (profilometry) of atomically flat samples. In that case, the probe tip is scanned over the sample and only the apical region of the probe interacts with the profiled surface. Therefore, AFM images will closely reproduce the topography regardless of distortions in the probe away from the apex. Accordingly, when stylus instruments are used in profilometry, the implicit assumption is that only the very apex of the stylus touches the surface at all points.
  • the AFM and stylus images can be quite different from the real topography. That is also the case when the dimensions of the sample are comparable to those of the employed probe (AFM probe or stylus).
  • the reason for this deviation is that areas of the probe other than the apex (for instance, the probe sides) interact with the sample as well.
  • the image distortion caused by the interaction of the probe with the surface is typically referred to as image convolution.
  • the characteristic dimensions of the samples or the features of interest lie well within the micrometer and sub-micrometer ranges. These fields include but are not limited to nanotechnology, micro-electromechanical systems (MEMS), semiconductor devices and storage media, micro-sensors, and blade fabrication.
  • MEMS micro-electromechanical systems
  • the investigated features could be photo-resist trenches in silicon wafers, memory pillars in quantum magnetic media devices, roughness in smooth optical surfaces, the radius of curvature of field emission probes and parts of micro-machines.
  • images of such structures can be distorted by convolution errors.
  • the level of convolution is greatly dependent on the relative size and shape of the employed probe with respect to the sample. Different probes can interact with the sample generating different distorted images. Convolution seriously reduces metrology accuracy. The deviations in the measurement of radius of curvature of high aspect ratio samples is proportional to the radius of the probe. Therefore, samples with dimensions similar to those of the probe will generate images with close to a 100% radius measurement error. Width measurements taken from images of photo-resist trenches are embedded with deviations proportional to the square of the height of the probe (for nearly parabolic probes). Probe geometry and size can also affect texture parameter (e.g. roughness) measurements due to convolution.
  • texture parameter e.g. roughness
  • Deconvolution is necessary in applications in which the sample feature dimensions are in the same range as the size of the probe tip, therefore resulting in strong image convolution.
  • Prior art deconvolution techniques include scanning a “standard sample” (i.e. a sample having a known surface shape) with a probe to provide a mean curve for the probe shape, an inner curve for the probe shape, and an outer curve for the probe shape, thereby characterising the probe shape.
  • a “standard sample” i.e. a sample having a known surface shape
  • Prior art deconvolution techniques include scanning a “standard sample” (i.e. a sample having a known surface shape) with a probe to provide a mean curve for the probe shape, an inner curve for the probe shape, and an outer curve for the probe shape, thereby characterising the probe shape.
  • One problem with this approach is that errors in the initial probe calibration propagate to all other measurements.
  • Another problem with this approach is that the shape of the probe change and to problems over time because of physical changes in probe due to probe wear.
  • an apparatus utilizing contact atomic force microscopy includes a probe, a controller operative to move the probe into a first vantage point relative to the sample to produce a first image, and to move the probe into a second vantage point relative to the sample to produce a second image.
  • the apparatus further includes a deconvolution processor which deconvolves the first image and the second image to produce an image of the sample.
  • a deconvolution method includes moving a probe in a first scanning pattern to generate a first image of the sample, changing said probe's vantage point, moving the probe in a second scanning pattern to generate a second image of the feature, and deconvolving the image of the sample using the first image and the second image.
  • the deconvolution technique utilizes an iterative process. After generating multiple images of a sample from different vantage points, the process begins by generating two or more estimates of the probe shape using a first image. The estimates can be generated using a blind deconvolution method. Next using an erosion technique, each of the estimates of the probe shapes is used to obtain a corresponding estimate of the sample shape. Thus if two estimates of the probe shape are generated, then two estimates of the sample shape are obtained. At least two estimates of the probe shape are then combined to provide a new probe estimate. Similarly, at least two estimates of the sample shape are combined to provide a new sample shape estimate. The process of combining sample estimates (including newly computed sample estimates) to generate more sample estimates can be repeated any desired number of times.
  • the new sample shape estimates can be utilized to generate new probe shape estimates.
  • the process of generating new probe shape estimates can also be repeated any desired number of times. In one embodiment, the above iterative process can be repeated until the changes in the newly computed sample shape and probe shape estimates are below a predetermined threshold when compared with the estimates from previous iteration. With such an arrangement an accurate estimate to the sample shape is obtained without the need for probe characterization. If additional images used, the fidelity of the estimated sample is increased.
  • the deconvolution technique utilizes Legendre transforms in the processing of images to produce sample and probe measurements. After generating two images of a sample from different vantage points, the Legendre transform of the first and second image are obtained. These transforms are used in a system of equations that relate the transform of the sample with a first transform of the probe and the transform of the first image, the transform of the sample with a transform of the probe from the changed vantage point and the transform of the second image and relating the first transform of the probe, the transform of the probe from the changed vantage point and the angle between the first and second vantage point. After obtaining a Legendre transform of the first image and obtaining a Legendre transform of the second image, a parametric function is provided to describe the probe geometry and then by establishing how the upright and the rotated versions are geometrically related.
  • the parametric function provides a functional dependency between the probe transform in the first orientation and the probe transform from the second vantage point which allows a Legendre transform of the probe to be eliminated from the system of equations.
  • Solving the system of equations by a least squares algorithm provides a Legendre transform of the sample.
  • the sample is recovered from the Legendre transform of the sample.
  • the probe geometry can also be recovered from the solution for the sample shape. With such an arrangement a partial but exact reconstruction of the sample shape is obtained without the need for probe characterization.
  • FIG. 1 shows a schematic diagram of an atomic force microscope
  • FIGS. 2A-2C demonstrate image deconvolution according to the present invention
  • FIG. 3 is a flow diagram illustrating the erosion and dilation deconvolution technique according to an aspect of the invention.
  • FIG. 4 is a flow diagram illustrating the Legendre transform deconvolution technique according to an aspect of the invention.
  • a contact mode atomic force microscope (AFM) 10 includes a probe assembly 12 having a probe 14 mounted to the free end of a cantilever beam 16 .
  • the probe assembly 12 has a probe tip 18 adapted to contact a sample surface 20 of the sample 22 to be profiled.
  • a sample 22 is disposed on a sample holder 24 . It should be understood, however, that in other applications, sample 22 may be held by any fixture or surface. Any means of locating the sample 22 proximate the probe 14 may be used.
  • a typical radius of the probe tip 18 is between 5 nm and 100 nm, a typical length of the probe 14 is on the order of 4.0 ⁇ m, and a typical length of the cantilever beam 16 is between 100 and 200 ⁇ m. It should be appreciated, however, that the present invention can be used in AFM's or other devices having probes and probe tips of any dimension.
  • the probe assembly 12 and sample 22 are moved relative to one another, causing the probe tip 18 to move over the sample surface 20 . Movement of the probe tip 18 over features including high aspect ratio structures 26 on the sample surface 20 causes the cantilever beam 16 to deflect. The deflection of the cantilever beam 16 is measured to provide a measure of the surface topography. More particularly, a laser beam 30 generated by a laser source 32 is directed towards, and reflects off a surface 16 a of the cantilever beam 16 . The reflected laser signal 31 impinges upon a sensor, or detector 40 , which may, for example, be provided as a photodetector array. In response to the reflected laser signal 31 , the detector 40 generates an output signal 42 which is provided to a signal processor 44 .
  • Signal processor 44 processes the output signal 42 to provide processed signals 46 , portions of which correspond to feed back control signals 48 and portions of which correspond to the topographical sample image signals 50 .
  • the control signals 48 are fed to a piezoelectric actuator 49 and the image signals 50 are fed to a display 52 on which the sample topographical image and other images can be displayed.
  • a topographical sample view of the image 54 can be viewed on a display 52 .
  • the signal processor 44 can include a deconvolution processor 45 .
  • the display 52 is optional and the signal processor 44 operations can be automated.
  • the signal processor 44 may process the resulting image itself to perform certain sample measurements, such as measuring the radius of the tip of a feature of the sample.
  • the functions of the signal processor 44 and the deconvolution processor 45 can be provided by a separate processor, or by computer software instructions or groups of instructions executed on another processor which is part of the AFM.
  • Such processing may be performed by a single processing apparatus which may, for example, be provided as part of microscope or measurement system such as that to be described below in conjunction with methods described in FIGS. 3 and 4.
  • the processing steps of the signal processor 44 and the deconvolution processor 45 can be performed by functionally equivalent circuits such as a digital signal processor circuit or an application specific integrated circuit (ASIC).
  • ASIC application specific integrated circuit
  • AFM cantilever beams 16 are inclined with respect to the sample plane in order to reflect the laser beam 30 into the detector 40 .
  • the probe tip 18 will touch one side of the high aspect ratio structure 26 , but possibly not the other side.
  • the sample holder 24 is inclined to the same degree as the cantilever beam 16 so that the probe 14 and sample holder 24 are initially normal to one another.
  • the probe tip 18 should be normal to the sample surface 20 being measured.
  • the method for image deconvolution uses multiple imaging or stereo microscope imaging.
  • the motivation for such an approach is drawn from the understanding that if the scanning is performed with the probe 14 always normal to the sample surface 20 , convolution errors are minimized. Since this is not always practical given the available AFM configurations, one approximation to the ideal configuration is obtained by obtaining multiple images in which the angle between the probe tip 18 and the sample 22 is different. This may be accomplished, for example by utilizing different sample tilt angles. The images are then combined to form the deconvolved sample image.
  • the angle between an a central longitudinal axis 56 of the probe 14 and a plane 58 parallel with the plane of the sample holder 24 corresponds to an initial angle ⁇ i .
  • the initial angle ⁇ i is measured with respect to a central longitudinal axis 56 of probe 14 and an axis formed by the intersection of plane 58 and a plane normal to the cantilever beam 16 and containing axis 56 .
  • a first image, image 1 61 is obtained by scanning the probe tip 18 over the sample surface 20 .
  • a second image, image 2 62 is obtained by tilting sample holder 24 causing the sample 22 ′ (tilted sample 22 ) to be imaged by the probe 14 at relative angle 74 T1 . Portions of the image 1 61 and image 2 62 can be presented on the display 52 as topographic images of the sample.
  • the sample images image 1 61 and image 2 62 are processed in accordance with techniques to be described below in conjunction with FIGS. 2-4.
  • FIGS. 2A-2C images obtained from measuring the sample 22 (FIG. 1) in a single scan in the X direction are shown.
  • the Z dimension represents heights and repetitive scans are in the X direction.
  • Curve 64 represents an image of a portion of the sample 22 which would be generated in the absence of convolution errors.
  • the portion of the sample 22 shown in FIG. 2A has a relatively high aspect ratio.
  • image 1 61 is obtained by scanning the sample 22 with the probe 14 in a first position.
  • the first position corresponds to the initial normal position as discussed above in conjunction with FIG. 1 .
  • the first position can be any position as long as subsequent positions (probe vantage points) differ from it (i.e. as long as the relative angle between the probe and the sample holder are different).
  • Curve probe 1 14 ′ represents a theoretical image of the probe 14 (FIG. 1) if the probe 14 could be perfectly characterized. When the probe 14 physically contacts the sample at point 14 b ′ the point that appears in image 1 61 is point 14 a ′.
  • the convolution error is the difference between the true probe-sample contact point 14 b ′ and the imaged point 14 a ′.
  • additional imaged points 14 e ′- 14 g ′ are generated to form image 1 61 with similar convolution errors. It should be appreciated that the number of points generated in the images depends upon the selected resolution of the microscope and the maximum resolution of the microscope scanner.
  • FIG. 2B shows a second image, image 2 62 , obtained by scanning the sample 22 (FIG. 1) with the probe 14 (FIG. 1) in a second position relative to the sample 22 (FIG. 1) which is different that that used to generate image 1 61 .
  • image 2 62 is generated when the sample holder 24 (FIG. 1) is tilted to an angle ⁇ T1 relative to a central longitudinal axis 56 of probe 14 and an axis 58 ′ formed by the tilted sample holder 24 ′ on which the tilted sample 22 ′ is disposed as shown in FIG. 1 .
  • the same result can be achieved by titling the probe 14 (FIG.
  • FIG. 2C shows a reconstructed image 66 using the deconvolution techniques described below.
  • the degree to which reconstructed image 66 approximates the true sample topography 64 depends upon a variety of factors including but not limited to the number of images processed, the deconvolution method, the extension of areas of the sample that were actually touched by the probe during scanning and the accuracy of any geometrical estimates of probe shape, including parametric geometric functions.
  • Two deconvolution techniques are described below in embodiments which can be used to recover the true sample shape using multiple images.
  • a first embodiment uses a minimum envelope method and a second embodiment uses a Legendre transformation-based reconstruction method.
  • a minimum envelope approach based on the concept of blind deconvolution is used to recover a deconvolved image.
  • the minimum envelope approach two or more images are taken with different vantage points (obtained by changing the relative angle between the probe and sample as described above).
  • the true sample profile for example curve 64 (FIG. 2A) is always bounded by the image sample, image 1 61 (FIG. 2 A).
  • a blind deconvolution approach is used to find initial estimates for probe and sample topographies based on the measured first image, image 1 61 , and the second image, image 2 62 .
  • the probe topography is bounded by the intersection of all estimates.
  • the sample topography is bounded by the intersections of all sample estimates obtained with the blind deconvolution approach.
  • step 70 a first image of a sample, image 1 61 , is obtained with a first relative angle between the probe and the sample.
  • step 74 a second image of the sample, image 2 62 is obtained with a different relative angle between the sample 22 and the probe 14 that is charged relative to the configuration used in step 72 .
  • the images, image 1 61 and image 2 62 are obtained from different vantage points. It should be appreciated that either the sample can be tilted, the probe can be tilted, or both the probe and sample can be tilted by different amounts to produce different angles between the probe and the sample.
  • a first blind estimate of the probe EstP 1 based on the first image 61 is obtained using Villarubia's blind deconvolution method based on the probe shape 120 as is known in the art.
  • the details of the blind estimation technique are explained in J. S. Villarubia, “Morphological Estimation of Tip Geometry for Scanned Probe Microscopy,” Surf. Sci. 321, 287-300 (1994) and G. S. Pingali, R. Jain, “Restoration of Scanning Probe Microscope Images,” IEEE Workshop on Applications of Computer Vision 282 (1992).
  • a second blind estimate of the probe EstP 2 is obtained using Villarubia's blind deconvolution method based on image 2 62 .
  • These blind estimates EstP 1 and EstP 2 are blunter than the real probe 14 and vary from the real shape, surface sample 20 , because the images do not offer enough high aspect ratio features for a good estimation. The lack of high aspect ratio features leads to poor shape estimates.
  • Step 80 uses erosion and probe estimate EstP 1 to find new sample estimate, EstS 1 (estimate 1 ) from the first sample image and step 82 uses erosion and probe estimate EstP 2 to find new sample estimate EstS 2 (estimate 2 ) from the second sample image.
  • the details of the erosion process are provided in J. S. Villarubia, “Morphological Estimation of Tip Geometry for Scanned Probe Microscopy,” Surf. Sci. 321, 287-300 (1994) and G. S. Pingali, R. Jain, “Restoration of Scanning Probe Microscope Images,” IEEE Workshop on Applications of Computer Vision 282 (1992).
  • step 86 the two estimates of the sample are combined to generate a better estimate of the sample surface topography.
  • the new sample estimate is used to “sharpen” the probe estimate, which is then used to iteratively create better sample estimates.
  • Steps 88 to 102 iteratively repeat the process until no change in the estimation above a predetermined threshold is calculated.
  • step 88 the new probe estimate from step 84 is used to obtain a new sample estimate 1 by erosion.
  • step 90 the new probe estimate from step 84 is used to obtain a new sample estimate 2 by erosion.
  • step 92 the two new sample estimates are combined to generate another more accurate new sample estimate.
  • the sample estimate for the tilted case must be rotated by minus ⁇ degrees rendering a rotated estimate denoted as EstS 1 *.
  • EstS 1 * The intersection is then found and will be the new estimate EstS 0new .
  • EstS 1 * are obtained by the erosion method. Since the estimate for the tilted case EstS 1 was obtained for the tilted sample, in order to combine the estimations, EstS 1 must be rotated back to render an estimate at 0 degrees (EstS 1 *).
  • the new sample estimates are used to obtain two new probe estimates.
  • step 100 the two new Probe estimates EStP 2new and EStP 1new are combined, generating a new sharper probe estimate EstP new .
  • Decision block 102 implements a loop in which the process of generating new probe and sample estimates in steps 88 to 102 are repeated.
  • new probe estimate EstP new new sample estimates EStS 0new′ , EstS new , can be generated by erosion. This is also described in B. D. Aumond and K. Youcef-Toumi, “Experimental High Precision Profilometry of High Aspect Ratio Samples,” Proceedings of the 1998 IEEE International Conference on Systems, Man, and Cybernetics, San Diego, Calif. (November 1998), Y. Yeo, B. D. Aumond, and K. Youcef-Toumi, “Precision Atomic Force Microscope Imaging,” Proceedings of the 2000 IEEE International Conference on Signal Processing, Beijing, China (August 2000), and B. D. Aumond, K. Youcef-Toumi, “High Precision Stereo Profilometry based on Atomic Force Microscopy Technology,” Mechatronics 2000 Conference, Atlanta, Ga. (September 2000).
  • the new sample estimates are combined to generate a better estimate of the real surfaces for the upright and tilted measurements.
  • This new surface estimate is in turn used to sharpen the probe estimates using both the first and second images (e.g. the upright sample image and tilted sample image).
  • the sharpened probes estimates are combined generating a new probe estimate EstP new′ .
  • a test is performed to determine whether the further new sample estimate or the still further new probe estimate changed from the last iteration. If the change is lower than a predetermined threshold, then the deconvolution process is complete as shown in step 104 . Otherwise the erosion process is repeated beginning at step 88 .
  • the sample combination steps 88 , 90 and 92 and probe sharpening steps 94 , 96 , 100 will eventually render better estimates of probe geometry and sample topography. Thus, when compared to single image blind deconvolution, the results are improved.
  • an underscanning technique (scanning the underside of the image with the probe reflected about the origin) as described in Y. Yeo, B. D. Aumond, and K. Youcef-Toumi, “ Precision Atomic Force Microscope Imaging, ” Proceedings of the 2000 IEEE International Conference on Signal Processing, Beijing, China (March 2000).
  • the reflected probe geometry (the geometry of the probe is reflected horizontally and then vertically) scans of the underside of the image instead of the real sample.
  • the reflected probe to simulate the scanning of the underside of the image very much in the same way that a probe scans a real sample.
  • Record the reflected probe apex positions as it scans the image underside very much in the same way one records the position of the real probe apex as it scans a real sample.
  • the record of all those reflected probe apex positions renders a sample estimate, very much in the same way that the record of the real probe apex positions, as it scans a real sample, will render the image.
  • This technique can be used to provide an initial probe estimate in steps 76 and 78 instead of using Villarubia's blind deconvolution method.
  • a Legendre based stereo profilometry can be used to deconvolve the sample image.
  • the Legendre Transformation-based method relies on two images only. It makes use of the fact that when the sample 22 is tilted relative to the probe 14 by a certain angle (as described above in FIG. 1 ), the Legendre transform of its shape changes in a predictable way relatively to the Legendre transform of the sample in the original orientation.
  • the Legendre transform of an AFM image is a linear combination of the Legendre transform of the probe shape and the Legendre transform of the sample shape
  • two images are obtained at steps 110 and 112 by scanning twice with different sample tilt angles ( ⁇ i and ⁇ T1 relative to the probe 14 ).
  • steps 114 - 138 solve a system of equations that renders the true sample denoted below as “sample” and probe shapes denoted below as “probe 1 ”, “probe 2 ” without a need for exact probe characterization.
  • Image 1 and image 2 refer to the measured images image 1 61 and image 2 62 as shown in FIG. 2 .
  • Probe 1 and probe 2 refer to the curves probe 1 14 ′ and 14 ′′ as shown in FIG. 2 and are derived mathematically in the equations below.
  • is the angle of rotation ( ⁇ T1 - ⁇ 1 ) between image 2 obtained in step 112 and in step 114 . It should be appreciated that some of these processing steps can occur in any order and that some of the steps can be combined mathematically.
  • the Legendre Transform for image 1 is calculated by calculating the slope m at each point in image 1 .
  • m r which is the slope of the rotated function
  • is the inclination angle of the line that is tangent to the surface at x′.
  • the Legendre Transform for the second image, image 2 62 is computed by calculating the slope m at each point in image 2 .
  • Equation 10 shows that given the angle of rotation ⁇ and the original Legendre Transform L(m), the rotated version L r (m) can be obtained. This is a necessary condition for the principle established by equation 2 to work.
  • a parametric function is defined to describe the probe.
  • the transformations are obtained for probe 1 in a first position and for probe 2 in a rotated position in steps 120 and 122 .
  • the Legendre transform for the probe in a first position is calculated.
  • a valve (L 0 /(cos( ⁇ ) ⁇ m 0 ⁇ sin( ⁇ )),m 0 ) is computed.
  • Step 122 produces the Legendre transform for the rotated probe 14 ′.
  • the deconvolution algorithm is derived as follows. For a certain sample with Legendre Transform S(m) and a probe with Legendre Transform P(m) an image can be generated.
  • the Legendre Transform of image 1 referred to as I 1 (m) should satisfy equation 1. That is:
  • the Legendre Transform for the sample, S(m), is defined as a function of the rotated slopes by a parametric polynomial.
  • the second image 2 referred to as I 2 (m r ) should be described by:
  • Equation 14 provides P trans ((m r ⁇ m′)/(1+m r ⁇ m′)) which is equivalent to P r (m r ) or P(m r ⁇ m′)/(1+m r ⁇ m′)/(cos( ⁇ ) ⁇ (m r ⁇ m′)/(1+m r ⁇ m′ sin( ⁇ ))).
  • step 132 the Legendre transform of the first image is converted to a function of m r from a function of m using m t and B t .
  • a least squares algorithm is used to obtain S(m r ) from which the sample measurement can be calculated.
  • the least squares solution is accomplished in Matlab (The Mathworks, Inc., Natick, Ma.). It should be appreciated that other commercial packages or software routines can provide the least squares solution.
  • both the sample shape and the probe shape are recovered.
  • a parametric function is provided to describe the probe by defining a parametric function to describe the probe geometry and then by establishing how the first (e.g. upright) and the rotated versions are geometrically related. It should be appreciated that the method described here can be further adapted to work with a non-symmetric probe by allowing different left and right-hand segment functions to describe the probe.
  • a parabolic function is chosen to describe the probe geometry. It should be appreciated that any parametric formula could be used.
  • I 1trans ( m r ⁇ m ′/1+ m r m ′) I 1r ( m r ).
  • step 124 the following parametric polynomial is used to describe S(m r ) as shown in equation 16.
  • This polynomial is used to reconstruct the sample shape. It should be appreciated that the number of terms can be varied and that other functions can be used to describe S(m r ) including sinusoidal, conical and exponential functions.
  • the Legendre Transform of the sample can be obtained using the parametric polynomial defined in step 124 , and converting the Image transforms to the same transform space in steps 126 - 130 .
  • m t is computed according to equation 17a.
  • the expression m t represents the derivative at the same point P when the curve is rotated ( ⁇ ) radians.
  • ⁇ t is computed for all the values m r contained in the domain of I 2 (m r ) ⁇ I 1r (m r ), according to equation 17b.
  • step 130 the Legendre transform of image 1 is converted from a function of m to a function of m r using m t and ⁇ t .
  • step 132 the subtraction of the converted first (e.g. upright) sample transform S trans from the rotated sample transform S(m r ) can be described by:
  • n is the power of the polynomial describing the Legendre transform of S(m r ), the Legendre transform of the sample evaluated at m r where m r are the values of the slopes of image 2 and m t1 n
  • step 134 the over constrained system of equation 2 represented by equation 20 with out the probe is solved with a least square method to determine the coefficients s i .
  • the sample shape is reconstructed by applying the reverse transform of equation 21 as follows:
  • the sample 22 shape can be recovered without any probe 14 characterization since the only information necessary is given by the two images and the rotation angle. Therefore, the sample 22 shape can be recovered by the Legendre Transform based method.
  • equation 2 the system established in equation 2 is solved by eliminating the probe shape and solving for the sample shape according to equations 15 and 20. This is done by assuming a general parametric shape for the sample, as the polynomial description shown in equation 16.
  • the reconstruction could be accomplished by eliminating the sample shape from the equations and then assuming a general parametric shape for the probe, instead of the sample.
  • a general parametric shape could be a sine wave of unknown frequency or amplitude, or a conic curve, or again, a polynomial. Instead of finding the shape of the sample first, one would find the probe shape. And then, by using equation 2, one can find the sample shape.
  • An example of this method is described in B. D. Aumond, K. Youcef-Toumi, “High Precision Stereo Profilometry based on Atomic Force Microscopy Technology”, Mechatronics 2000 Conference, Atlanta, Ga. (September 2000).
  • the least squares method offers a “best fit” to the points S(m r ). With a high order polynomial, the best fit curve will go through all points S(m r ) but in between those points, the curve will have some error representing S(m r ). With a low order polynomial, the curve may not exactly go through the points S(m r ) but will be a smoother curve which may provide more desirable results.
  • the techniques and apparatus of the present invention thus permit the successful profiling of high aspect ratio nanoscopic features using an atomic force microscope operated in the contact mode.
  • the difficulties of using an AFM with these types of structures stem from the fact that AFM technology was developed for mostly flat samples or features with small relative heights and/or smooth topography. While measurements of line width of micro trenches and micro steps in the semiconductor industry pushed the technology further, high aspect ratio structures present a new challenge.
  • the present invention addresses and overcomes these challenges, including new deconvolution techniques based on the use of multiple images of the sample without a continuing requirement for exact probe calibration.

Abstract

An apparatus and method are provided for processing the images obtained from an atomic force microscopy when profiling high aspect ratio features. A deconvolution technique for deconvolving the sample image includes the use of multiple images but does not require exact calibration of the scanning probe. In one embodiment, erosion and dilation techniques are used to obtain an undistorted image of the sample being measured. In another embodiment, Legendre transforms are used to obtain an undistorted image of the sample being measured. Also described is a technique for measuring the tip radius of the scanning probe.

Description

CROSS-REFERENCE TO RELATED APPLICATIONS
This application claims priority under 35 U.S.C. §119(e) to U.S. Provisional Patent Application No. 60/184,641 entitled “Image Deconvolution Techniques For Probe Scanning Apparatus” filed on Feb. 24, 2000, which is incorporated herein by reference in its entirety.
STATEMENTS REGARDING FEDERALLY SPONSORED RESEARCH
Not applicable.
FIELD OF THE INVENTION
This invention relates generally to profilometry and more particularly, to image deconvolution techniques for probe scanning devices.
BACKGROUND OF THE INVENTION
As is known in the art, a structure having one characteristic dimension (e.g. height) which is more pronounced than the others (e.g. width and length) is referred to as a high aspect ratio structure. Examples of these type of structures include probes of atomic force and scanning tunnelling microscopes, field emission probes, micro-indenters and Micro Electro-Mechanical systems (MEM's) structures. Other high aspect ratio structures are found in quantum magnetic media for high-density data storage, compact disk stampers, crystalline structures, blades and biological systems such as virus particles. Such high aspect ratio structures have typical heights on the order of tens of micrometers and tip radii in the range of tens of nanometers. Further, these structures may or may not be conductive.
Obtaining images of high aspect ratio structures poses certain challenges. In imaging such structures, it is sometimes desirable to image the side walls of the structure and to measure the radius of the tip of the structure in a non-destructive manner. Desired image resolutions can be on the order of 1 nm in the vertical direction (i.e. a direction which is normal to a measurement surface) and 10 nm in the lateral direction (i.e. a direction which is parallel to the measurement surface). Some profilometry methods are unable to provide such resolutions and thus such imaging criteria prevent the use of certain types of profilometry methods.
In one type of conventional stylus profilometer, a stylus having a sharp tip and coupled to a hinged arm is mechanically dragged along the sample surface. The deflection of the hinged stylus arm is measured and recorded as the surface profile. The use of a hinged stylus arm allows measurement of very rough surfaces, for example those with peak-to-peak heights greater than 1 mm. Probe-to-surface contact forces range from 10−3 N to 10−6 N. However, since the hinged stylus arm is partially supported by the stylus itself, physical rigidity limits the minimum stylus tip radius and hence the lateral resolution to about 0.1 mm.
In optical profilometry, many different optical phenomena (such as interference and internal reflection) can be utilized. The most popular technique is based on phase-measuring interferometry, in which a light beam reflecting off the sample surface is interfered with a phase-varied reference beam. The surface profile is deduced from the resulting fringe patterns. With a collimated light beam and a large photodetector array, the entire surface can be profiled simultaneously. This and other conventional optical profilometry methods are limited in lateral resolution by the minimum focussing spot size of about 0.5 μm (for visible light). In addition, measurement values are dependent upon the surface reflectivity of the material being profiled.
In the current state of the art, only scanning probe microscopes can meet a 10 nm lateral resolution requirement. In these microscopes, an atomically sharp (or nearly so) tip at a very close spacing to the sample surface is moved over the surface using a piezoactuator. One type of scanning probe microscope is the atomic force microscope (AFM), which measures the topography of a surface with a probe that has a very sharp tip. A probe assembly includes a cantilever beam from which the probe, or microstylus extends. The probe terminates at the probe tip having a typical tip radius of less than 0.1 μm. The probe typically has a length on the order of a couple of micrometers and the cantilever beam typically has a length between 100 μm and 200 μm.
In a contact mode atomic force microscope, the probe is moved relative to the surface of a sample and deflection of the cantilever is measured to provide a measure of the surface topography. More particularly, a laser beam is directed toward, and reflects off the back surface of the cantilever to impinge upon a sensor, such as a photodetector array. The electrical output signals of the photodetector array provide a topographical image of the sample surface and, further, provide feedback signals to a fine motion actuator, sometimes provided in the form of a piezoelectric actuator. In a constant force contact AFM, the fine motion actuator is responsive to the feedback signals for maintaining a substantially constant force between the probe tip and the sample, such as forces on the order of 10−8 N to 10−11 N.
Contact atomic force microscopy offers high lateral and vertical resolutions, such as less than 1 nm vertical resolution and less than 50 nm lateral resolution. Further, since the contact AFM relies on contact forces rather than on magnetic or electric surface effects, advantageously the contact AFM can be used to profile conductive and non-conductive samples. However, the maximum surface roughness that can be profiled is much less than that of conventional stylus profilometers which use a linear variable differential transducer (LVDT).
In the non-contact atomic force microscope, long range van der Waals forces are measured by vibrating the cantilever near its resonance frequency and detecting the change in the vibrational amplitude of a laser beam reflected off the cantilever due to a change in the force gradient caused by changes in the surface profile. The non-contact atomic force microscope offers non-invasive profiling. However, the technique has some disadvantages when compared to contact atomic force microscopy. First, van der Waals forces are hard-to-measure weak forces, rendering the microscope more susceptible to noise. Secondly, the probe tip must be maintained at a fixed height above the sample, typically on the order of a few nanometers, and the feedback control necessary to maintain this spacing must operate slowly to avoid crashing the probe tip on the sample. Thirdly, since the tip is always floating above the surface, the effective tip radius is increased and hence the achievable lateral resolution is decreased.
AFM was primarily developed for high-resolution 3-D imaging (profilometry) of atomically flat samples. In that case, the probe tip is scanned over the sample and only the apical region of the probe interacts with the profiled surface. Therefore, AFM images will closely reproduce the topography regardless of distortions in the probe away from the apex. Accordingly, when stylus instruments are used in profilometry, the implicit assumption is that only the very apex of the stylus touches the surface at all points.
However, when structures having relatively high aspect ratio features are imaged, the AFM and stylus images can be quite different from the real topography. That is also the case when the dimensions of the sample are comparable to those of the employed probe (AFM probe or stylus). The reason for this deviation is that areas of the probe other than the apex (for instance, the probe sides) interact with the sample as well. The image distortion caused by the interaction of the probe with the surface is typically referred to as image convolution. These two conditions for significant convolution distortion—reduced sample dimensions and high aspect ratio occur frequently.
In many engineering fields, the characteristic dimensions of the samples or the features of interest lie well within the micrometer and sub-micrometer ranges. These fields include but are not limited to nanotechnology, micro-electromechanical systems (MEMS), semiconductor devices and storage media, micro-sensors, and blade fabrication. The investigated features could be photo-resist trenches in silicon wafers, memory pillars in quantum magnetic media devices, roughness in smooth optical surfaces, the radius of curvature of field emission probes and parts of micro-machines. Thus, images of such structures can be distorted by convolution errors.
The level of convolution is greatly dependent on the relative size and shape of the employed probe with respect to the sample. Different probes can interact with the sample generating different distorted images. Convolution seriously reduces metrology accuracy. The deviations in the measurement of radius of curvature of high aspect ratio samples is proportional to the radius of the probe. Therefore, samples with dimensions similar to those of the probe will generate images with close to a 100% radius measurement error. Width measurements taken from images of photo-resist trenches are embedded with deviations proportional to the square of the height of the probe (for nearly parabolic probes). Probe geometry and size can also affect texture parameter (e.g. roughness) measurements due to convolution.
In order to achieve high metrology accuracy in micrometer and sub-micrometer measurements, convolution effects must be minimized or eliminated from stylus and AFM images. Techniques that correct such effects are known as deconvolution methods. Deconvolution is necessary in applications in which the sample feature dimensions are in the same range as the size of the probe tip, therefore resulting in strong image convolution. Prior art deconvolution techniques include scanning a “standard sample” (i.e. a sample having a known surface shape) with a probe to provide a mean curve for the probe shape, an inner curve for the probe shape, and an outer curve for the probe shape, thereby characterising the probe shape. One problem with this approach, however, is that errors in the initial probe calibration propagate to all other measurements. Another problem with this approach is that the shape of the probe change and to problems over time because of physical changes in probe due to probe wear.
It is, therefore, desirable to provide a deconvolution technique that does not rely on probe characterization. It is also desirable to provide a deconvolution technique which is not degraded by probe wear.
SUMMARY OF THE INVENTION
In view of the above problems and limitations of existing AFM probe calibration and deconvolution techniques, the existence of convolution errors when measuring high aspect ratio features and when sample feature dimensions are in the same range as the size of the probe tip, and in accordance with the present invention, it has been recognized that multiple images can be used to reduce convolution errors in the measured image of a sample without the need for probe characterization. It would therefore be desirable to provide an apparatus and method to deconvolve a sample image using multiple measured images without having to accurately characterize the measuring probe or to recalibrate the probe as it wears.
In accordance with the present invention, an apparatus utilizing contact atomic force microscopy (AFM) includes a probe, a controller operative to move the probe into a first vantage point relative to the sample to produce a first image, and to move the probe into a second vantage point relative to the sample to produce a second image. The apparatus further includes a deconvolution processor which deconvolves the first image and the second image to produce an image of the sample. With such an arrangement, a sample can be measured by obtaining and processing multiple images without utilizing probe characterization or correcting for probe wear. The apparatus thus utilizes multiple images at different vantage points to remove the convolution errors caused by the probe and sample shape. Furthermore, since the apparatus does not require a conventional probe characterization process, changes in the shape of the probe due to wear from use or from other causes do not effect the accuracy of the convolution.
In accordance with a further aspect of the invention, a deconvolution method includes moving a probe in a first scanning pattern to generate a first image of the sample, changing said probe's vantage point, moving the probe in a second scanning pattern to generate a second image of the feature, and deconvolving the image of the sample using the first image and the second image. With this particular technique, a process for deconvolving an image without utilizing probe characterization is provided. The process is also unaffected by changes in probe shape due to wear or other causes. This technique is thus advantageous when using a microscope to measure relatively high aspect ratio samples. Furthermore, since the microscope images the sample at different angles, problems which arise due to regions where there is no contact between the sample and probe tip (e.g. so-called shadow zones) are reduced.
In a first embodiment, the deconvolution technique utilizes an iterative process. After generating multiple images of a sample from different vantage points, the process begins by generating two or more estimates of the probe shape using a first image. The estimates can be generated using a blind deconvolution method. Next using an erosion technique, each of the estimates of the probe shapes is used to obtain a corresponding estimate of the sample shape. Thus if two estimates of the probe shape are generated, then two estimates of the sample shape are obtained. At least two estimates of the probe shape are then combined to provide a new probe estimate. Similarly, at least two estimates of the sample shape are combined to provide a new sample shape estimate. The process of combining sample estimates (including newly computed sample estimates) to generate more sample estimates can be repeated any desired number of times. Similarly, the new sample shape estimates can be utilized to generate new probe shape estimates. The process of generating new probe shape estimates can also be repeated any desired number of times. In one embodiment, the above iterative process can be repeated until the changes in the newly computed sample shape and probe shape estimates are below a predetermined threshold when compared with the estimates from previous iteration. With such an arrangement an accurate estimate to the sample shape is obtained without the need for probe characterization. If additional images used, the fidelity of the estimated sample is increased.
In another embodiment, the deconvolution technique utilizes Legendre transforms in the processing of images to produce sample and probe measurements. After generating two images of a sample from different vantage points, the Legendre transform of the first and second image are obtained. These transforms are used in a system of equations that relate the transform of the sample with a first transform of the probe and the transform of the first image, the transform of the sample with a transform of the probe from the changed vantage point and the transform of the second image and relating the first transform of the probe, the transform of the probe from the changed vantage point and the angle between the first and second vantage point. After obtaining a Legendre transform of the first image and obtaining a Legendre transform of the second image, a parametric function is provided to describe the probe geometry and then by establishing how the upright and the rotated versions are geometrically related.
The parametric function provides a functional dependency between the probe transform in the first orientation and the probe transform from the second vantage point which allows a Legendre transform of the probe to be eliminated from the system of equations. Solving the system of equations by a least squares algorithm provides a Legendre transform of the sample. Finally, the sample is recovered from the Legendre transform of the sample. The probe geometry can also be recovered from the solution for the sample shape. With such an arrangement a partial but exact reconstruction of the sample shape is obtained without the need for probe characterization.
BRIEF DESCRIPTION OF THE DRAWINGS
The foregoing features of this invention, as well as the invention itself, may be more fully understood from the following description of the drawings in which:
FIG. 1 shows a schematic diagram of an atomic force microscope;
FIGS. 2A-2C demonstrate image deconvolution according to the present invention;
FIG. 3 is a flow diagram illustrating the erosion and dilation deconvolution technique according to an aspect of the invention; and
FIG. 4 is a flow diagram illustrating the Legendre transform deconvolution technique according to an aspect of the invention.
DETAILED DESCRIPTION OF THE INVENTION
Referring now to FIG. 1, a contact mode atomic force microscope (AFM) 10 includes a probe assembly 12 having a probe 14 mounted to the free end of a cantilever beam 16. The probe assembly 12 has a probe tip 18 adapted to contact a sample surface 20 of the sample 22 to be profiled. In this example, a sample 22 is disposed on a sample holder 24. It should be understood, however, that in other applications, sample 22 may be held by any fixture or surface. Any means of locating the sample 22 proximate the probe 14 may be used. A typical radius of the probe tip 18 is between 5 nm and 100 nm, a typical length of the probe 14 is on the order of 4.0 μm, and a typical length of the cantilever beam 16 is between 100 and 200 μm. It should be appreciated, however, that the present invention can be used in AFM's or other devices having probes and probe tips of any dimension.
In use, the probe assembly 12 and sample 22 are moved relative to one another, causing the probe tip 18 to move over the sample surface 20. Movement of the probe tip 18 over features including high aspect ratio structures 26 on the sample surface 20 causes the cantilever beam 16 to deflect. The deflection of the cantilever beam 16 is measured to provide a measure of the surface topography. More particularly, a laser beam 30 generated by a laser source 32 is directed towards, and reflects off a surface 16 a of the cantilever beam 16. The reflected laser signal 31 impinges upon a sensor, or detector 40, which may, for example, be provided as a photodetector array. In response to the reflected laser signal 31, the detector 40 generates an output signal 42 which is provided to a signal processor 44.
Signal processor 44 processes the output signal 42 to provide processed signals 46, portions of which correspond to feed back control signals 48 and portions of which correspond to the topographical sample image signals 50. The control signals 48 are fed to a piezoelectric actuator 49 and the image signals 50 are fed to a display 52 on which the sample topographical image and other images can be displayed. A topographical sample view of the image 54 can be viewed on a display 52. The signal processor 44 can include a deconvolution processor 45. The display 52 is optional and the signal processor 44 operations can be automated.
Various types of signal processing may be performed, such as deconvolution, as described below in conjunction with FIG. 3 below. Further, the signal processor 44 may process the resulting image itself to perform certain sample measurements, such as measuring the radius of the tip of a feature of the sample. It should be appreciated the functions of the signal processor 44 and the deconvolution processor 45 can be provided by a separate processor, or by computer software instructions or groups of instructions executed on another processor which is part of the AFM. Such processing may be performed by a single processing apparatus which may, for example, be provided as part of microscope or measurement system such as that to be described below in conjunction with methods described in FIGS. 3 and 4. Alternatively, the processing steps of the signal processor 44 and the deconvolution processor 45 can be performed by functionally equivalent circuits such as a digital signal processor circuit or an application specific integrated circuit (ASIC).
Conventionally, AFM cantilever beams 16 are inclined with respect to the sample plane in order to reflect the laser beam 30 into the detector 40. However, if the probe 14 is inclined relative to sample holder 24, the probe tip 18 will touch one side of the high aspect ratio structure 26, but possibly not the other side. For this reason, preferably, the sample holder 24 is inclined to the same degree as the cantilever beam 16 so that the probe 14 and sample holder 24 are initially normal to one another. Ideally, the probe tip 18 should be normal to the sample surface 20 being measured.
In accordance with the present invention, the method for image deconvolution uses multiple imaging or stereo microscope imaging. The motivation for such an approach is drawn from the understanding that if the scanning is performed with the probe 14 always normal to the sample surface 20, convolution errors are minimized. Since this is not always practical given the available AFM configurations, one approximation to the ideal configuration is obtained by obtaining multiple images in which the angle between the probe tip 18 and the sample 22 is different. This may be accomplished, for example by utilizing different sample tilt angles. The images are then combined to form the deconvolved sample image.
For example, with the probe 14 and the sample holder 24 disposed in an initial position, the angle between an a central longitudinal axis 56 of the probe 14 and a plane 58 parallel with the plane of the sample holder 24 corresponds to an initial angle θi. The initial angle θi is measured with respect to a central longitudinal axis 56 of probe 14 and an axis formed by the intersection of plane 58 and a plane normal to the cantilever beam 16 and containing axis 56. With such relative probe and sample positions at a first vantage point, a first image, image1 61 is obtained by scanning the probe tip 18 over the sample surface 20. A second image, image2 62 is obtained by tilting sample holder 24 causing the sample 22′ (tilted sample 22) to be imaged by the probe 14 at relative angle 74 T1. Portions of the image1 61 and image2 62 can be presented on the display 52 as topographic images of the sample. The sample images image1 61 and image2 62 are processed in accordance with techniques to be described below in conjunction with FIGS. 2-4.
Now referring to FIGS. 2A-2C, images obtained from measuring the sample 22 (FIG. 1) in a single scan in the X direction are shown. The Z dimension represents heights and repetitive scans are in the X direction. Curve 64 represents an image of a portion of the sample 22 which would be generated in the absence of convolution errors. The portion of the sample 22 shown in FIG. 2A has a relatively high aspect ratio.
As illustrated in FIG. 2A, image1 61 is obtained by scanning the sample 22 with the probe 14 in a first position. In this particular example, the first position corresponds to the initial normal position as discussed above in conjunction with FIG. 1. Those of ordinary skill in the art should appreciate, however, that this need not be so. The first position can be any position as long as subsequent positions (probe vantage points) differ from it (i.e. as long as the relative angle between the probe and the sample holder are different). Curve probe1 14′ represents a theoretical image of the probe 14 (FIG. 1) if the probe 14 could be perfectly characterized. When the probe 14 physically contacts the sample at point 14 b′ the point that appears in image1 61 is point 14 a′. The convolution error is the difference between the true probe-sample contact point 14 b′ and the imaged point 14 a′. As the probe 14 (FIG. 1) is scanned over the sample 22 (FIG. 1), additional imaged points 14 e′-14 g′ are generated to form image1 61 with similar convolution errors. It should be appreciated that the number of points generated in the images depends upon the selected resolution of the microscope and the maximum resolution of the microscope scanner.
FIG. 2B shows a second image, image2 62, obtained by scanning the sample 22 (FIG. 1) with the probe 14 (FIG. 1) in a second position relative to the sample 22 (FIG. 1) which is different that that used to generate image1 61. In this particular example, image2 62 is generated when the sample holder 24 (FIG. 1) is tilted to an angle θT1 relative to a central longitudinal axis 56 of probe 14 and an axis 58′ formed by the tilted sample holder 24′ on which the tilted sample 22′ is disposed as shown in FIG. 1. It should be appreciated that rather than tilting the sample holder 24 (FIG. 1), the same result can be achieved by titling the probe 14 (FIG. 1). Alternatively, still the same result could be achieved by tilting both the probe 14 and the sample 22 by different angles to achieve the desired relative angle between the probe 14 and the sample 22. Image2 62 is generated as described above with respect to image1 61 and thus the image2 62 has similar convolution errors. As similarly described in FIG. 2A, curve probe2 14″ represents a theoretical image of the probe 14 (FIG. 1) at relative angle θT1 between the probe 14 and the tilted sample 22′, if the probe 14 could be perfectly characterized.
FIG. 2C shows a reconstructed image 66 using the deconvolution techniques described below. The degree to which reconstructed image 66 approximates the true sample topography 64 depends upon a variety of factors including but not limited to the number of images processed, the deconvolution method, the extension of areas of the sample that were actually touched by the probe during scanning and the accuracy of any geometrical estimates of probe shape, including parametric geometric functions. Two deconvolution techniques are described below in embodiments which can be used to recover the true sample shape using multiple images. A first embodiment uses a minimum envelope method and a second embodiment uses a Legendre transformation-based reconstruction method.
In the first embodiment, a minimum envelope approach based on the concept of blind deconvolution is used to recover a deconvolved image. In the minimum envelope approach, two or more images are taken with different vantage points (obtained by changing the relative angle between the probe and sample as described above). The true sample profile for example curve 64 (FIG. 2A) is always bounded by the image sample, image1 61 (FIG. 2A). A blind deconvolution approach is used to find initial estimates for probe and sample topographies based on the measured first image, image1 61, and the second image, image2 62. The probe topography is bounded by the intersection of all estimates. The sample topography is bounded by the intersections of all sample estimates obtained with the blind deconvolution approach.
Referring now to FIG. 3, the deconvolution process begins in step 70. In step 72 a first image of a sample, image1 61, is obtained with a first relative angle between the probe and the sample. In step 74, a second image of the sample, image2 62 is obtained with a different relative angle between the sample 22 and the probe 14 that is charged relative to the configuration used in step 72. Thus, the images, image1 61 and image2 62 are obtained from different vantage points. It should be appreciated that either the sample can be tilted, the probe can be tilted, or both the probe and sample can be tilted by different amounts to produce different angles between the probe and the sample.
In step 76, a first blind estimate of the probe EstP1 based on the first image 61 is obtained using Villarubia's blind deconvolution method based on the probe shape 120 as is known in the art. The details of the blind estimation technique are explained in J. S. Villarubia, “Morphological Estimation of Tip Geometry for Scanned Probe Microscopy,” Surf. Sci. 321, 287-300 (1994) and G. S. Pingali, R. Jain, “Restoration of Scanning Probe Microscope Images,” IEEE Workshop on Applications of Computer Vision 282 (1992).
In step 78, a second blind estimate of the probe EstP2 is obtained using Villarubia's blind deconvolution method based on image2 62. These blind estimates EstP1 and EstP2 are blunter than the real probe 14 and vary from the real shape, surface sample 20, because the images do not offer enough high aspect ratio features for a good estimation. The lack of high aspect ratio features leads to poor shape estimates.
Step 80 uses erosion and probe estimate EstP1 to find new sample estimate, EstS1 (estimate 1) from the first sample image and step 82 uses erosion and probe estimate EstP2 to find new sample estimate EstS2 (estimate 2) from the second sample image. The details of the erosion process are provided in J. S. Villarubia, “Morphological Estimation of Tip Geometry for Scanned Probe Microscopy,” Surf. Sci. 321, 287-300 (1994) and G. S. Pingali, R. Jain, “Restoration of Scanning Probe Microscope Images,” IEEE Workshop on Applications of Computer Vision 282 (1992).
In step 84, the two probe estimates are combined which results in the generation of a sharper estimate. If the estimate obtained from the upright sample image is EstP1 and the estimate obtained from the tilted sample image is EstP2 then the new estimate may be computed as EstP=EstP1∩EstP2. Since P⊂EstP1 and that P⊂EstP2 then it is guaranteed that P ⊂EstP1∩EstP2. Additionally, since (EstP1∩EstP2 532 EstP2 and (EstP1∩EstP2)⊂EstP1 the new estimate is sharper or more accurate than the previous estimates.
In step 86 the two estimates of the sample are combined to generate a better estimate of the sample surface topography. The new sample estimate is used to “sharpen” the probe estimate, which is then used to iteratively create better sample estimates.
Steps 88 to 102 iteratively repeat the process until no change in the estimation above a predetermined threshold is calculated.
In step 88, the new probe estimate from step 84 is used to obtain a new sample estimate 1 by erosion. In 90, the new probe estimate from step 84 is used to obtain a new sample estimate 2 by erosion.
In step 92, the two new sample estimates are combined to generate another more accurate new sample estimate. In order to do that, the sample estimate for the tilted case must be rotated by minus θ degrees rendering a rotated estimate denoted as EstS1*. The intersection is then found and will be the new estimate EstS0new. These estimates are obtained by the erosion method. Since the estimate for the tilted case EstS1 was obtained for the tilted sample, in order to combine the estimations, EstS1 must be rotated back to render an estimate at 0 degrees (EstS1*).
As shown in steps 94, 96 the new sample estimates are used to obtain two new probe estimates.
In step 100, the two new Probe estimates EStP2new and EStP1new are combined, generating a new sharper probe estimate EstPnew.
Decision block 102 implements a loop in which the process of generating new probe and sample estimates in steps 88 to 102 are repeated. With the new probe estimate EstPnew, new sample estimates EStS0new′, EstSnew, can be generated by erosion. This is also described in B. D. Aumond and K. Youcef-Toumi, “Experimental High Precision Profilometry of High Aspect Ratio Samples,” Proceedings of the 1998 IEEE International Conference on Systems, Man, and Cybernetics, San Diego, Calif. (November 1998), Y. Yeo, B. D. Aumond, and K. Youcef-Toumi, “Precision Atomic Force Microscope Imaging,” Proceedings of the 2000 IEEE International Conference on Signal Processing, Beijing, China (August 2000), and B. D. Aumond, K. Youcef-Toumi, “High Precision Stereo Profilometry based on Atomic Force Microscopy Technology,” Mechatronics 2000 Conference, Atlanta, Ga. (September 2000).
Thus, the new sample estimates are combined to generate a better estimate of the real surfaces for the upright and tilted measurements. This new surface estimate is in turn used to sharpen the probe estimates using both the first and second images (e.g. the upright sample image and tilted sample image). The sharpened probes estimates are combined generating a new probe estimate EstPnew′.
In decision block 102, a test is performed to determine whether the further new sample estimate or the still further new probe estimate changed from the last iteration. If the change is lower than a predetermined threshold, then the deconvolution process is complete as shown in step 104. Otherwise the erosion process is repeated beginning at step 88. The sample combination steps 88, 90 and 92 and probe sharpening steps 94, 96, 100 will eventually render better estimates of probe geometry and sample topography. Thus, when compared to single image blind deconvolution, the results are improved.
It should be appreciated, the concepts and process of sample combination and probe sharpening are completely equivalent in two-dimensions (2-D) or three-dimensions (3-D).
In an alternate embodiment of the minimum envelope method, an underscanning technique (scanning the underside of the image with the probe reflected about the origin) as described in Y. Yeo, B. D. Aumond, and K. Youcef-Toumi, “Precision Atomic Force Microscope Imaging,” Proceedings of the 2000 IEEE International Conference on Signal Processing, Beijing, China (August 2000).
Using the underscanning technique, the reflected probe geometry (the geometry of the probe is reflected horizontally and then vertically) scans of the underside of the image instead of the real sample. The reflected probe to simulate the scanning of the underside of the image very much in the same way that a probe scans a real sample. Record the reflected probe apex positions as it scans the image underside, very much in the same way one records the position of the real probe apex as it scans a real sample. The record of all those reflected probe apex positions renders a sample estimate, very much in the same way that the record of the real probe apex positions, as it scans a real sample, will render the image. This technique can be used to provide an initial probe estimate in steps 76 and 78 instead of using Villarubia's blind deconvolution method.
In a second embodiment, a Legendre based stereo profilometry can be used to deconvolve the sample image. The Legendre Transformation-based method relies on two images only. It makes use of the fact that when the sample 22 is tilted relative to the probe 14 by a certain angle (as described above in FIG. 1), the Legendre transform of its shape changes in a predictable way relatively to the Legendre transform of the sample in the original orientation. The Legendre transform of an AFM image is a linear combination of the Legendre transform of the probe shape and the Legendre transform of the sample shape
L{sample}=L{probe1}+L{image1}  (equation 1)
Now referring to FIG. 4, two images are obtained at steps 110 and 112 by scanning twice with different sample tilt angles (θi and θT1 relative to the probe 14). In general overview, by using the constraint that relates the Legendre transforms of the probe in the two different orientations as expressed below in equation 2, and mathematically derived from equations 3-10, steps 114-138 solve a system of equations that renders the true sample denoted below as “sample” and probe shapes denoted below as “probe 1”, “probe2” without a need for exact probe characterization. Image1 and image2 refer to the measured images image1 61 and image2 62 as shown in FIG. 2. Probe1 and probe2 refer to the curves probe1 14′ and 14″ as shown in FIG. 2 and are derived mathematically in the equations below.
L{sample}=L{probe1}+L{image1}  {equation 2}
L{sample}=L{probe2}+L{image2}
L{probe1}=f(L{probe2};θ)θ
Where θ is the angle of rotation (θT11) between image2 obtained in step 112 and in step 114. It should be appreciated that some of these processing steps can occur in any order and that some of the steps can be combined mathematically.
At step 114, the Legendre Transform for image 1 is calculated by calculating the slope m at each point in image 1. The Legendre Transform is L=f(x)−m·x for a function y=f(x).
If the function is rotated by θ, the Legendre Transform undergoes a trigonometrical change. The new points (x′,y′) of the rotated function (x,y)can be obtained as follows:
x′=x·cos(θ)−y·sin(θ)  (equation 3)
y′=x·sin(θ)+y·cos(θ)
The new Legendre Transform of the rotated function (denoted Lr) can be written as:
L r =y′−m r ·x′  (equation 4)
in which mr corresponds to the slope of the rotated function.
Now substituting equation system 3 into equation 5 one can obtain:
Lr =x sin(θ)+y cos(θ)−m r(x cos(θ)−y sin(θ))  (equation 5)
Now, mr, which is the slope of the rotated function, can be obtained as follows:
m=tan(α)  (equation 6)
m′=tan(θ)
m r=tan(α+θ)=(m·cos(θ)+sin(θ))/(cos(θ)−m·sin(θ))=(m+m′)/(1−m·m′)
Where α is the inclination angle of the line that is tangent to the surface at x′.
Therefore, the rotated Legendre Transform is given by:
L r=(x sin(θ)+y cos(θ))−((m+m′)/(1−m·m′))(x cos(θ)−y sin(θ))  (equation 7a)
L r(mr)=(x sin(θ)+y cos(θ))−m r·(x cos(θ)−y sin(θ))  (equation 7b)
At step 116, the Legendre Transform for the second image, image2 62, is computed by calculating the slope m at each point in image2.
The relationship between the Legendre Transforms L and Lr is derived as follows:
L r(m r)=y′−m r x′ and m r=(m+m′)/1−m·m′
L r=(x·sin(θ)+y·cos(θ))−((m+m′)/(1−m·m′))(x·cos(θ)−y·sin(θ))
L r=(x·sin(θ)+y·cos(θ))−((m+tan(θ))/(1−m·tan(θ)))(x·cos(θ)−y·sin(θ))
L r=(x·sin(θ)+y·cos(θ))−((m·cos(θ)+sin(θ))/(cos(θ)−m·sin(θ))) (x·cos(θ)−y·sin(θ))  (equation 8)
L r=(1/(cos(θ)−m·sin(θ)))(y·cos2(θ)−m·y·cos(θ)sin(θ)+x·cos(θ)sin(θ)−(m·x·sin2(θ)−(m·x cos2(θ)−m·y·cos(θ)sin(θ)+x·cos(θ)sin(θ)−y·sin2(θ)))  (equation 9)
Simplifying equation 9 leads to equation 10. Equation 10 shows that given the angle of rotation θ and the original Legendre Transform L(m), the rotated version Lr(m) can be obtained. This is a necessary condition for the principle established by equation 2 to work.
L r =y−m·x·/(cos(θ)−m·sin(θ))=L(m)/(cos(θ)−m·sin(θ)
m r=(m+m′)/(1−m·m′)
m=(m r −m′)/(1+m r ·m′)  (equation 10)
At step 118 a parametric function is defined to describe the probe. The transformations are obtained for probe1 in a first position and for probe2 in a rotated position in steps 120 and 122.
At step 120 the Legendre transform for the probe in a first position (e.g. an upright position) is calculated. At step 122 for each pair (L0, M0) a valve (L0/(cos(θ)−m0·sin(θ)),m0) is computed. Next, the value mr is substituted for mo where mr corresponds to mr=(m0+m′)/(1−m0·m′) is substituted for m0 to obtain the pair (L0/cos(θ)−m0 sin(θ),mr). The new pair (Lr,mr) is equal to (L0(m0)/(cos(θ)−m0·sin(θ)),mr). Step 122 produces the Legendre transform for the rotated probe 14′.
Additionally if a curve has a Legendre Transform with value L at derivative m (this value of derivative happens at some point P along the curve), then the rotated curve will have a Legendre Transform of value Ltrans at the same point P which will now have a derivative mr=(m+tan(θ))/(1−m·tan(θ)).
L trans(m)=L(m)/(cos(θ)−m·sin(θ))
L r(m r)=L trans(m r −m′(/(1+m r ·m′)  (equation 11)
The deconvolution algorithm is derived as follows. For a certain sample with Legendre Transform S(m) and a probe with Legendre Transform P(m) an image can be generated. The Legendre Transform of image1 referred to as I1(m) should satisfy equation 1. That is:
S(m)=I 1(m)+P(m)  (equation 12)
At step 124, the Legendre Transform for the sample, S(m), is defined as a function of the rotated slopes by a parametric polynomial. Now, if the probe is rotated by θ, the second image 2 referred to as I2(mr) should be described by:
S(m r)=I 2(m r)+P r(m r)  (equation 13)
Where Pr(mr) is the Legendre Transform of the rotated probe, probe2. Equation 11 shows that:
P r(m r)=P trans(m r −m′)/(1+m r ·m′)  (equation 14)
Applying the transformation stated in equation 14 provides Ptrans((mr−m′)/(1+mr·m′)) which is equivalent to Pr(mr) or P(mr−m′)/(1+mr·m′)/(cos(θ)−(mr−m′)/(1+mr·m′ sin(θ))).
Therefore, the system described by equations 12, 13 and 14 provides the Legendre Transform of the sample S(m) with no need for probe characterization.
Because the Legendre Transform of f(x) is a function of the slope m
LT[f(x)]Δ L(m)
∴f(x)=m(x)x+L(m(x))
The above leads to the true, deconvolved sample profile. In step 126, the equation S(m)=I1(m)+P(m) is used to obtain S(m)/(cos(θ)−m·sin(θ))=I1(m)/(cos(θ)−m·sin(θ))+P(m)/(cos(θ)−m·sin(θ)) or Strans(m)=I1trans(m)+Ptrans(m) or Strans(mr−m′)/(1+mr·m′))=I1trans(mr·m′)/(1+mr−m′)+Ptrans(mr−m′)/1+mrm′). Note that only the transform of the first image converted and then the probe is eliminated by subtraction in step 126.
From equation 14, Pr(mr)=Ptrans(mr−m′/1+mr·m′). This provides two equations:
S trans(m r −m′)/(1+m r ·m′))=I 1trans(m r −m′)/(1+m r ·m′))+P trans(m r −m′)/(1+m r −m′))
and
S(m r)−I 2(m r)=P r(m r).
Finally in step 126, probe transform Pr(mr) can be eliminated by simple subtraction resulting in: S ( m r ) - S trans m r - m 1 + m r · m = I 2 ( mr ) - I 1 trans m r - m 1 + m r · m = I 2 ( m r ) - I 1 r ( m r ) . Equation  (15)
Figure US06661004-20031209-M00001
At steps 128 and 130, intermediate calculations are performed to facilitate the elimination of the probe transform from equation 2 and allow the set of equations to be solved.
As step 132, the Legendre transform of the first image is converted to a function of mr from a function of m using mt and Bt.
At step 136, a least squares algorithm is used to obtain S(mr) from which the sample measurement can be calculated. Preferably the least squares solution is accomplished in Matlab (The Mathworks, Inc., Natick, Ma.). It should be appreciated that other commercial packages or software routines can provide the least squares solution.
At steps 138 and 140, both the sample shape and the probe shape are recovered.
The following example used to further describe the steps shown in FIG. 4. Returning to step 118, a parametric function is provided to describe the probe by defining a parametric function to describe the probe geometry and then by establishing how the first (e.g. upright) and the rotated versions are geometrically related. It should be appreciated that the method described here can be further adapted to work with a non-symmetric probe by allowing different left and right-hand segment functions to describe the probe.
At step 118, a parabolic function is chosen to describe the probe geometry. It should be appreciated that any parametric formula could be used. At step 118, the probe is described by parabolic function p(x)=½Rpx2, and at step 120 its Legendre transform is calculated by P(m)=−½Rp·m2. Where m is the slope at each point.
The Legendre Transform of the rotated probe is given by P ( m ) cos ( θ ) - m sin ( θ ) m = m r - m 1 + m r · m
Figure US06661004-20031209-M00002
That is:
P r(m r)=(−½R p((m r −m′)/((1+m r ·m′)2))/(cos(θ)−(m r −m′/1+m r ·m′)sin(θ))=
−½(Rp(m r −m′)2)/((1+m r ·m′)2 cos(θ)−(m r −m′)(1+m r ·m′)sin(θ)).
Pr(mr) is obtained as a function of P(m) and is identical to Pr=yr−mr·xr where yr(xr) is the function that describes the tilted probe with mr=dyr/dxr.
At step 122, the Legendre Transform, Pr(mr), for the tilted probe is calculated. Applying the transformation described in equation 11 to I1(m) will render
I1trans(m r −m′/1+m r m′)=I 1r(m r).
At step 124, the following parametric polynomial is used to describe S(mr) as shown in equation 16.
S(m r)=s n m r n +s n−1 m r n−1 + . . . +s 2 m r 2 +s 1 m r +s 0.  (equation 16)
This polynomial is used to reconstruct the sample shape. It should be appreciated that the number of terms can be varied and that other functions can be used to describe S(mr) including sinusoidal, conical and exponential functions.
The Legendre Transform of the sample can be obtained using the parametric polynomial defined in step 124, and converting the Image transforms to the same transform space in steps 126-130. At step 126, for all the values mr contained in the domain of I2(mr)−I1r(mr), mt is computed according to equation 17a.
m t(m r)=(m r −m′)/(1+m r ·m′).  (equation 17a)
Given a tilted curve with derivative mr at some point P, the expression mt represents the derivative at the same point P when the curve is rotated (−θ) radians.
In step 128, βt is computed for all the values mr contained in the domain of I2(mr)−I1r(mr), according to equation 17b.
βt(m r)=cos(θ)−m t sin(θ)=cos(θ)−(m r −m′/1+m r ·m′)sin(θ)  (equation 17b)
In step 130, the Legendre transform of image1 is converted from a function of m to a function of mr using mt and βt.
At step 132, subtracting I1r(mr) from I2(mr) results in S(mr)−Strans((mr−m′)/(1+mrm′)). Given that S(mr)−Strans((mr−m′)/(1+mr·m′))=I2 (mr)−I1r(mr) eliminates the probe from the system of equation 2.
In step 132, the subtraction of the converted first (e.g. upright) sample transform Strans from the rotated sample transform S(mr) can be described by:
S(m r)−S trans(m r −m′/1+m r ·m′)=s n m r n +s n−1 m r n−1 + . . . +s 2 m r 2 +s 1 m r
 +s 0−(s n m t n +s n−1 m t n−1 + . . . +s 2 m t 2 +s 1 m t +s 0)/βt  (equation 18)
or
s n(m r n −m t nt)+s n−1(m r n−1 −m t n−1t)+ . . . +s 2(m r 2 −m t 2t)+
s 1(m r −m tt)+s 0(1−1/βt)=I 2(m r)−I 1r(m r)  (equation 19)
The system above can then be re-written as follows: [ m r1 n - m t1 n / β t1 m r1 n - 1 - m t1 n - 1 / β t1 m r1 2 - m t1 2 / β t1 m r1 - m t1 / β t1 1 - 1 / β t1 m r2 n - m t2 n / β t2 m r2 n - m t2 n / β t2 m r2 n - m t2 n / β t2 m r2 n - m t2 n / β t2 1 - 1 / β t1 m rk n - m tk n / β tk m rk n - m tk n / β tk m rk n - m tk n / β tk m rk n - m tk n / β tk 1 - 1 / β t1 ] · { s n s 2 s 1 s 0 } = { I 2 ( m r1 ) - I 1 r ( m r1 ) I 2 ( m r2 ) - I 1 r ( m r2 ) I 2 ( m rk ) - I 1 r ( m rk ) } equation (20)
Figure US06661004-20031209-M00003
Where m r1 n
Figure US06661004-20031209-M00004
is the slope of point r1 . . . rk in image two, and n is the power of the polynomial describing the Legendre transform of S(mr), the Legendre transform of the sample evaluated at mr where mr are the values of the slopes of image2 and m t1 n
Figure US06661004-20031209-M00005
is a function of m r1 n
Figure US06661004-20031209-M00006
At step 134, the over constrained system of equation 2 represented by equation 20 with out the probe is solved with a least square method to determine the coefficients si.
At step 136, the sample shape is reconstructed by applying the reverse transform of equation 21 as follows:
S(m r)=y−m r ·x
x=−dS(m r)/dm r =−s n·(nm r n−1 −s n−1·(n−1) ·m r n−2 − . . . −s 2·(2)·m r −s 1
y=S(m r)+mr ·x=S(m r)−(dS(m r)/dm rm r  (equation 21)
The sample 22 shape can be recovered without any probe 14 characterization since the only information necessary is given by the two images and the rotation angle. Therefore, the sample 22 shape can be recovered by the Legendre Transform based method.
With this method, the system established in equation 2 is solved by eliminating the probe shape and solving for the sample shape according to equations 15 and 20. This is done by assuming a general parametric shape for the sample, as the polynomial description shown in equation 16.
It should be appreciated that the reconstruction could be accomplished by eliminating the sample shape from the equations and then assuming a general parametric shape for the probe, instead of the sample. A general parametric shape could be a sine wave of unknown frequency or amplitude, or a conic curve, or again, a polynomial. Instead of finding the shape of the sample first, one would find the probe shape. And then, by using equation 2, one can find the sample shape. An example of this method is described in B. D. Aumond, K. Youcef-Toumi, “High Precision Stereo Profilometry based on Atomic Force Microscopy Technology”, Mechatronics 2000 Conference, Atlanta, Ga. (September 2000). The parametric function used to describe the probe could also be chosen to be a conic shape, for example: y ( m ) = k ( 1 - 1 - α ( m 2 k 2 + α m 2 ) )
Figure US06661004-20031209-M00007
Using the method as described in FIG. 4 only part of the sample can be recovered. The reason is that the probe recovery can only happen at the probe point that touched the surface during the two different orientation scans. The sample recovery can only happen at the sample surface 20, with points that were touched by that set of probe points. In order to enlarge the surface recovery area, additional images taken at different angles can be used to obtain a better reconstruction for most surfaces. Even though the reconstruction is partial, the reconstruction is exact at all points S(mr).
The least squares method offers a “best fit” to the points S(mr). With a high order polynomial, the best fit curve will go through all points S(mr) but in between those points, the curve will have some error representing S(mr). With a low order polynomial, the curve may not exactly go through the points S(mr) but will be a smoother curve which may provide more desirable results.
If for example a polynomial of order two (a parabola) is used and S(mr) has a step like shape, the resulting fit will probably not be optimal. But if a polynomial of order 20 is used, it will better fit S(mr). If nothing is known about S(mr), different order polynomials can be tested to find the best fit. The selection of the function to represent S(mr) will have a corresponding effect on the reconstructed shape.
Additionally at step 138, given that S(mr) and I2(mr) are known, the probe Pr(mr) can be reconstructed by solving the linear equation S(mr)=I2(mr)+Pr(mr). Therefore, the technique can be used for probe characterization as well.
It is also important to note that the image recovery is not equivalent to selecting the minimum area below the two superimposed image curves (e.g. curves image1 61 and image2 62 in FIG. 2B). Selecting the minimum area below the two image curves offers only an estimate of the surface. The Legendre Transform based method is exact.
The techniques and apparatus of the present invention thus permit the successful profiling of high aspect ratio nanoscopic features using an atomic force microscope operated in the contact mode. The difficulties of using an AFM with these types of structures stem from the fact that AFM technology was developed for mostly flat samples or features with small relative heights and/or smooth topography. While measurements of line width of micro trenches and micro steps in the semiconductor industry pushed the technology further, high aspect ratio structures present a new challenge. The present invention addresses and overcomes these challenges, including new deconvolution techniques based on the use of multiple images of the sample without a continuing requirement for exact probe calibration.
All publications and references cited herein are expressly incorporated herein by reference in their entirety.
Having described the preferred embodiments of the invention, it will now become apparent to one of ordinary skill in the art that other embodiments incorporating their concepts may be used. It is felt therefore that these embodiments should not be limited to disclosed embodiments but rather should be limited only by the spirit and scope of the appended claims.

Claims (22)

What is claimed is:
1. A method for imaging a sample, comprising:
moving a probe in a first scanning pattern to generate a first image of the sample;
changing a vantage point of said probe;
moving said probe in a second scanning pattern to generate a second image of the sample; and
reconstructing the sample by deconvolution using said first image and said second image.
2. The method for imaging a sample according to claim 1, wherein the changing a vantage point further comprises tilting said probe and said sample relative to each other.
3. The method for imaging a sample according to claim 1, wherein reconstructing the sample by deconvolution, comprises the unordered steps of:
(a) estimating a first probe shape based on said first image to obtain a first probe estimate;
(b) estimating a second probe shape based on said image to obtain a second probe estimate;
(c) obtaining a first sample estimate by erosion using said first probe estimate;
(d) obtaining a second sample estimate by erosion using said second probe estimate;
(e) combining the first probe shape and the second probe shape to obtain a new probe estimate;
(f) combining said first sample shape and said second sample shape to obtain a new sample estimate;
(g) obtaining a new sample first estimate by erosion using the new probe estimate;
(h) obtaining a new sample second estimate by erosion using the new probe estimate;
(i) combining new sample first estimate and the new sample second estimate to obtain a further new sample estimate;
(j) using further new sample estimate, to obtain a further new probe first estimate;
(k) using further new sample estimate, to obtain a further new probe second estimate;
(l) combining further new probe first and second estimates to obtain a still further new probe estimate; and
(m) repeating steps until there is less than a predetermined threshold change in the further new sample estimate and the further new probe estimate from the previous iteration.
4. The method for imaging a sample according to claim 3, further comprising the step of generating additional images to increase image fidelity.
5. The method for imaging a sample according to claim 3, wherein estimating a first probe shape based on said first image to obtain a first probe estimate uses a blind deconvolution method; and
wherein estimating a second probe shape based on said image to obtain a second probe estimate uses a blind deconvolution method.
6. The method for imaging a sample according to claim 3, wherein estimating a first probe shape based on said first image to obtain a first probe estimate uses an underscanning method; and
wherein estimating a second probe shape based on said image to obtain a second probe estimate uses an underscanning method.
7. The method for imaging a sample according to claim 1, wherein the step of reconstructing the sample by deconvolution, comprises the steps of:
obtaining a Legendre transform of said first image;
obtaining a Legendre transform of said second image;
providing a parametric function to describe the probe;
providing a parametric function to describe the sample;
deriving the sample shape.
8. The method for imaging a sample according to claim 7, wherein the step of providing a parametric function to describe the probe further comprises modeling the probe to provide a probe model.
9. The method for imaging a sample according to claim 7, wherein the parametric function to describe the probe is derived from at least one of:
a conic function;
a sinusoidal function;
an exponential function; and
a low order polynomial function.
10. The method for imaging a sample according to claim 7, wherein the parametric function to describe the sample is derived from at least one of:
a conic function;
a sinusoidal function;
an exponential function; and
a low order polynomial function.
11. The method for imaging a sample according to claim 7, wherein the step of deriving the sample shape further comprises the steps of:
obtaining a Legendre transform of the probe;
obtaining a Legendre transform of the probe from said changed vantage point;
obtaining a Legendre transform of the sample from said changed vantage point; and
obtaining a Legendre transform of the sample.
12. The method for imaging a sample according to claim 11, wherein the step of obtaining a Legendre transform of the sample further comprises the steps of:
relating the Legendre transform of the sample from said changed vantage point to the Legendre transform of said probe, said Legendre transform of said first image, said Legendre transform of said second image and said Legendre transform of the probe from said changed vantage point; and
eliminating said Legendre transform of the probe and said Legendre transform of the probe from said changed vantage point from the relationship.
13. The method for imaging a sample according to claim 12, further comprising the step of using a least squares method to obtain a Legendre transform of the sample.
14. The method for imaging a sample according to claim 11, wherein the step of deriving the sample shape further comprises the step of using a reverse transform to reconstruct the sample shape.
15. The method for imaging a sample according to claim 11, further comprising the step of deriving the probe shape.
16. The method for imaging a sample according to claim 8, wherein the step of obtaining a Legendre transform of the probe further comprises calculating the Legendre transform of said probe model.
17. The method for imaging a sample according to claim 11, wherein the step of obtaining a Legendre transform of the probe further comprises the steps of:
relating the Legendre transform of the probe from said changed vantage point to said Legendre transform of the probe, said Legendre transform of said first image, said Legendre transform of said second image and said Legendre transform of the sample from said changed vantage point;
eliminating said Legendre transform of the sample and said Legendre transform of the sample from said changed vantage point from the relationship.
18. An apparatus for imaging a sample, comprising:
a probe;
a controller operative to move said probe into a first vantage point relative to the sample to produce a first image, and to move the probe into a second vantage point relative to the sample to produce a second image; and
a deconvolution processor which deconvolves said first image and said second image to reconstruct the sample shape.
19. The apparatus of claim 18 wherein the deconvolution processor uses blind deconvolution and erosion to reconstruct the sample shape.
20. The apparatus of claim 18 wherein the deconvolution processor uses relationships between a set Legendre transforms of the sample, the probe and said image from a first vantage point and a set of Legendre transforms of the sample, the probe and said image from a second vantage point to reconstruct the sample shape.
21. The apparatus of claim 18 wherein the deconvolution processor further reconstructs the probe shape.
22. The apparatus of claim 18 wherein the deconvolution processor uses underscanning and erosion to reconstruct the sample shape.
US09/789,992 2000-02-24 2001-02-21 Image deconvolution techniques for probe scanning apparatus Expired - Fee Related US6661004B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US09/789,992 US6661004B2 (en) 2000-02-24 2001-02-21 Image deconvolution techniques for probe scanning apparatus

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US18464100P 2000-02-24 2000-02-24
US09/789,992 US6661004B2 (en) 2000-02-24 2001-02-21 Image deconvolution techniques for probe scanning apparatus

Publications (2)

Publication Number Publication Date
US20010038072A1 US20010038072A1 (en) 2001-11-08
US6661004B2 true US6661004B2 (en) 2003-12-09

Family

ID=22677738

Family Applications (1)

Application Number Title Priority Date Filing Date
US09/789,992 Expired - Fee Related US6661004B2 (en) 2000-02-24 2001-02-21 Image deconvolution techniques for probe scanning apparatus

Country Status (2)

Country Link
US (1) US6661004B2 (en)
WO (1) WO2001063555A2 (en)

Cited By (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030218132A1 (en) * 2002-05-16 2003-11-27 Akihiko Honma Scanning probe microscope
US20050109925A1 (en) * 2003-11-26 2005-05-26 El Rifai Osamah M. Height calibration of scanning probe microscope actuators
US20060029253A1 (en) * 2004-07-30 2006-02-09 Euclid Discoveries Llc Apparatus and method for processing video data
US20060067585A1 (en) * 2004-09-21 2006-03-30 Euclid Discoveries Llc Apparatus and method for processing video data
US20060133681A1 (en) * 2004-11-17 2006-06-22 Euclid Discoveries Llc Apparatus and method for processing video data
US20060177140A1 (en) * 2005-01-28 2006-08-10 Euclid Discoveries Llc Apparatus and method for processing video data
US20060233448A1 (en) * 2005-03-31 2006-10-19 Euclid Discoveries Llc Apparatus and method for processing video data
WO2006055512A3 (en) * 2004-11-17 2007-03-15 Euclid Discoveries Llc Apparatus and method for processing video data
US20070176101A1 (en) * 2006-01-31 2007-08-02 Proksch Roger B Variable density scanning
US20070297645A1 (en) * 2004-07-30 2007-12-27 Pace Charles P Apparatus and method for processing video data
US20080161784A1 (en) * 2006-10-26 2008-07-03 Hogan Joseph M Method and system for remotely controlled MR-guided focused ultrasound ablation
US20080266575A1 (en) * 2007-04-27 2008-10-30 Angelo Gaitas Hybrid Contact Mode Scanning Cantilever System
US7472576B1 (en) 2004-11-17 2009-01-06 State Of Oregon Acting By And Through The State Board Of Higher Education On Behalf Of Portland State University Nanometrology device standards for scanning probe microscopes and processes for their fabrication and use
US20100031402A1 (en) * 2008-07-31 2010-02-04 Shigeru Wakiyama Probe aligning method for probe microscope and probe microscope operated by the same
US20100073458A1 (en) * 2007-01-23 2010-03-25 Pace Charles P Systems and methods for providing personal video services
US20110182352A1 (en) * 2005-03-31 2011-07-28 Pace Charles P Feature-Based Video Compression
US20120323513A1 (en) * 2011-06-20 2012-12-20 The University Of Sussex Apparatus and method for measuring charge density distribution
US8553782B2 (en) 2007-01-23 2013-10-08 Euclid Discoveries, Llc Object archival systems and methods
US8902971B2 (en) 2004-07-30 2014-12-02 Euclid Discoveries, Llc Video compression repository and model reuse
US8908766B2 (en) 2005-03-31 2014-12-09 Euclid Discoveries, Llc Computer method and apparatus for processing image data
US20150369838A1 (en) * 2012-12-12 2015-12-24 Universitat Basel Method and device for controlling a scanning probe microscope
US9366693B2 (en) 2006-01-31 2016-06-14 Oxford Instruments Plc Variable density scanning
WO2016118860A1 (en) * 2015-01-22 2016-07-28 The Board Of Trustees Of The Leland Stanford Junior University Methods and systems for determining proportions of distinct cell subsets
US9532069B2 (en) 2004-07-30 2016-12-27 Euclid Discoveries, Llc Video compression repository and model reuse
US9578345B2 (en) 2005-03-31 2017-02-21 Euclid Discoveries, Llc Model-based video encoding and decoding
US9621917B2 (en) 2014-03-10 2017-04-11 Euclid Discoveries, Llc Continuous block tracking for temporal prediction in video encoding
TWI580924B (en) * 2014-10-01 2017-05-01 耐諾股份有限公司 Deconvolution to reduce the effective spot size of a spectroscopic optical metrology device
US9743078B2 (en) 2004-07-30 2017-08-22 Euclid Discoveries, Llc Standards-compliant model-based video encoding and decoding
US10091507B2 (en) 2014-03-10 2018-10-02 Euclid Discoveries, Llc Perceptual optimization for model-based video encoding
US10097851B2 (en) 2014-03-10 2018-10-09 Euclid Discoveries, Llc Perceptual optimization for model-based video encoding

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8435373B2 (en) 2005-06-20 2013-05-07 Microcontinumm, Inc. Systems and methods for roll-to-roll patterning
US7833389B1 (en) * 2005-01-21 2010-11-16 Microcontinuum, Inc. Replication tools and related fabrication methods and apparatus
US9307648B2 (en) 2004-01-21 2016-04-05 Microcontinuum, Inc. Roll-to-roll patterning of transparent and metallic layers
US9039401B2 (en) 2006-02-27 2015-05-26 Microcontinuum, Inc. Formation of pattern replicating tools
EP1617293A1 (en) * 2004-07-14 2006-01-18 Universität Kassel A method of aligning a first article relative to a second article and an apparatus for aligning a first article relative to a second article
US8062495B2 (en) 2005-01-21 2011-11-22 Microcontinuum, Inc. Replication tools and related fabrication methods and apparatus
US7406860B2 (en) 2006-04-28 2008-08-05 Seagate Technology Llc Atomic force microscopy scanning and image processing
US7423264B2 (en) * 2006-09-08 2008-09-09 Kla-Tencor Technologies Corporation Atomic force microscope
US8310931B2 (en) * 2008-07-18 2012-11-13 International Business Machines Corporation Discovering network topology from routing information
US8296860B2 (en) * 2009-03-16 2012-10-23 Seagate Technology Llc Atomic force microscopy true shape measurement method
US8845912B2 (en) 2010-11-22 2014-09-30 Microcontinuum, Inc. Tools and methods for forming semi-transparent patterning masks
US9586817B2 (en) 2011-07-28 2017-03-07 Seagate Technology Llc Semi-auto scanning probe microscopy scanning
US9589797B2 (en) 2013-05-17 2017-03-07 Microcontinuum, Inc. Tools and methods for producing nanoantenna electronic devices
CN114236181B (en) * 2021-12-02 2023-10-20 中国电子科技集团公司第十三研究所 AFM probe measuring method, device, control equipment and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5591903A (en) * 1994-08-11 1997-01-07 Iowa State University Research Foundation Reconstructing the shape of an atomic force microscope probe
US5665905A (en) 1994-06-07 1997-09-09 International Business Machines Corporation Calibration standard for 2-D and 3-D profilometry in the sub-nanometer range and method of producing it
US5937111A (en) 1996-03-15 1999-08-10 Fuji Photo Film Co., Ltd. Image processing method and apparatus
US6489611B1 (en) * 1998-10-07 2002-12-03 Massachusetts Institute Of Technology Atomic force microscope for profiling high aspect ratio samples

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5665905A (en) 1994-06-07 1997-09-09 International Business Machines Corporation Calibration standard for 2-D and 3-D profilometry in the sub-nanometer range and method of producing it
US5591903A (en) * 1994-08-11 1997-01-07 Iowa State University Research Foundation Reconstructing the shape of an atomic force microscope probe
US5937111A (en) 1996-03-15 1999-08-10 Fuji Photo Film Co., Ltd. Image processing method and apparatus
US6489611B1 (en) * 1998-10-07 2002-12-03 Massachusetts Institute Of Technology Atomic force microscope for profiling high aspect ratio samples

Non-Patent Citations (43)

* Cited by examiner, † Cited by third party
Title
B.D. Aumond et al., "Experimental High Precision Profilometry of High Aspect Ratio Samples", Systems, Man, and Cybernetics, 1998, IEEE International Conference on San Diego, CA, pp. 4435-4440.
Bernardo D. Aumond, et al. "Experimental High Precision Profiliometry of High Aspect Ratio Samples" Mechanical Engineering Department Massachusetts Institute of Technology, Cambridge, MA 02139, 6 pages.
Bernardo D. Aumond, et al. "High Precision Stereo Profilometry Based on Atomic Force Microscopy" Mechanical Engineering Department Massachusetts Institute of Technology, Cambridge, MA 02139, 6 pages.
C. Durkan, et al. "40 nm resolution in reflection-mode SNOM with lambda = 685 nm" 3 pages.
C. Durkan, et al. "40 nm resolution in reflection-mode SNOM with λ = 685 nm" 3 pages.
C.C. Williams, et al. "Scanning capacitance microscopy on a 25 nm scale" Appl. Phys. Lett. 55 (2), Jul. 10, 1989, pp. 203-205.
C.C. Williams, et al. "Scanning thermal profiler" Appl. Phys. Lett. 49 (23), Dec. 8, 1986, pp. 1587-1589.
Copy of International Search Report of PCT/US01/05434.
D. Keller, "Reconstruction of STM and AFM Images Distorted by Finite-Size Tips", Surface Science, North-Holland Publishing Co., vol. 253, 1991, pp. 353-364.
David Keller, et al. "Reconstruction of STM and AFM images distorted by finite-size tips" Surface Science 253 (1991) 353-364 North Holland, pp. 353-364.
Derek G. Chietwynd, et al. "High Precision Surface Profilometry: From Stylus to Stm" Chapter 14, pp. 272-299.
Eiichi Furukawa, et al. "Development of a Flexure-Hinged Translation Mechanism Driven by Two Piezoelectric Stacks" JSME International Journal, Series C, vol. 38, No. 4 1995, pp. 743-748.
G. Binning, et al. "Atomic Force Microscope" vol. 56, No. 9, Mar. 3, 1986, pp. 930-933.
G. Binning, et al. "Surface Studies by Scanning Tunneling Microscopy" vol. 49, No. 1, Jul. 5, 1982 pp. 57-61.
G.L. Miller, et al. "A rocking beam electrostatic balance for the measurement of small forces" Rev. Sci. Instrum. 62 (3), Mar. 1991, pp. 705-709.
G.L. Miller, et al. "Resonant phase shift technique for the measurement of small changes in grounded capacitors" Rev. Sci. Instrum. 61 (4), Apr. 1990, pp. 1267-1272.
G.S. Pingali et al., "Estimation of Sample and Probe Tilts in Scanning Probe Microscopy", Proceeding of the Instrumentation and Measurement Technology Conference, 1993 New York, IEEE, pp. 327-332.
Gopal Sarma, et al. "Restoration of scanning probe microscope images" 0-8186-2840-May 1992 pp. 282-289.
GS Pingali et al., "Restoration of Scanning Probe Microscope Images", Applications of Computer Vision, Proceedings, 1992, pp. 282-289.
H. Bielefeldt, et al. "Reflection-scanning near-field optical microscopy and spectroscopy of opaque samples" Received Nov. 29, 1993 Accepted Mar. 22, 1994, pp. 1-3-108.
J.B. Pethica, et al. "Tip Surface Interactions in STM and AFM" Physica Scripta. vol. T19, 61 66. 1987, 6 pages.
J.F. Song, et al. "Stylus profiling at high resolution and low force" applied optics/vol. 30, No. 1/Jan. 1, 1991, pp. 42-50.
J.J. Saenz, et al. "Observation of magnetic forces by the atomic force microscope" J. Appl. Phys. 62 (10), Nov. 15, 1987, pp. 4293-4295.
J.M. Bennett, et al. "Stylus profiling instrument for measuring statistical properties of smooth optical surfaces" May 15, 1981/vol. 20, No. 10/Applied Optics, pp. 1785-1800.
J.S. Villarrubia "Scanned probe microscope tip characterization without calibrated tip characterizers" J. Vac. Sci. Technol. B. 14(2), Mar./Apr. 1996, 4 pages.
J.S. Villarrubia Morphological estimation of tip geometry for scanned probe microscopy; received May 25, 1994; accepted for publication Aug. 15, 1994, pp. 287-301.
J.S. Villarrubia, "Morphological Estimation of Tip Geometry for Scanned Probe Microscopy", Surface Science, vol. 321, 1994, pp. 287-300.
Joseph E. Griffith, et al "Dimensional metrology with scanning probe microscopes" J. Appl. Phys. 74 (9), Nov. 1, 1993, pp. R83-R109.
K.L. Westra "Effect of tip shape on surface roughness measurement from atomic force microscopy images of thin films" J. Vac. Sci. Technol. B 13(2), Mar./Apr. 1995, pp. 344-350.
Markiewicz et al. "Atomic Force MicroscopeTip Deconvolution using Calibration Arrays", Rev. Sci. Instrum. 66 (5), May 1995, pp. 3186-3190.* *
Masakazu Hayashi, et al. "Scanning Tunneling Microscope to Reveal Surface Roughness" Abstract, pp. 97-102.
Neil Goldfine, et al. "Near Surface Material Property Profiling for Determination of SCC Susceptibility" Fourth EPRI Balance-of-Plant Heat Exchanger NDE Symposiun, Jackson Hole, WY, Jun. 10-12, 1996, pp. 1-11.
P. Maivald, et al. "Using force modulation to image surface elasticities with the atomic force microscope" vol. 2, pp. 103-106 (1991).
P.M. Williams, et al. "Blind reconstruction of scanning probe image data" J. Vac. Sci. Technol. B 14(2), Mar./Apr. 1996, pp. 1557-1562.
Peter Markiewicz, et al. "Atomic force microscope tip deconvolution using calibration arrays" Rev. Sci. Instrum. 66 (5), May 1995, pp. 3186-3190.
Peter Markiewicz, et al. "Stimulation of atomic force microscope tip-sample/sample-tip" J. Vac. Sci. Technol. B 13(3), May/Jun. 1995, pp. 1115-1118.
S. Alexander, et al., "Atomic-resolution atomic-force microscope implemented using an optical lever" J. Appl. Phys. 65 (1), Jan. 1989, pp. 164-167.
Samuel Dongmo, et al. "Blind restoration method of scanning tunneling and atomic force microscopy images" J. Vac. Sci. Technol. B 14(2), Mar./Apr. 1996, pp. 1552-1556.
T. V. Vorburger "Methods for Characterizing Surface Topography" National Institute of Standards and Technology Gaithersburg, Maryland 29899, pp. 137-150.
W.P. Dyke, et al. "The Field Emitter: Fabrication, Electron Microscopy, and Electric Field Calculations" Journal of Applied Physics, vol. 24, May, 1953 pp. 570-575.
Y. Martin, et al. "Atomic force microscope-force mapping and profiling on a sub 100-A scale" J. Appl. Phys. 01 (10) May 15, 1987 pp. 4723-4729.
Y. Yeo, et al. "precision Atomic Force Microscope Imaging" Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge, MS 02139, 8 pages.
Yee Yeo "Image Processing for Precision Atomic force Microscopy" B.S., Purdue University (1998) 35 pages.

Cited By (61)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6904791B2 (en) * 2002-05-16 2005-06-14 Sii Nanotechnology Inc. Scanning probe microscope
US20050217354A1 (en) * 2002-05-16 2005-10-06 Akihiko Honma Scanning probe microscope
US7066015B2 (en) * 2002-05-16 2006-06-27 Sii Nanotechnology Inc. Scanning probe microscope
US20030218132A1 (en) * 2002-05-16 2003-11-27 Akihiko Honma Scanning probe microscope
US7041963B2 (en) 2003-11-26 2006-05-09 Massachusetts Institute Of Technology Height calibration of scanning probe microscope actuators
US20050109925A1 (en) * 2003-11-26 2005-05-26 El Rifai Osamah M. Height calibration of scanning probe microscope actuators
US20060029253A1 (en) * 2004-07-30 2006-02-09 Euclid Discoveries Llc Apparatus and method for processing video data
US9532069B2 (en) 2004-07-30 2016-12-27 Euclid Discoveries, Llc Video compression repository and model reuse
US8902971B2 (en) 2004-07-30 2014-12-02 Euclid Discoveries, Llc Video compression repository and model reuse
US7508990B2 (en) 2004-07-30 2009-03-24 Euclid Discoveries, Llc Apparatus and method for processing video data
US7158680B2 (en) 2004-07-30 2007-01-02 Euclid Discoveries, Llc Apparatus and method for processing video data
US9743078B2 (en) 2004-07-30 2017-08-22 Euclid Discoveries, Llc Standards-compliant model-based video encoding and decoding
US20070297645A1 (en) * 2004-07-30 2007-12-27 Pace Charles P Apparatus and method for processing video data
US7424157B2 (en) 2004-07-30 2008-09-09 Euclid Discoveries, Llc Apparatus and method for processing image data
US20060067585A1 (en) * 2004-09-21 2006-03-30 Euclid Discoveries Llc Apparatus and method for processing video data
US7426285B2 (en) 2004-09-21 2008-09-16 Euclid Discoveries, Llc Apparatus and method for processing video data
US20060133681A1 (en) * 2004-11-17 2006-06-22 Euclid Discoveries Llc Apparatus and method for processing video data
US7472576B1 (en) 2004-11-17 2009-01-06 State Of Oregon Acting By And Through The State Board Of Higher Education On Behalf Of Portland State University Nanometrology device standards for scanning probe microscopes and processes for their fabrication and use
WO2006055512A3 (en) * 2004-11-17 2007-03-15 Euclid Discoveries Llc Apparatus and method for processing video data
US7457435B2 (en) * 2004-11-17 2008-11-25 Euclid Discoveries, Llc Apparatus and method for processing video data
CN101103364B (en) * 2004-11-17 2010-05-12 欧几里得发现有限责任公司 Apparatus and method for processing video data
US7436981B2 (en) 2005-01-28 2008-10-14 Euclid Discoveries, Llc Apparatus and method for processing video data
US20060177140A1 (en) * 2005-01-28 2006-08-10 Euclid Discoveries Llc Apparatus and method for processing video data
US8964835B2 (en) 2005-03-31 2015-02-24 Euclid Discoveries, Llc Feature-based video compression
US20060233448A1 (en) * 2005-03-31 2006-10-19 Euclid Discoveries Llc Apparatus and method for processing video data
US8942283B2 (en) 2005-03-31 2015-01-27 Euclid Discoveries, Llc Feature-based hybrid video codec comparing compression efficiency of encodings
US8908766B2 (en) 2005-03-31 2014-12-09 Euclid Discoveries, Llc Computer method and apparatus for processing image data
US7457472B2 (en) 2005-03-31 2008-11-25 Euclid Discoveries, Llc Apparatus and method for processing video data
US9578345B2 (en) 2005-03-31 2017-02-21 Euclid Discoveries, Llc Model-based video encoding and decoding
US20110182352A1 (en) * 2005-03-31 2011-07-28 Pace Charles P Feature-Based Video Compression
US20110219479A1 (en) * 2006-01-31 2011-09-08 Asylum Research Corporation Variable Density Scanning
US20070176101A1 (en) * 2006-01-31 2007-08-02 Proksch Roger B Variable density scanning
WO2007089365A3 (en) * 2006-01-31 2008-04-03 Asylum Research Corp Variable density scanning
US9366693B2 (en) 2006-01-31 2016-06-14 Oxford Instruments Plc Variable density scanning
US8489356B2 (en) * 2006-01-31 2013-07-16 Oxford Instruments, PLC Variable density scanning
US7941286B2 (en) 2006-01-31 2011-05-10 Asylum Research Corporation Variable density scanning
US9106977B2 (en) 2006-06-08 2015-08-11 Euclid Discoveries, Llc Object archival systems and methods
US20080161784A1 (en) * 2006-10-26 2008-07-03 Hogan Joseph M Method and system for remotely controlled MR-guided focused ultrasound ablation
US8553782B2 (en) 2007-01-23 2013-10-08 Euclid Discoveries, Llc Object archival systems and methods
US8243118B2 (en) 2007-01-23 2012-08-14 Euclid Discoveries, Llc Systems and methods for providing personal video services
US8842154B2 (en) 2007-01-23 2014-09-23 Euclid Discoveries, Llc Systems and methods for providing personal video services
US20100073458A1 (en) * 2007-01-23 2010-03-25 Pace Charles P Systems and methods for providing personal video services
US20080266575A1 (en) * 2007-04-27 2008-10-30 Angelo Gaitas Hybrid Contact Mode Scanning Cantilever System
US8156568B2 (en) 2007-04-27 2012-04-10 Picocal, Inc. Hybrid contact mode scanning cantilever system
US20100031402A1 (en) * 2008-07-31 2010-02-04 Shigeru Wakiyama Probe aligning method for probe microscope and probe microscope operated by the same
US8495759B2 (en) * 2008-07-31 2013-07-23 Sii Nanotechnology Inc. Probe aligning method for probe microscope and probe microscope operated by the same
US20120323513A1 (en) * 2011-06-20 2012-12-20 The University Of Sussex Apparatus and method for measuring charge density distribution
US8798947B2 (en) * 2011-06-20 2014-08-05 The University of Sussex-Falmer Apparatus and method for measuring charge density distribution
US20150369838A1 (en) * 2012-12-12 2015-12-24 Universitat Basel Method and device for controlling a scanning probe microscope
US9500670B2 (en) * 2012-12-12 2016-11-22 Universitat Basel Method and device for controlling a scanning probe microscope
US10091507B2 (en) 2014-03-10 2018-10-02 Euclid Discoveries, Llc Perceptual optimization for model-based video encoding
US9621917B2 (en) 2014-03-10 2017-04-11 Euclid Discoveries, Llc Continuous block tracking for temporal prediction in video encoding
US10097851B2 (en) 2014-03-10 2018-10-09 Euclid Discoveries, Llc Perceptual optimization for model-based video encoding
TWI580924B (en) * 2014-10-01 2017-05-01 耐諾股份有限公司 Deconvolution to reduce the effective spot size of a spectroscopic optical metrology device
US9958327B2 (en) 2014-10-01 2018-05-01 Nanometrics Incorporated Deconvolution to reduce the effective spot size of a spectroscopic optical metrology device
US10274367B2 (en) 2014-10-01 2019-04-30 Nanometrics Incorporated Deconvolution to reduce the effective spot size of a spectroscopic optical metrology device
CN107430588A (en) * 2015-01-22 2017-12-01 斯坦福大学托管董事会 For the method and system for the ratio for determining different cell subsets
WO2016118860A1 (en) * 2015-01-22 2016-07-28 The Board Of Trustees Of The Leland Stanford Junior University Methods and systems for determining proportions of distinct cell subsets
US10167514B2 (en) 2015-01-22 2019-01-01 The Board Of Trustees Of The Leland Stanford Junior University Methods and systems for determining proportions of distinct cell subsets
CN107430588B (en) * 2015-01-22 2021-12-31 斯坦福大学托管董事会 Method and system for determining the proportion of different cell subsets
US11802314B2 (en) 2015-01-22 2023-10-31 The Board Of Trustees Of The Leland Stanford Junior University Methods and systems for determining proportions of distinct cell subsets

Also Published As

Publication number Publication date
US20010038072A1 (en) 2001-11-08
WO2001063555A2 (en) 2001-08-30
WO2001063555A3 (en) 2003-01-23

Similar Documents

Publication Publication Date Title
US6661004B2 (en) Image deconvolution techniques for probe scanning apparatus
US6489611B1 (en) Atomic force microscope for profiling high aspect ratio samples
US5898106A (en) Method and apparatus for obtaining improved vertical metrology measurements
US7684956B2 (en) Image reconstruction method
Hussain et al. Advances in the atomic force microscopy for critical dimension metrology
US9910064B2 (en) Force measurement with real-time baseline determination
EP1555676A2 (en) Method of operating a probe microscope
JPH0682248A (en) Method and apparatus for inspecting surface profile using interatomic-force scanning microscope
US8296860B2 (en) Atomic force microscopy true shape measurement method
Dixson et al. Spatial dimensions in atomic force microscopy: instruments, effects, and measurements
US6545273B1 (en) Use of multiple tips on AFM to deconvolve tip effects
US6427345B1 (en) Method and apparatus for a line based, two-dimensional characterization of a three-dimensional surface
US5591903A (en) Reconstructing the shape of an atomic force microscope probe
Wang et al. Large-area scanning probe nanolithography facilitated by automated alignment of probe arrays
JP2001194284A (en) Probe scanning method
Bugg Noncontact surface profiling using a novel capacitive technique: scanning capacitance microscopy
El Rifai et al. Imaging at the nano-scale
JP4497665B2 (en) Probe scanning control device, scanning probe microscope using the scanning control device, probe scanning control method, and measurement method using the scanning control method
EP3671224B1 (en) A method and apparatus for aligning a probe for scanning probe microscopy to the tip of a pointed sample
Kneedler et al. 3D metrology solution for the 65-nm node
JP3175342B2 (en) Surface shape measurement method
Aumond High precision profilometry
Aumond et al. High-precision metrology by means of a novel stereo imaging technique based on atomic force microscopy
Griffith et al. Metrology with scanning probe microscopes
Bhushan et al. Tip characterization for dimensional nanometrology

Legal Events

Date Code Title Description
AS Assignment

Owner name: MASSACHUSETTS INSTITUTE OF TECHNOLOGY, MASSACHUSET

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:AUMOND, BERNARDO D.;YOUCEF-TOUMI, KAMAL;REEL/FRAME:011676/0516

Effective date: 20010216

CC Certificate of correction
FPAY Fee payment

Year of fee payment: 4

FEPP Fee payment procedure

Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY

FPAY Fee payment

Year of fee payment: 8

REMI Maintenance fee reminder mailed
LAPS Lapse for failure to pay maintenance fees
STCH Information on status: patent discontinuation

Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362

FP Lapsed due to failure to pay maintenance fee

Effective date: 20151209