This application claims priority to Taiwan Application Serial Number 103111289, filed Mar. 26, 2014, which is herein incorporated by reference.
Field of Invention
The present invention relates to a metrology sampling method. More particularly, the present invention relates to a metrology sampling method for reducing a sampling rate.
Description of Related Art
Nowadays, most of the semiconductor and TFT-LCD plants adopt sampling test methods to monitor the quality of each product or workpiece (i.e. “wafer” in IC foundries or “glass” in TFT-LCD factories) processed by process tools. In general, after each N (for example, 25) workpieces are processed by the process tool, a manufacturing system designates the Nth workpiece in the each N workpieces as a workpiece expected (scheduled) to be measured, i.e. the sampling rate is 1/N. The workpiece expected to be measured then will be sent to a metrology tool for measurement, thereby inspecting the production quality of the process tool. This conventional sampling method is based on the presupposition that no abnormal circumstances will abruptly occur in the production processes, and thus the measurement result of the sampled product or workpiece can be used to determine the quality of the workpieces in the same lot with the sampled workpiece. Because the measurement time and tool required by the real workpiece metrology result in the increase of cycle time and production cost. Therefore, reducing the sampling rate to as low as possible is an important task for manufactures for reducing cycle time and production cost.
On the other hand, virtual metrology (VM) can be used to lower the frequency of actual measurement on workpiece to reduce the sampling rate. However, if a production variation occurs on the workpiece which is planned for measurement, no real metrology is available during this period for updating the VM models, thus resulting in poor VM prediction accuracy. Therefore, how to sample and obtain an appropriate workpiece in time affects the prediction accuracy of VM models.
Hence, there is a need to provide a metrology sampling method for overcoming the aforementioned shortcomings of the conventional skill.
An object of the present invention is to provide a metrology sampling method for reducing a metrology sampling for workpieces.
Another object of the present invention is to provide a metrology sampling method for providing an actual measurement value of a workpiece in time to retuning or retraining a VM model, thereby assuring VM accuracy.
According to the aforementioned objects, a metrology sampling method is provided, in this metrology sampling method, at first, a plurality of sets of historical process data used by a process tool for processing historical workpieces are collected. Then, a model-creation step is performed. In the model-creation step, the sets of historical process data are used to build a DQIX (Process Data Quality Index) model and a GSI (Global Similarity Index) model and compute a DQIX threshold and a GSI threshold. Thereafter, a metrology workpiece sampling step is performed. In the metrology workpiece sampling step, a workpiece is provided to the process tool, in which the process tool has a set of process data for processing the workpiece. Then, the set of process data is inputted into the DQIX model and the GSI model, thereby obtaining a DQIX value and a GSI value of the set of process data of the workpiece. When the DQIX value of the workpiece is greater than the DQIX threshold, a measurement of the workpiece is skipped. When the DQIX value of the workpiece is smaller than or equal to the DQIX threshold, a step is performed to check if the workpiece is a workpiece expected to be measured, thereby obtaining a first checking result. When the first checking result is true, metrology is performed on the workpiece. When the first checking result is false, a step is performed to check if the GSI value of the workpiece is smaller than or equal to the GSI threshold, thereby obtaining a second checking result. When the second checking result is true, the measurement of the workpiece is skipped.
According to the aforementioned objects, another metrology sampling method is provided. In the metrology sampling method, a predetermined workpiece sampling rate 1/N is assigned, in which the predetermined workpiece sampling rate is directed to selecting the Nth workpiece in each N workpieces processed by a process tool as a workpiece expected to be measured. In a metrology workpiece sampling step of the metrology sampling method, a workpiece count is increased by 1. Then, a first checking step is performed to check if the DQIX value is smaller than or equal to the DQIX threshold, thereby obtaining a first checking result. When the first checking result is false, a measurement of the workpiece is skipped. When the first checking result is true, a second checking step is performed to check if the workpiece count is greater than or equal to N, thereby obtaining a second checking result. When the second checking result is true, metrology is performed on the workpiece and the workpiece count is set to 0. When the second checking result is false, a third checking step is performed to check if the GSI value of the workpiece is smaller than or equal to the GSI threshold, thereby obtaining a third checking result. When the third checking result is true, the measurement of the workpiece is skipped.
According to the aforementioned objects, a computer program product stored on a non-transitory tangible computer readable recording medium is provided. When the computer program product is loaded by a computer and executed, the aforementioned metrology sampling methods are performed.
Thus, with the application of the embodiments of the present invention, the workpiece sampling rate can be greatly lowered, and the VM accuracy can be effectively assured.
It is to be understood that both the foregoing general description and the following detailed description are by examples, and are intended to provide further explanation of the invention as claimed.
The invention can be more fully understood by reading the following detailed description of the embodiment, with reference made to the accompanying drawings as follows:
Reference will now be made in detail to the present embodiments of the invention, examples of which are illustrated in the accompanying drawings. Wherever possible, the same reference numbers are used in the drawings and the description to refer to the same or like parts.
Embodiments of the present invention apply various indices to develop an Intelligent Sampling Decision (ISD) scheme, in which those indices can be used to detect various process-related status changes of a process tool during production (such as preventive-maintenance (PM) operation, process-tool malfunction, recipe change, etc.), or information of abnormality of the process tool (such as poor process data quality, deviation (shift or drift) of process data, poor metrology data quality, etc.). The indices includes a reliance Index (RI), a global similarity index (GSI), a process data quality index (DQIX) and a metrology data quality index (DQIy). The RI value, the GSI value, the DQIX value and the DQIy value used in the embodiments of the present invention can be referenced to U.S. Pat. No. 8,095,484 B2. Embodiments of the present invention can be combined with the VM system disclosed by U.S. Pat. No. 8,095,484 B2. U.S. Pat. No. 8,095,484 B2 is hereby incorporated by reference. The RI value is designed to gauge the reliance level of a virtual metrology value. The GSI value is used to assess the degree of similarity between the current set of input process data and all of the sets of process data used for building and training a conjecture model. The GSI value is provided to help the RI value gauge the reliance level of VM system. The DQIX value is used to evaluate whether a set of process data used for producing a workpiece is abnormal, and the DQIy value is used to evaluate whether the metrology data of the workpiece are abnormal.
Hereinafter, theoretical bases related to the RI value (RI model), the GSI value (GSI model), the DQIX value (DQIX model) and the DQIy value (DQIy model) are first explained.
The RI and GSI are used to learn in real time whether the VM value is reliable. The RI model is used to compute a RI value between 0 and 1 by analyzing the process data of the process tool, thereby determining whether the virtual metrology result can be trusted. The GSI model is used to compute the GSI value for the process. The GSI value is defined as the degree of similarity between the current set of input process data and all of the sets of the process data used for building or training the models.
Referring to Table 1, n sets of historical data are assumed to be collected, including process data (Xi, i=1, 2, . . . , n) and the corresponding actual measurement values (yi, i=1, 2, . . . , n), where each set of process data contains p individual parameters (from parameter 1 to parameter p), namely Xi=[xi,1, xi,2, . . . , xi,p]T. Additionally, (m-n) sets of process data in actual production were also collected, but no actual measurement values are available besides yn+1. That is, only the first among (m-n) pieces of the products is selected and actually measured. In the current manufacturing practice, the actual measurement value yn+1 obtained is used to infer and evaluate the quality of the (m−n−1) pieces of the products.
As shown in Table 1, y1, y2, . . . , yn are historical measurement values, and yn+1 is the actual measurement value of the first piece of the products being manufactured. Generally, a set of actual measurement values (yi, i=1, 2, . . . , n) is a normal distribution with mean μ and standard deviation σ, namely yi˜N(μ,σ2).
All the actual measurement values can be standardized in terms of the mean and standard deviation of the sample set (yi, i=1, 2, . . . , n). Their standardized values (also called z scores) Zy
wherein
The explanation herein adopts a neural-network (NN) algorithm as the conjecture algorithm for establishing the conjecture model performing virtual measurement, and uses such as a multi-regression (MR) algorithm to be the reference algorithm for establishing the reference model that serves as a comparison base for the conjecture model. However, the present invention can also apply other algorithms to be the conjecture algorithm or the reference algorithm, such as a back propagation neural network (BPNN) algorithm, a general regression neural network (GRNN) algorithm, a radial basis function neural network (RBFNN) algorithm, a simple recurrent network (SRN) algorithm, a support vector data description (SVDD) algorithm, a support vector machine (SVM) algorithm, a multiple regression (MR) algorithm, a partial least squares (PLS) algorithm, a nonlinear iterative partial least squares (NIPALS) algorithm or generalized linear models (GLMs), provided the reference algorithm differs from the conjecture algorithm, such as a SVM algorithm and other related algorithms, and thus the present invention is not limited thereto.
When the NN and MR algorithms are utilized, if their convergence conditions both are that SSE (Sum of Square Error) is minimized with n→∞, their standardized predictive measurement values (defined as Zy
with respect to the NN conjecture model differ from the standardized mean-estimating equation
and standard-deviation-estimating equation
with respect to the MR reference model.
The RI is designed to gauge the reliance level of the virtual metrology value. The RI thus should consider the degree of similarity between the statistical distribution of the virtual metrology value and the statistical distribution Zy
estimated by the reference algorithm which is such as the MR algorithm to replace Zy
Referring to of the reference prediction value from the reference model (built by such as the MR algorithm). As such, the RI equation is listed below:
and σ is set to be 1.
The RI increases with increasing overlap area A. This phenomenon indicates that the result obtained using the conjecture model is closer to that obtained from the reference model, and thus the corresponding virtual metrology value is more reliable. Otherwise, the reliability of the corresponding measurement value reduces with decreasing RI. When the distribution estimated from:
is fully overlapped with the distribution
estimated from
then according to the distribution theory of statistics, the RI value equals 1; and, when those two distributions are almost separate, the RI value approaches 0.
Hereinafter, the method for calculating the statistical distribution of the virtual metrology values ( and
) from the conjecture model is explained.
In the NN conjecture model, if the convergence condition is to minimize SSE, then it can be assumed that “for given is the distribution with mean
and standard deviation θZ
where the NN estimating equation of
and the NN estimating equation of σZ
Before the NN conjecture model is constructed, the process data must be standardized. The equations for standardizing the process data are presented below:
wherein
The n sets of standardized process data (Zx
Accordingly, the estimated value of μZy
and the estimated value of σZ
wherein
Hereinafter, the method for calculating the reference predication values
from the MR model is explained.
The basic assumption of the MR is that “for given
is the distribution with mean μZ
wherein the MR estimating equation of μZ
and the MR estimating equation of σZ
To obtain the MR relationship between the n sets of standardized process data (Zx
The least square method can obtain the estimating equation of βr, {circumflex over (β)}r=[{circumflex over (β)}r0, {circumflex over (β)}r1, . . . , {circumflex over (β)}rp]T as
βr=(ZxTZx)−1ZxTZy (14)
Therefore, the MR reference model can be obtained as
Zyr
Hence, during the conjecture phase, after inputting a set of process data, its MR estimating value Zŷ with
After obtaining the NN estimating equations
and the MR estimating equations
their normal distribution curves can be depicted, as illustrated in
After obtaining the RI, the RI threshold value (RIT) must be defined. If RI≧RIT, then the reliance level of the virtual metrology value is acceptable. A systematic approach for determining the RIT is described below.
Before determining the RIT, it is necessary to define a maximal tolerable error limit (EL). The error of the virtual metrology value is an absolute percentage of the difference between the actual measurement value yi and ŷNi obtained from the NN conjecture model divided by the mean of all the actual measurement values,
The EL can then be specified based on the error defined in equation (18) and the accuracy specification of virtual metrology (VM). Consequently, RIT is defined as the RI value corresponding to the EL, as shown in
with μ and σ defined in equation (4) and
Zcenter=Zŷ
where σy is specified in equation (3).
The following explains a method for building a GS model. As mentioned above, when virtual metrology is applied, no actual measurement value is available to verify the accuracy of the virtual metrology value. Therefore, instead of the standardized actual measurement value Zy
The GSI assesses the degree of similarity between any set of process data and the model set of process data. This model set is derived from all of the sets of historical process data used for building the conjecture model.
The present invention may utilize a statistical distance measure, such as Mahalanobis distance, to quantify the degree of similarity. Mahalanobis distance is a distance measure introduced by P.C. Mahalanobis in 1936. This measure is based on correlation between variables to identify and analyze different patterns of sample sets. Mahalanobis distance is a useful way of determining similarity of an unknown sample set to a known one. This method considers the correlation of the data set and is scale-invariant, namely it is not dependent on the scale of measurements. If the data set has high similarity, the calculated Mahalanobis distance calculated will be relatively small.
The present invention uses the calculated GSI (applying Mahalanobis distance) size to determine whether the newly input set of process data is similar to the model set of process data. If the calculated GSI is small, the newly input set is relatively similar to the model set. Thus the virtual metrology value of the newly input (high-similarity) set is relatively accurate. On the contrary, if the calculated GSI is too large, the newly input set is somewhat different from the model set. Consequently, the virtual metrology value estimated in accordance with the newly input (low-similarity) set has low reliance level in terms of accuracy.
The equations to calculate the standardized process data Zx
Assuming that the correlation coefficient between the s-th parameter and the t-th parameter is rst and that there are k sets of data, then
After calculating the correlation coefficients between the standardized model parameters, the matrix of correlation coefficients can be obtained as
Assuming that the inverse matrix (R−1) of R is defined as A, then
Hence, the equation for calculating the Mahalanobis distance (Dλ2) between the standardized λ-th set process data (Zλ) and the standardized model set process data (ZM) is as follows.
Finally, equation (25) is obtained.
The GSI of the standardized λ-th set process data is, then, equal to Dλ2/p.
After obtaining the GSI, the GSI threshold (GSIT) should be defined as the following:
GSIT=a*
In the so-called cross validation's leave-one-out (LOO) method, one sample data set is selected from all sets of process data used for building models as a simulated on-line test sample set, and then the remaining sets of process data are utilized to construct a GS model, and thereafter the GSI model is used to compute a GSI value for the test sample set, i.e. GSILOO. The above steps are repeated on all of the sample data sets (process data) used for building models, thus computing all of the GSILOO values of the respective sample data sets. Therefore,
A method for constructing a DQIX model is described as follows. Assume that n sets of historical process data are received for constructing the first DQIX model, wherein each set of historical process data is composed of p parameters. These n sets of historical process data are applied to generate p eigenvectors with p corresponding eigenvalues (λ1≧λ2≧ . . . ≧λp) in descending order by the principal component analysis (PCA). Then, a set of k significant eigenvalues (with λ≧1) is selected for constructing a feature-extraction matrix M, which is expressed as:
The method for computing the DQIX value is explained as follows.
At first, equation (28) is applied to transform the ith input set of process data Xi to k data feature variables Ai[a1, a2, . . . , ak].
Ai=M·Xi (28)
Then, these k data feature variables are transformed to ZA=[Za
where i represents the ith input set of process data;
a
: the mean of the j-th standardized variable of the training samples.
Theoretically, the value of
Meanwhile, the cross validation's leave-one-out (LOO) method is used to determine process data quality threshold (DQIX
DQIX
In the so-called cross validation's leave-one-out (LOO) method, one sample data set is selected from all sets of process data used for building models as a simulated on-line test sample set, and then the remaining sets of historical process data are utilized to construct a DQIX model, and thereafter the DQIX model newly built is used to compute a DQIX value for the simulated on-line test sample set, i.e. DQIX
It is noted that the feature-extraction matrix M and the DQIX compose a DQIX model, and the DQIX model will be updated (in the model-refreshing procedure) in accordance with a condition for retraining or tuning.
Thereafter, Z-score values of the historical process data are computed. Then, a DQIy model is created, in which the DQIy model is composed of m similar patterns
In the present embodiment, the m similar patterns {P1, P2, . . . , Pm} are sorted from Z-score values of those n sets of historical process data by applying adaptive resonance theory 2 (ART2) with ρ=0.98.
The method for computing the DQIy values is described as follows. At first, when a new actual measurement value yi is collected, the Z-score value ZX
where
The DQIy
yi=
where Rmax is the maximal-tolerable variance;
Rmax=max(RP
where RP
By adding yi into the similar pattern Pq, the DQIy
After obtaining the DQIy
The aforementioned PCA, LOO, ART2, Z-score and ED algorithms all are known to one having ordinary skill in the art, and thus the details thereof are not described herein.
Referring to
Scenario 1: a stable process is considered. When no status-changes of a process tool occur and all of the RI, GSI, DQIX, and DQIy values of a workpiece processed by the process tool are within their respective thresholds in a manufacturing process, i.e. RI≧RIT; GSI≦GSIT; DQIX≦DQIX
Scenario 2: the status of the process tool is changed. A possible status change of the process tool may occur when tool maintenance, repair, or recipe adjustment is performed. At this situation, embodiments of the present invention will request an actual measurement for the workpiece which is being processed when the status of the process tool is changed, and reset the workpiece count to 0. For example, if the workpiece originally expected to be measured is the 25th workpiece, and the status of the process tool is changed when the second workpiece is being processed, embodiments of the present invention performs actual metrology on the second workpiece, and the next workpiece expected to be measured is the 27th workpiece.
Scenario 3: the DQIX value of the workpiece is abnormal (i.e. DQIX>DQIX
Scenario 4: the GSI value or the RI value of the workpiece is abnormal (GSI>GSIT or RI<RIT). The purpose of the RI value is to gauge the reliance level of a VM value. If the RI value of the workpiece is abnormal (RI<RIT), it represents that the degree of reliability of the VM value of the workpiece is poor, and the actual measurement value of the workpiece is needed to retune or retrain the VM models. The purpose of the GSI value is to evaluate deviations of process data. A process-data deviation of the workpiece may result in a deviation of its corresponding actual measurement value. As such, the workpiece with the abnormal GSI value needs to be inspected. However, if the abnormal RI value or GSI value just happens once, then this may be a false alarm caused by noise or other factors. To confirm that a real deviation is detected, when at least a certain number of consecutive workpieces (such as two or four) have abnormal RI or GSI values, embodiments of the present invention then perform actual metrology on the second or fourth workpiece.
Scenario 5: the DQIy value of the workpiece is abnormal (i.e. DQIy>DQIy
Hereinafter, a flow process of the metrology workpiece sampling step 140 is explained, wherein the metrology workpiece sampling step 140 includes a sampling step 200 shown in
In the sampling step 200, at first, a workpiece is provided to a process tool (step 202), in which the process tool has a set of process data used to process the workpiece. The set of process data of the workpiece is inputted into the DQIX model and the GSI model built in the model-creation step 130, thereby obtaining a DQIX value, a GSI value and a RI value of the workpiece (step 204). In step 204, the set of process data of the workpiece is also inputted into the conjecture model built in the model-creation step 130, thereby computing a virtual metrology (VM) value of the workpiece; and the set of process data of the workpiece is also inputted into the reference model built in the model-creation step 130, thereby computing a reference prediction value of the workpiece. Thereafter, an overlap area between the statistical distribution of the virtual metrology value of the workpiece and the statistical distribution of the reference prediction value of the workpiece is calculated to generate a RI value of the workpiece, in which the RI increases with increasing overlap area, representing that the corresponding virtual metrology value is more reliable.
Thereafter, step 210 is performed to check if the process tool has been idled for a period of time (for example, idled for more than one day). When the checking result of step 210 is true (“yes”), the workpiece is regarded as the first workpiece after the idle period and step 292 has to be performed to perform metrology on the workpiece by using a metrology tool and set a workpiece count to 0, thereby confirming if the process tool is normal. When the checking result of step 210 is false (“no”), the workpiece count is increased by 1 (step 220). Then, step 230 is performed to check if the workpiece count is greater than or equal to N+p, where p is a positive integer (such as 5), and N is used to define the predetermined workpiece sampling rate 1/N. For preventing too many workpieces from skipping measurement, a safety threshold is set in step 230 to force a workpiece measurement to be performed after a certain amount of workpieces are skipped from measurement. Certainly, step 230 is optional, and embodiments of the present invention may omit this step. When the checking result of step 230 is true (“yes”), the metrology tool is used to perform metrology on the workpiece and the workpiece count is set to 0 (step 292). When the checking result of step 230 is false (“no”), step 240 is performed to check if a status change of the process tool occurs (for example, when tool maintenance, repair, or recipe adjustment etc. is performed.). When the checking result of step 240 is true (“yes”), the metrology tool is used to perform metrology on the workpiece by and the workpiece count is set to 0 (step 292), thereby confirming if the process tool is normal.
When the checking result of step 240 is false (“no”), step 250 is performed to check if the DQIX value of the workpiece is good or bad. When the DQIX value is greater than the DQIX threshold, it represents that the quality of the set of process data of the workpiece (the DQIX value) is not good (the checking result of step 250 is “bad”). Because the actual measurement value of the workpiece produced by using the set of process data with the abnormal DQIX value is not reliable, the measurement of the workpiece is skipped (step 290). When the DQIX value is smaller than or equal to the DQIX threshold, it represents that the quality of the set of process data of the workpiece is good (the checking result of step 250 is “good”), and step 260 is performed to check if the workpiece count is greater than or equal to N. When the checking result of step 260 is true (“yes”), it represents that the workpiece is the workpiece expected (scheduled) to be measured, metrology should be performed on the workpiece and the workpiece count is set to 0 (step 292).
When the checking result of step 260 is false (“no”), step 270 is performed to check if the GSI value and the RI value of the workpiece are good or bad. When the GSI value of the workpiece is smaller than or equal to the GSI threshold and the RI value of the workpiece is greater than or equal to the RI threshold, it represents that the virtual metrology value conjectured by using the set of process data of the workpiece is reliable (the checking result of step 270 is “good”), and thus the workpiece does not need to be measured (step 290). When the GSI value of the workpiece is greater than the GSI threshold or the RI value of the workpiece is smaller than the RI threshold, it represents that the virtual metrology value conjectured by using the set of process data of the workpiece is not reliable (the checking result of step 270 is “bad”), the workpiece may need to be measured. However, if the abnormal RI value or GSI just happens once, it may be a false alarm caused by noise or other factors, and thus, when at least a certain number of consecutive workpieces (such as two or four) have abnormal RI or GSI values, embodiments of the present invention then perform: actual metrology on the second or fourth workpiece. In other words, when the checking result of step 270 is “bad”, step 280 is performed to check if the GSI values of k workpieces (such as the previous one or three workpieces) processed before the workpiece all are greater than the GSI threshold or the RI values of k workpieces processed before the workpiece all are smaller than the RI threshold, where k is a positive integer. When the checking result of the step 280 is true (“yes”), the metrology tool is used to perform metrology on the workpiece and the workpiece count is set to 0 (step 292). When the checking result of the step 280 is false (“no”), the measurement of the workpiece is skipped (step 290). It is worthy to be noted that step 270 also may only check if the GSI value of the workpiece is good or bad. If the GSI value of the workpiece is too large, it represents that the set of process data of the workpiece have some differences from the sets of process data used for model building, such that the quality of the workpiece is likely to be abnormal, and needs actual metrology. It can be known from the above, with the applications of the embodiments of the present invention, a user may wait until the process tool has processed more workpieces to select a workpiece for measurement, i.e. N can be increased to lower the predetermined workpiece sampling rate 1/N without worrying to skip the measurement of the abnormal workpiece that ought to be measured. Thus, the embodiments of the present invention can effectively lower the predetermined workpiece sampling rate. However, the predetermined workpiece sampling rate also can be effectively lowered by only performing step 250 (checking the DQIX value of the workpiece), step 260 (checking if the workpiece is the workpiece expected to be measured) and step 270 (checking the GSI value and the RI value of the workpiece, or only checking the GSI value of the workpiece).
Besides, after the actual metrology is performed on the workpiece, metrology, the workpiece sampling step 140 is also directed to evaluation of a DQIy value of a workpiece, as shown in step 201 of
It is understood that the metrology sampling method of the present invention are performed by the aforementioned steps. A computer program of the present invention stored on a non-transitory tangible computer readable recording medium is used to perform the metrology sampling method described above. The aforementioned embodiments can be provided as a computer program product, which may include a machine-readable medium on which instructions are stored for programming a computer (or other electronic devices) to perform a process based on the embodiments of the present invention. The machine-readable medium can be, but is not limited to, a floppy diskette, an optical disk, a compact disk-read-only memory (CD-ROM), a magneto-optical disk, a read-only memory (ROM), a random access memory (RAM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), a magnetic or optical card, a flash memory, or another type of media/machine-readable medium suitable for storing electronic instructions. Moreover, the embodiments of the present invention also can be downloaded as a computer program product, which may be transferred from a remote computer to a requesting computer by using data signals via a communication link (such as a network connection or the like).
On the other hand, a metrology sampling method according to the embodiments of the present invention also can be combined with a method for automatic virtual metrology (AVM) disclosed by U.S. Pat. No. 8,095,484 B2. Referring to
It can be known from the aforementioned embodiments that, by using the ISD scheme constructed by various index values representing status changes or abnormal information of a process tool during a production process, the present invention can effectively assure the VM accuracy and greatly lower the workpiece sampling rate.
It will be apparent to those skilled in the art that various modifications and variations can be made to the structure of the present invention without departing from the scope or spirit of the invention, in view of the foregoing, it is intended that the present invention cover modifications and variations of this invention provided they fall within the scope of the following claims.
Number | Date | Country | Kind |
---|---|---|---|
103111289 A | Mar 2014 | TW | national |
Number | Name | Date | Kind |
---|---|---|---|
8095484 | Cheng et al. | Jan 2012 | B2 |
20090292386 | Cheng | Nov 2009 | A1 |
20120029662 | Cheng | Feb 2012 | A1 |
20130257461 | Heidmann | Oct 2013 | A1 |
Number | Date | Country |
---|---|---|
103513159 | Jan 2014 | CN |
I267012 | Nov 2006 | TW |
200745895 | Dec 2007 | TW |
Entry |
---|
D. Kurz, et. al., “Monitoring virtual metrology reliability in a sampling decision system,” Automation Science and Engineering (CASE), 2013 IEEE International Conference; Aug. 17-20, 2013; p. 20-25. |
Jae Yeon Claire Baek, Costas J. Spanos, “Optimization of blended virtual and actual metrology schemes,” Proc. SPIE 8324, Metrology, Inspection, and Process Control for Microlithography XXVI, 83241K (Apr. 5, 2012). |
D. Kurz, C. De Luca and J. Pilz, “A Sampling Decision System for Virtual Metrology in Semiconductor Manufacturing,” in IEEE Transactions on Automation Science and Engineering, vol. 12, No. 1, pp. 75-83, Jan. 2015. |
C. Yugma, J. Blue, S. Sauzere-Peres, A. Obeid, “Integration of scheduling and advanced process control in semiconductor manufacturing: review and outlook,” J. Sched. (2015) 18:195-205. |
D. Kurz, C. DeLuca and J. Pilz, “Sampling Decision System in semiconductor manufacturing using Virtual Metrology,” 2012 IEEE International Conference on Automation Science and Engineering (CASE), Seoul, 2012, pp. 74-79. |
C. F.Chen, F. T. Cheng, C. C. Wu and H. H. Huang, “Preliminary study of an intelligent sampling decision scheme for the AVM system,” 2014 IEEE International Conference on Robotics and Automation (ICRA), Hong Kong, 2014, pp. 3496-3501. |
Bousetta, A., and A. J. Cross., “Adaptive Sampling Methodology for In-Line Defect Inspection.” In Proceedings of IEEE, SEMI Advanced Semiconductor Manufacturing Conference, 25-31. Munich: IEEE, Apr. 11-12, 2005. |
Chien, C.-F., K.-H. Chang, and C.-P. Chen., “Design of a Sampling Strategy for Measuring and Compensating for Overlay Errors in Semiconductor Manufacturing.” International Journal of Production Research, vol. 41, No. 11, pp. 2547-2561, 2003. |
André Holfeld et al., “A Fab-wide APC Sampling Application.” IEEE Transactions on Semiconductor Manufacturing, vol. 20, No. 4, pp. 393-399, Nov. 2007. |
Daniel Kurz et al., “Sampling Decision System in semiconductor manufacturing using Virtual Metrology” In Proceedings of 2012 IEEE International Conference on Automation Science and Engineering (CASE 2012), 74-79. Seoul: IEEE, Aug. 20-24, 2012. |
Daniel Kurz et al., “Monitoring Virtual Metrology Reliability in a Sampling Decision System.” In Proceedings of 2013 IEEE International Conference on Automation Science and Engineering (CASE 2013), 20-25. Madison, WI: IEEE, Aug. 17-21, 2013. |
Lee, J. H., “Artificial intelligence-based sampling planning system for dynamic manufacturing process.” Expert System with Applications, 22:117-133, 2002. |
Justin Nduhura-Munga et al.,“Dynamic Management of Controls in Semiconductor Manufacturing.” In Proceedings of IEEE/SEMI Advanced Semiconductor Manufacturing Conference, 1-6. Saratoga Springs, NY: IEEE/SEMI, May 16-18, 2011. |
Justin Nduhura-Munga et al., “A Literature Review on Sampling Techniques in Semiconductor Manufacturing.” IEEE Transactions on Semiconductor Manufacturing, vol. 26, No. 2, pp. 188-195, May 2013. |
Alan Weber, “Virtual Metrology and Your Technology Watch List: Ten Things You Should Know about This Emerging Technology.” Future Fab International, Issue 22, Section 4, pp. 52-54, 2007. |
Number | Date | Country | |
---|---|---|---|
20150276558 A1 | Oct 2015 | US |