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