CN115469555B - Space image prediction and image quality optimization method for sensor chip projection lithography machine - Google Patents
Space image prediction and image quality optimization method for sensor chip projection lithography machine Download PDFInfo
- Publication number
- CN115469555B CN115469555B CN202211417846.6A CN202211417846A CN115469555B CN 115469555 B CN115469555 B CN 115469555B CN 202211417846 A CN202211417846 A CN 202211417846A CN 115469555 B CN115469555 B CN 115469555B
- Authority
- CN
- China
- Prior art keywords
- light source
- mask
- matrix
- model
- representing
- 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.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000001459 lithography Methods 0.000 title claims abstract description 33
- 238000003384 imaging method Methods 0.000 claims abstract description 74
- 238000009826 distribution Methods 0.000 claims abstract description 29
- 238000001259 photo etching Methods 0.000 claims abstract description 28
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 25
- 238000001228 spectrum Methods 0.000 claims abstract description 22
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 3
- 239000011159 matrix material Substances 0.000 claims description 62
- 230000001427 coherent effect Effects 0.000 claims description 25
- 238000005286 illumination Methods 0.000 claims description 21
- 238000013178 mathematical model Methods 0.000 claims description 16
- 210000001747 pupil Anatomy 0.000 claims description 16
- 230000003595 spectral effect Effects 0.000 claims description 16
- 230000003287 optical effect Effects 0.000 claims description 12
- 229920002120 photoresistant polymer Polymers 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 7
- 239000004576 sand Substances 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 230000006978 adaptation Effects 0.000 claims 1
- 230000002068 genetic effect Effects 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 239000013598 vector Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000002265 prevention Effects 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000000206 photolithography Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
Abstract
The invention discloses a space image prediction and image quality optimization method for a sensor chip projection lithography machine, which is a hybrid method of tabu search and genetic algorithm. The lithography system space image imaging model is based on an Abbe imaging method, and the calculation of space image intensity distribution is completed by extracting mask effective diffraction spectrum information. The light source S consists of four symmetrical subregions, and effective units of the subregions are extracted to be used as variables for optimizing the light source. The mask M is composed of symmetrical regular patterns, and cells near the edges of the sub-region patterns are extracted to be used as variables for mask optimization. Using the figure error and the edge placement error as an objective function F 1 And F 2 And performing Gaussian filtering operation on the optimized light source, and performing binary processing on the optimized gray mask by a threshold method. The invention provides a method for a photoetching system space image imaging model and provides a rapid light source mask optimization method, which provides effective help for improving the working efficiency of a sensor chip projection photoetching machine.
Description
Technical Field
The invention belongs to the field of projection lithography, and particularly relates to a space image prediction and image quality optimization method for a sensor chip projection lithography machine.
Background
Projection lithography is the dominant technology in the manufacture of very large scale integrated circuits. At present, the performance of very large scale integrated circuits is improved and miniaturized by reducing the feature size of a single transistor and increasing the number of transistors in the same area. According to Rayleigh criterion, the exposure wavelength is shortened, the numerical aperture of the objective lens is increased, and the photoetching imaging resolution ratio can be effectively improved. The continuous reduction of the critical dimension leads to the aggravation of the optical diffraction effect, so that in the imaging process, a more obvious optical proximity effect is introduced, and the quality of the photoetching imaging is reduced. Therefore, how to improve the lithography imaging performance has become a problem to be solved urgently.
The key of the lithography reverse optimization model based on the pixilated light source mask optimization lies in a forward imaging model and a reverse optimization model of lithography. The forward imaging model mainly comprises a light source model, a mask model and a pupil model. According to the Abbe imaging theory and the frequency spectrum relation among the Abbe imaging theory and the Abbe imaging theory, a forward forming model is established to obtain the intensity distribution of the photoetching space image. In the inverse lithography optimization model, in order to meet the optimization rule of the iterative algorithm, the sum of absolute values of differences of each element between the photoresist pattern and the ideal pattern is used as a cost function, namely a pattern error. Wherein, the photoresist pattern can approximately represent the photoresist effect by a sigmoid function. In the light source optimization model, light source effective units are marked and used as optimization variables of the iterative model. According to the characteristic that the illumination light source is symmetrical about the optical axis, only a quarter of optimization variables are needed, the complexity of an optimization model is reduced, and the optimization efficiency is improved. In the mask optimization model, an edge optimization strategy is adopted, namely, a unit near the edge of the feature pattern is used as an optimization variable, and the variable value is continuously updated through iteration.
In the inverse optimization model, the obtained effect is different according to different optimization algorithms. However, for the same feature pattern, the optimized light source intensity distribution has similar trends, and the optimization result of the mask pattern is different according to different optimization strategies. Prior art 1 (Yao Pen, jinyu Zhang, yan Wang, and Zhiping Yu, "Gradient-Based Source and Mask Optimization in Optical Lithography," IEEE trans. On Image Process. 20 (10), 2856-2864 (2011)) proposes a light Source Mask Optimization method of a Gradient descent method. The method adopts an Abbe imaging model to complete light source optimization, and adopts the sum of coherent systems (SOCS) to complete mask optimization. In the light source mask optimization model, the value range of an optimization variable is constrained by a cosine function, and meanwhile, a threshold function is adopted to carry out binary operation according to the value of a mask variable. Prior art 2 (x, ma, c, han, y, li, l, dong, and g.r. Area, "pixel source and mask optimization for imaging," j. Opt. Soc. Am. A30 (1), 112 (2013)), proposes a lithography light source mask optimization model based on a gradient method of pixel characterization. The model establishes a photoetching imaging model by a vector Abbe imaging method, and provides a synchronous light source mask optimization model and a sequence light source mask optimization model to optimize light source intensity distribution and mask pattern layout. In prior art 3 (c, yang, s, li, and x, wang, "Efficient source mask optimization using multiple source reconstruction," j, micro/nanolith, MEMS MOEMS 13 (4), 043001 (2014)), a genetic algorithm-based light source mask optimization method is proposed. The method adopts a polar coordinate mode to mark the distribution of the variable of the effective unit of the light source.
In the prior art 1 and the prior art 2, a model for optimizing a light source mask is established by a gradient-based method, although the convergence rate of the method is high, the derivation process of the cost function is complex, and the optimization model is prone to be in a locally optimal condition. In the prior art 3, a genetic algorithm is used as a kernel for optimizing the intensity distribution of the light source and the layout of the mask pattern, and the genetic algorithm is used as a heuristic algorithm for global optimization, so that the method has the advantages of simpler optimization structure, high search speed and the like. However, the photoetching imaging model and the photoetching model are complex, so that the dimension of the variable matrix is large, the method is easy to generate a state of premature convergence, and the convergence efficiency of the algorithm is reduced.
Disclosure of Invention
In order to solve the technical problem, the invention provides a space image prediction and image quality optimization method for a sensor chip projection lithography machine. The method realizes the space image imaging process of the photoetching system by adopting a mode of extracting effective characteristic pattern frequency spectrum information, reduces the calculation complexity of the process, and simultaneously improves the accuracy of an imaging model. The global light source mask optimization method is suitable for various objective functions F such as graphic errors, edge placement errors and the like, and effectively improves the imaging performance of a photoetching system while ensuring low complexity of a mask optimization result.
In order to achieve the purpose, the technical solution adopted by the invention is as follows:
a space image prediction and image quality optimization method for a sensor chip projection lithography machine is a hybrid method of tabu search-genetic algorithm, and comprises the following steps:
step 1: initializing space image imaging model parameters of photoetching system, inputting light sourceSAnd mask patternM;
And 2, step: coding the light source and the mask effective unit as target function variables;
and step 3: randomly generating a group variable whole matrix to be optimized according to a heuristic optimization algorithm structureP(ii) a Therefore, in the light source mask optimization model, a light source and mask initial population-light source variable population matrix is randomly generatedP S And mask variable population matrixP M The light source variable group matrixP S And mask variable population matrixP M Are respectively asP S =[S 1 ,S 2 ,S 3 ,…,S P ],P M =[M 1 ,M 2 ,M 3 ,…,M P ];
And 4, step 4: according to the initial light source variable matrixP S Calculating the intensity distribution of the space image generated by each individual through the space image imaging model of the photoetching systemI S The spatial image prediction under different imaging conditions is realized; and according to an objective function F 1 And F 2 Calculating adaptive value, and selecting current optimum initial light source variable individual by comparisonAnd an adapted value->(ii) a Wherein it is present>Expressed in the space image imaging model of the lithography system, the group matrix is formed according to different light source variablesP S Of (2)S P The resulting aerial image intensity distribution matrixI S I.e. based on>,/>;
And 5: the current optimal initial light source variable obtained in the step 4 is taken as an individualAs an input condition of a tabu search-genetic algorithm optimization model, iteratively updating the light source variable individuals until iteration is stopped;
and 6: according to the light source variable individuals after iteration is stopped, matrix inversion mirror image and Gaussian filtering manipulation are executed, and the appearance of the light source is recovered;
and 7: using the obtained light source as the light source condition of the mask optimization model; population matrix according to mask variablesP M Calculating and recording current best initial mask variable individualAnd an adapted value->;
And 8: the current best initial mask variable obtained in the step 7 is taken as an individualPerforming iterative update on mask variable individuals until iteration stops as an input condition of a tabu search-genetic algorithm optimization model;
and step 9: outputting the optimal light source intensity distribution and the mask pattern layout.
Further, in step 1, the mask pattern M and the input light source S are gridded, effective characteristic information of the mask pattern is extracted by an effective spectrum extraction method, and an aerial image generated by the lithography optical system is represented by the following mathematical model in combination with an Abbe imaging method:
wherein,Irepresenting the spatial image intensity distribution, i.e. the result of partially coherent imaging;N ' S representing the number of effective point light sources in the pixelated light source;CCI i the coherent imaging results produced by a single point source are represented by the following mathematical model:
wherein, (ii) (x i ,y i ) Representing a coherent image plane spatial coordinate system; (f,g) Representing a pupil plane spectral coordinate system; (f',g') denotes the mask spectral coordinate system;Hrepresents the optical transfer function, i.e. the pupil function, of the imaging system;H i representing the position of the light source according to points: (S x ,S y ) A pupil after translation;Mrepresenting a mask spectrum;
the coherent imaging model was discretized and represented by the following mathematical model:
wherein,representing a discretized coherent image plane spatial coordinate system; (j,k) Representing a discretized pupil plane spectral coordinate system; (j',k') denotes the discretized mask spectral coordinate system;N ext representing the number of extracted valid spectral samples,;
according to Abbe's method, aerial image obtained by partially coherent imaging processI ext Is represented by the following mathematical model:
wherein, in the followingm, nAnd establishing a photoetching system illumination light source model for the index-mode coordinate system.S(m,n) Is shown in the light source model coordinate system at (m,n) Point light source at the position, the size of the illumination light source model matrix of the photoetching system isN S ×N S ,m,n=1,2,3,…,N S 。
Further, in the step 4, the sum of absolute values of the differences between the photoresist pattern matrix and each unit of the ideal pattern is defined as a pattern error, and the pattern error is taken as an objective function F 1 (ii) a Marking the areas of the inner and outer pixel points near the boundary of the mask pattern, and calculating the edge placement error, taking the edge prevention error as an objective function F 2 (ii) a Objective function F 1 And F 2 Are respectively provided withIs represented as the following mathematical model:
wherein,RP(x i ,y i ) A photoresist pattern is represented and,M * (x i ,y i ) Representing an ideal graph;ICC Edge representing a matrix formed by the extracted edge pixel points;S ext representing the extracted light source matrix;M ' ext representing a matrix formed by extracting ideal graph edge pixel points according to the index positions;minimizethe method is characterized in that in a tabu search-genetic iteration optimization-based model, the numerical value of an objective function is guaranteed to be minimum when iteration stops.
Compared with the prior art, the invention has the following advantages:
1. the invention provides a hybrid genetic algorithm, improves the convergence efficiency of the genetic algorithm and ensures that a photoetching system has better imaging performance.
2. According to the optical characteristics of the photoetching illumination system, the light source appearance is divided into four identical parts, effective light source variables are encoded, and the light source intensity distribution gradient is improved by a Gaussian filtering method.
3. The optimization method provided by the invention is suitable for various objective functions, such as graphic errors, edge prevention errors and the like.
Drawings
FIG. 1 is a schematic flow chart of a method for spatial image prediction and image quality optimization for a sensor chip projection lithography machine according to the present invention;
FIG. 2 is a schematic diagram of a mask used in the present invention;
FIG. 3 is a schematic diagram of a mask effective feature spectrum extraction method according to the present invention;
FIG. 4 is a schematic view of a lithography system aerial image imaging method employed in the present invention;
FIG. 5 is a schematic diagram of a light source encoding and Gaussian filtering method according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
In order to improve the imaging performance of the aerial image of the lithography system, the invention discloses an aerial image prediction and image quality optimization method for a sensor chip projection lithography machine through the following embodiments.
As shown in FIG. 1, the space image prediction and image quality optimization method for sensor chip projection lithography machine of the present invention is a hybrid tabu search-genetic algorithm, wherein the space image prediction is performed in step 4, and the objective function is F 1 And F 2 The method comprises the following steps:
step 1: initializing space image imaging model parameters of photoetching system, inputting light sourceSAnd mask patternM;
Step 2: taking the coded light source and the mask effective unit as target function variables;
and step 3: according to the heuristic optimization algorithm structure, all matrixes of group variables to be optimized need to be randomly generatedP. Therefore, in the light source mask optimization model, a light source and mask initial population-light source variable population matrix is randomly generatedP S And mask variable population momentsMatrix ofP M The light source variable group matrixP S And mask variable population matrixP M Are respectively asP S =[S 1 ,S 2 ,S 3 ,…,S P ],P M =[M 1 ,M 2 ,M 3 ,…,M P ];
And 4, step 4: according to the initial light source variable matrixP S Calculating the intensity distribution of the space image generated by each individual through the space image imaging model of the photoetching systemI S And space image prediction under different imaging conditions is realized. And according to an objective function F 1 And F 2 Calculating adaptive value, and selecting current optimum initial light source variable individual by comparisonAnd an adaptive value>(ii) a Wherein it is present>Expressed in the space image imaging model of the lithography system, the group matrix is formed according to different light source variablesP S Of (2)S P The resulting aerial image intensity distribution matrixI S I.e. is->,/>;
And 5: the current optimal initial light source variable obtained in the step 4 is taken as an individualAs input conditions of a tabu search-genetic algorithm optimization modelIteratively updating the light source variable individual until iteration is stopped;
step 6: according to the light source variable individuals after iteration is stopped, matrix inversion mirror image and Gaussian filtering manipulation are executed, and the appearance of the light source is recovered;
and 7: the light source obtained above was used as the light source condition for the mask optimization model. Population matrix according to mask variablesP M Calculating and recording current best initial mask variable individualAnd an adapted value->;
And 8: the current best initial mask variable obtained in the step 7 is taken as an individualPerforming iterative update on mask variable individuals until iteration stops as an input condition of a tabu search-genetic algorithm optimization model;
and step 9: outputting the aerial image intensity distribution, the optimal light source intensity distribution and the mask pattern layout.
In this embodiment, the number of horizontal sampling points of the aerial image imaging model of the lithography system is setN=521, physical size of pixel point is 5.625nm×5.625nm. Then, the number of horizontal sampling points of the light sourceN S =41, coherence factor of annular illumination sourceσ inner =0.68 andσ outer =0.95. The mask feature pattern has a line width of 45 nm and a spacing of 45 nm as shown in fig. 2. The lithography system aerial image imaging model adopts an effective spectrum extraction method, as shown in FIG. 3. Wherein, step 301 represents the raster array pattern as the input pattern of the schematic diagram. Step 302 represents representing the feature information in the form of a frequency spectrum by using a fast fourier transform operation for the feature pattern, step 303 represents the feature information in the form of a three-dimensional frequency spectrum representation, step 304 represents an effective spectrum extraction operation, i.e. a spectrum information for the mask pattern according to a cut-off frequency rangeExtracting, wherein step 305 shows that in a lithography system space image imaging model, an effective three-dimensional spectrum of a characteristic graph can be obtained according to the relation between an effective spectrum and a cut-off spectrum which participate in imaging, and step 306 shows that the effective spectrum of the characteristic graph is represented by a two-dimensional distribution graph, wherein DC represents the radiusR f The circular region of (a), which is the effective cut-off spectral range.
As shown in the schematic diagram of the process of imaging the lithography system aerial image in fig. 4, assuming that the light source in the model of imaging the lithography system aerial image is the annular illumination light source S, step 401 represents representing the light source as the annular illumination light source S by using the discretization matrix, i.e. the intensity distribution of the illumination light source is represented by pixelization. Step 402 represents the extraction of the effective light source points and the matrix transformation operation. The matrix conversion operation converts the effective light source matrix from a multi-dimensional matrix to a column vectorS ext . Step 403 represents representing the active light sources by a one-dimensional column vector. According to the coherent imaging theory, in the discretized annular illumination light source S, each point light source can generate a coherent imageCCI i . Step 404 represents each partially coherent image to be extractedCCI i The edge pixels form an illumination cross coefficient matrixICC ext Each column of the matrix represents the intensity distribution of the coherent image produced by each point source. Step 405 represents a matrix transformation operation, which represents the transformation of the result of the product of the multi-dimensional matrix into a column vector and an illumination cross coefficient matrix, and the intensity of the spatial image of the lithography system can be obtainedI ext . Step 406 represents segmenting the aerial image intensityI ext And (4) discretizing expression. According to the size ofN M ×N M The mask pattern of (2) can be known as the size of the aerial image intensity distribution matrix obtained by the aerial image imaging model of the lithography system isN M ×N M Step 408 is required to operate. Said step 408 represents the acquisition of the complete aerial image intensity distribution by interpolationI. Step 409 represents the resulting standard aerial image intensity distributionI. Step 407 and step 410 represent the aerial image intensity portions to be extractedI ext AndIthe matrix of (2) is represented by local information amplification. According to the Abbe imaging method, the imaging model can be represented by the following mathematical model:
wherein,Irepresenting the spatial image intensity distribution, i.e. the result of partially coherent imaging.N ' S Representing the number of effective point sources in the pixelated light source.CCI i Represents the coherent imaging results produced by a single point source, which can be represented by the following mathematical model:
wherein, (ii) (x i ,y i ) Representing a coherent image plane spatial coordinate system; (f,g) Representing a pupil plane spectral coordinate system; (f',g') denotes a mask spectral coordinate system;Hrepresents the optical transfer function, i.e. the pupil function, of the imaging system;H i representing the position of the light source according to points: (S x ,S y ) A pupil after translation;Mrepresenting the mask spectrum.
According to the theory of lithographic imaging, the lithographic imaging system is a strictly optical diffraction limited system with a cut-off frequency approximately equal to the numerical apertureNAAnd wavelength of illuminationλRatio of (i) to (ii)NA/λ. In the lithography imaging model, the pupil function corresponds to a circular low-pass filter with a radius of the cutoff frequency of the optical system.
The lithographic imaging process being classical partially coherent imaging, the light source thereofCan be defined as the radius of the illumination sourcer S Radius of pupilr H Ratio of (i) to (ii)σ=r S /r H . Thus, in the partially coherent imaging model, the light source has a radius ofσ NA/λThe circle of (c). According to the Hopkins imaging method, the pupil matrix center is shifted according to the spot light source position. Thus, the maximum shift in the center of the pupil matrix is at a radius ofr S +r H Is circular. In the photoetching imaging process, the mask frequency spectrum participating in the process is the effective frequency spectrum, and the coverage area of the mask frequency spectrum is a radiusr S +2r H The circular area DC (as shown in fig. 3), only this part is also involved in the imaging process. Discretizing the coherent imaging model can be represented by the following mathematical model:
wherein,representing a discretized coherent image plane spatial coordinate system; (j,k) Representing a discretized pupil plane spectral coordinate system; (j',k') represents a discretized mask spectral coordinate system;N ext indicating the number of extracted active spectral sample points,;
according to Abbe's method, aerial image obtained by partially coherent imaging processI ext Is represented by the following mathematical model:
wherein, in the followingm, nAnd establishing a photoetching system illumination light source model for the index-mode coordinate system.S(m,n) Is shown in the light source model coordinate system at (m,n) Point light source at the position, the size of the illumination light source model matrix of the photoetching system isN S ×N S ,m,n=1,2,3,…,N S 。
Amplifying the mathematical model by interpolation method to obtain standard spatial image intensity distributionI。
The objective function F 1 And F 2 Defining the sum of absolute values of differences between each unit of photoresist pattern matrix and ideal pattern as pattern error, and using the pattern error as target function F 1 . Marking the areas of the inner and outer pixel points near the boundary of the mask pattern, and calculating the edge placement error, taking the edge prevention error as an objective function F 2 . Objective function F 1 And F 2 Are represented as the following mathematical models, respectively:
wherein,RP(x i ,y i ) A photoresist pattern is represented and,M * (x i ,y i ) Representing an ideal figure.ICC Edge And representing a matrix formed by the extracted edge pixel points.S ext Representing the extracted light source matrix.M ' ext And representing a matrix formed by extracting the edge pixel points of the ideal graph according to the index position.minimizeThe method is characterized in that in a tabu search-genetic iteration optimization-based model, the numerical value of an objective function is guaranteed to be minimum when iteration stops.
FIG. 5 is a schematic diagram of a light source encoding and Gaussian filtering method used in an imaging model of a lithography system. The light source can be divided into four identical parts according to the special optical symmetry of the illumination light source, namely a ring-shaped illumination light sourceSThe first part of (1)J 1 A second partJ 2 Third partJ 3 Fourth sectionJ 4 . In the optimization model, only one quarter of the light source points are required as optimization variables. Step 501 represents extracting the light source S quarter active area as an optimized area, marking the position of each active point by encoding, and performing a matrix transformation operation. Encoding the active source points into column vectors by encoding and matrix conversion operationsS ' Of matrix size ofmX 1. Step 502 represents randomly generating an optimized variable group in an initialization stage according to the number of variables in the optimized areaS r I.e. light source variable population matrixP S . Step 503 represents selecting the current best light source variable individual according to the objective function valueS opt . Step 504 represents obtaining the complete intensity distribution of the discrete light sources by matrix conversion and the like. Since the light source intensity distribution in the photolithography system is in continuous gradient change, the discrete light source is subjected to gaussian filtering operation in step 505 to obtain the illumination light source with continuously changing gray scaleS * 。
The application is described in detail above with reference to specific embodiments and the description of the drawings. The method carries out local extraction processing on the mask characteristic diffraction spectrum information, and reduces the calculation complexity of the space image imaging model. Quarter coding and optimization are carried out on the illumination light source and the mask, and the optimizing speed of the optimization model is improved. The taboo search algorithm is combined with the genetic algorithm, the global convergence expressive property of the genetic algorithm is improved, the target function graph error and the edge prevention error are combined, and the imaging expressive property of the photoetching system is effectively improved.
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (1)
1. A space image prediction and image quality optimization method for a sensor chip projection lithography machine is a hybrid method of tabu search-genetic algorithm, and is characterized by comprising the following steps:
step 1: initializing photoetching system space image imaging model parameters in photoetching light source mask optimization model based on tabu search-genetic algorithm, inputting light sourceSAnd mask patternM(ii) a Carrying out gridding processing on the mask pattern M and the input light source S, extracting effective characteristic information of the mask pattern by adopting an effective frequency spectrum extraction method, and representing a space image generated by a photoetching optical system by using the following mathematical model by combining an Abbe imaging method:
wherein,Irepresenting the spatial image intensity distribution, i.e. the result of partially coherent imaging;N ' S representing the number of effective point light sources in the pixelated light source;CCI i the coherent imaging results produced by a single point source are represented by the following mathematical model:
wherein, (ii) (x i ,y i ) Representing a coherent image plane spatial coordinate system; (f,g) Representing a pupil plane spectral coordinate system; (f',g') denotes a mask spectral coordinate system;Hrepresents the optical transfer function, i.e. the pupil function, of the imaging system;H i representing the position of the light source according to points: (S x ,S y ) A pupil after translation;Mrepresenting a mask spectrum;
the coherent imaging model was discretized and represented by the following mathematical model:
wherein,representing a discretized coherent image plane spatial coordinate system; (j,k) Representing a discretized pupil plane spectral coordinate system; (j',k') represents a discretized mask spectral coordinate system;N ext representing the number of extracted valid spectral samples,;
according to Abbe's method, aerial images obtained by partially coherent imaging processI ext Is represented by the following mathematical model:
wherein, in the followingm, nEstablishing a photoetching system illumination light source model for an index-mode coordinate system,S(m,n) Is shown in the light source model coordinate system at (m,n) Point light source at the position, the size of the illumination light source model matrix of the photoetching system isN S ×N S ,m,n=1,2,3,…,N S ;
Step 2: coding the light source and the mask effective unit as target function variables;
and step 3: randomly generating a group variable whole matrix to be optimized according to a heuristic optimization algorithm structureP(ii) a Therefore, in the light source mask optimization model, a light source and mask initial population-light source variable population matrix is randomly generatedP S And mask variable population matrixP M The light source variable group matrixP S And mask variable population matrixP M Are respectively asP S =[S 1 ,S 2 ,S 3 ,…,S P ],P M =[M 1 ,M 2 ,M 3 ,…,M P ];
And 4, step 4: according to the initial light source variable matrixP S Calculating the intensity distribution of the space image generated by each individual through the space image imaging model of the photoetching systemThe spatial image prediction under different imaging conditions is realized,i=1,2, …,p(ii) a And according to an objective function F 1 And F 2 Calculating an adaptation value and selecting the currently best initial light source variable individual by comparison>And an adapted value->(ii) a Wherein it is present>A group matrix formed in the space image imaging model of the lithography system according to different light source variablesP S Of (2)S P The resulting aerial image intensity distribution matrixI S I.e. is->,/>;/>
Defining the sum of absolute values of the difference values of each unit of the photoresist pattern matrix and the ideal pattern as a pattern error, and taking the pattern error as a target function F 1 (ii) a Marking areas of inner and outer pixel points near the boundary of the mask pattern, calculating edge placement error, and taking the edge placement error as an objective function F 2 (ii) a Objective function F 1 And F 2 Are represented as the following mathematical models, respectively:
wherein,RP(x i ,y i ) A photoresist pattern is represented and,M * (x i ,y i ) Representing an ideal graph;ICC Edge representing a matrix formed by the extracted edge pixel points;S ext representing the extracted light source matrix;M ' ext representing a matrix formed by extracting ideal graph edge pixel points according to the index positions;minimizeexpressed in a tabu search-genetic iteration-based optimization modelWhen iteration stops, the numerical value of the objective function is ensured to be minimum;
and 5: the current optimal initial light source variable obtained in the step 4 is used as an individualPerforming iterative updating on the light source variable individuals until iteration stops as an input condition of a taboo search-genetic algorithm optimization model;
step 6: according to the light source variable individual after iteration is stopped, matrix inversion mirror image and Gaussian filtering operation are executed, and the appearance of the light source is recovered;
and 7: using the obtained light source as the light source condition of the mask optimization model; population matrix according to mask variablesP M Calculating and recording current best initial mask variable individualAnd an adapted value->;
And step 8: the current best initial mask variable obtained in the step 7 is used as an individualPerforming iterative update on mask variable individuals until iteration stops as an input condition of a tabu search-genetic algorithm optimization model;
and step 9: outputting the optimal light source intensity distribution and the mask pattern layout.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211417846.6A CN115469555B (en) | 2022-11-14 | 2022-11-14 | Space image prediction and image quality optimization method for sensor chip projection lithography machine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211417846.6A CN115469555B (en) | 2022-11-14 | 2022-11-14 | Space image prediction and image quality optimization method for sensor chip projection lithography machine |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115469555A CN115469555A (en) | 2022-12-13 |
CN115469555B true CN115469555B (en) | 2023-03-31 |
Family
ID=84338115
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211417846.6A Active CN115469555B (en) | 2022-11-14 | 2022-11-14 | Space image prediction and image quality optimization method for sensor chip projection lithography machine |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115469555B (en) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101118421A (en) * | 2007-09-13 | 2008-02-06 | 北京航空航天大学 | Intelligent non-linearity PID controlling parameter tuning based on self-adapting ant colony |
CN101349871A (en) * | 2008-09-05 | 2009-01-21 | 上海微电子装备有限公司 | Photo-etching illuminating apparatus |
CN102193937A (en) * | 2010-03-12 | 2011-09-21 | 同济大学 | Image search method based on quick intelligent taboo search |
CN102692814A (en) * | 2012-06-18 | 2012-09-26 | 北京理工大学 | Light source-mask mixed optimizing method based on Abbe vector imaging model |
CN104111592A (en) * | 2014-08-06 | 2014-10-22 | 中国科学院光电技术研究所 | Method for realizing variable free illumination pupil based on micro-mirror array |
CN105051737A (en) * | 2013-02-22 | 2015-11-11 | 新思科技有限公司 | Hybrid evolutionary algorithm for triple-patterning |
CN105069194A (en) * | 2015-07-20 | 2015-11-18 | 中国科学院长春光学精密机械与物理研究所 | Genetic algorithm based optimization method for photoetching attenuation type mask |
TWI545519B (en) * | 2011-04-12 | 2016-08-11 | 尼康股份有限公司 | Imaging apparatus and program |
WO2020108861A1 (en) * | 2018-11-26 | 2020-06-04 | Asml Netherlands B.V. | Determining a mark layout across a patterning device or substrate |
CN111290281A (en) * | 2020-03-23 | 2020-06-16 | 中国科学院光电技术研究所 | Wavefront control method based on ADRC-Smith algorithm |
CN113242997A (en) * | 2018-12-19 | 2021-08-10 | Asml荷兰有限公司 | Method for sample plan generation and optimization |
CN113911728A (en) * | 2021-09-30 | 2022-01-11 | 江苏白王口腔护理用品有限公司 | Electric toothbrush brush head dynamic feeding system and feeding method based on vision |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002071738A1 (en) * | 2001-03-02 | 2002-09-12 | Dai Nippon Printing Co., Ltd. | Dither mask creating method and creating device |
-
2022
- 2022-11-14 CN CN202211417846.6A patent/CN115469555B/en active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101118421A (en) * | 2007-09-13 | 2008-02-06 | 北京航空航天大学 | Intelligent non-linearity PID controlling parameter tuning based on self-adapting ant colony |
CN101349871A (en) * | 2008-09-05 | 2009-01-21 | 上海微电子装备有限公司 | Photo-etching illuminating apparatus |
CN102193937A (en) * | 2010-03-12 | 2011-09-21 | 同济大学 | Image search method based on quick intelligent taboo search |
TWI545519B (en) * | 2011-04-12 | 2016-08-11 | 尼康股份有限公司 | Imaging apparatus and program |
CN102692814A (en) * | 2012-06-18 | 2012-09-26 | 北京理工大学 | Light source-mask mixed optimizing method based on Abbe vector imaging model |
CN105051737A (en) * | 2013-02-22 | 2015-11-11 | 新思科技有限公司 | Hybrid evolutionary algorithm for triple-patterning |
CN104111592A (en) * | 2014-08-06 | 2014-10-22 | 中国科学院光电技术研究所 | Method for realizing variable free illumination pupil based on micro-mirror array |
CN105069194A (en) * | 2015-07-20 | 2015-11-18 | 中国科学院长春光学精密机械与物理研究所 | Genetic algorithm based optimization method for photoetching attenuation type mask |
WO2020108861A1 (en) * | 2018-11-26 | 2020-06-04 | Asml Netherlands B.V. | Determining a mark layout across a patterning device or substrate |
CN113168110A (en) * | 2018-11-26 | 2021-07-23 | Asml荷兰有限公司 | Determining the layout of marks on the entire patterning device or substrate |
CN113242997A (en) * | 2018-12-19 | 2021-08-10 | Asml荷兰有限公司 | Method for sample plan generation and optimization |
CN111290281A (en) * | 2020-03-23 | 2020-06-16 | 中国科学院光电技术研究所 | Wavefront control method based on ADRC-Smith algorithm |
CN113911728A (en) * | 2021-09-30 | 2022-01-11 | 江苏白王口腔护理用品有限公司 | Electric toothbrush brush head dynamic feeding system and feeding method based on vision |
Non-Patent Citations (3)
Title |
---|
Global optimization of source and mask in inverse lithography via tabu search combined with genetic algorithm;haifeng Sun etal.;《Optics Express》;20220731;第24166-24185页 * |
基于遗传禁忌混合搜索算法的设备布局研究;齐继阳等;《2004 年"安徽制造业发展"博士科技论坛》;20041231;第1-5页 * |
基于遗传算法的光刻机性能匹配方法;茅言杰等;《第十七届全国光学测试学术交流会摘要集》;20181231;第143页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115469555A (en) | 2022-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yifan et al. | Patch-based progressive 3d point set upsampling | |
CN110678961B (en) | Simulating near field images in optical lithography | |
CN108133456A (en) | Face super-resolution reconstruction method, reconstructing apparatus and computer system | |
CN104133348B (en) | A kind of adaptive optical etching system light source optimization method | |
CN111028265B (en) | Target tracking method for constructing correlation filtering response based on iteration method | |
CN103631096B (en) | Source mask polarization optimization method based on Abbe vector imaging model | |
CN111968138A (en) | Medical image segmentation method based on 3D dynamic edge insensitivity loss function | |
Schulz et al. | Non-linear hypothesis testing of geometric object properties of shapes applied to hippocampi | |
CN102998896B (en) | Basic module-based mask main body graph optimization method | |
CN106910215B (en) | Super-resolution method based on fractional order gradient interpolation | |
KR20230173649A (en) | Training method, apparatus, device, and storage medium for photolithography mask generation model | |
CN113269689B (en) | Depth image complement method and system based on normal vector and Gaussian weight constraint | |
CN113158487B (en) | Wavefront phase difference detection method based on long-short term memory depth network | |
GB2581524A (en) | Facial shape representation and generation system and method | |
Lichtenstein et al. | Deep eikonal solvers | |
Mohamed et al. | A data and compute efficient design for limited-resources deep learning | |
JP2022045893A (en) | Mesh noise removing method based on graph convolution network | |
Knigge et al. | Modelling Long Range Dependencies in $ N $ D: From Task-Specific to a General Purpose CNN | |
CN114723884A (en) | Three-dimensional face reconstruction method and device, computer equipment and storage medium | |
CN115469555B (en) | Space image prediction and image quality optimization method for sensor chip projection lithography machine | |
Ouasfi et al. | Few ‘zero level set’-shot learning of shape signed distance functions in feature space | |
Durasov et al. | Double refinement network for efficient monocular depth estimation | |
Zhou et al. | MaskNet++: Inlier/outlier identification for two point clouds | |
CN103901713B (en) | Self-adaption optical proximity effect correction method adopting kernel regression technology | |
CN110766695B (en) | Image sparse representation-based matting method |
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 |