CN101192009B - Method for establishing OPC model - Google Patents

Method for establishing OPC model Download PDF

Info

Publication number
CN101192009B
CN101192009B CN2006101188210A CN200610118821A CN101192009B CN 101192009 B CN101192009 B CN 101192009B CN 2006101188210 A CN2006101188210 A CN 2006101188210A CN 200610118821 A CN200610118821 A CN 200610118821A CN 101192009 B CN101192009 B CN 101192009B
Authority
CN
China
Prior art keywords
value
additive
initial
initial value
opc model
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
Application number
CN2006101188210A
Other languages
Chinese (zh)
Other versions
CN101192009A (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.)
Semiconductor Manufacturing International Shanghai Corp
Original Assignee
Semiconductor Manufacturing International Shanghai Corp
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 Semiconductor Manufacturing International Shanghai Corp filed Critical Semiconductor Manufacturing International Shanghai Corp
Priority to CN2006101188210A priority Critical patent/CN101192009B/en
Publication of CN101192009A publication Critical patent/CN101192009A/en
Application granted granted Critical
Publication of CN101192009B publication Critical patent/CN101192009B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)

Abstract

A method of building optical proximity correction model comprises the following steps of collecting optical proximity correction data to a database from processes, determining an initial value X0, an end value Xn and a fixed step size h, calculating initial value slope coefficient according to the initial value X0, obtaining an initial optimal step size Y0 according to the initial value slope coefficient, obtaining a first additive value X1 according to the initial value X0 and the initial optimal step size Y0, repeating the steps, working out a second additive valueX2, a third additive value X3,..., till the end value Xn, fitting the initial value X0, the first additive value X1, the second additive value X2, the third additive value X3,..., till an n-1 additive value and the end value Xn. Through the above steps, as step size values can be regulated according to the slope coefficient worked out, regulation step number can be reduced and simulation time is shortened.

Description

