CN112394615A - Extreme ultraviolet lithography light source mask optimization method - Google Patents

Extreme ultraviolet lithography light source mask optimization method Download PDF

Info

Publication number
CN112394615A
CN112394615A CN202011279713.8A CN202011279713A CN112394615A CN 112394615 A CN112394615 A CN 112394615A CN 202011279713 A CN202011279713 A CN 202011279713A CN 112394615 A CN112394615 A CN 112394615A
Authority
CN
China
Prior art keywords
mask
matrix
light source
calculating
optimization
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011279713.8A
Other languages
Chinese (zh)
Other versions
CN112394615B (en
Inventor
张子南
李思坤
王向朝
成维
胡少博
刘珍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Optics and Fine Mechanics of CAS
Original Assignee
Shanghai Institute of Optics and Fine Mechanics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Optics and Fine Mechanics of CAS filed Critical Shanghai Institute of Optics and Fine Mechanics of CAS
Priority to CN202011279713.8A priority Critical patent/CN112394615B/en
Publication of CN112394615A publication Critical patent/CN112394615A/en
Application granted granted Critical
Publication of CN112394615B publication Critical patent/CN112394615B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/76Patterning of masks by 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
    • 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/80Etching
    • 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/70425Imaging strategies, e.g. for increasing throughput or resolution, printing product fields larger than the image field or compensating lithography- or non-lithography errors, e.g. proximity correction, mix-and-match, stitching or double patterning
    • G03F7/70433Layout for increasing efficiency or for compensating imaging errors, e.g. layout of exposure fields for reducing focus errors; Use of mask features for increasing efficiency or for compensating imaging errors

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
  • Exposure Of Semiconductors, Excluding Electron Or Ion Beam Exposure (AREA)

Abstract

A mask optimization method for an extreme ultraviolet lithography light source adopts pixels to represent a light source and a mask pattern, calculates a photoresist image through a thick mask imaging model corrected by point pulses, takes a pattern error between the photoresist pattern and a target pattern as an evaluation function, and optimizes light source distribution and the mask pattern by adopting a particle swarm optimization (SLPSO) based on a social learning mechanism. The invention utilizes the extreme ultraviolet lithography thick mask imaging model to establish the evaluation function, and improves the accuracy of imaging calculation in the optimization process. In addition, the optimization efficiency and the optimization capability are effectively improved by utilizing the SLPSO algorithm, and the imaging quality is improved.

Description

Extreme ultraviolet lithography light source mask optimization method
Technical Field
The invention relates to the technical field of extreme ultraviolet lithography resolution enhancement, in particular to an extreme ultraviolet lithography light source mask optimization method.
Background
The extreme ultraviolet lithography (EUVL) technology is a lithography technology which is extremely important in the large-scale mass production manufacturing of chips with technology nodes of 7nm and below. Similar to deep ultraviolet lithography (DUVL), in extreme ultraviolet lithography, the diffraction effects of the mask cause a degradation in the quality of the lithographic image, and resolution enhancement techniques are therefore required. The light Source Mask Optimization (SMO) technique is an important resolution enhancement technique, and can improve the imaging quality and enlarge the process window by jointly optimizing the light source distribution and the mask pattern. The mask of several hundreds of nanometers makes EUV imaging calculation complicated, and furthermore, the reflective mask and optical system cause shadow effect and defocus effect peculiar to EUV lithography, which are problems to be considered when applying SMO technology to EUV.
Prior art 1 (see prior art 1, Ma X, Wang Z, Chen X, et al, gradient-Based Source Mask Optimization for Extreme Ultraviolet Lithography [ J ]. IEEE Transactions on Computational Imaging,2018,5(1):120-135.) proposes a gradient-Based Extreme Ultraviolet Lithography light Source Mask Optimization technique, which optimizes an objective function using a gradient descent method. However, in the technology, a kirchhoff thin mask model is still adopted in the process of deducing the analytical expression of the target function, and the EUV thick mask effect cannot be reflected. In addition, gradient-based SMO techniques are limited in their application as lithographic imaging models and photoresist models become increasingly complex and even impossible to interpret. Prior art 2 (see prior art 2, Lei Wang, Sikun Li, Xiangzhao Wang, Guanyong Yang, charging Yang, Pixed source optimization for optical lithography via partial simulation, J.Micro/Nanolith. MEMS MOEMS,2016,15(1),013506) and prior art 3 (see prior art 3, Lei Wang, Sikun Li, Xiangzhao Wang, charging Yang, Feng Tang, Pixel based mask optimization video partial simulation, Proc. SPIE 2016,9780:97801V) propose a Particle Swarm Optimization (PSO) based light source and mask optimization technique, which is different from the intelligent method, and is applicable to any prior lithography modeling, without the need for knowledge-based optimization of a prior art, and without the need for a priori knowledge of a lithographic objective model. The technology can be transferred and applied to EUV lithography easily, but when the sampling points of a light source and a mask are gradually increased, the optimization dimension is rapidly increased, and the optimization efficiency and the optimization capability of a PSO algorithm are greatly reduced. Furthermore, the EUV complex lithography imaging model also causes the optimization speed to be greatly reduced, so that an optimization algorithm with better performance needs to be adopted.
Disclosure of Invention
The invention aims to provide an extreme ultraviolet lithography light source mask optimization method, which is based on an EUV thick mask model and a vector imaging model, optimizes a light source and a mask pattern by utilizing an SLPSO algorithm and improves imaging quality.
The technical solution of the invention is as follows:
an extreme ultraviolet lithography light source mask optimization method mainly comprises the following steps:
step 1, setting a light source coding and decoding mode and a mask coding and decoding mode:
the light source is characterized by a two-dimensional matrix S with the size of NS×NS,NSIs an odd number. Each element of the matrix represents a light source point, the value of the element is real number, the value of the element represents the intensity of the light source point, and the value range is [0,1 ]]. The position of the light source point in the matrix is denoted as S (m, n) and in frequency domain coordinates as S (f)s,gs) Wherein m is more than or equal to 1 and less than or equal to NS,1≤n≤NS, -1≤fs≤1,-1≤gsLess than or equal to 1. The incoherence and incoherence factors of the light source are expressed as sigmainAnd σoutIn order to ensure the symmetry of the light source, N which is the total number of the internal and external coherence factors in the first quadrant in the frequency domain coordinate system is selectedScThe individual light source points are coded and their intensity values are recorded as
Figure RE-GDA0002875501180000021
The position in the matrix is recorded as
Figure RE-GDA0002875501180000022
Mixing J and PSAs a code for the light source.
In the case of a lithographic imaging calculation,mixing J and PSDecoding into a complete light source matrix S: firstly resetting S to be a full 0 matrix, and then sequentially carrying out i is more than or equal to 1 and less than or equal to N on the ith light source point in the light source codeScFrom PSRead its position index (m)i,ni) Reading the intensity value J from JiSelecting the position S (m) in Si,ni),S(NS+1-mi,ni), S(mi,NS+1-ni),S(NS+1-mi,NS+1-ni) Is set to have an element value of JiAnd finally obtaining the decoded complete light source matrix S.
The mask is characterized by a two-dimensional matrix M with the size of NM×NM,NMIs an odd number. Each element of the matrix represents a pixel of the mask pattern, the value of the element is 0 or 1, the element of 0 represents that the pixel is a background region, and the element of 1 represents that the pixel is a pattern region. The position of a mask pixel in the matrix is denoted M (p, q), where 1 ≦ p ≦ NM,1≤q≤NM. Will be the (N) th of the matrixM+1)/2 line and (N)MUsing +1)/2 columns as a reference, and selecting all matrix elements as pixels to be coded when the shape of the mask pattern is asymmetric; when the shape of the mask graph is symmetrical up and down, selecting elements of a first quadrant and a second quadrant as pixels to be coded; when the shape of the mask graph is symmetrical left and right, selecting elements of a first quadrant and a fourth quadrant as pixels to be coded; when the shape of the mask pattern is symmetrical up, down, left and right, the element of the first quadrant is selected as the pixel to be encoded. The number of pixels to be coded is recorded as NMcRecord their element values as
Figure RE-GDA0002875501180000031
The position in the matrix is recorded as
Figure RE-GDA0002875501180000032
Mixing T and PMAs the coding of the mask pattern.
When photoetching imaging calculation is carried out, T and P are addedMDecoding into a complete mask matrix M: first, M is reset to all 0Matrix, then sequentially for the jth pixel in the mask code, 1 ≦ j ≦ NMcFrom PMRead its position index (p)j,qj) Reading the intensity value T from TjSelecting M (p) as the position of M in accordance with the symmetry of the original mask patternj,qj) And its corresponding symmetric element, if TjGreater than a threshold value TmrIts element value is set to 1, otherwise its element value is set to 0, and finally the decoded complete mask matrix M is obtained.
Step 2, initializing a light source matrix, a mask matrix and photoresist parameters:
initializing a light source matrix to S according to light source parametersinit(ii) a Initializing a mask matrix to M according to a mask patterninitWill MinitAs a matrix M of target patternstarget(ii) a Initializing a photoresist threshold trSensitivity of photoresist α, developing threshold t of photoresistrd. Initializing the normalized intensity I of the aerial imagenorm
Step 3, mask defocus optimization:
since the pattern on the silicon wafer is shifted due to the oblique incidence of the illumination light source on the mask in the extreme ultraviolet lithography, it is necessary to correct the shift of the pattern by setting the mask to be out of focus. Setting the mask defocusing amount as a parameter to be optimized, and taking the value as [ delta z [ ]min,Δzmax]Calculating photoresist images under different mask defocusing amounts by using a thick mask model and a vector imaging model according to real numbers in the range, taking Edge Placement Errors (EPE) of the photoresist images and a target image as evaluation functions, and optimizing by adopting a demarcation method to obtain the optimal mask defocusing amount delta zbestThe method comprises the following specific steps:
(1) initializing a maximum number of iterations NiterLower bound of variable zl=ΔzminUpper bound of variable zr=ΔzmaxIntermediate values z of the upper and lower bounds of the variablemid=(zl+zr) /2, number of initialization iterations niter=0。
(2) Will zl,zmidAnd zrInputting a thick mask model and vector imaging as mask defocus amounts, respectivelyThe model calculates the photoresist image and then calculates the corresponding evaluation function value fl,fmidAnd fr
(3) Calculating a temporary variable ztemp=(zl+zmid) And/2, inputting the mask defocus quantity into a thick mask model and a vector imaging model to calculate a photoresist image, and then calculating corresponding evaluation function values ftemp. If ftemp>fmidLet zl=ztempOn the contrary, z ismid=ztemp
(4) Calculating a temporary variable ztemp=(zr+zmid) And/2, inputting the mask defocus quantity into a thick mask model and a vector imaging model to calculate a photoresist image, and then calculating corresponding evaluation function values ftemp. If ftemp>fmidLet zr=ztempOn the contrary, z ismid=ztemp
(5) Let n beiter=niter+1, judgment of niter=NiterWhether or not this is true. If not, the optimization is continued to (2). If yes, the procedure is terminated, and the optimal mask defocus amount Delta z is obtainedbest=zmid
Step 4, calculating the thick mask parameters:
calibration of absorber transmittance t for mask background region in thick mask model by rigorous electromagnetic field simulationaTransmittance t of absorption layer in pattern regionbAnd a dot pulse correction amount δ at the pattern boundary position, which are all complex numbers.
According to the incident angle theta of the central point of the illumination light sourceincAnd azimuth angle
Figure RE-GDA0002875501180000041
Calculating the phase transmission factor phi of the incident light from the upper surface of the absorption layer to the middle position of the absorption layeraAnd a phase transmission factor phi of incident light from the middle position of the absorption layer to the lower surface of the absorption layerb。φaIs a complex number phibIs a complex matrix of size NO×NOIn which N isOFor diffraction orders of the mask in the x-directionOr the total number in the y-direction. Calculating the reflection matrix R of the multilayer film of the extreme ultraviolet lithography mask according to the equivalent film layer methodmlSize is NO×NO
The optimal mask defocus amount Delta z obtained according to the step 3bestCalculating the mask defocus matrix phidSize is NO×NO
Step 5, optimizing a light source:
(1) according to the initialized mask matrix MinitAn illumination cross-coefficient (ICC) matrix of the lithography imaging system is calculated from the thick mask model and the vector imaging model.
(2) And constructing an evaluation function for optimizing the light source. Calculating a normalized aerial image I according to the normalized intensity of the aerial image, the light source matrix and the ICC matrixAThen, calculating a photoresist image I through a Sigmoid photoresist modelR=1/(1+exp(-α*(IA-tr) ) is used in the process) and finally according to the photoresist development threshold t)rdCalculating the photoresist pattern
Figure RE-GDA0002875501180000042
The photoresist pattern has a pass size of NM×NMThe binary matrix of 0 and 1 indicates that the element values have the same meaning as the element values of the mask matrix. The evaluation function FSO for light source optimization is expressed as
Figure RE-GDA0002875501180000043
(3) Obtaining an initial light source code J according to the method in the step 1initThe method comprises the following steps of optimizing a light source through an SLPSO algorithm:
firstly, initializing the maximum iteration number NiterLet iteration number n iter0. Initializing particle swarm size MpopDetermining the particle dimension D according to the number of the coded light source pointspop=NScSetting a social influence factor epsilon as kappa multiplied by Dpop/MpopWhere κ is a coefficient.
② setting the first particle position of the particle group as JinitThe positions of the rest particles are randomly generated. Particle position through a size and JinitThe same row vector characterization with dimension DpopAnd each dimension has a value of 0-1, and the random numbers are generated uniformly. The position of the ith particle is denoted as piAll particles having a position pass size of Mpop×DpopMatrix P ofpopAnd (4) showing. The velocity of each particle is characterized by a row vector of the same magnitude as the particle position, and the velocities of all particles are initialized to 0. The velocity of the ith particle is denoted viAll particles have a velocity passage size of Mpop×DpopMatrix V ofpopAnd (4) showing.
Calculating learning probability factor P for all particles in the particle swarmLThe learning probability factor of the ith particle can be expressed as
Figure RE-GDA0002875501180000051
Fourthly, calculating the evaluation function value of all the particles in the particle group, wherein the evaluation function value of the ith particle is expressed as fi. Then, the evaluation function values are sorted in descending order, and the evaluation function values of all the particles are sorted into
Figure RE-GDA0002875501180000052
And adjusting the speed of the particle swarm according to the corresponding sequence. The particle position at which the evaluation function value is minimum is defined as an optimum particle position pbest
Calculating the average position of the particle swarm
Figure RE-GDA0002875501180000053
M (M) after sortingpopThe particles are the optimal particles for the 1 st to M thpop-1 particle is updated in turn. For the ith particle (in this case 1 < i < M)pop-1) with a learning probability factor Pi L. Generating a random number p between 0 and 1rIf p isr>Pi LIf the particle is not changed, the velocity and position of the particle are updated:
the j dimension v of the particle velocityi,jThe updating is as follows:
vi,j=r1×vi,j+r2×(pk,j-pi,j)+r3×ε×(pmean,j-pi,j),
wherein r is1,r2,r3Is a random real number between 0 and 1, and k is from i +1 to Mpop-1, random integer. p is a radical ofk,jDimension j, p, representing the position of the kth particlei,jDimension j, p, representing the position of the ith particlemean,jRepresents the j-th dimension of the mean position of the particle population.
The j dimension p of the particle positioni,jThe updating is as follows:
pi,j=pi,j+vi,j
sixthly, niter=niter+1, judgment of niter=NiterWhether or not this is true. If not, the optimization is continued to the fourth step. If yes, the procedure is terminated, and the optimal particle position p is outputbest
To pbestDecoding according to the light source decoding method in the step 1 to obtain the optimal light source
Figure RE-GDA0002875501180000054
Step 6, optimizing a mask:
(1) according to the optimized light source matrix
Figure RE-GDA0002875501180000055
A cross-transmission coefficient (TCC) matrix of the lithography imaging system is calculated through a vector imaging model.
(2) An evaluation function for mask optimization is constructed. Calculating a mask diffraction spectrum matrix B with a size of N according to the mask matrix and the thick mask parametersO×NOThen combining the TCC matrix and the normalized intensity of the space image to calculate a normalized space image IAThen passing through a Sigmoid photoresist dieCalculating the photoresist image IR=1/(1+exp(-α*(IA-tr) ) is used in the process) and finally according to the photoresist development threshold t)rdCalculating the photoresist pattern
Figure RE-GDA0002875501180000061
Evaluation function F for mask optimizationMOIs shown as
Figure RE-GDA0002875501180000062
The specific process of calculating the mask diffraction spectrum matrix B in the process of constructing the evaluation function is as follows:
generating an absorption layer transmittance matrix:
initializing the absorption layer transmittance matrix TabsorberIs of size NM×NMLet the position of the pixel in the background area with 0 element in the mask matrix M be denoted as Pbg(P, q), the position of the pixel of the graphic region having the element 1 is denoted as Pft(P, q) the position of the boundary pixel of the pattern detected by the difference operator is expressed as Pedge(p, q). At TabsorberIn the matrix, the position is equal to PbgThe element value of (p, q) is set to taWill be equal to PftThe element value of (p, q) is set to tbWill be equal to PedgeThe element value of (p, q) is set to δ. Will TabsorberFourier transform is carried out, and N is selected by taking 0 level as a centerO×NORange as TabsorberFourier coefficient matrix B ofT
Calculating a diffraction spectrum of the absorption layer after the first diffraction:
the diffraction spectrum of incident light propagating from the upper surface to the lower surface of the absorbing layer has a transit size NO×NOMatrix B ofabsorberExpressed, the calculation formula is:
Babsorber=φa⊙BT⊙φb
wherein |, indicates that the corresponding elements of the matrix are multiplied.
Calculating the diffraction spectrum after the reflection of the multilayer film:
light diffracted by the absorption layer passes through the multiple layersThe layer film is reflected, and the diffraction spectrum after reflection has a passing size of NO×NOMatrix B ofmlExpressed, the calculation formula is:
Bml=Babsorber⊙Rml
calculating the mask diffraction spectrum after the second diffraction of the absorption layer:
after the multilayer film is reflected, light in all directions is subjected to secondary diffraction through the absorption layer, and the light is mutually coherent and superposed to obtain a final mask diffraction spectrum B, wherein the calculation formula is as follows:
Figure RE-GDA0002875501180000063
wherein
Figure RE-GDA0002875501180000064
Representing a matrix convolution. Obtaining the optimal defocus amount Delta z of the mask according to optimizationbestAnd a mask defocus matrix phi obtained by pre-pressing calculation in advancedUpdating the mask diffraction spectrum:
B=B⊙Φd
(3) obtaining an initial mask code M according to the method in step 1initThe method comprises the following steps of optimizing a light source through an SLPSO algorithm: initializing a maximum number of iterations NiterLet iteration number n iter0. Initializing particle swarm size MpopDetermining the particle dimension D according to the number of the coded light source pointspop=NMcSetting a social influence factor epsilon as kappa multiplied by Dpop/MpopWhere κ is a coefficient. Then, optimizing according to the methods from step (5) to step (3), in particular, setting the position of the first particle to be the initial mask code M when the position of the particle swarm is initializedinitOptimizing the output pbestDecoding according to the mask decoding method in the step 1 to obtain an optimal mask
Figure RE-GDA0002875501180000071
Compared with the prior art, the invention has the following advantages:
1. the invention adopts the thick mask model and the vector imaging model to calculate the space image of the mask, has higher imaging simulation precision, can embody the thick mask effect, the shadow effect and the like inherently, and then improves the imaging quality through optimization without correcting the shadow effect of the optimized mask image.
2. The invention adopts the SLPSO algorithm as the optimization algorithm, and can effectively improve the optimization efficiency and the optimization capability by guiding the optimization direction by introducing a social learning mechanism.
Drawings
FIG. 1 is a schematic diagram showing the comparison of an initial light source, an initial mask, an optimized light source and an optimized mask adopted by the present invention.
FIG. 2 is a schematic diagram showing the comparison of the photoresist pattern in the initial state, the mask defocus state and the optimized state with the target pattern profile in the optimization process by using the method of the present invention.
FIG. 3 is a graph of the convergence of the merit function in the optimization process of the present invention.
Detailed Description
The present invention is further illustrated by the following examples, which should not be construed as limiting the scope of the invention.
The method comprises the following specific steps:
step 1, setting a light source coding and decoding mode and a mask coding and decoding mode:
the light source is characterized by a two-dimensional matrix S with the size of NS×NS,NS51. Each element of the matrix represents a light source point, the value of the element is real number, the value of the element represents the intensity of the light source point, and the value range is [0,1 ]]. The position of the light source point in the matrix is denoted as S (m, n) and in frequency domain coordinates as S (f)s,gs) Wherein m is more than or equal to 1 and less than or equal to NS,1≤n≤NS, -1≤fs≤1,-1≤gsLess than or equal to 1. The incoherence and incoherence factors of the light source are expressed as sigmain0.4 and σout0.9, in order to ensure the symmetry of the light source, selecting the total N in the first quadrant and between the internal and external coherence factors in the frequency domain coordinate systemSc255 light source points are coded and their intensity values are recorded as
Figure RE-GDA0002875501180000072
The position in the matrix is recorded as
Figure RE-GDA0002875501180000073
Mixing J and PSAs a code for the light source.
When photoetching imaging calculation is carried out, J and P are calculatedSDecoding into a complete light source matrix S: firstly resetting S to be a full 0 matrix, and then sequentially carrying out i is more than or equal to 1 and less than or equal to N on the ith light source point in the light source codeScFrom PSRead its position index (m)i,ni) Reading the intensity value J from JiSelecting the position S (m) in Si,ni),S(NS+1-mi,ni), S(mi,NS+1-ni),S(NS+1-mi,NS+1-ni) Is set to have an element value of JiAnd finally obtaining the decoded complete light source matrix S.
The mask pattern is shown in FIG. 2 with a periodicity of 130nm and a feature size CD of 18 nm. The mask is characterized by a two-dimensional matrix M with the size of NM×NM,NM65. Each element of the matrix represents a pixel of the mask pattern, the value of the element is 0 or 1, the element of 0 represents that the pixel is a background region, and the element of 1 represents that the pixel is a pattern region. The position of a mask pixel in the matrix is denoted M (p, q), where 1 ≦ p ≦ NM, 1≤q≤NM. Will be the (N) th of the matrixM+1)/2 line and (N)MUsing +1)/2 columns as a reference, and selecting all matrix elements as pixels to be coded when the shape of the mask pattern is asymmetric; when the shape of the mask graph is symmetrical up and down, selecting elements of a first quadrant and a second quadrant as pixels to be coded; when the shape of the mask graph is symmetrical left and right, selecting elements of a first quadrant and a fourth quadrant as pixels to be coded; when the shape of the mask pattern is symmetrical up, down, left and right, the element of the first quadrant is selected as the pixel to be encoded. Number of pixels to be encodedNumber is NMc,NMc1089, their element values are recorded as
Figure RE-GDA0002875501180000081
The position in the matrix is recorded as
Figure RE-GDA0002875501180000082
Mixing T and PMAs the coding of the mask pattern.
When photoetching imaging calculation is carried out, T and P are addedMDecoding into a complete mask matrix M: firstly, M is reset to be a full 0 matrix, and then j is more than or equal to 1 and less than or equal to N is sequentially carried out on the jth pixel in the mask codeMcFrom PMRead its position index (p)j,qj) Reading the intensity value T from TjSelecting M (p) as the position of M in accordance with the symmetry of the original mask patternj,qj) And its corresponding symmetric element, setting a threshold value Tmr0.5, if TjGreater than a threshold value TmrIts element value is set to 1, otherwise its element value is set to 0, and finally the decoded complete mask matrix M is obtained.
Step 2, initializing a light source matrix, a mask matrix and photoresist parameters:
initializing a light source matrix to S according to light source parametersinit(ii) a Initializing a mask matrix to M according to a mask patterninitWill MinitAs a matrix M of target patternstarget(ii) a Initializing a photoresist threshold tr0.5, resist sensitivity α 200, resist development threshold trd. Initializing the normalized intensity I of the aerial imagenorm0.2541. The illumination wavelength of the lithography machine is λ 13.5nm, the light source polarization is y-direction polarization, the numerical aperture of the projection objective lens is NA 0.33, and the zoom magnification is Rreduction=4。
Step 3, mask defocus optimization:
since the pattern on the silicon wafer is shifted due to the oblique incidence of the illumination light source on the mask in the extreme ultraviolet lithography, it is necessary to correct the shift of the pattern by setting the mask to be out of focus. Setting mask defocus to be optimizedParameter, value of [ Delta z ]min,Δzmax]Real number in the range, Δ zmin=0nm,Δz max200 nm. Calculating photoresist images under different mask defocusing amounts through a thick mask model and a vector imaging model, taking Edge Placement Errors (EPE) of a photoresist graph and a target graph as evaluation functions, and optimizing by adopting a demarcation method to obtain the optimal mask defocusing amount delta zbestThe method comprises the following specific steps:
(1) initializing a maximum number of iterations Niter20, lower bound of the variable zl=ΔzminUpper bound of variable zr=ΔzmaxIntermediate values z of the upper and lower bounds of the variablemid=(zl+zr) /2, number of initialization iterations niter=0。
(2) Will zl,zmidAnd zrRespectively inputting the mask defocusing amount into a thick mask model and a vector imaging model to calculate a photoresist image, and then calculating corresponding evaluation function values fl,fmidAnd fr
(3) Calculating a temporary variable ztemp=(zl+zmid) And/2, inputting the mask defocus quantity into a thick mask model and a vector imaging model to calculate a photoresist image, and then calculating corresponding evaluation function values ftemp. If ftemp>fmidLet zl=ztempOn the contrary, z ismid=ztemp
(4) Calculating a temporary variable ztemp=(zr+zmid) And/2, inputting the mask defocus quantity into a thick mask model and a vector imaging model to calculate a photoresist image, and then calculating corresponding evaluation function values ftemp. If ftemp>fmidLet zr=ztempOn the contrary, z ismid=ztemp
(5) Let n beiter=niter+1, judgment of niter=NiterWhether or not this is true. If not, the optimization is continued to (2). If yes, the procedure is terminated, and the optimal mask defocus amount Delta z is obtainedbest=zmid
Step 4, calculating the thick mask parameters:
calibration of absorber transmittance t for mask background region in thick mask model by rigorous electromagnetic field simulationaTransmittance t of absorption layer in pattern regionbAnd a correction quantity delta of point pulse at the boundary position of the graph, which are respectively ta=0.0762+0.2283j,tb=0.397-0.918j,δ=4.307+0.7426j。
Incident angle theta of central point of illumination light sourceincAt 6 deg. azimuth angle
Figure RE-GDA0002875501180000091
Calculating the phase transmission factor phi of the incident light from the upper surface of the absorption layer to the middle position of the absorption layeraAnd a phase transmission factor phi of incident light from the middle position of the absorption layer to the lower surface of the absorption layerb。φaIs a complex number phibIs a complex matrix of size NO×NOIn which N isO31 is the total number of mask diffraction orders in either the x-direction or the y-direction. Calculating the reflection matrix R of the multilayer film of the extreme ultraviolet lithography mask according to the equivalent film layer methodmlSize is NO×NO
The optimal mask defocus amount Delta z obtained according to the step 3bestCalculating the mask defocus matrix phidSize is NO×NO
Step 5, optimizing a light source:
(1) according to the initialized mask matrix MinitAn illumination cross-coefficient (ICC) matrix of the lithography imaging system is calculated from the thick mask model and the vector imaging model.
(2) And constructing an evaluation function for optimizing the light source. Calculating a normalized aerial image I according to the normalized intensity of the aerial image, the light source matrix and the ICC matrixAThen, calculating a photoresist image I through a Sigmoid photoresist modelR=1/(1+exp(-α*(IA-tr) ) is used in the process) and finally according to the photoresist development threshold t)rdCalculating the photoresist pattern
Figure RE-GDA0002875501180000101
Photoresist pattern passingSize NM×NMThe binary matrix of 0-1 indicates that the element values have the same meaning as the element values of the mask matrix. Evaluation function F for light source optimizationSOIs shown as
Figure RE-GDA0002875501180000102
(3) Obtaining an initial light source code J according to the method in the step 1initThe method comprises the following steps of optimizing a light source through an SLPSO algorithm:
firstly, initializing the maximum iteration number NiterLet n be 200iter0. Initializing particle swarm size MpopDetermining the particle dimension D according to the number of the coded light source points (100)pop=NSc255, setting a social influence factor epsilon as k x Dpop/MpopWherein κ ═ 0.01 is a coefficient.
② setting the first particle position of the particle group as JinitThe remaining particle positions are randomly generated. Particle position through a size and JinitThe same row vector characterization with dimension DpopAnd each dimension has a value of 0-1, and the random numbers are generated uniformly. The position of the ith particle is denoted as piAll particles having a position pass size of Mpop×DpopMatrix P ofpopAnd (4) showing. The velocity of each particle is characterized by a row vector of the same magnitude as the particle position, and the velocities of all particles are initialized to 0. The velocity of the ith particle is denoted viAll particles have a velocity passage size of Mpop×DpopMatrix V ofpopAnd (4) showing.
Calculating learning probability factor P for all particles in the particle swarmLThe learning probability factor of the ith particle can be expressed as
Figure RE-GDA0002875501180000103
Fourthly, calculating the evaluation function value of all the particles in the particle group, wherein the evaluation function value of the ith particle is expressed as fi. However, the device is not suitable for use in a kitchenThen, the evaluation function values are sorted in descending order, and the evaluation function values of all the particles are sorted into
Figure RE-GDA0002875501180000104
And adjusting the speed of the particle swarm according to the corresponding sequence. The particle position at which the evaluation function value is minimum is defined as an optimum particle position pbest
Calculating the average position of the particle swarm
Figure RE-GDA0002875501180000105
M (M) after sortingpopThe particles are the optimal particles for the 1 st to M thpop-1 particle is updated in turn. For the ith particle (in this case 1 < i < M)pop-1) with a learning probability factor Pi L. Generating a random number p between 0 and 1rIf p isr>Pi LIf the particle is not changed, the velocity and position of the particle are updated:
the j dimension v of the particle velocityi,jThe updating is as follows:
vi,j=r1×vi,j+r2×(pk,j-pi,j)+r3×ε×(pmean,j-pi,j),
wherein r is1,r2,r3Is a random real number between 0 and 1, and k is from i +1 to Mpop-1, random integer. p is a radical ofk,jDimension j, p, representing the position of the kth particlei,jDimension j, p, representing the position of the ith particlemean,jRepresents the j-th dimension of the mean position of the particle population.
The j dimension p of the particle positioni,jThe updating is as follows:
pi,j=pi,j+vi,j
sixthly, niter=niter+1, judgment of niter=NiterWhether or not this is true. If not, the optimization is continued to the fourth step. If the above-mentioned conditions are met, the computer system can,terminating the procedure and outputting the optimal particle position pbest
To pbestDecoding according to the light source decoding method in the step 1 to obtain the optimal light source
Figure RE-GDA0002875501180000111
Step 6, optimizing a mask:
(1) according to the optimized light source matrix
Figure RE-GDA0002875501180000112
A cross-transmission coefficient (TCC) matrix of the lithography imaging system is calculated through a vector imaging model.
(2) An evaluation function for mask optimization is constructed. Calculating a mask diffraction spectrum matrix B with a size of N according to the mask matrix and the thick mask parametersO×NOThen combining TCC matrix and the normalized intensity of the space image to calculate a normalized space image IAThen, calculating a photoresist image I through a Sigmoid photoresist modelR=1/(1+exp(-α*(IA-tr) ) is used in the process) and finally according to the photoresist development threshold t)rdCalculating the photoresist pattern
Figure RE-GDA0002875501180000113
Evaluation function F for mask optimizationMOIs shown as
Figure RE-GDA0002875501180000114
The specific process of calculating the mask diffraction spectrum matrix B in the process of constructing the evaluation function is as follows:
generating an absorption layer transmittance matrix:
initializing the absorption layer transmittance matrix TabsorberIs of size NM×NMLet the position of the pixel in the background area with 0 element in the mask matrix M be denoted as Pbg(P, q), the position of the pixel of the graphic region having the element 1 is denoted as Pft(P, q) the position of the boundary pixel of the pattern detected by the difference operator is expressed as Pedge(p, q). At TabsorberIn the matrix, the position is equal to PbgThe element value of (p, q) is set to taWill be equal to PftThe element value of (p, q) is set to tbWill be equal to PedgeThe element value of (p, q) is set to δ. Will TabsorberFourier transform is carried out, and N is selected by taking 0 level as a centerO×NORange as TabsorberFourier coefficient matrix B ofT
Calculating a diffraction spectrum of the absorption layer after the first diffraction:
the diffraction spectrum of incident light propagating from the upper surface to the lower surface of the absorbing layer has a transit size NO×NOMatrix B ofabsorberExpressed, the calculation formula is:
Babsorber=φa⊙BT⊙φb
wherein |, indicates that the corresponding elements of the matrix are multiplied.
Calculating the diffraction spectrum after the reflection of the multilayer film:
the light diffracted by the absorption layer is reflected by the multilayer film, and the passing size of the diffraction spectrum after reflection is NO×NOMatrix B ofmlExpressed, the calculation formula is:
Bml=Babsorber⊙Rml
calculating the mask diffraction spectrum after the second diffraction of the absorption layer:
after the multilayer film is reflected, light in all directions is subjected to secondary diffraction through the absorption layer, and the light is mutually coherent and superposed to obtain a final mask diffraction spectrum B, wherein the calculation formula is as follows:
Figure RE-GDA0002875501180000121
wherein
Figure RE-GDA0002875501180000122
Representing a matrix convolution. Obtaining the optimal defocus amount Delta z of the mask according to optimizationbestAnd a mask defocus matrix phi obtained by pre-pressing calculation in advancedUpdating the mask diffraction spectrum:
B=B⊙Φd
(3) obtaining an initial mask code M according to the method in step 1initThe method comprises the following steps of optimizing a light source through an SLPSO algorithm: initializing a maximum number of iterations Niter400, let iteration number n iter0. Initializing particle swarm size MpopDetermining the particle dimension D according to the number of the coded light source points (100)pop=NMc1089, the social influence factor epsilon is set to κ × Dpop/MpopWherein κ ═ 0.01 is a coefficient. Then, optimizing according to the methods from step (5) to step (3), in particular, setting the position of the first particle to be the initial mask code M when the position of the particle swarm is initializedinitOptimizing the output pbestDecoding according to the mask decoding method in the step 1 to obtain an optimal mask
Figure RE-GDA0002875501180000123
In this example, the Pattern Error (PE) between the photoresist pattern and the target pattern before optimization was 603, and the Edge Placement Error (EPE) was 2.142 nm. After the light source mask is optimized by the method, the Pattern Error (PE) of the photoresist pattern and the target pattern is 35 and is reduced by 94.2 percent, and the Edge Placement Error (EPE) is 0.2037nm and is reduced by 90.49 percent. The entire optimization process takes 303.8 s. The method has the advantages of high optimization efficiency and high optimization speed.
The above description is only one specific embodiment of the present invention, and the embodiment is only used to illustrate the technical solution of the present invention and not to limit the present invention. The technical solutions available to those skilled in the art through logical analysis, reasoning or limited experiments according to the concepts of the present invention are all within the scope of the present invention.

