Simulating anisotropic etching of silicon in any etchant: evolutionary algorithm for the calibration of the continuous cellular automaton

An evolutionary algorithm is presented for the automated calibration of the continuous cellular automaton for the simulation of isotropic and anisotropic wet chemical etching of silicon in as many as 31 widely different and technologically relevant etchants, including KOH, KOH+IPA, TMAH and TMAH+Triton, in various concentrations and temperatures. Based on state-of-the-art evolutionary operators, we implement a robust algorithm for the simultaneous optimization of roughly 150 microscopic removal rates based on the minimization of a cost function with four quantitative error measures, including (i) the error between simulated and experimental macroscopic etch rates for numerous surface orientations all over the unit sphere, (ii) the error due to underetching asymmetries and floor corrugation features observed in simulated silicon samples masked using a circular pattern, (iii) the error associated with departures from a step-flow-based hierarchy in the values of the microscopic removal rates, and (iv) the error associated with deviations from a step-flow-based clustering of the microscopic removal rates. For the first time, we present the calibration and successful simulation of two technologically relevant CMOS compatible etchants, namely TMAH and, especially, TMAH+Triton, providing several comparisons between simulated and experimental MEMS structures based on multi-step etching in these etchants.

[1]  Prem Pal,et al.  Complex three-dimensional structures in Si{1 0 0} using wet bulk micromachining , 2009 .

[2]  James P. Crutchfield,et al.  Revisiting the Edge of Chaos: Evolving Cellular Automata to Perform Computations , 1993, Complex Syst..

[3]  Lawrence. Davis,et al.  Handbook Of Genetic Algorithms , 1990 .

[4]  Bastien Chopard,et al.  Cellular automata modeling , 1998 .

[5]  Kazuo Sato,et al.  Silicon Micromachining Based on Surfactant-Added Tetramethyl Ammonium Hydroxide: Etching Mechanism and Advanced Applications , 2010 .

[6]  Zbigniew Michalewicz,et al.  Evolutionary Computation 2 : Advanced Algorithms and Operators , 2000 .

[7]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[8]  Kazuo Asaumi,et al.  Anisotropic-etching process simulation system MICROCAD analyzing complete 3D etching profiles of single crystal silicon , 1997, Proceedings IEEE The Tenth Annual International Workshop on Micro Electro Mechanical Systems. An Investigation of Micro Structures, Sensors, Actuators, Machines and Robots.

[9]  Zbigniew Michalewicz,et al.  Genetic algorithms + data structures = evolution programs (2nd, extended ed.) , 1994 .

[10]  Miguel A. Gosálvez,et al.  Reconstructing the 3D etch rate distribution of silicon in anisotropic etchants using data from vicinal {1 0 0}, {1 1 0} and {1 1 1} surfaces , 2011 .

[11]  Adam S. Foster,et al.  An atomistic introduction to anisotropic etching , 2007 .

[12]  Miguel A. Gosálvez,et al.  Discrete and continuous cellular automata for the simulation of propagating surfaces , 2009 .

[13]  Yan Xing,et al.  Atomistic methods for the simulation of evolving surfaces , 2008 .

[14]  Branislav Radjenović,et al.  Nonconvex Hamiltonians in three dimensional level set simulations of the wet etching of silicon , 2006 .

[15]  J. Fruhauf,et al.  A simulation tool for orientation dependent etching , 1993 .

[16]  David B. Fogel,et al.  An introduction to simulated evolutionary optimization , 1994, IEEE Trans. Neural Networks.

[17]  Yasuroh Iriye,et al.  Characterization of orientation-dependent etching properties of single-crystal silicon: effects of KOH concentration , 1998 .

[18]  K. Sato,et al.  Step flow-based cellular automaton for the simulation of anisotropic etching of complex MEMS structures , 2007 .

[19]  Qing‐An Huang,et al.  A cellular automaton-based simulator for silicon anisotropic etching processes considering high index planes , 2007 .

[20]  Miguel A. Gosálvez,et al.  Fabrication of novel microstructures based on orientation-dependent adsorption of surfactant molecules in a TMAH solution , 2010 .

[21]  Qing‐An Huang,et al.  Modeling and Simulations of Anisotropic Etching of Silicon in Alkaline Solutions with Experimental Verification , 2009 .

[22]  Heinz Mühlenbein,et al.  Predictive Models for the Breeder Genetic Algorithm I. Continuous Parameter Optimization , 1993, Evolutionary Computation.

[23]  Zhenjun Zhu,et al.  Micromachining process simulation using a continuous cellular automata method , 2000, Journal of Microelectromechanical Systems.

[24]  M. Gosálvez,et al.  Chapter Twenty Four – Wet Etching of Silicon , 2010 .

[25]  K. Sato,et al.  Analytical Solution of the Continuous Cellular Automaton for Anisotropic Etching , 2008, Journal of Microelectromechanical Systems.

[26]  Marija Radmilovic-Radjenovic,et al.  Level Set Approach to Anisotropic Wet Etching of Silicon , 2010, Sensors.

[27]  R. Nieminen,et al.  Atomistic simulations of surface coverage effects in anisotropic wet chemical etching of crystalline silicon , 2002 .

[28]  Kazuo Sato,et al.  Octree-based, GPU implementation of a continuous cellular automaton for the simulation of complex, evolving surfaces , 2011, Comput. Phys. Commun..