Structures, Energetics and Reactions of Hydrocarbons on Nickel

Author: Mueller, Jonathan Edward

Year: 2010

Degree: Dissertation (Ph.D.)

Advisor: Goddard, William A., III

Committee Members: Heath, James R.; Bercaw, John E.; Wang, Zhen-Gang; Goddard, William A., III

Option: Chemistry

DOI: 10.7907/RVXX-Z341

Abstract

To better understand and improve reactive processes on nickel surfaces such as the catalytic steam reforming of hydrocarbons, the decomposition of hydrocarbons at fuel cell anodes, and the growth of carbon nanotubes, we have performed atomistic studies of hydrocarbon adsorption and decomposition on low index nickel surfaces and nickel catalyst nanoparticles. Quantum mechanics (QM) calculations utilizing the PBE flavor of density functional theory (DFT) were performed on all CHx and C2Hy species to determine their structures and energies on Ni(111). In good agreement with experiments, we find that CH is the most stable form of CHx on Ni(111). It is a stable intermediate in both methane dehydrogenation and CO methanation, while CH(2,ad) is only stable during methanation. We also find that nickel surface atoms play an important catalytic role in C-H bond formation and cleavage. For the C2Hy species we find a low surface coverage decomposition pathway proceeding through CHCHad, the most stable intermediate, and a high surface coverage pathway which proceeds through CCH3,ad, the next most stable intermediate. Our enthalpies along these pathways are consistent with experimental observations.

To extend our study to larger systems and longer time scales, we have developed the ReaxFF reactive force field to describe hydrocarbon decomposition and reformation on nickel catalyst surfaces. The ReaxFF parameters were fit to geometries and energy surfaces from DFT calculations involving a large number of reaction pathways and equations of state for nickel, nickel carbides, and various hydrocarbon species chemisorbed on Ni(111), Ni(110) and Ni(100). The resulting ReaxFF description was validated against additional DFT data to demonstrate its accuracy, and used to perform reaction dynamics (RD) simulations on methyl decomposition for comparison with experiment. Finally ReaxFF RD simulations were applied to the chemisorption and decomposition of six different hydrocarbons (methane, acetylene, ethylene, benzene, cyclohexane and propylene) on a 468 atom nickel nanoparticle. These simulations realistically model hydrocarbon feedstock decomposition and provide reaction pathways relevant to this part of the carbon nanotube growth process. They show that C-C π bonds provide a low barrier pathway for chemisorption, and that the low energy of subsurface C is an important driving force in breaking C-C bonds.

Files