A SYSTEM FOR MONITORING AND ANALYZING MANUFACTURING PROCESSES USING STATISTICAL SIMULATION WITH SINGLE STEP FEEDBACK

Abstract: A statistical simulation of a semiconductor fabrication process is performed in parallel with the actual process. Input parameters derived from a probability density function are applied to the simulator which, in turn, simulates an actual fabrication process which is modeled as a probability density function. Each simulation step is repeated with a random seed value using a Monte Carlo technique, a trial-and-error method using repeated calculations to determine a best solution to a problem. The simulator generates an output in the form of a probability distribution. The statistical simulation uses single-step feedback in which a simulation run uses input parameters that are supplied or derived from actual in-line measured data. Output data generated by the simulator, both intermediate output structure data and WET data, are matched to actual in-line measured data in circumstances for which measured data is available. The probability density structure of the simulator is adjusted after each simulation step so that simulated data more closely matches in-line measured data.

Description: A SYSTEM FOR MONITORING AND ANALYZING MANUFACTURING PROCESSES USING STATISTICAL SIMULATION WITH SINGLE STEP FEEDBACK Field of the Invention The present invention relates to manufactuπng processes such as integrated circuit fabncation processes More specifically, the present invention relates to a feedback system for analyzing and monitormg a manufactuπng process Background of the Invention Process simulation is the usage of processing experiments, typically using a computer, as directed by mathematical models created to describe a process phenomena Many simulation and analysis tools (for example, Pisces, Medici, Suprem3, Suprem4 and PdFab) have been developed to assist process integration and device development These tools have not been as widely employed for integrated circuit manufactuπng Generally, these tools are developed pπmaπly for research and development purposes and do not adequately address vanous difficulties that aπse in the manufacturing environment Several characteristics are generally applicable to the manufacturing environment and distinguish the manufacturing environment from a process integration and device development environment One characteristic of the manufactuπng environment is that measurable aspects of processes are have a fundamentally statistical nature, rather than a deterministic nature Process variations and measurement errors are inherent to manufactuπng processes so that substantially all data measured in a manufacturing environment is statistical Exact measurement values are generally not available for each device at each stage of a manufacturing process so that a single data point is insufficient to justify a decision relating to the process For example, if it is known that application of input parameters A and B to a fabncation process to yield an output vaπable C, what is truly known is that input parameters A and B each have a statistical profile that, when combined in the fabncation process, yield a statistical profile C shown in Figure 1 The most useful information available in a manufacturing environment is the form of the statistical profile which results from the process Unfortunately, nearly all information that is utilized m the manufactunng environment and all input parameters to a fabrication tool are expressed m the form of single-valued parameters, rather than in statistical profiles While conventional simulation and analysis tools do not suitable address the statistical nature of manufactuπng processes, these tools are also deficient m taking advantage of the extensive process variables, in-line measurements and Wafer Electrical Testing (WET) data measurements that are available WET testing includes testing for various device electrical parameters including threshold voltage and drive current that are measured at a wafer level, before bonding. Conventional simulation and analysis tools designed for research and development generally presume that a fabrication process is not yet operational for actual manufacturing. Therefore, these conventional simulation and analysis tools are not sufficiently flexible for a manufacturing engineer to make optimizations of the process. Specifically, these conventional tools do not allow the manufacturing engineer to utilize the extensive statistical data that is available in a manufacturing process to optimize the simulation model. Furthermore, for simulation tools that use either an empirical approach or an analytical approach, the model fitting parameters do not have sufficient degrees of freedom to match the extensive data that is available from an actual manufacturing process. Another characteristic of the manufacturing environment is that monitoring of manufacturing processes and improvement of these processes is a fine-tuning process. Each tuning step includes a measurement of small differences in process variables with these differences being attributable at least in part to statistical fluctuations and also to complicated interactions between multiple reactions of the process as various process parameters are modified. Process results are typically difficult to measure with accuracy. A large number of highly variable factors influence process results. Modification of a single factor in isolation from other factors is difficult. This difficulty arises not only from a limited understanding of a factor's influence on the process but also because die various factors cause complex inter-related cross effects and interactions. Thus, a simulation and analysis tool that a production engineer confidently uses needs to supply a much higher order of measured precision of the data. Mere indications of data trends are insufficient. What is sought in the development of manufacturing tools and techniques is not a drastic change in a fabricated structure, but rather a small adjustment in characteristics. For example, what typically produces an improvement in an integrated circuit structure, such as an LDD structure in a transistor, is a change in dopant dosage of about ten percent or a change in applied temperature of 100°C. The combination of the small size of the adjustments which are achieved by process modifications and the difficulty in measuring results of the modifications accurately manifest a disadvantageous characteristic of the manufacturing environment akin to a poor signal-to-noise ratio in signal processing. Manufacturing simulation tools are calibrated prior to performing a simulation test. Calibration is typically accomplished by entering calibrated input parameters that are generated either experimentally or by previous simulation. In conventional manufacturing process calibration, a specified value of a parameter is fitted to produce a specified process output value. Realistic simulation results are rarely achieved using conventional simulation and calibration techniques since these techniques do not capture the true nature and complexity of the manufacturing process. Furthermore, conventional calibration processes require the intense study of device engineers before a simulation tool becomes useful. Thus, the calibration processes cause significant delay in process qualification and improvement. This problem is worsened by the fact that calibration procedures are repeated continuously as the environment in the manufacturing area changes over time. These characteristics of the manufacturing environment are applicable to an analysis of manufacturing monitoring as well as process improvements. For example, it is often desirable to know how processes change over time to track changes in fabrication results of a small amount, such as 3 percent, over time, for example 3 weeks to yield an ultimate result. Furthermore, these small differences are typically measured in an environment of statistical fluctuation and measurement error. Although many characteristics of the manufacturing environment are disadvantageous, several are advantageous. One advantage is that specification of the structure resulting from the manufacturing process is well defined. Another advantage is that actual in-line measurement data acquired at various stages of the manufacturing process and actual WET data are plentiful and easily available. These data include statistical profile data that are highly informative regarding the manufacturing process. What is sought is a technique for tracking and analyzing manufacturing processes such that inaccuracies arising from statistical fluctuations, complicated interactions, and measurement errors are avoided or compensated so that process modifications that produce even small differences can be measured, monitored and analyzed. Summary of the Invention In accordance with the present invention, a statistical simulation of a semiconductor fabrication process is performed in parallel with the actual process. Input parameters extracted from actual fabrication data and expressed in the form of a probability density function are applied to the simulator which, in turn, simulates an actual fabrication process which is modeled as a probability density function. Each simulation step is repeated with a random seed value using a Monte Carlo technique, a trial-and-error method using repeated calculations to determine a best solution to a problem. The simulator generates an ouφut in the form of a probability distribution. The statistical simulation uses single-step feedback in which a simulation run uses input parameters that are supplied or derived from actual in-line measured data. Ouφut data generated by the simulator, both intermediate output structure data and WET data, are matched to actual in-line measured data in circumstances for which measured data is available. The probability density structure of the simulator is adjusted after each simulation step so that simulated data more closely matches in-line measured data. The essence of the single-step feedback is that one dedicated series of Monte Carlo simulations exists for the calibration of each individual process step, using calibration results from all of the previous steps. For example, if a process mcludes fifty process steps and one hundred Monte Carlo simulation steps are needed to build the calibration model, then the simulation process will include 5000 simulation steps in the calibration process, each starting from the first step to the calibrating step. The statistical simulation includes two types of simulation runs including calibration runs and prediction runs. Both types of simulation runs are statistical simulations which employ a probability density function as a basic elementary unit. The calibration run calibrates simulated output data to match actual in-line measured data. A simulation operator specifies the actual in-line measured data to be employed for the calibration run. For example, in a typical fabrication process, paπicular process parameters including fabrication tool-specific parameters, temperatures, pressures, deposition material concentrations, and the like, are specified independently. Results of calibration runs are compared and parameters are updated as a result of these comparisons so that the process is specialized for selected fabrication characteristics. The described process has many advantages. One advantage is that the described apparatus and method combine statistical analysis of extensive available measurements with feedback to an in-line manufacturing process line. In addition, the automatic calibration feature of the statistical simulation method advantageously complements manual calibration of a fabrication process. A simulation operator maintains the capability to manually calibrate fabrication parameters by simply changing parameter values in the simulator. The automatic calibration feature facilitates manual calibration by furnishing a calibration framework in which specific parameters are modified independently and the results of these modifications are easily tracked. Thus, a whole range of calibration information is acquired that independently tracks modifications in specified parameters using single-step feedback. Another advantage of the disclosed method is that the defect rate of fabricated devices is very low due to the reduction in handling. A further advantage of the method that results from the reduced handling is a substantially reduced manufacturing cost. Brief Description of the Drawings The features of the invention believed to be novel are specifically set forth in the appended claims. However, the invention itself, both as to its structure and method of operation, may best be understood by referring to the following description and accompanying drawings. Figure 1 , labeled Prior Art, is a sequence of graphs which illustrate a black-box analysis of a manufacturing process, showing statistical profiles of two process input parameters which, in combination, yield a statistical profile of an output variable. Figure 2 is a block diagram showing components of a simulation system in accordance with an embodiment of the present invention. Figure 3 is a block diagram illustrating a two-run process of the simulation system shown in Figure 2. Figure 4 is a flow chart of a calibration process of a statistical simulation method in accordance with an embodiment of the present invention. Figure 5 is a flow chart of a measured data highpass filtering operation of the calibration process shown in Figure 4. Figures 6(A), 6(B) and 6(C) are a sequence of graphs which illustrate intermediate results of a measured data highpass filtering operation depicted in Figure 5. Figures 7(A), 7(B), 7(C), and 7(D) are a sequence of graphs which illustrate intermediate results of a specified data highpass filtering operation. Figure 8 is a flow chart which illustrates steps of a matching operation for matching actual data to simulated data. The matching operation is a step of the calibration process described with respect to Figure 4. Figure 9 is a flow chart which illustrates steps of a prediction process of a statistical simulation method in accordance with an embodiment of the present invention. Figure 10 is a flow chart showing elements of an overall structure of a fabrication system in accordance with an embodiment of the present invention. Detailed Description Figure 2 is a block diagram showing components of a simulation system 200, including various fabrication equipment 210, various test equipment 212 for measuring process parameters and a simulation computer system 214. In an actual fabrication process, typically a plurality of fabrication equipment 210 is employed to carry out different steps of the fabrication process. At several stages of the fabrication process, process parameters are acquired using the test equipment 212. Acquired process parameter data are conveyed to the simulation computer system over a transmission path 216. The transmission path 216 is any suitable communication system such as transmission wires or network lines. Other data transmission methods may include manual techniques such as conveyance by magnetic disk. Examples of fabrication equipment 210 include CVD reactors, vacuum pumps, ion implantation equipment, deposition equipment, photolithography equipment, optical aligners, photomasks, wafer cleaning equipment, wet etching equipment and the like. Examples of test equipment 212 include residual gas analyzers, spectral analyzers, electrical probes, optical dosimetry measurement equipment, optical microscopes, laser reflectometry equipment, spectroscopes and similar data acquisition tools that are known in the art of manufacturing. The simulation computer system 214 is any suitable computer system which includes common computational functionality and common communication operations for receiving process data. In one exemplary embodiment, a SUN SPARC 20™ application server with a network of x86 PC interface clients are employed as the computer system 214. A customized PdFab™ simulation tool and conventional software, including Microsoft Excel™ and Microsoft Access™ are used to process data. The simulation system 200 is a simulation and analysis tool for a manufacturing environment that incorporates several objective properties. First, the simulation system 200 provides for simulation and analysis in which data is input, processed and output with an arbitrary statistical distribution profile since actual fabrication data statistics may not be distributed in a Gaussian profile. Second, the simulation system 200 furnishes a capability to utilize any measured data, if available, to improve the accuracy of simulation results. The simulation system 200 also allows data supplementation to achieve a best possible simulation even if some statistical data is not measurable or is otherwise unavailable. Third, the simulation system 200 is flexible and allows models underlying the simulation to be upgraded and improved. Fourth, the simulation system 200 performs substantially all operations promptly and automatically without conflicts with existing data collection and analyzing systems in the fabrication. Referring to Figure 3, a two-run process of the simulation system 200 includes two different run modes of operation, specifically a calibration run 310 and a prediction run 320. In the calibration run 310, available data are processed to generate fitting functions for matching simulated and measured data. The calibration run 310 includes a single simulation step 312 and actual measurements 314. Results of the simulation step 312 and the actual measurements 314 are processed according to a fitting function 316. The calibration run 310 includes a series of single simulation steps 312, each with corresponding actual measurements (both process parameters and in- line/WET data) 314. Results of a plurality of the simulation steps 312 and the actual measurements 314 are processed to produce a fitting function 316. Each simulation step 312 processes data accumulated using all fitting functions previously calculated during the calibration run 310. In the prediction run 320, measured data and the fitting functions calculated in the calibration run 310 are used to predict unknown data. During the prediction run 320, if any process parameter is missing from the existing database, the statistical distribution function used in the corresponding calibration run 310 replaces the missing parameter. If any WET data or in-line parameter is contained in the existing dataset, these WET data or in-line parameters replace computed values without changing fitting functions. Referring to Figure 4, a flow chart of a calibration process 400 of the statistical simulation method is shown. The calibration process 400 includes two parallel processes including an actual in-line process 420 and a simulation process 450. The actual in-line process 420 includes with a wafer start step 422 for initializing process parameters takes place, a single process step 424 for performing a single selected fabrication process step, a measurement step 426 for measuring a result of the selected fabrication process, a process end step 428 for terminating the selected fabrication process and the measure WET data step 430 for measuring WET data parameters. Adjustment of different process steps, which are suitable for simulation analysis, of a plurality of suitable process steps is performed independently in this manner. A process step which is suitable for simulation is a process step that: (1) produces a fabrication result that is important to the functionality or structure of an integrated circuit, (2) produces a measurable result, (3) is capable of simulation, and (4) responds to changes in simulation parameters with differences in function or structure. Examples of suitable process steps include pad oxidation before channel implant, gate oxidation, gate etch and spacer oxide steps. Typically, two types of output data are measured that result from a fabrication process step, including vertical dimension data and horizontal dimension data. Examples of vertical dimension data include oxide growth, gate oxide thickness, polysilicon thickness data and the like. Horizontal dimension data include polysilicon gate length and LDD spacer width, for example. Doping profile data results from processes which are assumed to be stable so that a doping profile is typically known, allowmg simulation on the basis of theoretical estimations of parameter values, rather than measured data In one example of a suitable selected fabrication process, the single process step 424 performs a gate oxidation process and the in-line measurement step 426 measures a data profile of oxide thickness Corresponding steps of the simulation process 450 are performed in parallel with steps of the actual in-line process 420 A simulation start step 452 begins the simulation process 450 in response to initializing data from the actual m-l ine process 420 The wafer start step 422 generates initial data, such as oπentation data, that is measured and transferred to the simulation process 450, typically through a manufactunng control system, such as Workstream™, the remote access channel of the manufactuπng control system, such as Remote Workstream™, and a network connection to the application server, such as TCP/IP A simulation start step 452 initializes parameters of the simulation process 450 to arbitrary, used-defined values Followmg the simulation start step 452, a simulation step 454 simulates the actual process step performed in single process step 424, first usmg arbitrary, user-defined parameters and later adapting the parameter values on the basis of actual in-line measurements Vanous miscellaneous mput parameters such as processing time are designated by the test operator These input parameters are applied to the single process step 424 and the simulation step 454 Input data may be applied in several formats However, the input data is converted into a statistical distnbution function before actual processing begins For an array of input data points, data is sorted and the probability of a data value being between any two consecutive data points is assumed to be the same For data presented in statistical form, such as data with a mean, standard deviation and range limits, the data is modeled in a statistical distribution function as a truncated Gaussian profile for usage as a statistical distribution function For data presented in a statistical form, such as a mean and range limits, the data is modeled m a statistical distnbution function as a truncated Gaussian profile with each specified limit being presumed to deviate from the mean value by three standard deviations If the mean is not centered between the range limits, the function is modeled as an asymmetric profile and is considered the combination of two half-Gaussian profiles that have the same population and different standard deviations. For data presented in a statistical form, such as a mean and standard deviation, the data is modeled in a statistical distnbution function as a Gaussian profile Data presented as a single data pomt is used only for parameters that are insignificant when no additional information is unavailable. Each actual or simulation result, including intermediate results, is applied to the simulation and process as a statistical distnbution function, rather than a single data pomt Thus, a statistical distribution function is the elementary data type m the simulation system 200 Simulation step 454 operates on the basis of statistical simulation in which a simulation step is repeated many times Each single simulation run computes one value from each input statistical distribution function, processes the input values in accordance with a corresponding physical or device model, and obtains a value for each output statistical distnbution function A higher number of simulation repetitions results in improved accuracy in the output distnbution One example of an input processmg parameter is a processing time parameter In addition, vanous actual in-line measurements are acquired duπng the smgle process step Results of these in-line measurements are communicated to the simulation step 454 by a highpass filter step 456 The highpass filter step 456 filters input parameters and in-line measured data from the single process step 424 and performs digitized processmg to generate an output of a discrete probability density function Simulation results from the simulation step 454 are applied to a generate simulated profile step 458 and to a random seed step 460 The random seed step 460 updates a seed value for Monte Carlo technique simulation and applies the updated seed value to the simulation step 454 The generate simulated profile step 458 constructs a simulated profile which is used for matching analysis of actual and simulated results The simulation takes place in a plurality of Monte Carlo steps, each step mcludmg one pass through the simulation and random seed steps 454 and 460 Multiple-pass feedback steps that are performed in the simulation and random seed steps 454 and 460 are illustrative of a single step feedback concept for performing calibration of the simulation system 200 Each single simulation step is evaluated independently of other simulation steps In one simulation step, data are measured and applied to a simulation step and results of the simulation are fed back to the process Simulation results are used only to compute the fittmg function for data measured within the simulation step Using the single-step feedback concept, measured data is flexibly used to improve simulation accuracy while the interaction between different simulation steps and propagation error are reduced The single-step feedback concept is substantially different from conventional simulations in which an entire simulation is completed in each simulation run with intermediate output values of all simulation steps bemg assembled and the output fitting function being calculated all at one time Subsequent to the smgle process step 424, the in-line measurement step 426 measures the data profile of oxide thickness and communicates the data profile result to an actual data to simulated data compaπson step 462 The actual data to simulated data comparison step 462 receives simulated data from the generate simulated profile step 458 and actual in-line data from the in line measurement step 426, and compares these data profiles digital sample by digital sample, for example by calculating the point-by-point difference of the profiles The difference profile is stored in store difference profile step 464 A determine matching process to convert step 466 determines whether to convert simulated profile data usmg previous simulated data or actual data on the basis of the difference profile determmed by the data comparison step 462 A conversion operation is performed on the simulated profile in matching process conversion step 468 which iteratively converts the simulated profile received from the generate simulated profile step 458 on the basis of the difference profile received from the data compaπson step 462 The matching process conversion step 468 combines the difference profile and the simulated profile and produces a conversion function. The conversion function is then used to match any arbitrary simulated point to a converted simulated point for latter simulation One way to obtain the conversion function uses Newton's interpolation Thus, for a data range defined by actual measured data, Newton's interpolation is used to define a continuous function based on percentile matching For a data range not defmed by actual measured data, a third order polynomial is used to define the continuous function This matching method furnishes a suitable balance between flexibility and accuracy After each single process step and each simulation step is complete and the simulated profile is generated, additional process steps may be performed and additional simulation process steps including various simulation, highpass filter, generate simulated profile, random seed, actual data to simulated data comparison, store difference profile and determine matching process to convert and matching process conversion steps may be performed for the application of additional fabrication processes. For each single process step of actual in-line fabrication and each coπesponding simulation, a full actual in-line process and a full simulation are performed, using random input parameter data, so that independent processing stages are performed with actual fabrication steps proceeding in parallel with simulated steps. Following all actual process steps, actual processing terminates with an actual process end step 28 and the measure WET data step 430 measures the final output data of the fabrication process and generates a WET data profile. The WET data is processed in a highpass filter step 432 to generate filtered WET data. Following all simulated process steps, simulation terminates with a simulated process end step 470 and a full device simulation is performed by device simulation step 472. The device simulation step 472 incorporates profile information derived in previous simulation steps and utilizes Monte Carlo analysis using random seed simulations. A generate device simulation profile step 473 generates a device profile that corresponds to the actual calculated WET data. A device profile matching step 474 compares the WET data to the device data generated by the device simulation step 472 digital sample by digital sample, for example by calculating the point-by-point difference of the profiles. This device difference profile is stored in store device difference profile step 476. A determine matching process to convert WET data step 478 determines whether to convert simulated profile data using previous simulated WET data or actual WET data on the basis of the difference profile determined by the device profile matching step 474. A conversion operation is performed on the simulated WET data in WET statistical profile conversion step 480 to generate the final output profile of the calibration process 400. Output data that is generated by the calibration process 400 includes all of the matching profiles and all simulated profiles. In some embodiments of the calibration process 400, the actual in-line process 420 and the simulation process 450 are completely automatic, being performed in parallel, simultaneously and in real time. In other embodiments, the actual in-line process 420 and the simulation process 450 are performed in parallel but not in real time so that the various actual process steps and simulation steps are performed generally independently but with data acquired from measurements on the actual in-line process 420, filtered off-line in steps such as the highpass filter step 456, and communicated to the simulation process 450, off-line. In some embodiments of the calibration process 400, in-line data extraction from the actual in-line process 420 is automated and the highpass filtering in steps such as the highpass filter step 456 is operated for pre-selected samples of actual in-line data automatically and periodically. Referring to Figures 6(A), 6(B) and 6(C) in conjunction with a flow chart shown in Figure 5, a highpass filtenng operation for highpass filtering measured data is descnbed graphically The highpass filtering operation is performed in the highpass filter steps 432 and 456 shown Figure 4, for example Figure 5 is a flow chart depicting steps of a highpass filter operation 500 for filtering measured data Data is measured m measurement step 510 and measured data is counted in each of a plurality of quantized enor ranges δx m fill histogram step 512 Data is in the form of process parameters, in-line measurements and WET data Generally, process parameters mclude any suitable and meaningful process parameters that may be measured, observed or specified duπng a fabncation process Suitable process parameters include furnace temperatures and deposition times, for example In-line measurements mclude any measurements of process results such as oxide thickness for a gate oxidation process step, polysilicon thickness for deposition steps, and the like WET data includes electncal parameter measurements such as threshold voltage (Vth), saturation threshold voltage (Vsat), maximum transconductance (Gm) and saturation current (Idsat) Process steps of pad oxidation before channel implant, gate oxidation, gate etch and spacer oxide steps are considered most influential for controlling these WET data parameters A process parameter designates a specified parametnc value which serves as an input parameter to a process The in-line measurements and WET data are measured as results or applied as input parameters to the process Figure 6(A) shows a histogram of actual in-line data measurement values shown on the honzontal axis in a quantized eπor range δx and the number of data points falling withm each δx range of measurement values on the vertical axis The highpass filteπng operation described with respect to Figures 6(A), 6(B) and 6(C) operates on measured and observed data Referring again to Figure 5, enor quantization step 512 modifies quantization of the error term όx mto a step size όy in a predetermined manner so that each όy is the largest value which is less than όx and equal to the difference between the maximum and mmimum value x divided by an arbitrary integer n Figure 6(B) shows a histogram of actual in-line parameter values shown on the horizontal axis in terms of the n groupmgs of quantized step size όy These values are normalized in data normalize step 514 so that each histogram value indicates a percentage of the total number of data points falling within each step size range Denve probability density function step 516 calculates a probability density function of the actual in-line parameter values An example of the probability density function is shown in Figure 6(C) Generate probability density table step 518 assembles a table of n probability density function values in ascending order of steps όy Refernng to Figures 7(A), 7(B), 7(C) and 7(D), a highpass filtering operation for highpass filtering specified data is descnbed graphically For specified data, a probability density function is derived theoretically rather than expenmentally so that no data is processed Specified data is specified by vanous statistical parameters In one example, specified data is specified by a mean parameter value μ, a standard deviation value p, a maximum parameter value and a minimum parameter value so that a probability density function is fitted to a Gaussian distnbution as is shown in Figure 7(A) In another example, specified data is specified by a mean parameter value μ and a standard deviation value p alone so that a probability density function is fitted to a Gaussian distnbution as is shown in Figure 7(B) In a further example, specified data is specified by a mean parameter value μ, a maximum parameter value and a mmimum parameter value so that a probability density function is fitted to an asymmetric Gaussian distnbution as is shown in Figure 7(C) In another example, specified data is specified by a mean parameter value μ alone so that a probability density function is estimated to a Gaussian distnbution as is shown in Figure 7(D) Generally, a more suitable probability density function is calculated when more parameters are supplied as opposed to fewer parameters Referring to Figure 8, a matching operation 800 for usage in the actual data to simulated data companson step 462 and the device profile matching step 474, both shown m Figure 4, is illustrated by a flow chart First, actual data is highpass filtered in the manner descnbed with respect to Figure 5 Specifically, the matchmg operation 800 includes the fill histogram step 512, the error quantization step 512, the data normalize step 514, the denve probability density function step 516 and the generate probability density table step 518 In a match profile step 820, the probability density table for simulated data generated in a step which generates simulated data such as generate simulated profile step 458 and the generate device simulation profile step 473 is fitted to the probability density table for the actual data generated by step 518 Profile fittmg is performed by forcing the step sizes of the simulated profile δy' so that the percentage of the total number of data points falling withm each vaπable step size range of the simulated profile matches the percentage of data points m each fixed step size range of the actual data profile In this manner, a difference profile is generated which includes a 4 times n array of step samples. For each of the n step size samples, four elements are stored including an amplitude number indicative of the percentage of total data pomts for each step element, the ascending fixed-size δy steps for the actual data, the ascending vaπable-size όy' steps for the fitted simulated data and elements indicative of the difference step sizes δy - δy * In a determine best fit step 822, a best fit function for the amplitude number is derived using a selected fitting function, as is known m the art, such as a polynomial fit function, a splme function or the like In one example of a best fit function, a constant function is applied and residue (difference) calculated, then a linear polynomial fit is applied and residue calculated If the difference is residue is less than a designated amount, a constant function furnishes a best fit If neither the constant nor the linear polynomial functions supply a best fit, a second order polynomial function is applied Similarly, higher order polynomial functions are applied until a best fit is found Generally, the best simulator is a simulator in which the best fit is a constant function or, if a constant function does not yield a best fit, a lower order polynomial function One advantage that is achieved by the matching operation 800 is a high simulation accuracy based on the automatically-obtained precise fitting function between a calibration output statistical distribution and actual data This operation statistically fits two data populations for which the entire percentile distribution is known so that the fitting function for a range of parameters and data measurements is simply calculated by a direct percentile correlation This means, for example, that 10%, 35 % , and 75% of the output statistical distribution function is mapped to 10%, 35 %, and 75% of the actual data Data values between two data points is accurately measured by interpolation A further advantage of the matching operation 800 is that the process is flexible and applicable for substantially any profile for any number of data samples Another advantage is that the fitting function applies to out-of-range data points. A further advantage is that low pass and high pass filters may be applied to the resulting data to compensate for measurement inaccuracy. The matching operation 800 also handles out-of-range data and predicts possible outcomes that result from application of experimental process parameters that have values outside conventional limits. Therefore, the matching operation 800 extends beyond the region defined by calibration data. From numerical analysis theory, the higher the extrapolation function order, the larger the possible error. To compromise me accuracy and flexibility of the fitting function, a third-order polynomial is used for the fitting function in a region not defined by calibration data. Following the calibration process, a prediction process is performed. The operator specifies input parameters for the prediction process including any operator-selected input parameters or profiles. Any input parameters that are not specified default to parameters or profiles that are derived in the calibration process. Referring to Figure 9, a prediction process 900 begins with a simulation start step 910 in which operator-specified input parameters and profiles are entered and default input parameters and profiles are installed. In user-defined simulation step 912, a simulation is performed and profile derived using user-specified data and data from a processed table. Following the simulation step 912, a simulated profile is generated in simulated profile step 914. A matching conversion step 916 combines the difference profile and the simulated profile and produces a conversion function. A converted structure is derived in step 918. Various process parameters are simulated in this manner until a process end step 920 terminates the prediction simulation. A device simulation is then performed in device simulation step 922. A profile generated from the device simulation is then matched in matching to convert step 924, thus generating simulated WET data. The WET data derived from device simulation is determined in precisely the same manner as the in-line data after process simulation, using the same matching and conversion techniques. One advantage of the simulated process method is that simulated and actual data are all treated uniformly and in a statistical, ra er than deterministic, manner. Another advantage is that all actual in-line data and WET data are available for usage in improving the accuracy of the simulation. Furthermore, for an improved-stability fabrication process, data fitting is substantially automatic. In addition, physical calibration of the process is not necessary so that automatic calibration is possible by simulation is achieved before a fabrication tool is used in manufacturing. An additional advantage is that simulated data achieve improved parameter accuracy. Referring to Figure 10, an overall structure of a fabrication system 1000 is shown in block diagram form. In this system 1000, a fabrication process 1010, as circuits are fabricated, generates actual in-line data that is assembled in in-line data block 1012 and information for derivation of WET data measurement step 1016. This WET data is assembled in WET data block 1018. The in-line data and WET data are entered into a database 1020. Data in the database 1020 are processed in automatic data extraction step 1022 at scheduled intervals such as weeks. These data are assembled in data process block 1024 and applied to a calibration process 1026. Matching procedures are implemented in matching function block 1028 to derive profiles and data which are entered into the database 1020 The profiles and data are applied to a prediction process 1030. Selected mfotmation in the database 1020 are also manually extracted in manual extraction block 1032 and this data along with operator-supplied data 1034 are applied to a manual data processmg block 1036 Data denved in the data processmg block 1036 are also applied to the prediction process 1030 A calibration process 1038 is also performed using the data from manual data processmg block 1036 Information from the calibration process 1038 is used m a manual matching function block 1040 to derive profiles and data resulting from the manual extraction of selected data, which are also entered for performance of the prediction process 1030 The prediction process 1030 generates prediction information for various applications in block 1042 One practical use of a calibration run is to simulate a fabrication process, holding all parameters constant, but taking actual in-line data measurements over time. For example, m some embodiments of the method, data may be measured weekly over multiple weeks to detect shifting in process properties over time Other applications for statistical simulation include momtonng of manufactuπng processes, process development and improvement and comparison of data for different manufactuπng processes The description of certain embodiments of this invention is intended to be illustrative and not limiting Numerous other embodiments will be apparent to those skilled in the art, all of which are included withm the broad scope of this invention What is claimed is*

Claims: CLAIMS 1 A manufactuπng monitoπng system comprising a fabrication equipment for performing a fabrication process step applied to a workpiece; a test equipment for acquiring a sample of a process test parameter sensed from the workpiece; a computer coupled to the test equipment to receive the process test parameter sample; 5 a software program operable upon the computer system, the software program having a plurality of routines including a routme for receiving the process test parameter sample, a routine for simulating a simulated fabrication process step which coπesponds to the fabrication process step to generate a plurality of simulated fabrication process results; and 10 a routine for generating a simulated profile indicative of the simulated fabncation process results 2 A system accordmg to Claim 1, wherein the software program further includes a Monte Carlo simulation routine for iteratively updating simulation results usmg a random seed value 3. A system according to Claim 1 , wherein the software program further includes a routine for formatting a plurality of process test parameter samples mto a process data profile; and a routme for matching the process data profile to the simulated profile 4 A system accordmg to Claim 1 , wherein the software program simulating routine is a statistical simulation of a semiconductor fabrication process further compnsmg. a subroutine for deπving an input parameter to the simulation from an actual fabrication process parameter, the input parameter being derived from a probability density function; and 5 a subroutine for simulating an actual fabrication process, the actual fabrication process being modeled as a probability density function 5 A system accordmg to Claim 4, wherein the software program simulating routme is performed in parallel wim the fabrication process step. 6 A system accordmg to Claim 4, wherein the software program simulating routme is a statistical simulation of a semiconductor fabrication process further compnsmg a subroutine for repeating a plurality of simulation steps, the repeated steps applying a varying random seed value that is vaned using a Monte Carlo technique 7 A system accordmg to Claim 4, wherein the software program simulating routme is a statistical simulation of a semiconductor fabncation process further compnsmg a simulation loop routine usmg smgle-step feedback in which a simulation step of a plurality of simulation steps uses input parameters that are supplied from actual in-line measured data and generates 5 output data, a subroutine for matching output data generated by the simulation step against actual in-line measured data, and a subroutine for adjusting a probability density structure of the simulator after each simulator step so that the simulated data more closely matches the actual in-line measured data 8 A method of momtonng a fabrication process comprising the steps of performing a fabrication process step applied to a workpiece, acquiring a sample of a process test parameter sensed from the workpiece, simulating a simulated fabncation process step which conesponds to the fabncation process step to 5 generate a plurality of simulated fabrication process results; and generating a simulated profile indicative of the simulated fabncation process results 9 A method according to Claim 8, wherein the simulating step further includes the step of performing a Monte Carlo simulation routine for iteratively updating simulation results using a random seed value 10 A method according to Claim 8, further comprising the steps of formatting a plurality of process test parameter samples into a process data profile, and matching the process data profile to the simulated profile 11 A method according to Claim 8, wherein the simulating step is a statistical simulation step of a semiconductor fabncation process, the simulating step further compnsmg the steps of denvmg an input parameter to the simulation from an actual fabrication process parameter, the input parameter being derived from a probability density function; and 5 simulating an actual fabrication process, the actual fabrication process bemg modeled as a probability density function 12 A method according to Claim 11 , wherein the simulating step is performed in parallel with the fabncation process step 13 A method according to Claim 11, wherein the simulating step is a statistical simulation of a semiconductor fabncation process, the simulating step further compnsmg the steps of repeating a plurality of simulation steps, the repeated steps applying a varying random seed value that is vaned using a Monte Carlo technique 14 A method accordmg to Claim 11 , wherem the simulating step is a statistical simulation of a semiconductor fabncation process, the simulating step further compnsmg the steps of looping through a simulation loop usmg smgle-step feedback in which a simulation step of a plurality of simulation steps uses input parameters that are supplied from actual in-line measured data and generates output data, matching output data generated by the simulation step against actual m-line measured data, and adjusting a probability density structure of the simulator after each simulation step so that the simulated data more closely matches the actual in-line measured data

Inventors: CHEN MING CHUN, CHEN, MING, CHUN
Applicants: ADVANCED MICRO DEVICES INC, ADVANCED MICRO DEVICES, INC.
Application year: 1996
Main classification: H01L21
IPC codes: H01L21