Abstract. Anisotropic wet chemical etching of quartz is a bulk micromachining process for the fabrication of Micro-Electro-Mechanical Systems (MEMS), such as resonators and temperature sensors. Despite the success of the Continuous Cellular Automaton (CCA) for the simulation of wet etching of silicon, the simulation of the same process for quartz has received little attention -especially from an atomistic perspective-resulting in a lack of accurate modeling tools. This paper analyzes the crystallographic structure of the main surface orientations of quartz and proposes a novel classification of the surface atoms as well as an Evolutionary Algorithm (EA) to determine suitable values for the corresponding atomistic removal rates. Not only the presented Evolutionary Continuous Cellular Automaton (ECCA) reproduces the correct macroscopic etch rate distribution for quartz hemispheres but it is also capable of performing fast and accurate 3D simulations of MEMS structures. This is shown by several comparisons between simulated and experimental results and, in particular, by a detailed, quantitative comparison for an extensive collection of trench profiles.