Claims (2)

1. An extreme ultraviolet lithography light source mask optimization method is characterized by comprising the following steps:
step 1, encoding and decoding of light source and mask patterns:
(1) encoding and decoding of light sources
The light source is characterized by a two-dimensional matrix S with the size of NS×NS,NSEach element of the matrix represents a light source point, the value of the element is real number and represents the intensity of the light source point, and the value range is [0,1 ]](ii) a The position of the light source point in the matrix is denoted as S (m, n) and in frequency domain coordinates as S (f)s,gs) Wherein m is more than or equal to 1 and less than or equal to NS,1≤n≤NS,-1≤fs≤1,-1≤gs≤1;
The incoherence and incoherence factors of the light source are expressed as sigmainAnd σoutSelecting N in the first quadrant of the frequency domain coordinate system and between the internal and external coherence factorsScThe light source points are coded, and the intensity value is recorded as
Figure RE-FDA0002875501170000011
The position in the matrix is recorded as
Figure RE-FDA0002875501170000012
Mixing J and PSA code as a light source;
when photoetching imaging calculation is carried out, J and P are calculatedSAnd (3) decoding: firstly resetting a light source matrix S to be a full 0 matrix, and then sequentially carrying out i is more than or equal to 1 and less than or equal to N for the ith light source point in the light source codeScFrom PSRead its position index (m)i,ni) Reading the intensity value J from JiSelecting the position S (m) in Si,ni),S(NS+1-mi,ni),S(mi,NS+1-ni),S(NS+1-mi,NS+1-ni) Is set to have an element value of JiFinally, obtaining a decoded complete light source matrix S;
(2) encoding and decoding of masks
The mask is characterized by a two-dimensional matrix M with the size of NM×NM,NMIs odd, each element of the matrix represents a pixel of the mask pattern, the value of the element is 0 or 1, the element is 0 representing that the pixel is a background area, the element is a background areaA value of 1 indicates that the pixel is a pattern region, and the position of the mask pixel in the matrix is represented as M (p, q), where 1. ltoreq. p.ltoreq.NM,1≤q≤NM
Will be the (N) th of the matrixM+1)/2 line and (N)MUsing +1)/2 columns as a reference, and selecting all matrix elements as pixels to be coded when the shape of the mask pattern is asymmetric; when the shape of the mask graph is symmetrical up, down, left and right, the element of the first quadrant is selected as the pixel to be coded; when the shape of the mask graph is symmetrical up and down, selecting elements of a first quadrant and a second quadrant as pixels to be coded; when the shape of the mask graph is symmetrical left and right, selecting elements of a first quadrant and a fourth quadrant as pixels to be coded;
the number of pixels to be coded is recorded as NMcThe element value is recorded as
Figure RE-FDA0002875501170000013
The position in the matrix is recorded as
Figure RE-FDA0002875501170000021
Mixing T and PMA code as a mask pattern;
when photoetching imaging calculation is carried out, T and P are addedMDecoding into a complete mask matrix M: firstly resetting a mask matrix M to be a full 0 matrix, and then sequentially carrying out j is more than or equal to 1 and less than or equal to N on the jth pixel in the mask codeMcFrom PMRead its position index (p)j,qj) Reading the intensity value T from TjSelecting M (p) as the position of M in accordance with the symmetry of the original mask patternj,qj) And its corresponding symmetric element, if TjGreater than a threshold value TmrSetting the element value to 1, otherwise, setting the element value to 0, and finally obtaining a decoded complete mask matrix M;
step 2, initializing a light source matrix, a mask matrix and photoresist parameters:
initializing a light source matrix J according to light source parametersinitMask pattern initialization mask matrix MinitWill MinitAs a matrix M of target patternstargetInitializing the photoresist threshold trSensitivity of photoresist α, developing threshold t of photoresistrdInitializing the normalized intensity I of the aerial imagenorm
Step 3, mask defocus optimization:
setting the mask defocusing amount as a parameter to be optimized, and taking the value as [ delta z [ ]min,Δzmax]Calculating photoresist images under different mask defocusing amounts by using a thick mask model and a vector imaging model according to real numbers in the range, taking Edge Placement Errors (EPE) of the photoresist images and a target image as evaluation functions, and optimizing by adopting a binary search method to obtain the optimal mask defocusing amount delta zbest
Step 4, calculating the thick mask parameters:
calibration of absorber transmittance t for mask background region in thick mask model by rigorous electromagnetic field simulationaTransmittance t of absorption layer in pattern regionbAnd a point pulse correction δ at the pattern boundary position, all of which are complex;
according to the incident angle theta of the central point of the illumination light sourceincAnd azimuth angle
Figure RE-FDA0002875501170000022
Calculating the phase transmission factor phi of the incident light from the upper surface of the absorption layer to the middle position of the absorption layeraAnd a phase transmission factor phi of incident light from the middle position of the absorption layer to the lower surface of the absorption layerb,φaIs a complex number phibIs a complex matrix of size NO×NOIn which N isOIs the total number of mask diffraction orders in the x-direction or y-direction;
calculating the reflection matrix R of the multilayer film of the extreme ultraviolet lithography mask according to the equivalent film layer methodmlSize is NO×NOAccording to the optimal mask defocus amount Delta z obtained in step 3bestCalculating the mask defocus matrix phidSize is NO×NO
Step 5, optimizing a light source:
(1) according to the initialized mask matrix MinitThrough a thick mask dieThe type and vector imaging models calculate an Illumination Cross Coefficient (ICC) matrix of the photoetching imaging system;
(2) constructing an evaluation function of light source optimization: calculating a normalized aerial image I according to the normalized intensity of the aerial image, the light source matrix and the ICC matrixAThen, calculating a photoresist image I through a Sigmoid photoresist modelR=1/(1+exp(-α*(IA-tr) ) is used in the process) and finally according to the photoresist development threshold t)rdCalculating the photoresist pattern
Figure RE-FDA0002875501170000031
The photoresist pattern has a pass size of NM×NMIs represented by a binary matrix of 0,1, the element values having the same meaning as those of the mask matrix, and the evaluation function F for light source optimizationSOIs shown as
Figure RE-FDA0002875501170000032
(3) Carrying out light source coding according to the method in the step 1, and obtaining an optimal light source matrix through optimization of an SLPSO algorithm
Figure RE-FDA0002875501170000033
Step 6, mask optimization:
(1) according to the optimized light source matrix
Figure RE-FDA0002875501170000034
Calculating a cross transmission coefficient (TCC) matrix of the photoetching imaging system through a vector imaging model;
(2) constructing an evaluation function of mask optimization: calculating a mask diffraction spectrum matrix B of size N from the mask matrix and the thick mask parametersO×NOThen combining the TCC matrix and the normalized intensity of the space image to calculate a normalized space image IAThen, calculating a photoresist image I through a Sigmoid photoresist modelR=1/(1+exp(-α*(IA-tr) ) is used in the process) and finally according to the photoresist development threshold t)rdCalculating the photoresist pattern
Figure RE-FDA0002875501170000035
Evaluation function F for mask optimizationMOIs shown as
Figure RE-FDA0002875501170000036
(3) Carrying out mask coding according to the method in the step 1, and obtaining an optimal mask matrix through optimization of an SLPSO algorithm
Figure RE-FDA0002875501170000037
2. The euv lithography light source mask optimization method according to claim 1, wherein in the mask optimization process, a mask diffraction spectrum matrix B is calculated according to the mask matrix M and the thick mask parameters, specifically comprising the steps of:
step 1, generating an absorption layer transmittance matrix:
initializing the absorption layer transmittance matrix TabsorberIs of size NM×NMLet the position of the pixel in the background area with 0 element in the mask matrix M be denoted as Pbg(P, q), the position of the pixel of the graphic region having the element 1 is denoted as Pft(P, q) the position of the boundary pixel of the pattern detected by the difference operator is expressed as Pedge(p, q); at TabsorberIn the matrix, the position is equal to PbgThe element value of (p, q) is set to taWill be equal to PftThe element value of (p, q) is set to tbWill be equal to PedgeSetting the element value of (p, q) to delta, and setting TabsorberFourier transform is carried out, and N is selected by taking 0 level as a centerO×NORange as TabsorberFourier coefficient matrix B ofT
Step 2, calculating a diffraction spectrum of the absorption layer after the first diffraction:
the diffraction spectrum of incident light propagating from the upper surface to the lower surface of the absorbing layer has a transit size NO×NOMatrix B ofabsorberExpressed by the formula of:
Babsorber=φa⊙BT⊙φb
Wherein [ ] indicates a multiplication of corresponding elements of the matrix;
step 3, calculating the diffraction spectrum after the reflection of the multilayer film:
the light diffracted by the absorption layer is reflected by the multilayer film, and the passing size of the diffraction spectrum after reflection is NO×NOMatrix B ofmlExpressed, the calculation formula is:
Bml=Babsorber⊙Rml
step 4, calculating a mask diffraction spectrum of the absorption layer after the second diffraction:
after the multilayer film is reflected, light in all directions is subjected to secondary diffraction through the absorption layer, and the light is mutually coherent and superposed to obtain a final mask diffraction spectrum B, wherein the calculation formula is as follows:
Figure RE-FDA0002875501170000041
wherein
Figure RE-FDA0002875501170000042
Representing a matrix convolution. Obtaining the optimal defocus amount Delta z of the mask according to optimizationbestAnd a mask defocus matrix phi obtained by pre-pressing calculation in advancedUpdating the mask diffraction spectrum:
B=B⊙Φd
CN202011279713.8A 2020-11-16 2020-11-16 Extreme ultraviolet lithography light source mask optimization method Active CN112394615B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011279713.8A CN112394615B (en) 2020-11-16 2020-11-16 Extreme ultraviolet lithography light source mask optimization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011279713.8A CN112394615B (en) 2020-11-16 2020-11-16 Extreme ultraviolet lithography light source mask optimization method