Set up the method for OPC model
Technical field
The present invention relates to semiconductor fabrication process, particularly the photomask data processing technique is especially set up the optical proximity correction model.
Background technology
The current realization of photomask data deviation is usually based on some models, and described model is calibrated in concrete photoetching process, is commonly referred to as OPC model.For example, the correction of optical proximity effect (OPE) often need be attempted " calibration " photoetching process so that compensation OPE.Prior art comprises so-called calibration parameter is associated with OPC model that this need carry out one group of detailed optical parameter measurement at different feature locations.It is many more that measurement data is collected, and the accuracy of calibration parameter is just good more.Yet, for reliable OPC model parametric calibration different under environment, optical parameter measurement more than the common needs in different critical characteristic positions are hundreds of, this is labour intensive and time consuming work, what is worse owing to the reason of experience level, take what kind of optical parameter measurement often to become and depend on the operator, thus this obvious negative effect parametric calibration limited the overall validity of optical nearing correction technique.
U.S. Patent number is 6,563,566 disclose and carry out illumination by a series of calculating of attempting linearization mask transmission optimization and optimize and photomask optimization, wherein disclose maximization minimum NILS (normalized image log slope) and selection and will be used in various conditions in this calculating.Because rely on the symmetry of photomask, this calculating is restricted, yet the linearization of the photomask transmission of use need be used several approximate in calculating, rather than actual imaging equation itself, this use photomask formation institute important plan as the time can produce error.The linearization of photomask transmission also needs to use a lot of variablees, and this needs a large amount of computing times and carries out this calculating.Therefore along with the reducing of logical implication size, be necessary to provide the photomask implementation of using the minimum of computation time accurately to form desired image.
The method for making of existing OPC model, for example this optical parametric of logarithm value aperture is optimized, and as shown in Figure 1, S101 at first gathers the numerical aperture data to database from technology; S102 determines initial value X 0, end value X nWith fixed step size h; S103 is with initial value X 0Carry out addition with fixed step size h, draw the first additive value X 1S104 is with the first additive value X 2Carry out addition with fixed step size h, draw the second additive value X 2S105 continues above-mentioned steps, and cycle calculations goes out the value added X of third phase 3, the 4th additive value X 4To end value X nS106 is with initial value X 0, the first additive value X 1, the second additive value X 2, the value added X of third phase 3, the 4th additive value X 4End value X nCarry out match.
OPC model is a function, has characterized the difference of the optical field distribution that the optical field distribution that calculates by optical theory and actual amount record.This function is relevant with some optical parametrics, such as numerical aperture.When that was some occurrences when numerical aperture, function had minimum value, and it is best just to illustrate that the light field of model description and actual light field are coincide.
In general our function F (X of investigating 0) all be the multivariate function usually, its independent variable parameter has a plurality of, and what Fig. 2 represented is when independent variable X is two (two-dimensional function), F (X 0) possible distribution shape.
As shown in Figure 3, suppose the function F (X of investigation 0) be two-dimensional function, that at first determines the initial value X of two parameters before beginning modeling with the existing method of setting up OPC model 0, end value X NAnd fixed step size value h, the combination of the value of used parameter just can be regarded as on the lattice point that is distributed in a two-dimensional mesh (among the figure in the ellipse) in optimizing process so, and the parameter value combination outside the lattice point is not to be used in calculating.
The method for making of existing OPC model because step-length is a fixed value, if intermediate data is a lot of between initial value and the end value, can cause the adjusting step number too many, and simulated time is elongated.
Summary of the invention
The problem that the present invention solves provides a kind of method of setting up OPC model, prevents because step-length is a fixed value, if intermediate data is a lot of between initial value and the end value, can cause the adjusting step number too many, and simulated time is elongated.
For addressing the above problem, the invention provides a kind of method of setting up OPC model, comprise the following steps: from technology, to gather optical nearing correction data to database; Determine initial value X 0, end value X NWith fixed step size h; According to initial value X 0Calculate the initial value slope
Figure G061B8821020061218D000031
According to the initial value slope
Figure G061B8821020061218D000032
Draw initial optimization step-length Y with fixed step size h 0According to initial value X 0With initial optimization step-length Y 0Draw the first additive value X 1Repeat above-mentioned steps, calculate the second additive value X 2, the value added X of third phase 3The n-1 additive value is to end value X nWith initial value X 0, the first additive value X 1, the second additive value X 2, the value added X of third phase 3... n-1 additive value and end value X nCarry out match.The formula that calculates additive value is: X n=X N-1+ Y N-1, Y wherein N-1For optimizing step-length, n is a natural number.
Described calculation optimization step-length Y N-1Formula be: Y n - 1 = ∂ F ∂ X | X = X n - 1 × h , Wherein
Figure G061B8821020061218D000034
Be slope.
The formula of described slope calculations is: ∂ F ∂ X | X = X n - 1 = F ( X n - 1 + ΔX ) - F ( X n - 1 ) ΔX , Wherein Δ X is the increment of independent variable X, F (X 0) equal X for optical nearing correction data 0The time, the degree of agreement when OPC model and optical nearing correction actual correction.
Compared with prior art, the present invention has the following advantages: because step-length can be regulated according to the slope that calculates, can realize regulating step number like this and reduce, simulated time shortens.
Description of drawings
Fig. 1 is the process flow diagram that prior art is set up the OPC model method;
Fig. 2 is prior art F (X when working as independent variable X and being two 0) possible distribution shape;
Fig. 3 is the grid chart that prior art is set up OPC model;
Fig. 4 is the process flow diagram that the present invention sets up the OPC model method;
Fig. 5 is F (X of the present invention 0) curve map of function;
Fig. 6 is the effect comparison diagram of the OPC model method of prior art and the present invention's foundation.
Embodiment
The linearization of photomask transmission need be used a lot of variablees, this needs a large amount of computing times carries out this calculating, and prior art is because step-length is a fixed value, if intermediate data is a lot of between initial value and the end value, can cause the adjusting step number too many, simulated time is elongated.So reducing along with the logical implication size, be necessary to provide the photomask implementation of using the minimum of computation time accurately to form desired image, the present invention can realize regulating step number like this and reduce because step-length can be regulated according to the slope that calculates, and simulated time shortens.For above-mentioned purpose of the present invention, feature and advantage can be become apparent more, the specific embodiment of the present invention is described in detail below in conjunction with accompanying drawing.
Fig. 4 is the process flow diagram that the present invention sets up the OPC model method.As shown in Figure 4, execution in step S201 gathers optical nearing correction data to database from technology; S202 determines initial value X 0, end value X nWith fixed step size h; S203 is according to initial value X 0Calculate the initial value slope
Figure G061B8821020061218D000041
S204 is according to the initial value slope Draw initial optimization step-length Y with fixed step size h 0S205 is according to initial value X 0With initial optimization step-length Y 0Draw the first additive value X 1S206 repeats above-mentioned steps, calculates the second additive value X 2, the value added X of third phase 3The n-1 additive value is to end value X nS207 is with initial value X 0, the first additive value X 1, the second additive value X 2, the value added X of third phase 3... n-1 additive value and end value X nCarry out match.In the present embodiment, the formula that calculates additive value is: X n=X N-1+ Y N-1, Y wherein N-1For optimizing step-length, n is a natural number.The concrete first additive value X 1=X 0+ Y 0, the second additive value X 2=X 1+ Y 1, the value added X of third phase 3=X 2+ Y 2N-1 additive value X N-1=X N-2+ Y N-2
Described optimization step-length Y N-1Computing formula be: Y n - 1 = ∂ F ∂ X | X = X n - 1 × h , Wherein
Figure G061B8821020061218D000052
Be slope.In the present embodiment, concrete initial optimization step-length Y 0 = ∂ F ∂ X | X = X 0 × h , Second optimizes step-length Y 1 = ∂ F ∂ X | X = X 1 × h , The 3rd optimizes step-length Y 2 = ∂ F ∂ X | X = X 2 × h . . . . . . N-1 optimizes step-length
Y n - 2 = ∂ F ∂ X | X = X n - 2 × h .
In the present embodiment, the formula of described slope calculations is: ∂ F ∂ X | X = X n - 1 = F ( X n - 1 + ΔX ) - F ( X n - 1 ) ΔX , Wherein Δ X is the increment of independent variable X, F (X 0) equal X for optical nearing correction data 0The time, the degree of agreement when OPC model and optical nearing correction actual correction.Concrete initial value slope is ∂ F ∂ X | X = X 0 = F ( X 0 + ΔX ) - F ( X 0 ) ΔX , The first additive value slope is ∂ F ∂ X | X = X 1 = F ( X 1 + ΔX ) - F ( X 1 ) ΔX , The second additive value slope is ∂ F ∂ X | X = X 2 = F ( X 2 + ΔX ) - F ( X 2 ) ΔX , The value added slope of third phase is ∂ F ∂ X | X = X 3 = F ( X 3 + ΔX ) - F ( X 3 ) ΔX , N-1 additive value slope is
∂ F ∂ X | X = X n - 1 = F ( X n - 1 + ΔX ) - F ( X n - 1 ) ΔX .
Fig. 5 is F (X of the present invention 0) curve map of function.Suppose that we optimize F (X 0) time only consider a variable parameter, function is with respect to the funtcional relationship of variable element such as the curve representation among Fig. 5 so, OPC model is from initial value X 0Set out, find optimized parameter value (the curve the lowest point part among the figure), after the optimum value that finds all parameters, these the most optimized parameter values are carried out match, finish optimization OPC model.
Fig. 6 is the effect comparison diagram of the OPC model method of prior art and the present invention's foundation.As shown in Figure 6, when parameters optimization is 4,6 or 8 respectively, set up the method for OPC model and method that the present invention sets up OPC model is finished the needed time of last optimization with existing.When parameters optimization was 4, existing to set up the required time of OPC model be 12468s, and the present invention to set up the required time of OPC model be 4513s; When parameters optimization was 6, existing to set up the required time of OPC model be 30465s, and the present invention to set up the required time of OPC model be 10145s; When parameters optimization is 4, having now and setting up the required time of OPC model is 48627s, and the present invention sets up the required time of OPC model is 16031s, the present invention is owing to carry out slope calculating with each value, and then slope be multiply by fixed step size as optimizing step-length, make and regulate the step number minimizing, simulated time shortens.
Though the present invention with preferred embodiment openly as above; but it is not to be used for limiting the present invention; any those skilled in the art without departing from the spirit and scope of the present invention; can make possible change and modification, so protection scope of the present invention should be as the criterion with the scope that claim of the present invention was defined.

