We predict mechanochemical formation of heterogeneous diamond structures from rapid uniaxial compression in graphite using quantum molecular dynamics simulations. Ensembles of simulations reveal the formation of different diamondlike products starting from thermal graphite crystal configurations. We identify distinct classes of final products with characteristic probabilities of formation, stress states, and electrical properties and show through simulations of rapid quenching that these products are nominally stable and can be recovered at room temperature and pressure. Some of the diamond products exhibit significant disorder and partial closure of the energy gap between the highest-occupied and lowest-unoccupied molecular orbitals (i.e., the HOMO-LUMO gap). Seeding atomic vacancies in graphite significantly biases toward forming products with small HOMO-LUMO gap. We show that a strong correlation between the HOMO-LUMO gap and disorder in tetrahedral bonding configurations informs which kinds of structural defects are associated with gap closure. The rapid diffusionless transformation of graphite is found to lock vacancy defects into the final diamond structure, resulting in configurations that prevent $s{p}^{3}$ bonding and lead to localized HOMO and LUMO states with a small gap.
Read full abstract