Potential energy surface and quasiclassical trajectory studies of the CN+H2 reaction

Marc A. Ter Horst, George C Schatz, Lawrence B. Harding

Research output: Contribution to journalArticle

103 Citations (Scopus)

Abstract

We present dynamical studies of the CN+H2 reaction based on an empirical potential energy surface that is derived from high quality ab initio calculations. The ab initio calculations, which use a multireference configuration interaction method with large correlation consistent basis sets, indicate that the linear HHCN barrier is about 4.3 kcal/mol above CN+H2, and that there is no reaction path which connects CN+H2 to the stable intermediate H2CN, although there is a path for dissociation of H2CN to H+HCN. The empirical surface is written as a sum of two-, three-, and four-body terms, with the two- and three-body terms for HCN based on an accurate global surface that describes both the HCN and HNC force fields. The four-body terms are developed so as to describe the HHCN linear saddle point and the H2CN minimum accurately, as well as dissociation of H2CN into HCN+H, and the ridge which separates the abstraction and H2CN dissociation pathways. Other features of the potential surface, such as the HCNH cis and trans minima, and the pathways leading to the formation of HNC+H are also described, though less accurately. Three different choices for the HHCN saddle point properties are considered. We find that the surface which matches the ab initio barrier energy most accurately gives rate constants that are too low. Much better agreement is obtained using a 3.2 kcal/mol barrier. The trajectory results show typical dependence of the CN+H2 reactive cross sections on initial translational energy and initial vibration/rotation state, with CN behaving as a spectator and H2 playing an active role in the reaction dynamics. Analysis of the H+HCN products indicates that both the C-H stretch and bend modes are significantly excited, with bend excitation showing strong sensitivity to the saddle point properties and to reagent translational energy. At translational energies below 20 kcal/mol, direct H abstraction is strongly favored over addition elimination.

Original languageEnglish
Pages (from-to)558-571
Number of pages14
JournalJournal of Chemical Physics
Volume105
Issue number2
Publication statusPublished - 1996

Fingerprint

Potential energy surfaces
potential energy
Trajectories
trajectories
saddle points
dissociation
Energy barriers
Rate constants
energy
field theory (physics)
configuration interaction
reagents
ridges
elimination
vibration
sensitivity
cross sections
products
excitation

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Potential energy surface and quasiclassical trajectory studies of the CN+H2 reaction. / Ter Horst, Marc A.; Schatz, George C; Harding, Lawrence B.

In: Journal of Chemical Physics, Vol. 105, No. 2, 1996, p. 558-571.

Research output: Contribution to journalArticle

Ter Horst, Marc A. ; Schatz, George C ; Harding, Lawrence B. / Potential energy surface and quasiclassical trajectory studies of the CN+H2 reaction. In: Journal of Chemical Physics. 1996 ; Vol. 105, No. 2. pp. 558-571.
@article{d898f1b96b95414cb786c6fa106bffa5,
title = "Potential energy surface and quasiclassical trajectory studies of the CN+H2 reaction",
abstract = "We present dynamical studies of the CN+H2 reaction based on an empirical potential energy surface that is derived from high quality ab initio calculations. The ab initio calculations, which use a multireference configuration interaction method with large correlation consistent basis sets, indicate that the linear HHCN barrier is about 4.3 kcal/mol above CN+H2, and that there is no reaction path which connects CN+H2 to the stable intermediate H2CN, although there is a path for dissociation of H2CN to H+HCN. The empirical surface is written as a sum of two-, three-, and four-body terms, with the two- and three-body terms for HCN based on an accurate global surface that describes both the HCN and HNC force fields. The four-body terms are developed so as to describe the HHCN linear saddle point and the H2CN minimum accurately, as well as dissociation of H2CN into HCN+H, and the ridge which separates the abstraction and H2CN dissociation pathways. Other features of the potential surface, such as the HCNH cis and trans minima, and the pathways leading to the formation of HNC+H are also described, though less accurately. Three different choices for the HHCN saddle point properties are considered. We find that the surface which matches the ab initio barrier energy most accurately gives rate constants that are too low. Much better agreement is obtained using a 3.2 kcal/mol barrier. The trajectory results show typical dependence of the CN+H2 reactive cross sections on initial translational energy and initial vibration/rotation state, with CN behaving as a spectator and H2 playing an active role in the reaction dynamics. Analysis of the H+HCN products indicates that both the C-H stretch and bend modes are significantly excited, with bend excitation showing strong sensitivity to the saddle point properties and to reagent translational energy. At translational energies below 20 kcal/mol, direct H abstraction is strongly favored over addition elimination.",
author = "{Ter Horst}, {Marc A.} and Schatz, {George C} and Harding, {Lawrence B.}",
year = "1996",
language = "English",
volume = "105",
pages = "558--571",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "2",

}