Publications (2)

Publication Number Publication Date
CN112394615A true CN112394615A (en) 2021-02-23
CN112394615B CN112394615B (en) 2024-02-09

Family

ID=74599477

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011279713.8A Active CN112394615B (en) 2020-11-16 2020-11-16 Extreme ultraviolet lithography light source mask optimization method

Country Status (1)

Country Link
CN (1) CN112394615B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113311669A (en) * 2021-04-14 2021-08-27 长春理工大学 Photoetching image obtaining method capable of improving imaging quality
CN113589643A (en) * 2021-07-09 2021-11-02 中国科学院上海光学精密机械研究所 Curve type reverse photoetching method based on modulation block superposition synthesis mask pattern
CN115469512A (en) * 2022-09-13 2022-12-13 武汉宇微光学软件有限公司 Light source calibration method and system in light source mask joint optimization

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070038972A1 (en) * 2005-08-09 2007-02-15 Christoph Noelscher Method for optimizing a photolithographic mask
CN104155852A (en) * 2014-08-26 2014-11-19 中国科学院上海光学精密机械研究所 Light source optimization method of photolithography machine
CN104635434A (en) * 2015-01-29 2015-05-20 北京理工大学 Extreme ultraviolet lithography mask optimization method
CN104914684A (en) * 2015-06-24 2015-09-16 北京理工大学 Extreme ultraviolet lithography light source-mask joint optimization method
CN106125511A (en) * 2016-06-03 2016-11-16 北京理工大学 Low error suseptibility multiple target source mask optimization method based on vector imaging model
CN106164777A (en) * 2014-04-14 2016-11-23 Asml荷兰有限公司 The Optimizing Flow of photoetching process
CN109426083A (en) * 2017-08-24 2019-03-05 中芯国际集成电路制造(上海)有限公司 The optimization method and its optimization system and photolithography method of photoetching process
JP2019117418A (en) * 2019-04-25 2019-07-18 株式会社ニコン Apparatus for controlling lithography device, program for the apparatus, exposure system, and method of manufacturing device
CN110554580A (en) * 2019-09-27 2019-12-10 中国科学院上海光学精密机械研究所 multi-parameter joint optimization method for photoetching machine
CN111781804A (en) * 2020-07-16 2020-10-16 中国科学院上海光学精密机械研究所 Light source mask optimization method based on covariance matrix adaptive evolution strategy algorithm
CN111897188A (en) * 2020-06-28 2020-11-06 北京理工大学 Method for detecting polarization aberration of ultra-large numerical aperture strict vector imaging system

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070038972A1 (en) * 2005-08-09 2007-02-15 Christoph Noelscher Method for optimizing a photolithographic mask
CN106164777A (en) * 2014-04-14 2016-11-23 Asml荷兰有限公司 The Optimizing Flow of photoetching process
CN104155852A (en) * 2014-08-26 2014-11-19 中国科学院上海光学精密机械研究所 Light source optimization method of photolithography machine
CN104635434A (en) * 2015-01-29 2015-05-20 北京理工大学 Extreme ultraviolet lithography mask optimization method
CN104914684A (en) * 2015-06-24 2015-09-16 北京理工大学 Extreme ultraviolet lithography light source-mask joint optimization method
CN106125511A (en) * 2016-06-03 2016-11-16 北京理工大学 Low error suseptibility multiple target source mask optimization method based on vector imaging model
CN109426083A (en) * 2017-08-24 2019-03-05 中芯国际集成电路制造(上海)有限公司 The optimization method and its optimization system and photolithography method of photoetching process
JP2019117418A (en) * 2019-04-25 2019-07-18 株式会社ニコン Apparatus for controlling lithography device, program for the apparatus, exposure system, and method of manufacturing device
CN110554580A (en) * 2019-09-27 2019-12-10 中国科学院上海光学精密机械研究所 multi-parameter joint optimization method for photoetching machine
CN111897188A (en) * 2020-06-28 2020-11-06 北京理工大学 Method for detecting polarization aberration of ultra-large numerical aperture strict vector imaging system
CN111781804A (en) * 2020-07-16 2020-10-16 中国科学院上海光学精密机械研究所 Light source mask optimization method based on covariance matrix adaptive evolution strategy algorithm

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HUNG-FEI KUO等: "Forming Freeform Source Shapes by Utilizing Particle Swarm Optimization to Enhance Resolution in Extreme UV Nanolithography", IEEE TRANSACTIONS ON NANOTECHNOLOGY, vol. 14, no. 2, pages 322 - 329, XP011574876, DOI: 10.1109/TNANO.2015.2393916 *
LEI WANG等: "Pixelated source optimization for optical lithography via particle swarm optimization", JOURNAL OF MICRO-NANOLITHOGRAPHY MEMS AND MOEMS, vol. 15, no. 1, pages 1 - 11 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113311669A (en) * 2021-04-14 2021-08-27 长春理工大学 Photoetching image obtaining method capable of improving imaging quality
CN113311669B (en) * 2021-04-14 2023-02-10 长春理工大学 Photoetching image obtaining method capable of improving imaging quality
CN113589643A (en) * 2021-07-09 2021-11-02 中国科学院上海光学精密机械研究所 Curve type reverse photoetching method based on modulation block superposition synthesis mask pattern
CN113589643B (en) * 2021-07-09 2024-05-17 中国科学院上海光学精密机械研究所 Curve type reverse photoetching method for synthesizing mask pattern based on superposition of modulation blocks
CN115469512A (en) * 2022-09-13 2022-12-13 武汉宇微光学软件有限公司 Light source calibration method and system in light source mask joint optimization

