Exploratory Direct Dynamics Simulations of 3O2 Reaction with Graphene at High Temperatures

Direct chemical dynamics simulations at high temperatures of reaction between 3O2 and graphene containing varied number of defects were performed using the VENUS-MOPAC code. Graphene was modeled using (5a,6z)-periacene, a poly aromatic hydrocarbon with 5 and 6 benzene rings in the armchair and zigzag directions, respectively. Up to six defects were introduced by removing carbon atoms from the basal plane. Usage of the PM7/unrestricted Hartree–Fock (UHF) method, for the simulations, was validated by benchmarking singlet-triplet gaps of n-acenes and (5a,nz) periacenes with high-level theoretical calculations. PM7/UHF calculations showed that graphene with different number of vacancies has different ground electronic states. Dynamics simulations were performed for two 3O2 collision energies Ei of 0.4 and 0.7 eV, with the incident angle normal to the graphene plane at 1375 K. Collisions on graphene with one, two, three, and four vacancies (1C-, 2C-, 3C-, and 4C-vacant graphene) showed no reactive trajectories...