Three-Dimensional Current Density Distribution Simulations for a Resistive Patterned Wafer

A combined boundary element method (BEM)-finite element method (FEM) numerical approach is used for the simulation of current density and layer thickness distributions in a wafer plating reactor. The current and potential distribution effects due to the electrolyte resistivity are modeled with BEM, while the transient internal resistive 'terminal' effect of the wafer is modeled using FEM. A nonlinear Butler-Volmer type overpotential relation is considered to describe cathode kinetics. The declining internal wafer resistivity that is due to the growth of the initial copper seed layer, is modeled over a number of discrete time steps. Different contacting methods (4 or 8 contact points, ring contact) are investigated, and their terminal effect time is compared. The wafer consists of a ring-shaped current thief, with adjacent photoresist area's and a central electroactive patterned zone.

[1]  R. Wijsman,et al.  Theory of the Effect of Electrode Resistance on Current Density Distribution in Electrolytic Cells , 1953 .

[2]  S. Boyer,et al.  Study of sacrificial anode cathodic protection of buried tanks: Numerical modelling , 1999 .

[3]  M. Matlosz,et al.  Secondary Current Distribution in a Hull Cell Boundary Element and Finite Element Simulation and Experimental Verification , 1987 .

[4]  J. Thonstad,et al.  Mathematical modelling of aluminium cells with prebaked anodes: Part II: Current distribution and influence of sideledge , 1997 .

[5]  J. Deconinck,et al.  Optimisation of a cupplater reactor for gold deposition on wafers , 2001 .

[6]  R. M. Degerstedt,et al.  Mathematical Models for Cathodic Protection of an Underground Pipeline with Coating Holidays: Part 1 — Theoretical Development , 1995 .

[7]  H. Y. Cheh,et al.  Optimization of Electrodeposit Uniformity by the Use of Auxiliary Electrodes , 1990 .

[8]  Mark E. Orazem,et al.  Application of Boundary Element Models to Predict Effectiveness of Coupons for Accessing Cathodic Protection of Buried Structures , 2000 .

[9]  Richard C. Alkire,et al.  Anisotropic chemical pattern etching of copper foil. III: Mathematical model , 1994 .

[10]  A. West,et al.  Nonuniform Current Distribution and Thickness during Electrodeposition onto Resistive Substrates , 1992 .

[11]  Herman Deconinck,et al.  Object‐oriented three‐dimensional hybrid grid generation , 2003 .

[12]  Gregory J. Wilson,et al.  Two- and three-dimensional numerical modeling of copper electroplating for advanced ULSI metallization , 2000 .

[13]  Jack J. Dongarra,et al.  Software Libraries for Linear Algebra Computations on High Performance Computers , 1995, SIAM Rev..

[14]  R. White,et al.  Predicted Secondary Current Distributions for Linear Kinetics in a Modified Three‐Dimensional Hull Cell , 1990 .

[15]  C. Hussey,et al.  Aluminum Bromide‐1‐Methyl‐3‐Ethylimidazolium Bromide Ionic Liquids II . Transport Numbers , 1986 .

[16]  G. Maggetto,et al.  Calculation of Current Distribution and Electrode Shape Change by the Boundary Element Method , 1985 .

[17]  Young S. Choi,et al.  Three-dimensional calculation of current distribution in electrodeposition on patterned cathode with auxiliary electrode , 1996 .

[18]  Ralph E. White,et al.  A Semianalytical Method for Predicting Primary and Secondary Current Density Distributions: Linear and Nonlinear Boundary Conditions , 2000 .

[19]  L. Bortels,et al.  The multi-dimensional upwinding method as a new simulation tool for the analysis of multi-ion electrolytes controlled by diffusion, convection and migration. Part 1. Steady state analysis of a parallel plane flow channel , 1996 .

[20]  R. M. Degerstedt,et al.  Mathematical Models for Cathodic Protection of an Underground Pipeline with Coating Holidays: Part 2 Case Studies of Parallel Anode Cathodic Protection Systems underground pipelines by parallel-ribbon anodes. , 1997 .

[21]  Ralph E. White,et al.  Simulating Shape Changes during Electrodeposition Primary and Secondary Current Distribution , 2002 .

[22]  L. Bortels,et al.  Numerical 3-D Simulation of a Cathodic Protection System for a Buried Pipe Segment Surrounded by a Load Relieving U-Shaped Vault , 2003 .

[23]  C. Hussey,et al.  Aluminum Bromide‐1‐Methyl‐3‐Ethylimidazolium Bromide Ionic Liquids I . Densities, Viscosities, Electrical Conductivities, and Phase Transitions , 1986 .

[24]  R. Coifman,et al.  Fast wavelet transforms and numerical algorithms I , 1991 .