Skip to content

Reaction and relaxation at surface hotspots: using molecular dynamics and the energy-grained master equation to describe diamond etching

Research output: Contribution to journalArticle

Original languageEnglish
JournalPhilosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
Volume375
Issue number2092
DOIs
DateAccepted/In press - 11 Jan 2017
DatePublished (current) - 28 Apr 2017

Abstract

The extent to which vibrational energy transfer dynamics can impact reaction outcomes beyond the gas phase remains an active research question. Molecular dynamics (MD) simulations are the method of choice for investigating such questions; however, they can be extremely expensive, and therefore it is worth developing cheaper models that are capable of furnishing reasonable results. This paper has two primary aims. First, we investigate the competition between energy relaxation and reaction at 'hotspots' that form on the surface of diamond during the chemical vapour deposition process. To explore this, we developed an efficient reactive potential energy surface by fitting an empirical valence bond model to higher-levelab initioelectronic structure theory. We then ran 160 000 NVE trajectories on a large slab of diamond, and the results are in reasonable agreement with experiment: they suggest that energy dissipation from surface hotspots is complete within a few hundred femtoseconds, but that a small fraction of CH3does in fact undergo dissociation prior to the onset of thermal equilibrium. Second, we developed and tested a general procedure to formulate and solve the energy-grained master equation (EGME) for surface chemistry problems. The procedure we outline splits the diamond slab into system and bath components, and then evaluates microcanonical transition-state theory rate coefficients in the configuration space of the system atoms. Energy transfer from the system to the bath is estimated using linear response theory from a single long MD trajectory, and used to parametrize an energy transfer function which can be input into the EGME. Despite the number of approximations involved, the surface EGME results are in reasonable agreement with the NVE MD simulations, but considerably cheaper. The results are encouraging, because they offer a computationally tractable strategy for investigating non-equilibrium reaction dynamics at surfaces for a broader range of systems.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

    Research areas

  • Journal Article

Download statistics

No data available

Documents

Documents

  • Full-text PDF (final published version)

    Rights statement: This is the final published version of the article (version of record). It first appeared online via The Royal Society at https://doi.org/10.1098/rsta.2016.0206 . Please refer to any applicable terms of use of the publisher.

    Final published version, 1 MB, PDF-document

    Licence: CC BY

DOI

View research connections

Related faculties, schools or groups