TY - JOUR

T1 - Potential energy surface and quasiclassical trajectory studies of the CN+H2 reaction

AU - Ter Horst, Marc A.

AU - Schatz, George C

AU - Harding, Lawrence B.

PY - 1996

Y1 - 1996

N2 - We present dynamical studies of the CN+H2 reaction based on an empirical potential energy surface that is derived from high quality ab initio calculations. The ab initio calculations, which use a multireference configuration interaction method with large correlation consistent basis sets, indicate that the linear HHCN barrier is about 4.3 kcal/mol above CN+H2, and that there is no reaction path which connects CN+H2 to the stable intermediate H2CN, although there is a path for dissociation of H2CN to H+HCN. The empirical surface is written as a sum of two-, three-, and four-body terms, with the two- and three-body terms for HCN based on an accurate global surface that describes both the HCN and HNC force fields. The four-body terms are developed so as to describe the HHCN linear saddle point and the H2CN minimum accurately, as well as dissociation of H2CN into HCN+H, and the ridge which separates the abstraction and H2CN dissociation pathways. Other features of the potential surface, such as the HCNH cis and trans minima, and the pathways leading to the formation of HNC+H are also described, though less accurately. Three different choices for the HHCN saddle point properties are considered. We find that the surface which matches the ab initio barrier energy most accurately gives rate constants that are too low. Much better agreement is obtained using a 3.2 kcal/mol barrier. The trajectory results show typical dependence of the CN+H2 reactive cross sections on initial translational energy and initial vibration/rotation state, with CN behaving as a spectator and H2 playing an active role in the reaction dynamics. Analysis of the H+HCN products indicates that both the C-H stretch and bend modes are significantly excited, with bend excitation showing strong sensitivity to the saddle point properties and to reagent translational energy. At translational energies below 20 kcal/mol, direct H abstraction is strongly favored over addition elimination.

AB - We present dynamical studies of the CN+H2 reaction based on an empirical potential energy surface that is derived from high quality ab initio calculations. The ab initio calculations, which use a multireference configuration interaction method with large correlation consistent basis sets, indicate that the linear HHCN barrier is about 4.3 kcal/mol above CN+H2, and that there is no reaction path which connects CN+H2 to the stable intermediate H2CN, although there is a path for dissociation of H2CN to H+HCN. The empirical surface is written as a sum of two-, three-, and four-body terms, with the two- and three-body terms for HCN based on an accurate global surface that describes both the HCN and HNC force fields. The four-body terms are developed so as to describe the HHCN linear saddle point and the H2CN minimum accurately, as well as dissociation of H2CN into HCN+H, and the ridge which separates the abstraction and H2CN dissociation pathways. Other features of the potential surface, such as the HCNH cis and trans minima, and the pathways leading to the formation of HNC+H are also described, though less accurately. Three different choices for the HHCN saddle point properties are considered. We find that the surface which matches the ab initio barrier energy most accurately gives rate constants that are too low. Much better agreement is obtained using a 3.2 kcal/mol barrier. The trajectory results show typical dependence of the CN+H2 reactive cross sections on initial translational energy and initial vibration/rotation state, with CN behaving as a spectator and H2 playing an active role in the reaction dynamics. Analysis of the H+HCN products indicates that both the C-H stretch and bend modes are significantly excited, with bend excitation showing strong sensitivity to the saddle point properties and to reagent translational energy. At translational energies below 20 kcal/mol, direct H abstraction is strongly favored over addition elimination.

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

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

M3 - Article

VL - 105

SP - 558

EP - 571

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 2

ER -