The invention relates to a method for simulating illumination and imaging properties of an optical production system during the illumination and imaging of an object by use of an optical measurement system. Further, the invention relates to a metrology system for carrying out such a method.
Such a method and a metrology system to this end are known from DE 10 2019 208 552 A1 and DE 10 2019 215 800 A1. A metrology system for measuring an aerial image of a lithography mask in three dimensions is known from WO 2016/012426 A1. DE 10 2013 219 524 A1 describes a device and a method for determining an imaging quality of an optical system, and an optical system. DE 10 2013 219 524 A1 has described a phase retrieval method for determining a wavefront on the basis of the imaging of a pinhole. The specialist article “A new system for a wafer level CD metrology on photomasks,” proceedings of SPIE—The International Society for Optical Engineering, 2009, 7272, by Martin et al. has disclosed a metrology system for determining a wafer level critical dimension (CD).
It is an aspect of the present invention to improve a method for simulating illumination and imaging properties of an optical production system when illuminating and imaging an object by use of an optical measurement system.
This aspect is achieved according to the invention by a simulation method having the features specified in claim 1.
According to the invention, it was recognized that the recording of measurement aerial images by use of the plurality of pupil stops, in particular the recording of measurement aerial images in a plurality of measurement positions of a pupil stop, which was selected in advance for the best possible simulation of the illumination setting of the optical production system, provides the possibility of improving the accuracy of the simulation method overall and in particular provides the option of reducing illumination angle-dependent artefacts, in particular, in the reconstructed complex mask transfer function, that is to say in the transfer function of the imaged object. It is then possible to correctly take account of 3-D mask effects. This can be taken into account when examining lithography masks, especially when examining masks used for EUV lithography.
Within the scope of the simulation method, it is possible to select precisely one pupil stop from the plurality of pupil stops provided, which may differ in terms of their stop boundary shape and/or in terms of their stop boundary orientation. Alternatively, it is possible to select and use a plurality of different pupil stops for the purpose of specifying different measurement positions. The pupil stops provided may specify at least one of the following illumination settings in particular: Quadrupole, C-quad, dipole, annular, conventional. A person skilled in the art finds examples for such settings in, inter alia, WO 2012/028303 A1. First of all, there may be an initial determination of a best focal plane (defocus value zm=0) within the preparation of the imaging method. z-increments when determining the 3-D aerial image in the last step of the simulation method, i.e. when determining the aerial image from the reconstructed mask transfer function and the optical production system illumination setting, may differ from defocus values that may at first be specified in the simulation method. Pixel sizes of the recorded measurement aerial images may be sampled for the purpose of an adaptation to the desired pixel resolution.
The target pupil stop, which can be specified, and the target stop boundary shape thereof may relate to a plurality of, or else a multiplicity of, individual illumination or pupil spots, i.e. a plurality of stop apertures for example arranged in a grid-like manner. Such illumination or pupil spots may yield an illumination setting used within the scope of the production illumination, the said illumination setting for example being able to be set by way of an illumination optical unit having a field facet mirror and a pupil facet mirror.
A displacement drive as claimed in claim 2 has proven its worth for the reproducible specification of pupil stop measurement positions. This correspondingly applies to the object holder that is displaceable perpendicular to the object plane.
Different stop boundary shapes and/or stop boundary orientations of the provided pupil stops as claimed in claim 3 increase a flexibility when carrying out the simulation method.
A simulation method as claimed in claim 4 with the aid of a qualification algorithm has particularly proven its worth. The defocus values and/or the pupil stop measurement positions can be specified with the aid of an object holder that is displaceable by actuator, and with the aid of a displacement drive for pupil stop displacement. A piezo drive and/or a stepper motor drive can be used as displacement drive or displacement actuator.
Recording the measurement aerial images as claimed in claim 5 has proven its worth in practice.
A central measurement position and a plurality of offset measurement positions as claimed in claim 6 have proven their worth within the scope of the practical implementation of the simulation method. In the central measurement position, the pupil stop is arranged in the center of the utilized pupil of the optical measurement system. Provision can be made of two to ten offset measurement positions, in particular two to five, for example three or four offset measurement positions. The offset measurement positions can be arranged uniformly distributed in the circumferential direction about the central measurement position. The offset measurement positions may be displaced relative to the central measurement position in the Cartesian directions or else in the directions of the quadrants. The measurement positions may be arranged randomly in the circumferential direction and may be arranged on one or more radii, in particular on two or three different radii. A completely random arrangement of the measurement positions within the measurement pupil or else partially outside of the measurement pupil is also possible. Should a random arrangement be mentioned above, the latter can be determined by the use of an algorithmic random function.
Defocus value/measurement position combinations as claimed in claim 7 have proven their worth in practice. It was found that it is not necessary to home in on all pupil stop measurement positions, specified within the scope of the method, for each defocus value. This reduces the measurement time.
A selection method for the pupil stop as claimed in claim 8 ensures the best possible simulation of the target pupil stop using the selected pupil stop. In the respective pupil spots, illumination light is present in the illumination pupil. A distance qualification of assigned pupil spots of the target stop boundary shape of the respective pupil stop can be carried out within the scope of the selection method. A merit function can be defined and minimized within the scope of the selection method.
Modeling the illumination direction-dependent mask spectrum as claimed in claim 9 has proven its worth during the reconstruction as this helps reduce the number of degrees of freedom present during the optimization within the scope of the reconstruction.
Displacing the imaging pupil stop as claimed in claim 10 extends the modelling possibilities within the scope of the simulation method.
A reconstruction as claimed in claim 11 leads to a particularly good simulation.
The result of the simulation method as claimed in claim 12 lies in the option of an aerial image description also dependent on a chief ray angle of an illumination by the optical production system. Therefore, a different illumination chief ray angle of the production system may also be taken into account when determining the 3-D aerial image. This increases the power of the simulation method.
The advantages of the metrology system as claimed in claims 13 to 15 correspond to those that have already been explained above with reference to the method claims.
A selection apparatus having a stop storage unit as claimed in claim 16 advantageously enables the pupil stop selection step of the simulation method. In particular, the selection can be made with the aid of a robotic actuation system which takes the respective selected pupil stop from the stop storage unit and moves it to its use location in the pupil plane. The selection apparatus moreover ensures a substitution of a last-used pupil stop for a newly selected pupil stop. In particular, the last-used pupil stop can be transferred from the use location back to the stop storage unit by use of the robotic actuation system in that case.
An aperture of the stop, i.e. of the illumination pupil stop and/or the imaging pupil stop, may also be variably specifiable, for example in the style of an iris diaphragm.
The metrology system may comprise a light source for the illumination light. A light source of this type may be configured as an EUV light source.
An EUV wavelength of the light source may range between 5 nm and 30 nm. A light source in the DUV wavelength range, for example of the order of 193 nm, is also possible.
Exemplary embodiments of the invention are explained in greater detail below with reference to the drawing, in which:
In order to facilitate the representation of positional relationships, a Cartesian xyz-coordinate system will be used hereinafter. In
In a view that corresponds to a meridional section,
An example of the test structure 5 is depicted in a plan view in
The metrology system 2 is used to analyze a three-dimensional (3-D) aerial image (aerial image metrology system). One application is found in the simulation of an aerial image of a lithography mask, in the way that the aerial image would also appear in an optical production system of a producing projection exposure apparatus, for example in a scanner. To this end, an imaging quality of the metrology system 2 itself, in particular, can be measured and optionally adjusted. Consequently, the analysis of the aerial image can serve to determine the imaging quality of a projection optical unit of the metrology system 2, or else to determine the imaging quality of, in particular, projection optical units within a projection exposure apparatus. Metrology systems are known from WO 2016/012426 A1, from US 2013/0063716 A1 (cf.
The illumination light 1 is reflected and diffracted at the test structure 5. A plane of incidence of the illumination light 1 is parallel to the yz-plane in the case of the central, initial illumination.
The EUV illumination light 1 is produced by an EUV light source 8. The light source 8 can be a laser plasma source (LPP; laser produced plasma) or a discharge source (DPP; discharge produced plasma). In principle, a synchrotron-based light source can also be used, e.g. a free electron laser (FEL). A used wavelength of the EUV light source can range between 5 nm and 30 nm. In principle, in one variant of the metrology system 2, a light source for another used light wavelength can also be used instead of the light source 8, for example a light source for a used wavelength of 193 nm.
An illumination optical unit 9 of the metrology system 2 is arranged between the light source 8 and the test structure 5. The illumination optical unit 9 serves for the illumination of the test structure 5 to be examined, with a defined illumination intensity distribution over the object field 3 and at the same time with a defined illumination angle distribution with which the field points of the object field 3 are illuminated. Such an illumination angle distribution is also referred to as illumination setting.
The respective illumination angle distribution of the illumination light 1 is specified by way of a pupil stop 10, which is arranged in an illumination optical unit pupil plane 11. The pupil stop 10 is also referred to as a sigma stop.
Further variants of pupil stops 10 with a central passage pole I of increasingly larger radius are shown in
Corresponding annular illumination settings can be realized using the embodiments of the pupil stops 10 according to
Measured from the x-coordinate of the pupil stop 10 of
The pupil stop 10 of the illumination optical unit 9 is embodied as a stop which is displaceable in driven fashion and which is arranged in front of the object plane 4 in an illumination light beam path 15 of the illumination light 1. A drive unit used for the driven displacement of the pupil stop 10 is depicted at 16 in
With the aid of the displacement drive 16, it is possible to displace the selected pupil stop 10 along the pupil coordinates kx and ky in the pupil plane 11.
The displacement drive 16 may also include a stop interchange unit, by means of which a specific pupil stop 10 is replaced with another, specific pupil stop 10. To this end, the stop interchange unit may take the respective selected pupil stop from a stop storage unit and return the replaced stop to this stop storage unit.
The test structure 5 is held by an object holder 17 of the metrology system 2. The object holder 17 cooperates with an object displacement drive 18 for displacing the test structure 5, in particular along the z-coordinate.
Following reflection at the test structure 5, the electromagnetic field of the illumination light 1 has a distribution 19 which is depicted in
The illumination light 1 reflected by the test structure 5 enters an imaging optical unit or projection optical unit 20 of the metrology system 2.
A diffraction spectrum 21 arises in a pupil plane of the projection optical unit 20 on account of the periodicity of the test structure 5 (cf.
In the diffraction spectrum 21, the 0th order of diffraction of the test structure 5 is present centrally. In addition,
The orders of diffraction of the diffraction spectrum 21 depicted in
The imaging pupil stop 23 is operatively connected to a displacement drive 25, the function of which corresponds to that of the displacement drive 16 for the sigma stop 10.
The pupils 24 (cf.
Firstly, the images of the −1st, 0th and +1st orders of diffraction and, secondly, an imaging contribution of the optical system, specifically the projection optical unit 20, contribute to the intensity distribution in the exit pupil 26. This imaging contribution which is elucidated in
The projection optical unit 20 images the test structure 5 towards a spatially resolving detection device 27 of the metrology system 2. The detection device 27 is in the form of a camera, in particular a CCD camera or CMOS camera.
The projection optical unit 20 is embodied as magnifying optical unit. A magnification factor of the projection optical unit 20 may be greater than 10, may be greater than 50, may be greater than 100 and may even be greater still. As a rule, this magnification factor is less than 1000.
In a manner corresponding to
The following procedure is carried out to simulate the illumination and imaging properties of the optical production system when illuminating and imaging the object, using the example of the test structure 5, by use of the optical measurement system 1 of the metrology system 2:
Initially, a plurality of pupil stops 10 with in each case different stop boundary shapes are provided for the purposes of specifying correspondingly different measurement illumination settings. This is implemented by providing pupil stops 10, for example in the style of the pupil stops 10 of
Then, a target pupil stop with a target stop boundary shape is specified proceeding from an illumination setting of the optical production system to be simulated. The target pupil stop can be an arrangement of a plurality or multiplicity of individual pupil spots or stop spots. In this case, the intensity of individual illumination spots or pupil spots generally differs between the individual spots.
The target pupil stop 36 can be specified by way of a definition of appropriate stop aperture contours, especially continuous stop aperture contours. Such stop aperture contours can be described by polygonal chains, for example.
These continuous openings are then approximated by a finite number of pupil spots 37 within the openings. These spots are depicted in
For the specific example in
The finer the grid of illumination spots, the more accurately the actual stop shape can be approximated.
Proceeding from this target pupil stop 36, at least one pupil stop 10 is then selected from the provided plurality of pupil stops 10 by use of an algorithm which qualifies deviations between the respective stop boundary shape of the provided pupil stops 10 and the target stop boundary shape of the target pupil stop 36. To this end, the pupil stop 10 currently under examination during the selection (also referred to as pupil stop to be qualified below) can in turn be decomposed within its stop boundary into a plurality of pupil spots 38 arranged in grid-like fashion and represented by circles in
The scope of qualification comprises determining the similarity between the target illumination pupil (also denoted “T” below) and the possible measurement stops 10 (also denoted “M” below). For instance, this can be implemented by calculating an overlap function Q.
Here, A is a function for (approximately) calculating the area. The first term corresponds to the normalized area of the overlap between measurement stop and target illumination pupil. The second and third terms correspond to the normalized difference area between the measurement stop and the target illumination pupil, and vice versa. The difference area is intended to refer to the area contained only in the first pupil and not in the second.
The operators “∩”, “∪” and “\” correspond to the intersection (∩), union (∪) and relative complement (\) operators from set theory. In this case, the intersection M1∩M2 of the sets/areas M1 and M2 is intended to mean the set/area which is contained both in M1 and in M2, i.e. corresponds to the overlap area of M1 and M2. The union M1∪ M2 of the sets/areas M1 and M2 describes the set/area which is contained in M1 or M2, i.e. corresponds to the overall area covered by M1 or M2. The relative complement M1\M2 of the sets/areas M1 and M2 describes the set/area which is covered by M1 but not contained in M2.
For instance, the area function A can be implemented as counting illumination spots in the pupil. To this end, the target illumination pupil and measurement pupil are equipped with the same grid. Typically, the grid corresponds to the pupil facet grid in the scanner on which the target illumination pupil is sampled (cf.
Thus, the selection of the pupil stop 10 encompasses a comparison between the poses of pupil spots 37 of the target stop boundary shape and the poses of pupil spots 38 of the provided pupil stops 10.
Moreover, a plurality of defocus values zm (cf.
Moreover, a plurality of measurement positions (kx, ky) of the selected pupil stop 10 are specified within the scope of the simulation method.
Now, measurement aerial images I(x, y) in the style of the intensity distributions 31 according to
The sequence in
In comparison with the imaging pupil stop 23,
In comparison with the centered position according to
An alternative sequence of measurement positions (kx, ky) of the pupil stop 10 is depicted in
Relative to the imaging pupil stop 23,
Relative to the imaging pupil stop 23,
Relative to the imaging pupil stop 23,
The completed sequence of measurement positions (kx, ky) is shown in
The selection of the respective measurement position sequence, or optionally a subset therefrom, is implemented on the basis of the arrangement of individual structures of the test structure 5 and/or on the basis of the illumination setting of the optical production system to be simulated. For instance, the measurement position sequence can be selected in a manner analogous to the stop selection algorithm (see above), with all stop positions of a sequence being taken into account and the sequence being selected for which the overlap of the measurement sequence with the target illumination pupil is maximal.
The poses of the pupil stop 10 which differ from the center position in terms of the relative pose with respect to the imaging pupil stop 23 are also referred to as offset measurement positions. Within the scope of a measurement position sequence, two to ten such offset measurement positions can be homed in on, this typically being two to five offset measurement positions, for example three or four offset measurement positions. The offset measurement positions can be arranged uniformly distributed in the circumferential direction. To reduce the measurement time, it is also possible to use only a subset, e.g. every second measurement position, from the measurement schemes (
The specified defocus values zm are all measured with the aid of the respective measurement position sequence. In an alternative, it is possible that the entire respective measurement position sequence is used only for one defocus value or for individual defocus values zm, with the measurement aerial images being recorded for fewer measurement positions of the pupil stop relative to the imaging pupil stop 23 in the case of other defocus values zm. In an extreme case, it is possible for example to home in on the full measurement position sequence and record a respective measurement aerial image there for only one defocus value zm, whereas the measurement aerial image Imeas(x, y) is only recorded at one respective measurement position, in particular for the case of the centered pupil stop 10, in the case of the other specified defocus values zm.
For instance, the following defocus value/measurement position combinations can be recorded: A central defocus value zm and a plurality of measurement positions (kx, ky) of the pupil stop 10, i.e., in particular, a centered measurement position and a plurality of offset measurement positions, and defocus values zmin, zmax maximally offset from the central defocus value on both sides, with exactly one central measurement position (kx, ky) of the pupil stop 10 being adopted at these positions zmin, zmax.
Then, a complex mask transfer function is reconstructed from the totality of measurement aerial images recorded with the selected pupil stop 10. A similar reconstruction step is also described in DE 10 2019 215 800 A1.
The reconstruction is implemented within the scope of a modelled description, within which the projection optical unit 20 of the metrology system 2 with the illumination setting specified by the pupil stop 10 is described by a function σ({right arrow over (p)}) which reproduces the illumination directions {right arrow over (p)} which are passed through the pupil stop 10. A displacement of the pupil stop 10 by a vector {right arrow over (q)} with absolute coordinate values kx and ky leads to a displaced illumination function σ({right arrow over (p)}−{right arrow over (q)}).
Each illumination direction generates a complex-valued field distribution m({right arrow over (r)},{right arrow over (p)}) (cf. the field distribution 19 in
is the curtailment by the numerical aperture of the imaging optical unit 20, that is to say by the imaging pupil stop 23, and
is a wavefront error caused by a defocus z (displacement by the object holder 17). The propagated spectrum (cf.
In this case, {right arrow over (r)} is the xy-position of the intensity measurement, i.e. the respective pixel of the camera 27.
The object now is to determine the mask spectrum M({right arrow over (k)},{right arrow over (p)}). In this case, k is the pupil coordinates in the entrance pupil 24 of the projection optical unit 20 and {right arrow over (p)} is the illumination direction. The Fourier transform of the respective mask spectrum is the associated mask transfer function.
The reconstructed spectra can then be used to calculate the aerial image for any other illumination setting σtarget({right arrow over (p)})) and any defocus ztarget.
The determination of M({right arrow over (k)},{right arrow over (p)}) can be formulated as an optimization problem: Sought are the spectra M({right arrow over (k)},{right arrow over (p)}) for which there is a minimum deviation F between the simulated aerial images and the aerial images Imeas measured at the defocus positions z1, z2 . . . zN and the illumination directions {right arrow over (q)}1, {right arrow over (q)}2 . . . {right arrow over (q)}M. The following optimization problem should be solved:
A separate spectrum needs to be reconstructed for each illumination direction 5. As a rule, the optimization problem is underdetermined. There are various options for handling this problem.
The simplest solution is a Hopkins approximation, which assumes that the spectrum is only displaced by the same value in the case of a displacement of the illumination direction, that is to say M({right arrow over (k)},{right arrow over (p)})=M0({right arrow over (k)}−{right arrow over (p)}). As a result, there now is only still one spectrum that it needs to be reconstructed. The angle dependence of the reflectivity, shadowing effects and mask-induced aberrations means that the dependence of the mask spectrum on the illumination direction is not completely negligible in the case of real EUV lithography masks as test structures 5. The Hopkins approximation is stretched to its limits.
To take account of the dependence of the spectrum and the illumination direction 5, the following ansatz can be considered for the angle-dependent spectrum M of the test structure 5:
In this case, M0({right arrow over (k)}) is a spectrum that is independent of the illumination direction, in a manner analogous to the Hopkins approximation. C({right arrow over (k)},{right arrow over (p)},{right arrow over (α)}) is any complex-valued function albeit defined prior to the reconstruction, which models the dependence of the amplitude and phase on the illumination direction. α1 . . . N are free parameters, which are determined within the scope of the optimization.
By way of example, the following the function C({right arrow over (k)}, {right arrow over (p)}, α1, α2, . . . , αN) could be used:
Within the scope of reconstructing the complex mask transfer function M, a mask spectrum M({right arrow over (p)}) that is dependent on the illumination direction is modelled as a product of a spectrum that is independent of the illumination direction and a correction function (C({right arrow over (k)}, {right arrow over (p)}, {right arrow over (α)}1, {right arrow over (α)}2, . . . , αN)).
Now, the mask spectrum M0({right arrow over (k)}) and the parameters α1 . . . N which minimize the difference between measured and simulated aerial images are sought after. The following optimization problem is solved:
Thus, the number of free parameters is only increased by N vis-a-vis the Hopkins approximation, with N typically being small.
Using the reconstructed, now directionally dependent spectrum, it is possible to calculate a simulated aerial image Isim for the target illumination setting Gtarget and the target defocus ztarget:
Equation (6) then allows comparison between the simulated aerial image Isim and the respectively measured aerial image Imeas, and this can be used to reconstruct the mask spectrum M and, accordingly, the complex mask transfer function.
From Equation (6), the 3-D aerial image can be calculated with the aid of the reconstructed mask transfer function M and the illumination setting Gtarget of the optical production system. In this way, it is possible to ascertain what the aerial image of the test structure 5 would look like if it were imaged by the optical production system.
In a variant of the simulation method, it is also possible to use a plurality of different pupil stops 10 to specify the various measurement positions (kx, ky).
To prepare the simulation method, it is possible to record an aerial image stack in order to make sure which z-pose of the object plane 4 supplies an optimally sharp image in the image plane 29 (zero of the z-pose).
z-increments which are used in Equation (6) when determining the aerial image Isim may differ from the defocus values zm that are specified within the scope of the simulation method.
Pixel sizes of the recorded measurement aerial images Imeas may be re-sampled for the purpose of matching to a desired pixel resolution.
A plurality of kx, ky positions of the imaging pupil stop 23 can also be set by way of the displacement drive 25 in a simulation method.
When reconstructing the mask transfer function, it is accordingly possible to take account of imaging aberrations of the optical measurement system, in particular imaging aberrations of the imaging optical unit 20 of the metrology system 2.
The determination of the 3-D aerial image Imeas and/or the calculation of the simulated aerial image Isim may be carried out using a different illumination chief ray angle to that of the reconstruction of the mask transfer function.
For selecting the respective pupil stop 10 from the provided plurality of pupil stops 10 with in each case different stop boundary shapes and/or stop boundary orientations, the metrology system 2 has a selection apparatus not depicted in detail in the drawing. This selection apparatus has a stop storage unit, in which the plurality of pupil stops 10 with in each case different stop boundary shapes and/or stop boundary orientations are stored for the purpose of specifying correspondingly different measurement illumination settings.
In the selection step of the simulation method, the last pupil stop inserted is firstly removed from its use location in the pupil plane 11 and supplied to the stop storage unit in the selection apparatus with the aid of an actuator system of the selection apparatus, in particular with the aid of a robotic actuator system. Subsequently, the pupil stop 10 selected according to the simulation method is selected from the stop storage unit and inserted in the use position in the pupil plane 11 with the aid of the robotic actuator system.
Number | Date | Country | Kind |
---|---|---|---|
102022200372.1 | Jan 2022 | DE | national |
The present patent application is a continuation of and claims benefit under 35 U.S.C. § 120 from PCT Application No. PCT/EP2023/050450, filed on Jan. 10, 2023, which claims priority to German patent application DE 10 2022 200 372.1, filed on Jan. 14, 2022. The entire contents of each of these earlier applications are incorporated by reference herein.
Number | Date | Country | |
---|---|---|---|
Parent | PCT/EP2023/050450 | Jan 2023 | WO |
Child | 18768476 | US |