Yield Modeling

Keywords: yield modeling,yield,defect density,poisson yield,negative binomial,murphy model,critical area,semiconductor yield,die yield,wafer yield

Yield Modeling: Mathematical Foundations Yield modeling in semiconductor manufacturing is the mathematical framework for predicting the fraction of functional dies on a wafer. Since fabrication involves hundreds of process steps where defects can occur, accurate yield prediction is critical for: - Cost estimation and financial planning - Process optimization and control - Manufacturing capacity decisions - Design-for-manufacturability feedback Fundamental Definitions Yield (Y) is defined as: Y = fractextNumber of good diestextTotal dies on wafer The mathematical challenge involves relating yield to: - Defect density (D) - Die area (A) - Defect clustering behavior (alpha) - Process variations (sigma) The Poisson Model (Baseline) The simplest model assumes defects are randomly and uniformly distributed across the wafer. Basic Equation Y = e^-AD Where: - A = die area (cm²) - D = average defect density (defects/cm²) Mathematical Derivation If defects follow a Poisson distribution with mean lambda = AD, the probability of zero defects (functional die) is: P(X = 0) = frace^-lambda lambda^00! = e^-AD Limitations - Problem: This model consistently underestimates real yields - Reason: Actual defects cluster—they don't distribute uniformly - Result: Some wafer regions have high defect density while others are nearly defect-free Defect Clustering Models Real defects cluster due to: - Particle contamination patterns - Equipment-related issues - Process variations across the wafer - Lithography and etch non-uniformities Murphy's Model (1964) Assumes defect density is uniformly distributed between 0 and 2D_0: Y = frac1 - e^-2AD_02AD_0 For large AD_0, this approximates to: Y approx frac12AD_0 Seeds' Model Assumes exponential distribution of defect density: Y = e^-sqrtAD Negative Binomial Model (Industry Standard) This is the most widely used model in semiconductor manufacturing. Main Equation Y = left(1 + fracADalpharight)^-alpha Where alpha is the clustering parameter: - alpha to infty: Reduces to Poisson (no clustering) - alpha to 0: Extreme clustering (highly non-uniform) - Typical values: alpha approx 0.5 to 5 Mathematical Origin The negative binomial arises from a compound Poisson process: 1. Let X sim textPoisson(lambda) be the defect count 2. Let lambda sim textGamma(alpha, beta) be the varying rate 3. Marginalizing over lambda gives X sim textNegative Binomial The probability mass function is: P(X = k) = binomk + alpha - 1k left(fracbetabeta + 1right)^alpha left(frac1beta + 1right)^k The yield (probability of zero defects) becomes: Y = P(X = 0) = left(fracbetabeta + 1right)^alpha = left(1 + fracADalpharight)^-alpha Model Comparison At AD = 1: | Model | Yield | |:------|------:| | Poisson | 36.8% | | Murphy | 43.2% | | Negative Binomial (alpha = 2) | 57.7% | | Negative Binomial (alpha = 1) | 50.0% | | Seeds | 36.8% | Critical Area Analysis Not all die area is equally sensitive to defects. Critical area (A_c) is the region where a defect of given size causes failure. Definition For a defect of radius r: - Short critical area: Region where defect center causes a short circuit - Open critical area: Region where defect causes an open circuit Stapper's Critical Area Model For parallel lines of width w, spacing s, and length l: A_c(r) = begincases 0 & textif r < fracs2 [8pt] 2lleft(r - fracs2right) & textif fracs2 leq r < fracw+s2 [8pt] lw & textif r geq fracw+s2 endcases Integration Over Defect Size Distribution The total critical area integrates over the defect size distribution f(r): A_c = int_0^infty A_c(r) cdot f(r) , dr Common distributions for f(r): - Log-normal: f(r) = frac1rsigmasqrt2pi expleft(-frac(ln r - mu)^22sigma^2right) - Power-law: f(r) propto r^-p for r_min leq r leq r_max Yield with Critical Area Y = expleft(-int_0^infty A_c(r) cdot D(r) , drright) Yield Decomposition Total yield is typically factored into independent components: Y_texttotal = Y_textgross times Y_textrandom times Y_textparametric Component Definitions | Component | Description | Typical Range | |:----------|:------------|:-------------:| | Y_textgross | Catastrophic defects, edge loss, handling damage | 95–99% | | Y_textrandom | Random particle defects (main focus of yield modeling) | 70–95% | | Y_textparametric | Process variation causing spec failures | 90–99% | Extended Decomposition For more detailed analysis: Y_texttotal = Y_textgross times prod_i=1^N_textlayers Y_textrandom,i times prod_j=1^M_textparams Y_textparam,j Parametric Yield Modeling Dies may function but fail to meet performance specifications due to process variation. Single Parameter Model If parameter X sim mathcalN(mu, sigma^2) with specification limits [L, U]: Y_p = Phileft(fracU - musigmaright) - Phileft(fracL - musigmaright) Where Phi(cdot) is the standard normal cumulative distribution function: Phi(z) = frac1sqrt2pi int_-infty^z e^-t^2/2 , dt Process Capability Indices Cp (Process Capability) C_p = fracUSL - LSL6sigma Cpk (Process Capability Index) C_pk = minleft(fracUSL - mu3sigma, fracmu - LSL3sigmaright) Cpk to Yield Conversion | C_pk | Sigma Level | Yield | DPMO | |:--------:|:-----------:|:-----:|-----:| | 0.33 | 1σ | 68.27% | 317,300 | | 0.67 | 2σ | 95.45% | 45,500 | | 1.00 | 3σ | 99.73% | 2,700 | | 1.33 | 4σ | 99.9937% | 63 | | 1.67 | 5σ | 99.999943% | 0.57 | | 2.00 | 6σ | 99.9999998% | 0.002 | Multiple Correlated Parameters For n parameters with mean vector boldsymbolmu and covariance matrix boldsymbolSigma: Y_p = int int cdot int_mathcalR frac1(2pi)^n/2|boldsymbolSigma|^1/2 expleft(-frac12(mathbfx-boldsymbolmu)^T boldsymbolSigma^-1(mathbfx-boldsymbolmu)right) dmathbfx Where mathcalR is the specification region. Computational Methods: - Monte Carlo integration - Gaussian quadrature - Importance sampling Spatial Yield Models Modern fabs analyze spatial patterns using wafer maps to identify systematic issues. Radial Defect Density Model Accounts for edge effects: D(r) = D_0 + D_1 r^2 Where: - r = distance from wafer center - D_0 = baseline defect density - D_1 = radial coefficient General Spatial Model D(x, y) = D_0 + sum_i beta_i phi_i(x, y) Where phi_i(x, y) are spatial basis functions (e.g., Zernike polynomials). Spatial Autocorrelation (Moran's I) I = fracn sum_i sum_j w_ij(Z_i - barZ)(Z_j - barZ)W sum_i (Z_i - barZ)^2 Where: - Z_i = pass/fail indicator for die i (1 = fail, 0 = pass) - w_ij = spatial weight between dies i and j - W = sum_i sum_j w_ij - barZ = mean failure rate Interpretation: - I > 0: Clustered failures (systematic issue) - I approx 0: Random failures - I < 0: Dispersed failures (rare) Variogram Analysis The semi-variogram gamma(h) measures spatial dependence: gamma(h) = frac12|N(h)| sum_(i,j) in N(h) (Z_i - Z_j)^2 Where N(h) is the set of die pairs separated by distance h. Multi-Layer Yield Modern ICs have many process layers, each contributing to yield loss. Independent Layers Y_texttotal = prod_i=1^N Y_i = prod_i=1^N left(1 + fracA_i D_ialpha_iright)^-alpha_i Simplified Model If defects are independent across layers with similar clustering: Y = left(1 + fracA cdot D_texttotalalpharight)^-alpha Where: D_texttotal = sum_i=1^N D_i Layer-Specific Critical Areas Y = prod_i=1^N expleft(-A_c,i cdot D_iright) For Poisson model, or: Y = prod_i=1^N left(1 + fracA_c,i D_ialpha_iright)^-alpha_i For negative binomial. Yield Learning Curves Yield improves over time as processes mature and defect sources are eliminated. Exponential Learning Model D(t) = D_infty + (D_0 - D_infty)e^-t/tau Where: - D_0 = initial defect density - D_infty = asymptotic (mature) defect density - tau = learning time constant Power Law (Wright's Learning Curve) D(n) = D_1 cdot n^-b Where: - n = cumulative production volume (wafers or lots) - D_1 = defect density after first unit - b = learning rate exponent (typically 0.2 leq b leq 0.4) Yield vs. Time Combining with yield model: Y(t) = left(1 + fracA cdot D(t)alpharight)^-alpha Yield-Redundancy Models (Memory) Memory arrays use redundant rows/columns for defect tolerance through laser repair or electrical fusing. Poisson Model with Redundancy If a memory has R spare elements and defects follow Poisson: Y_textrepaired = sum_k=0^R frac(AD)^k e^-ADk! This is the CDF of the Poisson distribution: Y_textrepaired = fracGamma(R+1, AD)Gamma(R+1) = fracgamma(R+1, AD)R! Where gamma(cdot, cdot) is the lower incomplete gamma function. Negative Binomial Model with Redundancy Y_textrepaired = sum_k=0^R binomk+alpha-1k left(fracalphaalpha + ADright)^alpha left(fracADalpha + ADright)^k Repair Coverage Factor Y_textrepaired = Y_textbase + (1 - Y_textbase) cdot RC Where RC is the repair coverage (fraction of defective dies that can be repaired). Statistical Estimation Maximum Likelihood Estimation for Negative Binomial Given wafer data with n_i dies and k_i failures per wafer i: Likelihood function: mathcalL(D, alpha) = prod_i=1^W binomn_ik_i (1-Y)^k_i Y^n_i - k_i Log-likelihood: ell(D, alpha) = sum_i=1^W left[ lnbinomn_ik_i + k_i ln(1-Y) + (n_i - k_i) ln Y right] Estimation: Requires iterative numerical methods: - Newton-Raphson - EM algorithm - Gradient descent Bayesian Estimation With prior distributions P(D) and P(alpha): P(D, alpha mid textdata) propto P(textdata mid D, alpha) cdot P(D) cdot P(alpha) Common priors: - D sim textGamma(a_D, b_D) - alpha sim textGamma(a_alpha, b_alpha) Model Selection Use information criteria to compare models: Akaike Information Criterion (AIC): AIC = -2ln(mathcalL) + 2k Bayesian Information Criterion (BIC): BIC = -2ln(mathcalL) + kln(n) Where k = number of parameters, n = sample size. Economic Model Die Cost textCost_textdie = fractextCost_textwaferN_textdies times Y Dies Per Wafer Accounting for edge exclusion (dies must fit entirely within usable area): N approx fracpi D_w^24A - fracpi D_wsqrt2A Where: - D_w = wafer diameter - A = die area More accurate formula: N = fracpi (D_w/2 - E)^2A cdot eta Where: - E = edge exclusion distance - eta = packing efficiency factor (approx 0.9) Cost Sensitivity Analysis Marginal cost impact of yield change: fracpartial textCost_textdiepartial Y = -fractextCost_textwaferN cdot Y^2 Break-Even Analysis Minimum yield for profitability: Y_textmin = fractextCost_textwaferN cdot textPrice_textdie Key Models Yield Models Comparison | Model | Formula | Best Application | |:------|:--------|:-----------------| | Poisson | Y = e^-AD | Lower bound estimate, theoretical baseline | | Murphy | Y = frac1-e^-2AD2AD | Moderate clustering | | Seeds | Y = e^-sqrtAD | Exponential clustering | | Negative Binomial | Y = left(1 + fracADalpharight)^-alpha | Industry standard, tunable clustering | | Critical Area | Y = e^-int A_c(r)D(r)dr | Layout-aware prediction | Parameters | Parameter | Symbol | Typical Range | Description | |:----------|:------:|:-------------:|:------------| | Defect Density | D | 0.01–1 /cm² | Defects per unit area | | Die Area | A | 10–800 mm² | Size of single chip | | Clustering Parameter | alpha | 0.5–5 | Degree of defect clustering | | Learning Rate | b | 0.2–0.4 | Yield improvement rate | Equations Basic yield: Y = e^-AD Industry standard: Y = left(1 + fracADalpharight)^-alpha Total yield: Y_texttotal = Y_textgross times Y_textrandom times Y_textparametric Die cost: textCost_textdie = fractextCost_textwaferN times Y Practical Implementation Workflow 1. Data Collection - Gather wafer test data (pass/fail maps) - Record lot/wafer identifiers and timestamps 2. Parameter Estimation - Estimate D and alpha via MLE or Bayesian methods - Validate with holdout data 3. Spatial Analysis - Generate wafer maps - Calculate Moran's I to detect clustering - Identify systematic defect patterns 4. Parametric Analysis - Model electrical parameter distributions - Calculate C_pk for key parameters - Estimate parametric yield losses 5. Model Integration - Combine: Y_texttotal = Y_textgross times Y_textrandom times Y_textparametric - Validate against actual production data 6. Trend Monitoring - Track D and alpha over time - Fit learning curve models - Project future yields 7. Cost Optimization - Calculate die cost at current yield - Identify highest-impact improvement opportunities - Optimize die size vs. yield trade-off

Source: ChipFoundryServices — Search this topic — Ask CFSGPT