Also Published As

Publication number Publication date
CN112394615B (en) 2024-02-09

Similar Documents

Publication Publication Date Title
TWI820070B (en) Inspection of reticles using machine learning
KR102459381B1 (en) A method for training a machine learning model for computational lithography.
CN112394615A (en) Extreme ultraviolet lithography light source mask optimization method
JP7362744B2 (en) Method and device for selecting layout patterns
TW202113633A (en) Method for training a deep learning model of a patterning process and related computer program product
CN110554580B (en) Multi-parameter joint optimization method for photoetching machine
US20060248495A1 (en) Method and apparatus of model-based photomask synthesis
CN114766012A (en) Method and system for predicting process information using parameterized models
US11989873B2 (en) Stochastic contour prediction system, method of providing the stochastic contour prediction system, and method of providing EUV mask using the stochastic contour prediction system
CN110244523B (en) Integrated photoetching method and photoetching system
Awad et al. Accurate prediction of EUV lithographic images and 3D mask effects using generative networks
Xiao et al. A novel optical proximity correction (OPC) system based on deep learning method for the extreme ultraviolet (EUV) lithography
WO2023117250A1 (en) Method and apparatus to determine overlay
CN116982006A (en) Method and system for predicting process information using parameterized models
CN112236722B (en) Method for determining electromagnetic fields associated with computing a lithographic mask model
TWI818397B (en) Modular autoencoder model for manufacturing process parameter estimation
EP3828632A1 (en) Method and system for predicting electric field images with a parameterized model
Yu et al. Mask 3D model based on complex-valued convolution neural network for EUV lithography
Woldeamanual et al. Accurate prediction of EUV lithographic images and 3D mask effects using generative networks
TW202240316A (en) Selecting patterns based on their representation characterized by lithographic apparatus or process characteristic
Moreau et al. A deep learning workflow to generate free-form masks for grayscale lithography
WO2023186446A1 (en) Methods related to an autoencoder model or similar for manufacturing process parameter estimation
WO2024017808A1 (en) Deep learning models for determining mask designs associated with semiconductor manufacturing
CN116802647A (en) Modular automatic encoder model for manufacturing process parameter estimation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant