Application of the ReaxFF reactive force field to reactive dynamics of hydrocarbon chemisorption and decomposition

Jonathan E. Mueller, Adri Van Duin, William A. Goddard

Research output: Contribution to journalArticle

99 Citations (Scopus)

Abstract

We report here reactive dynamics (RD) simulations of the adsorption and decomposition of a gas of 20-120 methane, ethyne, ethene, benzene, cyclohexane, or propene molecules interacting with a 21 Å diameter nickel nanoparticle (468 atoms). These RD simulations use the recently developed ReaxFF reactive force field to describe decomposition, reactivity, and desorption of hydrocarbons as they interact with nickel surfaces. We carried out 100 ps of RD as the temperature is ramped at a constant rate from 500 to 2500 K (temperature programmed reactions). We find that all four unsaturated hydrocarbon species chemisorb to the catalyst particle with essentially no activation energy (attaching to the surface through π electrons) and then proceed to decompose by breaking C-H bonds to form partially dehydrogenated species prior to decomposition to lower order hydrocarbons. The eventual breaking of C-C bonds usually involves a surface Ni atom inserting into the C-C bond to produce an atomic C that simultaneously with C-C cleavage moves into the subsurface layer of the particle. The greater stability of this subsurface atomic C (forming up to four Ni-C bonds) over adatom C on the particle surface (forming at most three Ni-C bonds) is critical for favorable cleaving of C-C bonds. For the two saturated hydrocarbon species (methane and cyclohexane), we observe an activation energy associated with dissociative chemisorption. These results are consistent with available experimental reactivity data and quantum mechanics (QM) energy surfaces, validating the accuracy of ReaxFF for studying hydrocarbon decomposition on nickel clusters.

Original languageEnglish (US)
Pages (from-to)5675-5685
Number of pages11
JournalJournal of Physical Chemistry C
Volume114
Issue number12
DOIs
StatePublished - Apr 1 2010

Fingerprint

Chemisorption
Hydrocarbons
chemisorption
field theory (physics)
hydrocarbons
Decomposition
decomposition
Nickel
nickel
cyclohexane
Methane
Cyclohexane
methane
reactivity
activation energy
Activation energy
Acetylene
Atoms
Adatoms
alkanes

All Science Journal Classification (ASJC) codes

  • Electronic, Optical and Magnetic Materials
  • Energy(all)
  • Physical and Theoretical Chemistry
  • Surfaces, Coatings and Films

Cite this

@article{10d2a00b3cc849c6a1b21ae544a13b91,
title = "Application of the ReaxFF reactive force field to reactive dynamics of hydrocarbon chemisorption and decomposition",
abstract = "We report here reactive dynamics (RD) simulations of the adsorption and decomposition of a gas of 20-120 methane, ethyne, ethene, benzene, cyclohexane, or propene molecules interacting with a 21 {\AA} diameter nickel nanoparticle (468 atoms). These RD simulations use the recently developed ReaxFF reactive force field to describe decomposition, reactivity, and desorption of hydrocarbons as they interact with nickel surfaces. We carried out 100 ps of RD as the temperature is ramped at a constant rate from 500 to 2500 K (temperature programmed reactions). We find that all four unsaturated hydrocarbon species chemisorb to the catalyst particle with essentially no activation energy (attaching to the surface through π electrons) and then proceed to decompose by breaking C-H bonds to form partially dehydrogenated species prior to decomposition to lower order hydrocarbons. The eventual breaking of C-C bonds usually involves a surface Ni atom inserting into the C-C bond to produce an atomic C that simultaneously with C-C cleavage moves into the subsurface layer of the particle. The greater stability of this subsurface atomic C (forming up to four Ni-C bonds) over adatom C on the particle surface (forming at most three Ni-C bonds) is critical for favorable cleaving of C-C bonds. For the two saturated hydrocarbon species (methane and cyclohexane), we observe an activation energy associated with dissociative chemisorption. These results are consistent with available experimental reactivity data and quantum mechanics (QM) energy surfaces, validating the accuracy of ReaxFF for studying hydrocarbon decomposition on nickel clusters.",
author = "Mueller, {Jonathan E.} and {Van Duin}, Adri and Goddard, {William A.}",
year = "2010",
month = "4",
day = "1",
doi = "10.1021/jp9089003",
language = "English (US)",
volume = "114",
pages = "5675--5685",
journal = "Journal of Physical Chemistry C",
issn = "1932-7447",
publisher = "American Chemical Society",
number = "12",

}

Application of the ReaxFF reactive force field to reactive dynamics of hydrocarbon chemisorption and decomposition. / Mueller, Jonathan E.; Van Duin, Adri; Goddard, William A.

In: Journal of Physical Chemistry C, Vol. 114, No. 12, 01.04.2010, p. 5675-5685.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Application of the ReaxFF reactive force field to reactive dynamics of hydrocarbon chemisorption and decomposition

AU - Mueller, Jonathan E.

AU - Van Duin, Adri

AU - Goddard, William A.

PY - 2010/4/1

Y1 - 2010/4/1

N2 - We report here reactive dynamics (RD) simulations of the adsorption and decomposition of a gas of 20-120 methane, ethyne, ethene, benzene, cyclohexane, or propene molecules interacting with a 21 Å diameter nickel nanoparticle (468 atoms). These RD simulations use the recently developed ReaxFF reactive force field to describe decomposition, reactivity, and desorption of hydrocarbons as they interact with nickel surfaces. We carried out 100 ps of RD as the temperature is ramped at a constant rate from 500 to 2500 K (temperature programmed reactions). We find that all four unsaturated hydrocarbon species chemisorb to the catalyst particle with essentially no activation energy (attaching to the surface through π electrons) and then proceed to decompose by breaking C-H bonds to form partially dehydrogenated species prior to decomposition to lower order hydrocarbons. The eventual breaking of C-C bonds usually involves a surface Ni atom inserting into the C-C bond to produce an atomic C that simultaneously with C-C cleavage moves into the subsurface layer of the particle. The greater stability of this subsurface atomic C (forming up to four Ni-C bonds) over adatom C on the particle surface (forming at most three Ni-C bonds) is critical for favorable cleaving of C-C bonds. For the two saturated hydrocarbon species (methane and cyclohexane), we observe an activation energy associated with dissociative chemisorption. These results are consistent with available experimental reactivity data and quantum mechanics (QM) energy surfaces, validating the accuracy of ReaxFF for studying hydrocarbon decomposition on nickel clusters.

AB - We report here reactive dynamics (RD) simulations of the adsorption and decomposition of a gas of 20-120 methane, ethyne, ethene, benzene, cyclohexane, or propene molecules interacting with a 21 Å diameter nickel nanoparticle (468 atoms). These RD simulations use the recently developed ReaxFF reactive force field to describe decomposition, reactivity, and desorption of hydrocarbons as they interact with nickel surfaces. We carried out 100 ps of RD as the temperature is ramped at a constant rate from 500 to 2500 K (temperature programmed reactions). We find that all four unsaturated hydrocarbon species chemisorb to the catalyst particle with essentially no activation energy (attaching to the surface through π electrons) and then proceed to decompose by breaking C-H bonds to form partially dehydrogenated species prior to decomposition to lower order hydrocarbons. The eventual breaking of C-C bonds usually involves a surface Ni atom inserting into the C-C bond to produce an atomic C that simultaneously with C-C cleavage moves into the subsurface layer of the particle. The greater stability of this subsurface atomic C (forming up to four Ni-C bonds) over adatom C on the particle surface (forming at most three Ni-C bonds) is critical for favorable cleaving of C-C bonds. For the two saturated hydrocarbon species (methane and cyclohexane), we observe an activation energy associated with dissociative chemisorption. These results are consistent with available experimental reactivity data and quantum mechanics (QM) energy surfaces, validating the accuracy of ReaxFF for studying hydrocarbon decomposition on nickel clusters.

UR - http://www.scopus.com/inward/record.url?scp=77950243995&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=77950243995&partnerID=8YFLogxK

U2 - 10.1021/jp9089003

DO - 10.1021/jp9089003

M3 - Article

AN - SCOPUS:77950243995

VL - 114

SP - 5675

EP - 5685

JO - Journal of Physical Chemistry C

JF - Journal of Physical Chemistry C

SN - 1932-7447

IS - 12

ER -