Claims (1)

1. a method of setting up OPC model is characterized in that, comprises the following steps:
From technology, gather optical nearing correction data to database;
Determine initial value X 0, end value X nWith fixed step size h;
According to initial value X 0Calculate the initial value slope
Figure FSB00000472174900011
According to the initial value slope
Figure FSB00000472174900012
Draw initial optimization step-length Y with fixed step size h 0
According to initial value X 0With initial optimization step-length Y 0Draw the first additive value X 1
Repeat above-mentioned steps, calculate the second additive value X 2, the value added X of third phase 3... n-1 additive value X N-1To end value X n, n is a natural number;
The described first additive value X 1, the second additive value X 2, the value added X of third phase 3... n-1 additive value X N-1To end value X n, obtain in the following way: X n=X N-1+ Y N-1, Y wherein N-1For optimizing step-length;
Described optimization step-length Y N-1Obtain in the following way:
Figure FSB00000472174900013
Wherein
Figure FSB00000472174900014
Be slope;
With initial value X 0, the first additive value X 1, the second additive value X 2, the value added X of third phase 3... n-1 additive value X N-1With end value X nCarry out match, set up OPC model photoetching process is calibrated.
CN2006101188210A 2006-11-28 2006-11-28 Method for establishing OPC model Active CN101192009B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2006101188210A CN101192009B (en) 2006-11-28 2006-11-28 Method for establishing OPC model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2006101188210A CN101192009B (en) 2006-11-28 2006-11-28 Method for establishing OPC model

Publications (2)

Publication Number Publication Date
CN101192009A CN101192009A (en) 2008-06-04
CN101192009B true CN101192009B (en) 2011-07-06

Family

ID=39487072

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2006101188210A Active CN101192009B (en) 2006-11-28 2006-11-28 Method for establishing OPC model

Country Status (1)

Country Link
CN (1) CN101192009B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103676490B (en) * 2012-09-20 2015-11-25 中芯国际集成电路制造(上海)有限公司 A kind of method monitoring weakness Crack cause

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103792785B (en) * 2012-10-29 2017-03-08 中芯国际集成电路制造(上海)有限公司 A kind of method that optical proximity correction is carried out to the figure with low picture contrast
JP6491974B2 (en) * 2015-07-17 2019-03-27 日立化成株式会社 EXPOSURE DATA CORRECTION DEVICE, WIRING PATTERN FORMING SYSTEM, AND WIRING BOARD MANUFACTURING METHOD

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1392453A (en) * 2001-06-20 2003-01-22 旺宏电子股份有限公司 Optical nearing correcting method
US6563566B2 (en) * 2001-01-29 2003-05-13 International Business Machines Corporation System and method for printing semiconductor patterns using an optimized illumination and reticle
CN1818790A (en) * 2005-02-07 2006-08-16 中芯国际集成电路制造(上海)有限公司 Optical adjacent correction for mask pattern during photoetching process

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6563566B2 (en) * 2001-01-29 2003-05-13 International Business Machines Corporation System and method for printing semiconductor patterns using an optimized illumination and reticle
CN1392453A (en) * 2001-06-20 2003-01-22 旺宏电子股份有限公司 Optical nearing correcting method
CN1818790A (en) * 2005-02-07 2006-08-16 中芯国际集成电路制造(上海)有限公司 Optical adjacent correction for mask pattern during photoetching process

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103676490B (en) * 2012-09-20 2015-11-25 中芯国际集成电路制造(上海)有限公司 A kind of method monitoring weakness Crack cause

Also Published As

Publication number Publication date
CN101192009A (en) 2008-06-04

Similar Documents

Publication Publication Date Title
KR20210047351A (en) Structural topology optimization method based on material field reduction and water distribution
AU2021233949A1 (en) SWMM and EFDC coupling model-based regulation and storage project environmental effect assessment method and device
CN103177289B (en) Modeling method for noise-uncertainty complicated nonlinear dynamic system
CN102968529B (en) Method for quantifying computed result non-determinacy interval of water supply pipe network model
CN112904266B (en) Method and device for predicting service life of electric energy meter
CN109977460A (en) A kind of multi-objective optimization design of power method based on vehicle body section parameter
CN101192009B (en) Method for establishing OPC model
CN107860889A (en) The Forecasting Methodology and equipment of the soil organism
CN110442911B (en) High-dimensional complex system uncertainty analysis method based on statistical machine learning
EP2312506A1 (en) Computer-based method and device for automatically providing a prediction on a future energy demand to an energy source
CN101339431B (en) Dynamic material consumption energy consumption real-time control method
CN115169818A (en) Digital modeling-based technical engineering measurement and calculation method and system
CN110377989B (en) Non-uniformity two-variable design flood calculation method based on hydrothermal coupling balance
CN109840368B (en) Irrigation area channel water level flow relation calibration method
CN115146929A (en) Three-stage DEA-Tobit port energy system efficiency evaluation method containing unexpected output
CN113361774B (en) Ecological space intensive utilization optimization method and system
CN116029037A (en) IFC-based assembly type construction engineering amount calculating method
CN117273200A (en) Runoff interval forecasting method based on convolution optimization algorithm and Pyraformer neural network
CN116628979A (en) Multiple regression total primary productivity pixel-by-pixel prediction method, device and medium
CN107025330B (en) Method for calculating transverse distribution coefficient of single-box multi-chamber wide beam bridge
CN115657166A (en) Method and device for determining weather forecast value across multiple climate areas
CN112100867B (en) Power grid load prediction method
US20220173821A1 (en) Machine learning apparatus
CN109858699B (en) Water quality quantitative simulation method and device, electronic equipment and storage medium
CN104251184A (en) Method and system for automatic recommendation of optimal wind power generation set

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant