Kinetic pathways of coesite densification from metadynamics

We study compression of coesite to pressures above 35 GPa, substantially beyond the equilibrium transition pressure to octahedral phases (8 GPa to stishovite). Experiments at room temperature showed that up to 30 GPa the metastable coesite structure develops only minor displacive changes (coesite-II and coesite-III) while the Si atoms remain 4-coordinated. Beyond 30 GPa, reconstructive transformations start, following different pathways from the complex structure of coesite. In addition to amorphization, two different crystalline outcomes were observed. One is the formation of defective high-pressure octahedral phases [Hu et al., Nat. Commun. 6, 6630 (2015)], and another one is the formation of unusual and complex dense phases coesite-IV and coesite-V with Si atoms in 4-fold, 5-fold, and 6-fold coordination [Bykova et al., Nat. Commun. 9, 4789 (2018)]. Capturing these structural transformations computationally represents a challenge. Here, we show that employing metadynamics with Si-O coordination number and volume as generic collective variables in combination with a machine-learning based ACE potential [Erhard et al., Nat. Commun. 15, 1927 (2024)], one naturally observes all three mentioned pathways, resulting in the phases observed experimentally. We describe the atomistic mechanisms along the transformation pathways. While the pathway to coesite-IV is simpler, the transformation to octahedral phases involves two steps: first, a hcp sublattice of O atoms is formed where Si atoms occupy octahedral positions, but the octahedron chains do not form a regular pattern. In the second step, the Si atoms order and the chains develop a more regular arrangement. We predict that the pathway to coesite-IV is preferred at room temperature while, at 600 K, the formation of octahedral phases is more likely.