NL2003719A - Delta tcc for fast sensitivity model computation. - Google Patents

Delta tcc for fast sensitivity model computation. Download PDF

Info

Publication number
NL2003719A
NL2003719A NL2003719A NL2003719A NL2003719A NL 2003719 A NL2003719 A NL 2003719A NL 2003719 A NL2003719 A NL 2003719A NL 2003719 A NL2003719 A NL 2003719A NL 2003719 A NL2003719 A NL 2003719A
Authority
NL
Netherlands
Prior art keywords
image
tcc
difference
kernels
optical
Prior art date
Application number
NL2003719A
Other languages
English (en)
Inventor
Yu Cao
Jun Ye
Wenjin Shao
Hua Cao
Original Assignee
Brion Tech Inc
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
Priority claimed from PCT/US2009/049792 external-priority patent/WO2010005957A1/en
Application filed by Brion Tech Inc filed Critical Brion Tech Inc
Publication of NL2003719A publication Critical patent/NL2003719A/nl

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/39Circuit design at the physical level
    • G06F30/398Design verification or optimisation, e.g. using design rule check [DRC], layout versus schematics [LVS] or finite element methods [FEM]
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F1/00Originals for photomechanical production of textured or patterned surfaces, e.g., masks, photo-masks, reticles; Mask blanks or pellicles therefor; Containers specially adapted therefor; Preparation thereof
    • G03F1/36Masks having proximity correction features; Preparation thereof, e.g. optical proximity correction [OPC] design processes
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F1/00Originals for photomechanical production of textured or patterned surfaces, e.g., masks, photo-masks, reticles; Mask blanks or pellicles therefor; Containers specially adapted therefor; Preparation thereof
    • G03F1/68Preparation processes not covered by groups G03F1/20 - G03F1/50
    • G03F1/70Adapting basic layout or design of masks to lithographic process requirements, e.g., second iteration correction of mask patterns for imaging
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70Microphotolithographic exposure; Apparatus therefor
    • G03F7/70483Information management; Active and passive control; Testing; Wafer monitoring, e.g. pattern monitoring
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70Microphotolithographic exposure; Apparatus therefor
    • G03F7/70483Information management; Active and passive control; Testing; Wafer monitoring, e.g. pattern monitoring
    • G03F7/70491Information management, e.g. software; Active and passive control, e.g. details of controlling exposure processes or exposure tool monitoring processes
    • G03F7/705Modelling or simulating from physical phenomena up to complete wafer processes or whole workflow in wafer productions
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70Microphotolithographic exposure; Apparatus therefor
    • G03F7/70483Information management; Active and passive control; Testing; Wafer monitoring, e.g. pattern monitoring
    • G03F7/70491Information management, e.g. software; Active and passive control, e.g. details of controlling exposure processes or exposure tool monitoring processes
    • G03F7/70525Controlling normal operating mode, e.g. matching different apparatus, remote control or prediction of failure
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70Microphotolithographic exposure; Apparatus therefor
    • G03F7/708Construction of apparatus, e.g. environment aspects, hygiene aspects or materials
    • G03F7/70991Connection with other apparatus, e.g. multiple exposure stations, particular arrangement of exposure apparatus and pre-exposure and/or post-exposure apparatus; Shared apparatus, e.g. having shared radiation source, shared mask or workpiece stage, shared base-plate; Utilities, e.g. cable, pipe or wireless arrangements for data, power, fluids or vacuum
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/39Circuit design at the physical level

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Environmental & Geological Engineering (AREA)
  • Epidemiology (AREA)
  • Public Health (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
  • Exposure Of Semiconductors, Excluding Electron Or Ion Beam Exposure (AREA)

Description

DELTA TCC FOR FAST SENSTTTVTTY MODEL COMPUTATION
FIELD
[0001] The present invention generally relates to a method and computer product for performing simulations of imaging results associated with a lithographic process.
BACKGROUND
[0002] Lithographic apparatuses can be used, for example, in the manufacture of integrated circuits (ICs). In such a case, the mask may contain a circuit pattern corresponding to an individual layer of the IC, and this pattern can be imaged onto a target portion (e.g. comprising one or more dies) on a substrate (silicon wafer) that has been coated with a layer of radiation-sensitive material (resist). In general, a single wafer will contain a whole network of adjacent target portions that are successively irradiated via the projection system, one at a time. In one type of lithographic projection apparatus, each target portion is irradiated by exposing the entire mask pattern onto the target portion in one go; such an apparatus is commonly referred to as a wafer stepper. In an alternative apparatus, commonly referred to as a step-and-scan apparatus, each target portion is in'adiated by progressively scanning the mask pattern under the projection beam in a given reference direction (the "scanning" direction) while synchronously scanning the substrate table parallel or anti-parallel to this direction. Since, in general, the projection system will have a magnification factor M (generally < 1), the speed V at which the substrate table is scanned will be a factor M times that at which the mask table is scanned. More information with regard to lithographic devices as described herein can be gleaned, for example, from US 6,046,792, incorporated herein by reference.
[0003] The manufacturing process of an integrated circuit involves imaging multiple layers of circuit patterns onto a semiconductor wafer. To address the increasingly demanding needs of the semiconductor industry where the dimensions of the circuit elements are continually being reduced, technologies such as customized illumination schemes, phase shifting masks, and optical proximity correction (OPC) may be used. Numerical methods known as "computational lithography" are often required to simulate the lithography process to ensure the printed patterns meet the design target and to keep experimental costs under control. It is desirable that these numerical methods be highly accurate in order to predict the changes in aerial images (AIs) and critical dimensions (CDs) caused by small perturbations, also known as "sensitivities", to the optical settings.
SUMMARY
[0004] It is desirable to provide a method that accurately and rapidly simulates the change in the imaging of arbitrary mask patterns resulting from variations in the optical setting of an optical system [0005] In accordance with an aspect of the invention, there is provided a method for determining a difference between a reference image and a further image of a pattern, the method including determining a reference imaging function; determining parameters of a difference function representative of a difference between the reference imaging function and a further imaging function; calculating a difference between the reference image and the further image of the pattern based on the difference function and the determined parameters.
[0006] In one embodiment of the invention, the further imaging function is determined based upon the reference imaging function and the difference function.
[0007] In one embodiment of the invention, the further imaging function is determined by concatenating the reference imaging function and the difference function. In a further embodiment of the invention, calculating the difference in the reference image and the further image of the pattern includes determining the difference in image intensity of the reference image and the further image of the pattern. In addition, calculating the difference in the reference image and the further image of the pattern further may include determining the difference between a measure for critical dimensions in the reference image and the further image.
[0008] In another embodiment of the invention, at least a part of the terms of the difference function are represented in eigenfunctions and corresponding eigenvalues, and a simplified difference function is determined by disregarding eigenfunctions and eigenvalues based upon the eigenvalues. In one embodiment of the invention, the simplified difference function includes less than 300 eigenvalues or less than 100 eigenvalues.
[0009] In one embodiment of the invention, the reference image models an image of the pattern formed by a reference optical system under reference conditions, and the further image models an image of the pattern formed by a further optical system, or formed by the reference optical system under further conditions, or formed by a modified reference optical system.
[0010] In one embodiment of the invention, the difference between the reference image and the further image is decreased by changing the conditions of the further optical system, the reference optical system or the modified reference optical system. The reference optical system may be a projection system of a lithographic exposure apparatus.
[0011] In another aspect of the invention, there is provided a computer product having machine executable instructions, the instructions to be executed by a machine to perform a method for determining a change in an image of a mask pattern. In one embodiment, the method includes determining a reference imaging function; determining parameters of a difference function representative of a difference between the reference imaging function and a further imaging function; calculating a difference between the reference image and the further image of the pattern based on the difference function and the determined parameters.
BRIEF DESCRIPTION OF THE DRAWINGS
[0012] These and other aspects and features of the present invention will become apparent to those ordinarily skilled in the art upon review of the following description of specific embodiments of the invention in conjunction with the accompanying figures, wherein: [0013] Figure 1 is an exemplary block diagram illustrating a typical lithographic projection system.
[0014] Figure 2 is an exemplary block diagram illustrating the functional modules of a lithographic simulation model.
[0015] Figure 3 is a flowchart representing a method for determining a difference between a reference image and a further image of a pattern , in accordance with an embodiment of the invention; [0016] Figures 4a-b show two different intensity profiles of an illumination system.
[0017] Figure 4c shows the image sensitivity to the illumination change as calculated using a conventional method.
[0018] Figure 4d shows the image sensitivity to the illumination change as calculated using the method of Figure 3.
[0019] Figures 5a-b show two different intensity profiles of an illumination system.
[0020] Figure 5c shows the image sensitivity to the illumination change as calculated using a conventional method.
[0021] Figure 5d shows the image sensitivity to the illumination change as calculated using the method of Figure 3.
[0022] Figure 6 is a block diagram that illustrates a computer system which can assist in the implementation of the simulation method of the present invention.
[0023] Figure 7 schematically depicts a lithographic projection apparatus suitable for use with the method of the present invention.
DETAILED DESCRIPTION
[0024] The present invention will now be described in detail with reference to the drawings, which are provided as illustrative examples of the invention so as to enable those skilled in the art to practice the invention. Notably, the figures and examples below are not meant to limit the scope of the present invention to a single embodiment, but other embodiments are possible by way of interchange of some or all of the described or illustrated elements. Moreover, where certain elements of the present invention can be partially or fully implemented using known components, only those portions of such known components that are necessary for an understanding of the present invention will be described, and detailed descriptions of other portions of such known components will be omitted so as not to obscure the invention. Embodiments described as being implemented in software should not be limited thereto, but can include embodiments implemented in hardware, or combinations of software and hardware, and vice-versa, as will be apparent to those skilled in the art, unless otherwise specified herein. In the present specification, an embodiment showing a singular component should not be considered limiting; rather, the invention is intended to encompass other embodiments including a plurality of the same component, and vice-versa, unless explicitly stated otherwise herein. Moreover, applicants do not intend for any term in the specification or clauses to be ascribed an uncommon or special meaning unless explicitly set forth as such. Further, the present invention encompasses present and future known equivalents to the known components referred to herein by way of illustration.
[0025] Prior to discussing the present invention, a brief discussion regarding the overall simulation and imaging process to be calibrated is provided. Figure 1 illustrates an exemplary lithographic projection system 10. The major components are a light or radiation source 12, which may be a deep-ultraviolet excimer laser source, illumination optics which define the partial coherence (denoted as sigma) and which may include specific source shaping optics 14, 16a and 16b; a patterning device (e.g. a mask or reticle) 18; and projection optics 16c that produce an image of the patterning device pattern onto the wafer plane 22. An adjustable filter or aperture 20 at the pupil plane may restrict the range of beam angles that impinge on the wafer plane 22, where the largest possible angle defines the numerical aperture of the projection optics NA-sm(0max).
[0026] In a lithography simulation system, these major system components can be described by separate functional modules, for example, as illustrated in Figure 2. Referring to Figure 2, the functional modules include the design layout module 26, which defines the target design; the patterning device layout module (e.g. mask layout module) 28, which defines the patterning device (e.g. mask) to be utilized in imaging process; the patterning device model module (e.g. mask model module) 30, which defines the model of the patterning device layout (e.g. mask layout) to be utilized during the simulation process; the optical model module 32, which defines the performance of the optical components of lithography system; and the resist model module 34, which defines the performance of the resist being utilized in the given process. As is known, the result of the simulation process produces, for example, predicted contours and CDs in the result module 36.
[0027] More specifically, it is noted that the properties of the illumination and projection optics are captured in the optical model 32 that includes, but not limited to, NA-sigma (σ) settings as well as any particular illumination source shape, where σ(θΓ sigma) is outer radial extent of the illuminator. NA and σ are known parameters, as will be appreciated by one skilled in the art. The optical properties of the photo-resist layer coated on a substrate - i.e. refractive index, film thickness, propagation and polarization effects- may also be captured as part of the optical model 32. The patterning device model (e.g. mask model) 30 captures the design features of the patterning device (e.g. reticle or mask) and may also include a representation of detailed physical properties of the patterning device (e.g. reticle or mask), as described, for example, in U.S. Patent Appln. No. 10/530,402. Finally, the resist model 34 describes the effects of chemical processes which occur during resist exposure, PEB and development, in order to predict, for example, contours of resist features formed on the substrate wafer. The objective of the simulation is to accurately predict, for example, edge placements and CDs, which can then be compared against the target design. The target design, is generally defined as the pre-OPC patterning device layout, and will be provided in a standardized digital file format such as GDS11 or OASIS.
[0028] In general, the connection between the optical and the resist model is a simulated aerial image within the resist layer, which arises from the projection of radiation or light onto the substrate, refraction at the resist interface and multiple reflections in the resist film stack. The light intensity distribution (aerial image) is turned into a latent "resist image" by absorption of photons, which is further modified by diffusion processes and various loading effects. Efficient simulation methods that are fast enough for full-chip applications approximate the realistic 3-dimensional intensity distribution in the resist stack by a 2-dimensional aerial (and resist) image.
[0029] As should be therefore apparent from the above, the model formulation describes all of the known physics and chemistry of the overall process, and each of the model parameters corresponds to a distinct physical or chemical effect. The model formulation thus sets an upper bound on how well the model can be used to simulate the overall lithography process. So precise calibration of the model parameters is extremely important.
[0030] There are various ways that the model parameters can be defined. One efficient implementation of a lithography model is possible using the following formalism, where the image (here in scalar form, which may be extended to include polarization vector effects) is expressed as a Fourier sum over signal amplitudes in the pupil plane. According to the standard Hopkins theory, the aerial image intensity may be defined by:
Figure NL2003719AD00071
(Eq.l) [0031] where, A I(x) is the aerial image intensity at point x within the image plane (for notational simplicity, a two-dimensional coordinate represented by a single variable is utilized), k represents a point on the source plane, S (k) is the source intensity from point k, k ' and k" are points on the pupil plane, M is the Fourier transform of the mask image, P is the pupil function, and
Figure NL2003719AD00081
[0032] An important aspect of the foregoing derivation is the change of summation order (moving the sum over k inside) and indices (replacing k' with k+k' and replacing k" with k+ k"), which results in the separation of the Transmission Cross Coefficients (TCCs), defined by the term inside the square brackets in the third line of the equation. These coefficients arc independent of the patterning device pattern (e.g. mask pattern) and therefore can be precomputed using knowledge of the optical elements or configuration only (e.g., NA and σ or the detailed illuminator profile). It is further noted that although in the given example (Eq.l) is derived from a scalar imaging model, this formalism can also be extended to a vector imaging model, where the polarized light components arc summed separately.
[0033] Furthermore, the approximate aerial image A1 can be calculated by using only a limited number of dominant TCC kernels, which can be determined by diagonalizing the TCC matrix to obtain a diagonalized TCC matrix, also named hereinafter "diagTCC," and retaining (by truncation) the kemels corresponding to its largest eigenvalues, i.e., [0034]
Figure NL2003719AD00082
where A. (/= 1,...,A) denotes the N largest eigenvalues and φ;(·) denotes the corresponding eigenvector of the TCC matrix. It is noted that (Eq.2) is exact when all kemels are retained in the Eigen series expansion, i.e., when N is equal to the rank of the TCC matrix. However, in actual applications, it is typical to truncate the series by selecting a smaller N (i.e. n) to increase the speed of the computation process. The more kernels kept in the Eigen series, the better accuracy retained, at the cost of more computation. In general, the TCC matrix prior to the diagonalization operation can be referred to as a "raw" TCC matrix, and will be denoted hereinafter "rawTCC." [0035] Thus, (Eq.l) can be rewritten as: (Eq.3)
Figure NL2003719AD00091
where Φ;(λ:) = 'Yui>i{k")M(k")QVÇ){-jk"x) and |·| denotes the magnitude of a complex number. k" [0036] Using a sufficiently large number of TCC kemels and a suitable model calibration methodology allows for an accurate description of the optical projection process and provides an accurate representation of the aerial image. Additional information regarding the calculation of an aerial image can be gleaned from U.S. Patent No. 6,871,337 issued on March 22, 2005 titled "Illumination Optimization For Specific Mask Patterns", the contents of which is incorporated herein in its entirety by reference.
[0037] As will be appreciated, the aerial image AI only depends on the patterning device (e.g. mask) image and the transmission cross-coefficients (TCC). The TCCs capture all the optical characteristics of the optical system. If two models, which correspond to two optical systems, have the same TCCs, then the aerial images from the two models match perfectly for the same patterning device (e.g. mask). If the resist parts of the models are also the same, then the printing results also match perfectly.
[0038] In a manufacturing environment, it is often desirable to match the performance of an optical system (e.g. a second scanner or a second projection system in a lithographic exposure apparatus) to a reference optical system (e.g. a first reference scanner or a first optical system in a lithographic exposure apparatus). Further, it is also desirable to predict the change in the aerial image resulting from a change of the optical setting in a given optical system. This involves tweaking the optical settings on the to-be-tuned optical system, such as the numerical aperture (NA) of the optical system, sigma, stage tilt, and exposure dosage. These changes might be very small. For instance, a change of 0.01 in NA, or 10 mo in sigma inner/outer values is typical. To predict the changes in aerial images (AIs) and critical dimensions (CDs) caused by such small perturbations, also known as "sensitivities" to the optical settings, highly accurate optical models are needed.
[0039] As optical settings are represented by truncated TCC eigen-series, the truncation-induced errors play a significant role in the model sensitivities. Simulation experiments showed that the conventional setting of 64 TCC kernels is inadequate to capture the perturbation with satisfactory accuracy. In fact, for some cases, up to 1024 TCC kemels are used.
[0040] One conventional approach that may be used to predict the changes in aerial images and critical dimensions (1) from a reference optical system to another optical system or (2) within an optical system involves calculating two sets of aerial images using two TCC series in the following manner: [0041] In the conventional approach, as shown below in Equation 4, first, the reference optical system having a first optical setting (denoted as "state 1") is modeled using a TCC matrix (also denoted, hereinafter, as "rawTCC_l"). The first TCC matrix is then diagonalized and truncated to select the kemels with the greatest eigenvalues. The results of the diagonalization and truncation operations provide a limited number of TCC kernels denoted, hereinafter, as diagTCC l |iarge. As mentioned previously, typically, a large number of kernels is needed when calculating the aerial image, in order to accurately account for small perturbations. The aerial image (AI(diagTCC_l)|iarge) can be determined using Eq. 3. The critical dimension (CD(diagTCC_l)|iarge) of a given pattern will then be obtained based on the profile of the aerial image (AI(diagTCC_l)|iarge).
Figure NL2003719AD00101
[0042] Similar calculations are performed for a second optical setting, as shown in Equation 5. The second optical setting may correspond to the optical setting of a different optical system or the same optical system but in a perturbed state, for instance cause by drift or heating of optical elements. Thus, the optical system having a second optical setting (denoted as "state_2") is modeled using a TCC matrix (also denoted, hereinafter, as "rawTCC_2"). The second TCC matrix is then diagonalized and truncated to select the kemels with the greatest eigenvalues. The results of the diagonalization and truncation operations provide a limited number of TCC kernels denoted, hereinafter, as diagTCC_2|iarge. The aerial image (AI(diagTCC_2)|iarge) can be determined using Eq. 3. The critical dimension (CD(diagTCC_2)|iarge) of a given pattern will then be obtained based on the profile of the aerial image.
Figure NL2003719AD00111
[0043] The critical dimension sensitivity resulting from the changes in optical setting between the reference optical system and the second optical system or between different optical settings within a given optical system can be calculated in the following manner:
Figure NL2003719AD00112
[0044] While the above method for calculating changes in critical dimensions resulting from small perturbations provide good results when a large number of kernels are selected (typically up to 1024 TCC kemels), this is done at the expense of the computation speed and, thus, the simulation throughput. It is therefore desirable to provide an approach that can accurately and rapidly determine variations in an aerial image resulting from small perturbations in the optical settings.
[0045] In one embodiment of the invention, instead of using two optical models with extremely large sets of TCC kemels to describe the corresponding optical settings and compute the sensitivities, two small sets of TCC kernels are calculated. The first set of TCC kernels describes a reference imaging function, which may correspond to one of the optical system states (denoted as "base state" below), and the second set, referred to as "ATCC", describes a difference imaging function representative of the change in optical settings. A further imaging function is calculated by creating a new set of TCC kernels using the reference imaging function (first set of TCC kernels) and the difference imaging function (second set of TCC kernels). The new set of TCC kernels may be representative of a second optical system state (denoted as "perturbed state" below). The two TCC sets, i.c. the new set of TCC kernels representative of the perturbed state and the set of TCC kernels representative of the base state are then used for sensitivity calculation.
[0046] In one embodiment of the invention, the base state and the perturbed state are representative of different optical settings of different optical systems (e.g. a first or reference lithographic apparatus and a second or further lithographic apparatus). Alternatively, the base state and the perturbed state are representative of different optical settings of a same optical system. The optical setting may be described by various parameters including, for example, the numerical aperture of the optical system, the exposure dose, and the profile of illumination in an illuminator (characterized, e.g., by riuincr and o0uter)· It will be appreciated that additional parameters can be used to characterize the change of optical settings in other embodiments of the invention.
[0047] Referring now to Figure 3, this figure shows a method 300 for determining a difference between a reference image and a further image of a pattern . The method 300 begins at procedure 310 where a reference imaging function is determined. The reference imaging function may be representative of a base state ("state_l") of a first or reference optical system. The reference imaging function can be constructed by (a) using a first matrix of transmission cross-coefficients (TCCk,kO defined by rawTCC _\k,k" = ^S(k)P(k + k')P*(k + /c"), where k S(k) is the source intensity of the first optical system and P(k) is the pupil function of the first optical system, and by (b) diagonalizing that first matrix of transmission cross-coefficients (rawTCC lk,k') to obtain a diagonalized first set of transmission cross-coefficient kernels (diagTCC_lk,kO defined by diagTCC_\vk, = ιλ!φ!(^)φ. (k"), where λ; (i=l,...,N) are the N largest eigenvalues and Φ; are the corresponding eigenvectors. The diagonalized matrix diagTCC_lk,k' is truncated to retain a small number of TCC kernels with the highest eigenvalues, which would have the largest contribution to imaging. The truncation operation provides a matrix of TCC kernels denoted as "diagTCC_lk,k'|smaii," which corresponds to the reference imaging function. Typically, the number n of selected TCC kernels will depend upon the properties of the optical system, the nature of the optical change, geometry of the mask patterns to be imaged, and the required level of accuracy. In one embodiment, the number n of TCC kernels may be less than about 300. In another embodiment, the number n of TCC kernels may be less than about 100.
[0048] After calculating the reference imaging function function, the method proceeds to procedure 320, where parameters of a difference function representative of a difference between the reference imaging function and a further imaging function are determined. The difference function is representative of a change in the optical settings from the base state to a perturbed state. As noted previously, the base state and the perturbed state may be representative of different optical settings of either (a) different optical systems or (b) a same optical system. In one embodiment, the difference function is calculated by determining the difference (ArawTCC) between the first matrix of transmission cross-coefficients (rawTCC_lk,k') representative of the base state of, for example, a first optical system and a second matrix of transmission cross-coefficients (rawTCC_2k,k') representative of the perturbed state. The perturbed state may be associated with a second or further optical system or the same optical system. The second matrix of transmission cross-coefficients (rawTCC_2k,k') is defined by rawTCC _2k,k" = ^S(k)P(k + k')P* (k + k"). Thus, the difference function k
ArawTCC=rawTCC_2k,k'- rawTCC_lk,k'· In practice, the difference function ArawTCC corresponds to a matrix of transmission cross-coefficients.
[0049] The difference function (ArawTCC) between the first matrix of transmission cross-coefficients (rawTCC_lk,kO representative of the base state of, for example, a first optical system and the second matrix of transmission cross-coefficients (rawTCC_2k,k') representative of the perturbed state provides a matrix of transmission cross-coefficients, which is diagonalized and truncated. Similarly to the reference imaging function, a truncation operation is performed to retain a small number n! of TCC kernels, denoted as "AdiagTCCk,k'|smaii", having the highest eigenvalues. The number n' of TCC kernels will depend upon the properties of the optical system, the nature of the optical change, geometry of the mask patterns to be imaged, and the required level of accuracy. In one embodiment, the number n' of TCC kernels may be less than about 300. In another embodiment, the number n' of TCC kemels may be less than about 100. In one embodiment, the number n of kemels selected at procedure 310 may be the same as the number n' of kernels selected at procedure 320. In other embodiment, the number n can be smaller or greater than n'. The truncation operation provides a simplified difference function.
[0050] The further imaging function is based upon the reference imaging function and the difference function. In one embodiment, the further imaging function is determined by concatenating the reference imaging function and the difference function. For example, the further imaging function, denoted as diagTCC_lk,k'|small : AdiagTCCycIsmaii, is created using the reference imaging function diagTCC lk,k'Ismaii and the difference function AdiagTCC k,k'|smaii, where the sign denotes a combination of the two sets of diagonalized TCC kernels. For example, in one embodiment, the further imaging function is created by concatenating the set of TCC kemels diagTCC_lk,k'|smaii and the set of TCC kerneis ÄdiagTCCk,k'|Smaii· The further imaging function is representative of the perturbed state and is used to calculate the aerial images AT(diagTCCJk,k'Ismail : AdiagTCC|smaii) and critical dimensions CD(diagTCCJk,k'|smaii : AdiagTCC'k.k-Ismail) of the mask pattern.
[0051] Referring back to Figure 3, the method proceeds to procedure 330 where a difference between the reference image and the further image of the pattern based on the difference function and the determined parameters are calculated. In one embodiment, calculating the difference in the reference image and the further image of the pattern includes determining the difference in image intensity of the reference image and the further image of the pattern. Further, calculating the difference in the reference image and the further image of the pattern may include determining the difference between a measure for critical dimensions in the reference image and the further image.
[0052] For example, the aerial image sensitivity ΔΑΙ is calculated by determining the difference between the aerial image determined using the further imaging function and the aerial image determined using the reference imaging function. Thus, ΔΑΙ = AI(diagTCC_lk,k'|smaii : AdiagTCC|smaii) - AI((diagTCC_lk,k'|smaii)· The critical dimension sensitivity can be deduced from the aerial image sensitivity by calculating the difference between the critical dimension determined using the further imaging function CD(diagTCC_lk,k'Ismail : AdiagTCCk,k'|smaii) and the critical dimension determined using the reference imaging function CD(diagTCC_lk,k'|smaii)· [0053] In one embodiment, the difference between the reference image and the further image is decreased by changing the conditions of the further optical system, the reference optical system or the modified reference optical system.
[0054] In the approach defined in the embodiment of Figure 3, the perturbation is treated "separately", as opposed to the conventional approach where the perturbation is coupled with the baseline aerial images. The reasoning is that the difference function AdiagTCCk,k'|smaii corresponding to the perturbation will typically have a smaller rank, and can therefore be sufficiently accurately approximated by a smaller number of TCC kernels, as opposed to the complete TCC matrices that take many more TCC kernels to capture the sensitivity effectively.
[0055] There are numerous benefits of the approach of Figure 3. Compared to conventional approaches, keeping less kernels during TCC diagonalization saves computation time significantly. It will be appreciated that, using the method of Figure 3, a reduction in computation time by at least a factor of 3, and in some embodiments, a factor of 8, can be obtained. In addition, the memory consumption is significantly reduced. Indeed, using less TCC kernels requires less memory in use (for example, 16x for 64 TCC kernels vs. 1024 TCC kernels) at run time, which in turn allows more simulation tasks running in parallel on a same computational platform. This is particularly important for applications where multiple optical conditions need to be simulated simultaneously, and therefore multiple TCC's need to be held in memory. There is more flexibility in simulation job set-up. Using the approach of the embodiment of Figure 3, users may be able to concatenate a difference function, ÀTCC, to an existing TCC and calculate aerial images AIs or criticial dimensions, CDs, without the full knowledge of the numerical settings for the existing TCC. For example, one may perturb a base model using the present invention and generate a new model, which has the required sensitivity accuracy, and converges to the existing model when the perturbation approaches zero. These benefits will be further illustrated in Figures 4a-d and Figures 5a-d.
[0056] Figures 4a-d show critical dimension sensitivity calculations (Figures 4c and 4d) for two different intensity profiles of an illuminator (Figures 4a and 4b). The mask pattern used in the simulation is an array of 50nm lines arranged in various pitches ranging from lOOnm to lOOOnm. The intensity profiles in Figures 4a-b differ from each other in that the aoulcr slope of the intensity profile of Figure 4b is 3 times as sharp as the o0uter slope of the intensity profile of Figure 4a.
[0057] Figure 4c shows the variations of critical dimensions induced by a change in the optical settings calculated using (a) an infinite or very large number of diagonalized transmission cross coefficient kemels, or the undiagonalized rawTCC and (b) a limited number of TCC kemels for different gauges using a conventional approach. Each point in the curve represents a single gauge. Each gauge corresponds to a pattern. The baseline, i.e. ACD (delta critical dimension) =0, corresponds to the ideal result obtained using an infinite or very large number of transmission cross-coefficient kernels when calculating the image sensitivity that results from the change in illuminator intensity profile using a conventional approach. Thus, when ACD=0 for each gauge, the aerial image sensitivity is accurately determined and there is no difference between an aerial image calculated with an infinite or very large number of transmission cross coefficient kernels and an aerial image calculated with a limited number of transmission cross coefficient kernels. As shown in Figure 4c, it is desirable to use at least 1024 TCC kernels to obtain satisfactory results for each gauge using the conventional approach. If the number of TCC kemels is reduced (from 1024 to 512, to 256, to 128 and to 64), the difference between the results obtained using a very large number of transmission cross-coefficient kemels (baseline) and a limited number of transmission cross-coefficient kernels increases significantly.
[0058] Figure 4d differs from Figure 4c in that the difference in sensitivity is calculated in accordance with the method of Figure 3. In Figure 4c, the number of TCC kernels corresponding to the base state (i.e. the intensity profile of Figure 4a) is 64. Each curve corresponds to a different number of TCC kernels in the difference fonction, i.e. Ad iagTCCk,k'| small- This number varies from 32 to 256. As can be seen in Figure 4d, it is possible to obtain results that match those obtained with an infinite number of TCC kernels by using only 64 TCC kernels for the difference function and 64 TCC kernels for the base state. Thus, instead of using 1024 TCC kemels with the conventional approach, it is possible to obtain similar results with only 64 TCC kernels for the base state and 64 + 64 = 128 kernels for the perturbed state. Thus, it will be appreciated that the computation time is significantly reduced with the method of Figure 3 and that less memory is needed.
[0059] Calculations similar to those of Figures 4c-d were performed for other changes in intensity profiles shown in Figures 5a-b. Sensitivity calculations are shown in Figures 5c-d. The intensity profile of Figure 5a corresponds to an annular shape having a o0uter of 0.96 and a of 0.76. The intensity profile of Figure 5b corresponds to an annular shape having a o0uter of 0.997 and a amncr of 0.797. Calculations were performed for the same mask pattern as in Figures 4a-d.
[0060] As can be seen in Figure 5c, it is desirable to use at least 512 TCC kemels to accurately determine the change in aerial image with the conventional approach. By contrast, as shown in Figure 5d, the perturbed state modeled with only about 64 + 256 = 320 TCC kernels shows reasonable convergence for all of the gauges. Thus, the same changes as in Figure 5c can be determined with only about 320 TCC kernels, as opposed to 512 TCC kemels.
[0061] Figure 6 is a block diagram that illustrates a computer system 100 which can assist in implementing the lithographic simulation methods disclosed herein. Computer system 100 includes a bus 102 or other communication mechanism for communicating information, and a processor 104 coupled with bus 102 for processing information. Computer system 100 also includes a main memory 106, such as a random access memory (RAM) or other dynamic storage device, coupled to bus 102 for storing information and instructions to be executed by processor 104. Main memory 106 also may be used for storing temporary variables or other intermediate information during execution of instructions to be executed by processor 104. Computer system 100 further includes a read only memory (ROM) 108 or other static storage device coupled to bus 102 for storing static information and instructions for processor 104. A storage device 110, such as a magnetic disk or optical disk, is provided and coupled to bus 102 for storing information and instructions.
[0062] Computer system 100 may be coupled via bus 102 to a display 112, such as a cathode ray tube (CRT) or flat panel or touch panel display for displaying information to a computer user. An input device 114, including alphanumeric and other keys, is coupled to bus 102 for communicating information and command selections to processor 104. Another type of user input device is cursor control 116, such as a mouse, a trackball, or cursor direction keys for communicating direction information and command selections to processor 104 and for controlling cursor movement on display 112. This input device typically has two degrees of freedom in two axes, a first axis (e.g., x) and a second axis (e.g., y), that allows the device to specify positions in a plane. A touch panel (screen) display may also be used as an input device.
[0063] According to one embodiment of the invention, portions of the simulation process may be performed by computer system 100 in response to processor 104 executing one or more sequences of one or more instructions contained in main memory 106. Such instructions may be read into main memory 106 from another computer-readable medium, such as storage device 110. Execution of the sequences of instructions contained in main memory 106 causes processor 104 to perform the process steps described herein. One or more processors in a multi-processing arrangement may also be employed to execute the sequences of instructions contained in main memory 106. In alternative embodiments, hard-wired circuitry may be used in place of or in combination with software instructions to implement the invention. Thus, embodiments of the invention are not limited to any specific combination of hardware circuitry and software.
[0064] The term "computer-readable medium" as used herein refers to any medium that participates in providing instructions to processor 104 for execution. Such a medium may take many forms, including but not limited to, non-volatile media, volatile media, and transmission media. Non-volatile media include, for example, optical or magnetic disks, such as storage device 110. Volatile media include dynamic memory, such as main memory 106. Transmission media include coaxial cables, copper wire and fiber optics, including the wires that comprise bus 102. Transmission media can also take the form of acoustic or light waves, such as those generated during radio frequency (RF) and infrared (IR) data communications. Common forms of computer-readable media include, for example, a floppy disk, a flexible disk, hard disk, magnetic tape, any other magnetic medium, a CD-ROM, DVD, any other optical medium, punch cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, a carrier wave as described hereinafter, or any other medium from which a computer can read.
[0065] Various forms of computer readable media may be involved in carrying one or more sequences of one or more instructions to processor 104 for execution. For example, the instructions may initially be borne on a magnetic disk of a remote computer. The remote computer can load the instructions into its dynamic memory and send the instructions over a telephone line using a modem. A modem local to computer system 100 can receive the data on the telephone line and use an infrared transmitter to convert the data to an infrared signal. An infrared detector coupled to bus 102 can receive the data carried in the infrared signal and place the data on bus 102. Bus 102 carries the data to main memory 106, from which processor 104 retrieves and executes the instructions. The instructions received by main memory 106 may optionally be stored on storage device 110 either before or after execution by processor 104.
[0066] Computer system 100 also preferably includes a communication interface 118 coupled to bus 102. Communication interface 118 provides a two-way data communication coupling to a network link 120 that is connected to a local network 122. For example, communication interface 118 may be an integrated services digital network (ISDN) card or a modem to provide a data communication connection to a corresponding type of telephone line. As another example, communication interface 118 may be a local area network (LAN) card to provide a data communication connection to a compatible LAN. Wireless links may also be implemented. In any such implementation, communication interface 118 sends and receives electrical, electromagnetic or optical signals that carry digital data streams representing various types of information.
[0067] Network link 120 typically provides data communication through one or more networks to other data devices. For example, network link 120 may provide a connection through local network 122 to a host computer 124 or to data equipment operated by an Internet Service Provider (ISP) 126. ISP 126 in turn provides data communication services through the worldwide packet data communication network, now commonly referred to as the "Internet" 128. Local network 122 and Internet 128 both use electrical, electromagnetic or optical signals that carry digital data streams. The signals through the various networks and the signals on network link 120 and through communication interface 118, which carry the digital data to and from computer system 100, are exemplary forms of carrier waves transporting the information.
[0068] Computer system 100 can send messages and receive data, including program code, through the network(s), network link 120, and communication interface 118. In the Internet example, a server 130 might transmit a requested code for an application program through Internet 128, ISP 126, local network 122 and communication interface 118. In accordance with the invention, one such downloaded application provides for the TCC computation of the embodiment, for example. The received code may be executed by processor 104 as it is received, and/or stored in storage device 110, or other non-volatile storage for later execution. In this manner, computer system 100 may obtain application code in the form of a carrier wave.
[0069] Figure 7 schematically depicts an exemplary lithographic projection apparatus for which lithographic processing can be simulated utilizing the process of present invention. The apparatus comprises: - a radiation system Ex, IL, for supplying a projection beam PB of radiation. In this particular case, the radiation system also comprises a radiation source LA; - a first object table (patterning device support or mask table) MT provided with a mask holder for holding a patterning device (e.g. mask) MA (e.g., a reticle), and connected to first positioner for accurately positioning the mask with respect to item PL; - a second object table (substrate table) WT provided with a substrate holder for holding a substrate W (e.g., a resist-coated silicon wafer), and connected to second positioner for accurately positioning the substrate with respect to item PL; - a projection system ("lens") PL (e.g., a refractive, catoptric or catadioptric optical system) for imaging an irradiated portion of the mask MA onto a target portion C (e.g., comprising one or more dies) of the substrate W.
[0070] As depicted in Figure 7, the apparatus is of a reflective type (i.e., has a reflective mask). However, in general, it may also be of a transmissive type, for example (with a transmissive mask). Alternatively, the apparatus may employ another kind of patterning means as an alternative to the use of a mask; examples include a programmable mirror array or LCD matrix.
[0071] The source LA (e.g., a mercury lamp or excimer laser) produces a beam of radiation. This beam is fed into an illumination system (illuminator) IL, either directly or after having traversed a conditioner, such as a beam expander Ex, for example. The illuminator IL may include an adjuster AM for setting the outer and/or inner radial extent (commonly referred to as σ-outer and σ-inner, respectively) of the intensity distribution in the beam. In addition, it will generally include various other components, such as an integrator IN and a condenser CO. In this way, the beam PB impinging on the patterning device (e.g. mask) MA has a desired uniformity and intensity distribution in its cross-section.
[0072] It should be noted with regard to Figure 7 that the source LA may be within the housing of the lithographic projection apparatus (as is often the case when the source LA is a mercury lamp, for example), but that it may also be remote from the lithographic projection apparatus, the radiation beam that it produces being led into the apparatus (e.g., with the aid of suitable directing mirrors); this latter scenario is often the case when the source LA is an excimer laser (e.g., based on KrF, ArF or F2 lasing). The current invention encompasses at least both of these scenarios.
[0073] The beam PB subsequently intercepts the patterning device (e.g. mask) MA, which is held on a patterning device support (e.g. mask table) MT. After being reflected by the patterning device (e.g. mask) MA, the beam PB passes through the lens PF, which focuses the beam PB onto a target portion C of the substrate W. With the aid of the second positioner (and interferometric measuring means IF), the substrate table WT can be moved accurately, e.g. so as to position different target portions C in the path of the beam PB. Similarly, the first positioner can be used to accurately position the patterning device (e.g. mask) MA with respect to the path of the beam PB, e.g., after mechanical retrieval of the patterning device (e.g. mask) MA from a patterning device library, or during a scan, in general, movement of the object tables MT, WT will be realized with the aid of a long-stroke module (coarse positioning) and a short-stroke module (fine positioning), which arc not explicitly depicted in Figure 7. However, in the case of a wafer stepper (as opposed to a step-and-scan tool) the patterning device support (e.g. mask table) MT may just be connected to a short stroke actuator, or may be fixed.
[0074] The depicted tool can be used in two different modes: - In step mode, the patterning device support (e.g. mask table) MT is kept essentially stationary, and an entire mask image is projected in one go (i.e., a single "flash") onto a target portion C. The substrate table WT is then shifted in the x and/or y directions so that a different target portion C can be irradiated by the beam PB; - In scan mode, essentially the same scenario applies, except that a given target portion C is not exposed in a single "flash". Instead, the patterning device support (e.g. mask table) MT is movable in a given direction (the so-called "scan direction", e.g., the y direction) with a speed v, so that the projection beam PB is caused to scan over a mask image; concurrently, the substrate table WT is simultaneously moved in the same or opposite direction at a speed V= Mv, in which M is the magnification of the lens PL (typically, M = 1/4 or 1/5). In this manner, a relatively large target portion C can be exposed, without having to compromise on resolution.
[0075] The concepts disclosed herein may simulate or mathematically model any generic imaging system for imaging sub wavelength features, and may be especially useful with emerging imaging technologies capable of producing wavelengths of an increasingly smaller size. Emerging technologies already in use include EUV (extreme ultra violet) lithography that is capable of producing a 193nm wavelength with the use of a ArF laser, and even a 157nm wavelength with the use of a Fluorine laser. Moreover, EUV lithography is capable of producing wavelengths within a range of 20-5nm by using a synchrotron or by hitting a material (either solid or a plasma) with high energy electrons in order to produce photons within this range. Because most materials are absorptive within this range, illumination may be produced by reflective mirrors with a multi-stack of Molybdenum and Silicon. The multi-stack mirror has a 40 layer pairs of Molybdenum and Silicon where the thickness of each layer is a quarter wavelength. Even smaller wavelengths may be produced with X-ray lithography. Typically, a synchrotron is used to produce an X-ray wavelength. Since most material is absorptive at x-ray wavelengths, a thin piece of absorbing material defines where features would print (positive resist) or not print (negative resist).
[0076] While the concepts disclosed herein may be used for imaging on a substrate such as a silicon wafer, it shall be understood that the disclosed concepts may be used with any type of lithographic imaging systems, e.g., those used for imaging on substrates other than silicon wafers.
[0077] The descriptions above are intended to be illustrative, not limiting. Thus, it will be apparent to one skilled in the art that modifications may be made to the invention as described without departing from the scope of the clauses set out below. Other aspects of the invention are set out as in the following numbered clauses: 1. A method for determining a difference between a reference image and a further image of a pattern, the method comprising: determining a reference imaging function; determining parameters of a difference function representative of a difference between the reference imaging function and a further imaging function; calculating a difference between the reference image and the further image of the pattern based on the difference function and the determined parameters. 2. The method according to clause 1, comprising determining the further imaging function based upon the reference imaging function and the difference function. 3. The method according to clause 2, wherein the further imaging function is determined by concatenating the reference imaging function and the difference function. 4. The method according to any of the clauses 1 to 3, wherein calculating the difference in the reference image and the further image of the pattern comprises determining the difference in image intensity of the reference image and the further image of the pattern. 5. The method according to clause 4, wherein calculating the difference in the reference image and the further image of the pattern further comprises determining the difference between a measure for critical dimensions in the reference image and the further image. 6. The method according to any of the clauses 1 to 5, comprising representing at least a part of the terms of the difference function in eigenfunctions and corresponding eigenvalues, and determining a simplified difference function by disregarding eigenfunctions and eigenvalues based upon the eigenvalues. 7. The method according to clause 6, wherein the simplified difference function comprises less than 300 eigenvalues or less than 100 eigenvalues. 8. The method according to any of the clauses 1 to 7, wherein the reference image models an image of the pattern formed by a reference optical system under reference conditions, and the further image models an image of the pattern formed by a further optical system, or formed by the reference optical system under further conditions, or formed by a modified reference optical system. 9. The method according to clause 8, comprising decreasing the difference between the reference image and the further image by changing the conditions of the further optical system, the reference optical system or the modified reference optical system. 10. The method according to clause 9, wherein the reference optical system is a projection system of a lithographic exposure apparatus. 11. A computer product having machine executable instructions, the instructions to be executed by a machine to perform a method for determining a change in an image of a mask pattern as in any one of the preceding clauses.

Claims (1)

1. Een lithografieinrichting omvattende: een belichtinginrichting ingericht voor het leveren van een stralingsbundel; een drager geconstrueerd voor het dragen van een patroneerinrichting, welke patroneerinrichting in staat is een patroon aan te brengen in een doorsnede van de stralingsbundel ter vorming van een gepatroneerde stralingsbundel; een substraattafel geconstrueerd om een substraat te dragen;en een projectieinrichting in gericht voor bet projecteren van de gepatroneerde stralingsbundel op een doelgebied van het substraat, met het kenmerk, dat de substraattafel is ingericht voor het positioneren van het doelgebied van het substraat in een brandpuntsvlak van de projectieinrichting.
NL2003719A 2008-11-10 2009-10-28 Delta tcc for fast sensitivity model computation. NL2003719A (en)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US11302408P 2008-11-10 2008-11-10
US11302408 2008-11-10
US2009049792 2009-07-07
PCT/US2009/049792 WO2010005957A1 (en) 2008-07-07 2009-07-07 Illumination optimization

Publications (1)

Publication Number Publication Date
NL2003719A true NL2003719A (en) 2010-05-11

Family

ID=42232499

Family Applications (3)

Application Number Title Priority Date Filing Date
NL2003719A NL2003719A (en) 2008-11-10 2009-10-28 Delta tcc for fast sensitivity model computation.
NL2003718A NL2003718A (en) 2008-11-10 2009-10-28 Methods and system for model-based generic matching and tuning.
NL2003729A NL2003729A (en) 2008-11-10 2009-10-30 Smart selection and/or weighting of parameters for lithographic process simulation.

Family Applications After (2)

Application Number Title Priority Date Filing Date
NL2003718A NL2003718A (en) 2008-11-10 2009-10-28 Methods and system for model-based generic matching and tuning.
NL2003729A NL2003729A (en) 2008-11-10 2009-10-30 Smart selection and/or weighting of parameters for lithographic process simulation.

Country Status (4)

Country Link
US (4) US8443307B2 (nl)
JP (2) JP5033859B2 (nl)
CN (1) CN101846886B (nl)
NL (3) NL2003719A (nl)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL2003702A (en) 2008-11-10 2010-05-11 Brion Tech Inc Pattern selection for lithographic model calibration.
NL2003719A (en) 2008-11-10 2010-05-11 Brion Tech Inc Delta tcc for fast sensitivity model computation.
NL2003716A (en) * 2008-11-24 2010-05-26 Brion Tech Inc Harmonic resist model for use in a lithographic apparatus and a device manufacturing method.
JP2010156866A (ja) * 2008-12-27 2010-07-15 Toshiba Corp 特徴量抽出方法、テストパターン選択方法、レジストモデル作成方法および設計回路パターン検証方法
US8805630B2 (en) * 2009-08-25 2014-08-12 Taiwan Semiconductor Manufacturing Company, Ltd. Method and system for modeling in semiconductor fabrication
US8355807B2 (en) * 2010-01-22 2013-01-15 Synopsys, Inc. Method and apparatus for using aerial image sensitivity to model mask errors
US8812145B2 (en) 2010-01-22 2014-08-19 Synopsys, Inc. Modeling mask errors using aerial image sensitivity
KR101096979B1 (ko) * 2010-05-07 2011-12-20 주식회사 하이닉스반도체 반도체 소자의 패턴 균일도 조절 방법
NL2007577A (en) 2010-11-10 2012-05-14 Asml Netherlands Bv Optimization of source, mask and projection optics.
NL2007579A (en) * 2010-11-10 2012-05-14 Asml Netherlands Bv Pattern-dependent proximity matching/tuning including light manipulation by projection optics.
NL2007578A (en) * 2010-11-17 2012-05-22 Asml Netherlands Bv Pattern-independent and hybrid matching/tuning including light manipulation by projection optics.
NL2008924A (en) 2011-06-22 2013-01-02 Asml Netherlands Bv System and method to ensure source and image stability.
NL2008957A (en) * 2011-07-08 2013-01-09 Asml Netherlands Bv Methods and systems for pattern design with tailored response to wavefront aberration.
US9466100B2 (en) 2012-06-06 2016-10-11 Kla-Tencor Corporation Focus monitoring method using asymmetry embedded imaging target
KR101396088B1 (ko) 2012-12-10 2014-05-27 한양대학교 산학협력단 3차원 cmos 전계효과 트랜지스터 및 이의 제조방법
US9519285B2 (en) * 2013-01-23 2016-12-13 Taiwan Semiconductor Manufacturing Company, Ltd. Systems and associated methods for tuning processing tools
US10338480B2 (en) * 2014-05-30 2019-07-02 Nikon Corporation Lithography system, simulation apparatus, and pattern forming method
US10691863B2 (en) 2015-10-19 2020-06-23 Asml Netherlands B.V. Method and apparatus to correct for patterning process error
WO2017067765A1 (en) 2015-10-19 2017-04-27 Asml Netherlands B.V. Method and apparatus to correct for patterning process error
KR102059018B1 (ko) 2015-10-19 2019-12-24 에이에스엠엘 네델란즈 비.브이. 패터닝 공정 오차를 보정하는 장치 및 방법
WO2017067748A1 (en) 2015-10-19 2017-04-27 Asml Netherlands B.V. Method and apparatus to reduce effects of nonlinear behavior
US11314172B2 (en) 2018-03-20 2022-04-26 Asml Netherlands B.V. Instant tuning method for accelerating resist and etch model calibration
US10867112B2 (en) * 2018-06-28 2020-12-15 Taiwan Semiconductor Manufacturing Co., Ltd. Methods of making mask using transmission cross coefficient (TCC) matrix of lithography process optical system
EP3588191A1 (en) * 2018-06-29 2020-01-01 ASML Netherlands B.V. Tuning patterning apparatus based on optical characteristic
US10990002B2 (en) * 2019-07-23 2021-04-27 Taiwan Semiconductor Manufacturing Co., Ltd. Sub-resolution assist features

Family Cites Families (63)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5523193A (en) 1988-05-31 1996-06-04 Texas Instruments Incorporated Method and apparatus for patterning and imaging member
JP2938568B2 (ja) 1990-05-02 1999-08-23 フラウンホファー・ゲゼルシャフト・ツール・フォルデルング・デル・アンゲバンテン・フォルシュング・アインゲトラーゲネル・フェライン 照明装置
US5229872A (en) 1992-01-21 1993-07-20 Hughes Aircraft Company Exposure device including an electrically aligned electronic mask for micropatterning
WO1997033205A1 (en) 1996-03-06 1997-09-12 Philips Electronics N.V. Differential interferometer system and lithographic step-and-scan apparatus provided with such a system
US6689519B2 (en) * 2000-05-04 2004-02-10 Kla-Tencor Technologies Corp. Methods and systems for lithography process control
TW500987B (en) * 2000-06-14 2002-09-01 Asm Lithography Bv Method of operating an optical imaging system, lithographic projection apparatus, device manufacturing method, and device manufactured thereby
TW552561B (en) * 2000-09-12 2003-09-11 Asml Masktools Bv Method and apparatus for fast aerial image simulation
US6809809B2 (en) 2000-11-15 2004-10-26 Real Time Metrology, Inc. Optical method and apparatus for inspecting large area planar objects
TWI285295B (en) * 2001-02-23 2007-08-11 Asml Netherlands Bv Illumination optimization in lithography
TWI285299B (en) 2001-04-04 2007-08-11 Asml Netherlands Bv Lithographic manufacturing process, lithographic projection apparatus, and device manufactured thereby
WO2003104921A2 (en) * 2002-06-07 2003-12-18 Praesagus, Inc. Characterization adn reduction of variation for integrated circuits
US7170604B2 (en) * 2002-07-03 2007-01-30 Tokyo Electron Limited Overlay metrology method and apparatus using more than one grating per measurement direction
JP2004103674A (ja) * 2002-09-06 2004-04-02 Renesas Technology Corp 半導体集積回路装置の製造方法
CA2501571A1 (en) 2002-10-07 2004-04-15 Mempile Inc. Tight focusing method and system
US7245356B2 (en) * 2003-02-11 2007-07-17 Asml Netherlands B.V. Lithographic apparatus and method for optimizing illumination using a photolithographic simulation
JP2004341160A (ja) * 2003-05-15 2004-12-02 Seiko Epson Corp 露光用マスク、光近接効果補正装置、光近接効果補正方法、半導体装置の製造方法および光近接効果補正プログラム
KR101056142B1 (ko) * 2004-01-29 2011-08-10 케이엘에이-텐코 코포레이션 레티클 설계 데이터의 결함을 검출하기 위한 컴퓨터로구현되는 방법
US7242459B2 (en) 2004-01-30 2007-07-10 Asml Masktools B.V. Method of predicting and minimizing model OPC deviation due to mix/match of exposure tools using a calibrated Eigen decomposition model
JP2007520892A (ja) 2004-02-03 2007-07-26 メンター・グラフィクス・コーポレーション イメージの忠実度およびスループットに対する光源の最適化
US7079223B2 (en) 2004-02-20 2006-07-18 International Business Machines Corporation Fast model-based optical proximity correction
US7653890B2 (en) * 2004-04-02 2010-01-26 Cadence Design Systems, Inc. Modeling resolution enhancement processes in integrated circuit fabrication
JP2007535135A (ja) * 2004-04-20 2007-11-29 ライテル・インストルメンツ リソグラフ投影ツールのエミュレーション法
US7403264B2 (en) * 2004-07-08 2008-07-22 Asml Netherlands B.V. Lithographic projection apparatus and a device manufacturing method using such lithographic projection apparatus
US7116411B2 (en) 2004-08-26 2006-10-03 Asml Masktools B.V. Method of performing resist process calibration/optimization and DOE optimization for providing OPE matching between different lithography systems
US7331033B2 (en) 2004-08-27 2008-02-12 Applied Materials, Israel, Ltd. Simulation of aerial images
US7251807B2 (en) 2005-02-24 2007-07-31 Synopsys, Inc. Method and apparatus for identifying a manufacturing problem area in a layout using a process-sensitivity model
US7262831B2 (en) * 2004-12-01 2007-08-28 Asml Netherlands B.V. Lithographic projection apparatus and device manufacturing method using such lithographic projection apparatus
US7297453B2 (en) * 2005-04-13 2007-11-20 Kla-Tencor Technologies Corporation Systems and methods for mitigating variances on a patterned wafer using a prediction model
US7853920B2 (en) * 2005-06-03 2010-12-14 Asml Netherlands B.V. Method for detecting, sampling, analyzing, and correcting marginal patterns in integrated circuit manufacturing
US7374957B2 (en) * 2005-07-11 2008-05-20 Asml Netherlands B.V. Method of calibrating or qualifying a lithographic apparatus or part thereof, and device manufacturing method
US7488933B2 (en) * 2005-08-05 2009-02-10 Brion Technologies, Inc. Method for lithography model calibration
KR100958714B1 (ko) * 2005-08-08 2010-05-18 브라이언 테크놀로지스, 인코포레이티드 리소그래피 공정의 포커스-노광 모델을 생성하는 시스템 및방법
CN102662309B (zh) * 2005-09-09 2014-10-01 Asml荷兰有限公司 采用独立掩模误差模型的掩模验证系统和方法
US7425397B2 (en) * 2005-09-12 2008-09-16 Asml Netherlands B.V. Method of determining an illumination profile and device manufacturing method
US7830493B2 (en) * 2005-10-04 2010-11-09 Asml Netherlands B.V. System and method for compensating for radiation induced thermal distortions in a substrate or projection system
US7921383B1 (en) * 2006-01-11 2011-04-05 Olambda, Inc Photolithographic process simulation including efficient result computation for multiple process variation values
US7433791B2 (en) * 2006-02-17 2008-10-07 Asml Masktools B.V. Method of performing multiple stage model calibration for optical imaging simulation models
WO2007109103A2 (en) * 2006-03-16 2007-09-27 Kla-Tencor Technologies Corporation Method and system for optimizing alignment performance in a fleet of exposure tools
JP2009535609A (ja) * 2006-04-28 2009-10-01 マイクロニック レーザー システムズ アクチボラゲット 画像の記録及び表面の調査のための方法及び装置
US7743357B2 (en) * 2006-05-31 2010-06-22 Synopsys, Inc. Method and apparatus for determining a process model that models the impact of CAR/PEB on the resist profile
US8102408B2 (en) * 2006-06-29 2012-01-24 Kla-Tencor Technologies Corp. Computer-implemented methods and systems for determining different process windows for a wafer printing process for different reticle designs
US8040573B2 (en) * 2006-08-14 2011-10-18 Asml Masktools B.V. Method, program product and apparatus for translating geometrical design rules into boundary conditions in the imaging space so as to define test patterns for use in optical model calibration
US8120753B2 (en) * 2006-11-08 2012-02-21 Asml Masktools B.V. Method, program product and apparatus for generating a calibrated pupil kernel and method of using the same in a lithography simulation process
US8611637B2 (en) * 2007-01-11 2013-12-17 Kla-Tencor Corporation Wafer plane detection of lithographically significant contamination photomask defects
US7873204B2 (en) 2007-01-11 2011-01-18 Kla-Tencor Corporation Method for detecting lithographically significant defects on reticles
WO2008089222A1 (en) * 2007-01-18 2008-07-24 Nikon Corporation Scanner based optical proximity correction system and method of use
US7962863B2 (en) * 2007-05-07 2011-06-14 Kla-Tencor Corp. Computer-implemented methods, systems, and computer-readable media for determining a model for predicting printability of reticle features on a wafer
US7999920B2 (en) * 2007-08-22 2011-08-16 Asml Netherlands B.V. Method of performing model-based scanner tuning
WO2009042671A1 (en) * 2007-09-24 2009-04-02 The Board Of Trustees Of The University Of Illinois Three-dimensional microfabricated bioreactors with embedded capillary network
US7673278B2 (en) * 2007-11-29 2010-03-02 Tokyo Electron Limited Enhanced process yield using a hot-spot library
US7925369B2 (en) * 2007-12-18 2011-04-12 Globalfoundries Inc. Method and apparatus for optimizing models for extracting dose and focus from critical dimension
NL1036750A1 (nl) * 2008-04-14 2009-10-15 Brion Tech Inc A Method Of Performing Mask-Writer Tuning and Optimization.
KR101749987B1 (ko) * 2008-06-03 2017-06-22 에이에스엠엘 네델란즈 비.브이. 모델-기반 공정 시뮬레이션 시스템들 및 방법들
US8006203B2 (en) * 2008-08-28 2011-08-23 Synopsys, Inc. Bulk image modeling for optical proximity correction
US8181128B2 (en) * 2008-10-13 2012-05-15 Synopsys, Inc. Method and apparatus for determining a photolithography process model which models the influence of topography variations
NL2003702A (en) * 2008-11-10 2010-05-11 Brion Tech Inc Pattern selection for lithographic model calibration.
NL2003719A (en) * 2008-11-10 2010-05-11 Brion Tech Inc Delta tcc for fast sensitivity model computation.
US8612903B2 (en) * 2010-09-14 2013-12-17 Luminescent Technologies, Inc. Technique for repairing a reflective photo-mask
NL2007579A (en) * 2010-11-10 2012-05-14 Asml Netherlands Bv Pattern-dependent proximity matching/tuning including light manipulation by projection optics.
NL2007578A (en) * 2010-11-17 2012-05-22 Asml Netherlands Bv Pattern-independent and hybrid matching/tuning including light manipulation by projection optics.
JP2013004672A (ja) * 2011-06-15 2013-01-07 Toshiba Corp シミュレーションモデル作成方法
NL2008924A (en) * 2011-06-22 2013-01-02 Asml Netherlands Bv System and method to ensure source and image stability.
NL2008957A (en) * 2011-07-08 2013-01-09 Asml Netherlands Bv Methods and systems for pattern design with tailored response to wavefront aberration.

Also Published As

Publication number Publication date
JP2010114443A (ja) 2010-05-20
JP5033859B2 (ja) 2012-09-26
US20100260427A1 (en) 2010-10-14
JP2010118655A (ja) 2010-05-27
US20100146475A1 (en) 2010-06-10
US20150074619A1 (en) 2015-03-12
NL2003718A (en) 2010-05-11
NL2003729A (en) 2011-01-10
CN101846886A (zh) 2010-09-29
CN101846886B (zh) 2013-01-09
US8893058B2 (en) 2014-11-18
US8379991B2 (en) 2013-02-19
US8443307B2 (en) 2013-05-14
US20130263064A1 (en) 2013-10-03
JP5033860B2 (ja) 2012-09-26
US10169522B2 (en) 2019-01-01

Similar Documents

Publication Publication Date Title
US8379991B2 (en) Delta TCC for fast sensitivity model computation
US9390206B2 (en) Methods and systems for lithography process window simulation
US8120753B2 (en) Method, program product and apparatus for generating a calibrated pupil kernel and method of using the same in a lithography simulation process
JP4890517B2 (ja) モデルベーススキャナ調整を実行する方法
US9053280B2 (en) Rule optimization in lithographic imaging based on correlation of functions representing mask and predefined optical conditions
EP1688795A2 (en) Method, computer program and apparatus for improving calibration of resist models used in critical dimension calculation
US8521481B2 (en) Method, program product and apparatus for modeling resist development of a lithography process
KR20160062141A (ko) 프로파일 인식 소스-마스크 최적화
US11506984B2 (en) Simulation of lithography using multiple-sampling of angular distribution of source radiation
US10025198B2 (en) Smart selection and/or weighting of parameters for lithographic process simulation
US8792147B2 (en) Method, program product and apparatus for creating optimal test patterns for optical model calibration and for selecting suitable calibration test patterns from an arbitrary layout
CN114514476A (zh) 用于确定视场设定的方法
EP4194950A1 (en) Systems and methods for reducing pattern shift in a lithographic apparatus

Legal Events

Date Code Title Description
WDAP Patent application withdrawn

Effective date: 20100824