Statistical Process Control and Computer Integrated Manufacturing Run to Run Control, Real-Time SPC, Computer Integrated Manufacturing. 1 The Equipment Controller Today, the operation of individual pieces of equipment can be streamlined with the help of external software applications. SPC is just one of them. CIM database Maintenance Workcell Coordinator Monitoring Statistical Process Control Equipment BCAM Fault Diagnosis Local Database(s) Modeling and Simulation Recipe Generation 2
The Workcell Controller Most process steps are so interrelated that must be controlled together using feed forward and feedback loops. Crucial pieces of equipment must by controlled by SPC throughout this operation. Adaptive Control Equipment Model Adaptive Control Equipment Model test Step test Step test 3 Model Based Control All actions are based on the comparison of response surface models to actual equipment behavior. Malfunction alarms are detected using a multivariate extension of the regression chart on the prediction residuals of the model. Control alarms are detected with a multivariate CUSUM chart of the prediction residuals. Control limits are based on experimental errors as well as on the model prediction errors due to regression. Hard limits on equipment inputs are also used. 4
The Idea of Statistically Based Feedback Control response 1. Original Operation 2. Process Shifts 3. Shift is detected by MBSPC 4. RSM is adapted 5. New Recipe is generated y n 2 4 3 y f y o 1 5 new RSM original RSM x o x f controlling input 5 Feedback Control Model-based, adaptive Feedback Control has been employed on several processes. Initial Settings Parameter Estimator Equipment Model yes Recipe Update Controller t Test Model Update? No Spin Coat & Bake M 6
Feedback Control in Resist Application Adaptive Controller Input Settings Incoming Wafer Spin-Coat & Bake Equipment Thickness & Reflectance Measurement Outgoing Wafer 7 Feedback Control in Resist Application (cont.) 12600 12400 12200 12000 Target Model Prediction Experimental Data 11800 0 10 20 30 40 50 60 44 40 36 32 28 24 20 Model Prediction Experimental Data 0 10 20 30 40 50 Wafer Number Target 60 8
Alarms in Resist Application Control 20 16 12 8 4 Alarm (a) Alarm (b) Alarm (c) 0 0 10 20 30 40 50 60 15 10 5 Alarm Alarm 0 0 10 20 30 40 50 60 Wafer Number Alarm UCL UCL 9 Adapting the Regression Model The regression model has many coefficients that may need adaptation. What can be adapted depends on what measurements are available. 2 3 1 4 y n = a o x 2 +b o x +c n y n = a n x 2 +b n x +c n y n = a o x 2 +b n x +c n y o = a o x 2 +b o x +c o 10
Adapting the Regression Model (cont) In multivariate situations it is often not clear which of the gain or higher order variables can be adapted in the original model. y x 2 A B C In these cases the model is rotated so that it has orthogonal coefficients, along the principal components of the available observations. These coefficients are updated one at a time. x 1 11 Model Adaptation in Resist Application Control -13500-14000 -14500-15000 0 10 20 30 40 50 60 2.00 1.95 1.90 1.85 135 134 133 132 131 0 10 20 30 40 50 60 Wafer Number 12
Feed-Forward Control Models can also be used to predict the outcome and correct ahead of time if necessary. Model Projected CD Spread Recipe Update No In Spec? Yes Standard Setting Exposure M Develop 13 Modified Charts for Feed Forward Control LCL LSL σ/ n (prediction spread) σ (equipment spread) UCL USL No FF Yield Loss False Alarm Probability 14
Feed Forward/Feedback Control Results (cont) 92 90 88 86 84 82 80 78 Open Loop FeedForward Target = 88.75% 0 5 10 15 20 25 30 Wafer No 15 Concurrent Control of Multiple Steps Equip. Model Equip. Model Equip. Model Adaptive Control Adaptive Control Adaptive Control Original Inputs Spin Coat & Bake T&R Meas Expose R Meas Develop CD Meas 16
Resist Thickness Example 12000 Target = 11906Å 11500 11000 0 5 10 15 20 Wafer Number Closed Loop Open Loop 17 Resist Thickness Input 100 80 60 Closed Loop Open Loop 40 0 5 10 15 20 Wafer Number 18
100 Latent Image output Closed Loop Control Alarm Open Loop Malfunction Alarm 90 Target = 87.75% 80 70 0 5 10 15 20 Wafer Number 19 Latent Image Input 1.00 Closed Loop Open Loop 0.90 0.80 0.70 0.60 0 5 10 15 20 Wafer Number 20
Developer Output Closed Loop Control Alarm Open Loop Malfunction Alarm 2.80 2.70 2.60 Target = 2.66 um 2.50 0 5 10 15 20 Wafer Number 21 Developer Input 65 Closed Loop Open Loop 60 55 50 0 5 10 15 20 Wafer Number 22
Process Improvement Due to Run to Run Control 5 4 3 2 1 Target Dev. Malf. 5 4 3 2 1 Target 2.55 2.65 2.75 2.85 2.55 2.65 2.75 2.85 CD in µm CD in µm Closed Loop Operation Open Loop Operation 23 Supervisory Control The complete controller must be able to perform feedback and feed-forward control, along with automated diagnosis. Feed-forward control must be performed in an optimum fashion over several pieces of the equipment that follow. max Cpk control equip 24
The Concept of Dynamic Specifications Specs are enforced by a cost function which is defined in terms of the parameters passed between equipment. A change is propagated upstream through the system by redefining specifications for all steps preceding the change. Model Model Constraint Mapping Control Constraint Mapping Control Step Step 25 Specs to step A must change if step B "ages" Out2 A New spec limits for A Out2 B Fixed spec limits for B Mapping PC PC 2 1 Model of B Out1 A Step Step Step A B C Out1 B 26
Specs are outlined automatically - Stepper example Thickness (Develop time varies between 55 and 65 seconds) 27 An Example of Supervisory Control 20 15 Baseline FF/FB only Supervisory Target Target 20 15 Target 20 15 10 5 Outliers 10 5 Outliers 10 5 Outliers 0 0.80.91.01.11.21.31.41.51.61.71.81.9 CD (µm) 0 0.80.91.01.11.21.31.41.51.61.71.81.9 CD (µm) 0 0.80.91.01.11.21.31.41.51.61.71.81.9 CD(µm) 28
Results from Supervisory Control Application Baseline FF/FB Only Supervisory Th PAC LI 29 Control Results - CD Baseline FF/FB Only Supervisory 30
Summary, so far Response surface models can be built based on designed experiments and regression analysis. Model-based Run-to-Run control is based on control alarms and on malfunction alarms. RSM models are being updated automatically as equipment age. Optimal, dynamic specifications can be used to guide a complex process sequence. Next stop: Real-time Statistical Process Control! 31 Autocorrelated Data One important and widespread assumption in SPC is that the samples take random values that are independently and identically distributed according to a normal distribution y t = µ + e t t = 1,2,... e t ~ N (0, σ 2 ) With automated readings and high sampling rates, each reading statistically depends on its previous values. This implies the presence of autocorrelation defined as: ρk = Σ (y i - y)(y i+k - y) Σ(y i - y) 2 = 0 k = 1,2,... The IIND property must be restored before we apply any traditional SPC procedures. 32
Time Series Modeling Various models have been used to describe and eliminate the autocorrelation from continuous data. A simple case exists when only one autocorrelation is present: y t = µ + ϕ y t-1 + e t e t ~ N (0, σ 2 ) µ' = µ / (1- ϕ), σ' = σ / 1-ϕ 2 e t = y t - y t The IIND property can be restored if we use this model to "forecast" each new value and then use the forecasting error (an IIND random number) in the SPC procedure. 33 Example: LPCVD Temperature Readings Temp Readings from LPCVD Tube 608 607 606 605 604 603 602 LPCVD Temp Autocorrelation Temp(t+1) = 758-0.253 Temp(t) 608 607 606 605 604 603 602 601 0 20 40 60 Time 601 80 100 601 603 605 Temp (t) Temps are not IIND since future readings can be predicted! 607 34
The Residuals of the Prediction can be Used for SPC.. Temp. Readings from LPCVD Tube 608 607 606 605 604 603 602 3 2 1 0-1 -2 Temperature Residuals 601 0 20 40 60 Time 80 100-3 0 20 40 60 Time 80 100 35 Estimated Time Series A "Time Series" is a collection of observations generated sequentially through time. Successive observations are (usually) dependent. Our objectives are to: Describe - features of a time series process Explain - relate observations to rules of behavior Forecast - see into the future Control - alter parameters of the model 36
Two Basic Flavors of Time Series Models Stationary data (i.e. time independent mean, variance and autocorrelation structure) can be modelled as: Autoregressive Moving Average z t = ϕ 1 z t-1 + e t e t ~ N (0, σ 2 ) z t = y t - µ z t = θ 1 e t-1 + e t e t ~ N (0, σ 2 ) z t = y t - µ Mixture (i.e. Autoregressive + Moving Average) models. y t = µ + ϕ 1 z t-1 + ϕ 2 z t-2 +...+ ϕ p z t-p + e t - θ 1 e t-1 - θ 2 e t-2 -...- θ q e t-q ϕ: autoregressive θ: moving average 37 First Order Autoregressive Model AR(1) This model assumes that the next reading can be predicted from the last reading according to a simple regression equation. z t = ϕ 1 z t-1 + e t e t ~ N (0, σ 2 ) z t = y t - µ 0.6 0.5 0.4 0.3 0.2 0.1 0.0-0.1-0.2-0.3-0.4-0.5-0.6 35 40 45 50 55 60 AR Forecast sample Error 38
Higher order AR(p) and ACF, PACF representation Higher order autoregressive models are common in engineering. Their structure can be inferredacf and pacf plots. Autocorrelation Function (acf): ρk = Σ (y i - y)(y i+k - y) Σ(y i - y) 2 = 0 k = 1,2,... k Partial Autocorrelation Function (pacf): z t = ϕ 1 z t-1 z t = ϕ 1 z t-1 + ϕ 2 z t-2 z t = ϕ 1 z t-1 + ϕ 2 z t-2 + ϕ 3 z t-3 z t = ϕ 1 z t-1 + ϕ 2 z t-2 + ϕ 3 z t-3 + ϕ 4 x t-4... k 39 First Order Moving Average Model MA(1) This model assumes that the next reading can be predicted from the last residual according to a simple regression equation. z t = θ 1 e t-1 + e t e t ~ N (0, σ 2 ) z t = y t - µ zt 0 time 40
Mixed AR & MA Models: ARMA(p,q) In general, each new value depends not only on past readings but on past residuals as well. The general form is: y t = µ + ϕ 1 z t-1 + ϕ 2 z t-2 +...+ ϕ p z t-p + e t - θ 1 e t-1 - θ 2 e t-2 -...- θ q e t-q ϕ: autoregressive θ: moving average This structure is called ARMA (autoregressive moving average). The particular model is an ARMA(p,q). If the data is differentiated to become stationary, we get an ARIMA (Autoregressive, Integrated Moving Average) model. Structures also exist that describe seasonal variations and multivariate processes. 41 Summary on Time Series Models Time Series Models are used to describe the "autocorrelation structure" within each real-time signal. After the autocorrelation structure has been described, it can be removed by means of time series filtering. The models we use are known as Box-Jenkins linear models. The generation of these models involves some statistical judgmenttically. 42
Shewhart and CUSUM time series residuals 43 Fitted ARIMA(0,1,1) Example Data ARIMA(0,1,1) Residuals 44
So, what is RTSPC? RTSPC reads real-time signals from processing tools. It automatically does ACF and PACF analysis to build and save time series models. During production, RTSPC filters the real-time signals. The filtered residuals are combined using T 2 statistics. This analysis is done simultaneously in several levels: Real-Time Signals Wafer Averages Lot Averages The multivariate T 2 chart provides a robust real-time summary of machine goodness. 45 Emerging Factory Control Structure Central SPC Alarm Handler Planning and Scheduling "CIM-Bus" Maintenance Workstream Recipe Management SECSII Server 46
The SPC Server Accesses data base and draws simple X-R charts. Disables machine upon alarm Benefits from automated data collection Performs arbitrary correlations across the process Can to build causal models across the process Monitors process capabilities of essential steps Maintains 2000+ charts across a typical fab Keeps track of alarm explanations given by operators and engineers. 47 Training for SPC Operators: understand and "own" basic charts. Process Engineers: be able to decide what to monitor and what chart to use (grouping, etc.) Equipment Engineers: be able to collect tool data. Understand how to control real-time tool data. Manufacturing Manager: understand process capabilities. Monitor several charts collectively. Fab Statistician: understand the technology and its limitations. Appreciate cost of measurements. 48
Summary of SPC topics Random variables and distributions. Sampling and hypothesis testing. The assignable cause. Control chart and operating characteristic. p, c and u charts. XR, XS charts and pattern analysis. Process capability. Acceptance charts. Maximum likelihood estimation, CUSUM. Multivariate control. Evolutionary operation. Regression chart. Time series modeling. 49 The 1997 Roadmap Year 1997 1999 2001 2003 2006 2009 2012 Feature nm 250 180 150 130 100 70 50 Area mm 2 300 340 385 430 520 620 750 Density cm -2 3.7M 6.2M 10M 18M 39M 84M 180M Cost µc/tr 3000 1735 1000 580 255 110 50 technology 248 248 193? 157? 14 14 14 wafer size 200 300 300 300 300 450 450 Function/milicent 20 15 10 5 0 1997 1999 2001 2003 2006 2009 2012 Function/milicent 50
Where will the Extra Productivity Come from? (Jim Owens, Sematech) 51 The Opportunities Year 1997 1999 2001 2003 2006 2009 2012 Feature nm 250 180 150 130 100 70 50 Yield 85% 95% 100? 100? 100? 100? Equipment 35% 50% utilization Test 5-15% wafers Setup 10% Speed 15% OEE 30% Test Wafers 8% Quality 2% No Oper 10% No Prod 7% Down Pl 3% Down Un 15% 52
Basic CD Economics (D. Gerold et al, Sematech AEC/APC, Sept 97, Lake Tahoe, NV) 53 Application of Run-to-Run Control at Motorola Dose n = Dose n-1 - β (CD n-1 - CD target ) (D. Gerold et al, Sematech AEC/APC, Sept 97, Lake Tahoe, NV) 54
CD Improvement at Motorola σ Leff reduced by 60% (D. Gerold et al, Sematech AEC/APC, Sept 97, Lake Tahoe, NV) 55 Why was this Improvement Important? 600k APC investment, recovered in two days... 56
Less Tangible Opportunities Reduce cost of second sourcing (facilitate technology transfer) Dramatically increase flexibility (beat competition with more customized options) Extend life span of older technologies Cut time to market (by linking manufacturing to design) 57 What is the current extend of control in our industry? Widespread inspection and SPC Systematic setup and calibration (DOE) Widespread use of RSM / Taguchi techniques Factory statistics is an established discipline 58
Advances for the Semiconductor Industry An old idea - ensure equipment integrity - automatically A new idea - perform feedback control on the workpiece Isolate performance from technology Isolate technology from equipment Create a process with truly interchangeable parts 59 How do you steer your Process? 1980: short distance, wide road, but no steering wheel, bad front end alignment 1998: longer distance, narrow road, no steering wheel, tricky front end alignment 2005? 60
Changing the do not touch my process attitude A stable process is one that is locally characterized and locked. SPC is used to make sure it stays there. An agile process is one that is characterized over a region of operation. Process data and control algorithms are used to obtain goals. Can we reach the 2010 process goals with a stable process? 61 In Conclusion SPC provides open loop control. In-situ data can be used for tighter run-to-run and supervisory control. Several technical (and some cultural) problems must be addressed before that happens: Need sensors that are simple, non-intrusive and robust. User interfaces suitable for the production floor. Next step in the evolution of manufacturing: From hand crafted products, to hand crafted lines to lines with interchangeable parts. 62