ReaxFF Study of Surface Chemical Reactions between α-Al2O3 Substrates and H2O/H2 Gas-Phase Molecules

We developed an Al/O/H ReaxFF force field to explore chemical reactions on α-Al 2 O 3 surfaces in H 2 O/H 2 gas-phase environments. This force field generates surface energy profiles of A-, C-, R-, and M-planes with various terminations (Al- or O-) and predicts the thermodynamic and kinetic behaviors of hydrolysis on Al-terminated α-Al 2 O 3 (0001), consistent with quantum chemical studies. Molecular dynamics (MD) simulations of H 2 O/α-Al 2 O 3 (0001) reveal that water autocatalysis plays a significant role in accelerating H 2 O dissociations on Al-terminated α-Al 2 O 3 (0001). Compared with the 50% Al-terminated surface, the 100% Al-terminated surface becomes more easily hydroxylated at temperatures as low as 350 K, relying more on an O x H y clustering mechanism than complete H 2 O dissociations, and desorbs significantly more H 2 O molecules once heated up to 500 K or higher. But heating cannot eliminate surface hydroxyls for either case, and achieving a Gibbsite-like surface by H 2 O exposure is unlikely. H 2 O dissociations on α-Al 2 O 3 (0001) terminated with randomly distributed surface Al species deviate from 1–2 and 1–4 pathways due to irregular vacancy defects, and a random surface appears to be more reactive to H 2 O than the ordered one with the same surface Al coverage. Simulations of H 2 /α-Al 2 O 3 suggest that the combination of a dense surface O coverage and a low thermodynamic surface stability leads to elevated H 2 dissociation kinetics. To accelerate the surface O removals of 100% O-terminated α-Al 2 O 3 (0001) in H 2 gas exposure, we reduced the H–H σ bond energy parameter, equivalent to lowering the H 2 dissociation barrier by ∼ 19.4 kcal/mol during the simulation. After ∼ 1.5 ns, the surface termination became comparable to the 100% Al-terminated one but retained a small quantity of hydroxyls. This force field reveals how the α-Al 2 O 3 crystallographic plane and the surface termination influence the dissociation behaviors of H 2 O/H 2 gas molecules and lays the foundation for future force field developments targeted at thin film epitaxy on sapphire.