Abstract:This study explores for the first time the possibility of using an evolutionary algorithm (EA) for the determination of atomistic rates for use in kinetic Monte Carlo (KMC) simulations of anisotropic etching of silicon for the manufacture of microelectromechanical systems (MEMS). Traditionally, KMC rates are determined based on (i) computationally expensive density functional theory (DFT) calculations or, when possible, (ii) a combination of physical insight and a labor-intensive, manual procedure where, e.g., experimental and simulated surface morphologies are visually matched for a collection of surface orientations. Compared to these approaches, the evolutionary KMC method proposed in this study provides a more flexible, autonomous procedure to describe correctly a wide variety of etching conditions. We focus on the use of a functional representation of the atomistic rates, referred to as the removal probability function (RPF). This simplifies the EA search space to just a few parameters and reduces the number of required experimental data points to just a few etch rates. By proposing two alternative RPFs with four and six parameters, respectively, we show that the ability to explain the orientation dependence of the etch rate for a wide variety of etchants increases with the number of parameters and conclude that the six-parameter RPF provides sufficiently good simulations for a wide range of etching conditions, including KOH, KOH+IPA, TMAH and TMAH+Triton at different concentrations and temperatures. By uncovering the relationships between the parameters and the concentration of the etchant, it is possible to extend the simulations to nonmeasured etching conditions. Although the use of an RPF effectively restricts the search to a subspace of the atomistic rates, the present results suggest that physically meaningful KMC rates can probably be determined in the near future by direct comparison of macroscopic experiments and simulations through the use of evolutionary strategies.

Monte Carlo study of kinetic smoothing during dissolution and etching of the Kossel (100) and silicon (111) surfaces

The Dissociated Amorphous Silica Surface: Model Development and Evaluation.

Monte Carlo Simulation for the Surface Morphology of Anisotropic Etching of Crystalline Silicon

Surface Morphology Of Ge(111) During Etching By Kev Ions

A new theory for the anisotropic etching of silicon and some underdeveloped chemical micromachining concepts

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

Surface Structure Dependence of Mechanochemical Etching: Scanning Probe-Based Nanolithography Study on Si(100), Si(110), and Si(111)

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

Surface termination and hydrogen bubble adhesion on Si(100) surfaces during anisotropic dissolution in aqueous KOH

Relaxation Of Nanometer-Scale Surface Morphology

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

Study of anisotropic etching of (1 0 0) Si with ultrasonic agitation

Coarsening and Slope Selection During Crystal Growth and Etching of Ge(001)

A kinetic Monte Carlo study of the C 3 S dissolution mechanism

A Metropolis Monte Carlo Simulation Approach for Anisotropic Wet Etching and Its Applications

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

Mechanism of Molten KOH Etching of SiC Single Crystals: Comparative Study with Thermal Oxidation

Numerical Simulations of Dynamic Scaling Behavior of the Etching Model on Fractal Substrates

Evolutionary determination of kinetic Monte Carlo rates for the simulation of evolving surfaces in anisotropic etching of silicon

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