Abstract:This study focuses on the simulation of wet etching of sapphire at different concentrations and temperatures. We use the Evolutionary Kinetic Monte Carlo (E-KMC) method, which combines (i) a removal probability function (RPF), to reduce the large collection of adjustable removal rates to a smaller number of energy parameters, and (ii) an evolutionary algorithm, to optimize those parameters by minimizing errors between simulated and experimental rates for a selected group of surface orientations. To improve the resemblance between simulated and experimental etch rate distributions at different concentrations, additional surface orientations located in the region with worst similarity are used. This is validated by a comparison of experimental and simulated trench profiles obtained on C-plane sapphire wafers with rectangular mask openings. Most importantly, this study addresses for the first time the temperature dependence considered in traditional formulations of the RPF, concluding that it is only partially correct. To overcome this deficiency, a physically meaningful modified removal probability function (M-RPF) is formulated as the product of a standard prefactor and a standard Boltzmann factor, in agreement with Transition State Theory. Effectively, the prefactor transforms the facet-specific orientation-dependence of the etch rate at a given temperature into the atom-specific site-dependence of removal rates on the local neighborhood, while the Boltzmann factor assigns a different activation energy to every surface atom type. It is shown that the M-RPF provides a suitable approach to correctly model the overall temperature dependence. The study extends the range of wet etching applications of the E-KMC method. [2021-0221]

Anisotropic-etching process simulation system MICROCAD analyzing complete 3D etching profiles of single crystal silicon

Development of an orientation-dependent anisotropic etching simulation system MICROCAD

Simulation of three-dimensional etch profile of silicon during orientation dependent anisotropic etching

Simulation of two-dimensional etch profile of silicon during orientation-dependent anisotropic etching

Isotropic Etch for SiO2 Microcantilever Release with ICP System

Development of quartz etching database and 3-D micromachining simulation system

Simulation of Anisotropic Chemical Etching of Single Crystalline Silicon using Cellular-Automata

Micromachining process simulation using a continuous cellular automata method

A cellular automaton-based simulator for silicon anisotropic etching processes considering high index planes

A 3-D Simulator for Silicon Anisotropic Wet Chemical Etching Process Based on Cellular Automata

Evolutionary Kinetic Monte Carlo method for the simulation of anisotropic etching of Z-cut, at-cut and BT-cut quartz

Simulation of Anisotropic Etching of Silicon using a Cellular Automata Model

Simulation of anisotropic chemical etching of crystalline silicon using a cellular automata model

Level Set Approach to Anisotropic Wet Etching of Silicon

Anisotropic Crystalline Etching Simulation Using a Continuous Cellular Automata Algorithm

Orientation-dependent surface morphology of crystalline silicon during anisotropic etching using a continuous cellular automaton

Chemical Etching Processes at the Dynamic GaAs/Electrolyte Interface in the Electrochemical Direct-Writing Micromachining

Application of the evolutionary kinetic Monte Carlo method for the simulation of anisotropic wet etching of sapphire

Evolutionary Kinetic Monte Carlo Simulation of Anisotropic Wet Etching of Sapphire at Different Concentrations and Temperatures

Anisotropic Etching of Three-Dimensional Shapes in Silicon - the Important Role of Galvanic Interaction

Three-Dimensional Simulation of DRIE Process Based on the Narrow Band Level Set and Monte Carlo Method