Accurate ab initio potential energy surface, thermochemistry, and dynamics of the Br(2P, 2P3/2) CH4 → HBr CH3 reaction

Research output: Contribution to journalArticle

23 Citations (Scopus)

Abstract

Chemically accurate full-dimensional non-spin-orbit and spin-orbit (SO) ground-state potential energy surfaces (PESs) are obtained for the Br + CH 4 → HBr CH3 reaction by fitting 21 574 composite ab initio energy points. The composite method considers electron correlation methods up to CCSD(T), basis sets up to aug-cc-pwCVTZ-PP, correlation of the core electrons, scalar relativistic effects via an effective core potential (ECP), and SO corrections, thereby achieving an accuracy better than 0.5 kcalmol. Benchmark structures and relative energies are computed for the stationary points using the ab initio focal-point analysis (FPA) scheme based on both ECP and Douglas-Kroll approaches providing all-electron relativistic CCSDT(Q)complete-basis-set quality energies. The PESs accurately describe the saddle point of the abstraction reaction and the van der Waals complexes in the entrance and product channels. The SO-corrected PES provides a classical barrier height of 7285(7232 ± 50) cm-1, De values of 867(799 ± 10) and 399(344 ± 10) cm-1 for the complexes CH3-HBr and CH3-BrH, respectively, and reaction endothermicity of 7867(7857 ± 50) cm-1, in excellent agreement with the new, FPA-based benchmark data shown in parentheses. The difference between the Br CH4 asymptotes of the non-SO and SO PESs is 1240 cm-1, in good agreement with the experiment (1228 cm-1). Quasiclassical trajectory calculations based on more than 13 million trajectories for the late-barrier Br CH4(vk= 0, 1) [k= 1, 2, 3, 4] reactions show that the vibrational energy, especially the excitation of the stretching modes, activates the reaction much more efficiently than translational energy, in agreement with the extended Polanyi rules. Angular distributions show dominant backward scattering for the ground-state reaction and forward scattering for the stretching-excited reactions. The reactivity on the non-SO PES is about 3-5 times larger than that on the SO PES in a wide collision energy range of 8000-16 000 cm-1.

Original languageEnglish
Article number134301
JournalThe Journal of Chemical Physics
Volume138
Issue number13
DOIs
Publication statusPublished - Apr 7 2013

Fingerprint

Thermochemistry
Potential energy surfaces
thermochemistry
Orbits
potential energy
orbits
Ground state
Stretching
energy
Trajectories
trajectories
Forward scattering
Electron correlations
asymptotes
electrons
composite materials
ground state
Electrons
Correlation methods
Angular distribution

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

@article{746d504ee28743d4b5eb279a963cac06,
title = "Accurate ab initio potential energy surface, thermochemistry, and dynamics of the Br(2P, 2P3/2) CH4 → HBr CH3 reaction",
abstract = "Chemically accurate full-dimensional non-spin-orbit and spin-orbit (SO) ground-state potential energy surfaces (PESs) are obtained for the Br + CH 4 → HBr CH3 reaction by fitting 21 574 composite ab initio energy points. The composite method considers electron correlation methods up to CCSD(T), basis sets up to aug-cc-pwCVTZ-PP, correlation of the core electrons, scalar relativistic effects via an effective core potential (ECP), and SO corrections, thereby achieving an accuracy better than 0.5 kcalmol. Benchmark structures and relative energies are computed for the stationary points using the ab initio focal-point analysis (FPA) scheme based on both ECP and Douglas-Kroll approaches providing all-electron relativistic CCSDT(Q)complete-basis-set quality energies. The PESs accurately describe the saddle point of the abstraction reaction and the van der Waals complexes in the entrance and product channels. The SO-corrected PES provides a classical barrier height of 7285(7232 ± 50) cm-1, De values of 867(799 ± 10) and 399(344 ± 10) cm-1 for the complexes CH3-HBr and CH3-BrH, respectively, and reaction endothermicity of 7867(7857 ± 50) cm-1, in excellent agreement with the new, FPA-based benchmark data shown in parentheses. The difference between the Br CH4 asymptotes of the non-SO and SO PESs is 1240 cm-1, in good agreement with the experiment (1228 cm-1). Quasiclassical trajectory calculations based on more than 13 million trajectories for the late-barrier Br CH4(vk= 0, 1) [k= 1, 2, 3, 4] reactions show that the vibrational energy, especially the excitation of the stretching modes, activates the reaction much more efficiently than translational energy, in agreement with the extended Polanyi rules. Angular distributions show dominant backward scattering for the ground-state reaction and forward scattering for the stretching-excited reactions. The reactivity on the non-SO PES is about 3-5 times larger than that on the SO PES in a wide collision energy range of 8000-16 000 cm-1.",
author = "G{\'a}bor Czak{\'o}",
year = "2013",
month = "4",
day = "7",
doi = "10.1063/1.4797467",
language = "English",
volume = "138",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "13",

}

TY - JOUR

T1 - Accurate ab initio potential energy surface, thermochemistry, and dynamics of the Br(2P, 2P3/2) CH4 → HBr CH3 reaction

AU - Czakó, Gábor

PY - 2013/4/7

Y1 - 2013/4/7

N2 - Chemically accurate full-dimensional non-spin-orbit and spin-orbit (SO) ground-state potential energy surfaces (PESs) are obtained for the Br + CH 4 → HBr CH3 reaction by fitting 21 574 composite ab initio energy points. The composite method considers electron correlation methods up to CCSD(T), basis sets up to aug-cc-pwCVTZ-PP, correlation of the core electrons, scalar relativistic effects via an effective core potential (ECP), and SO corrections, thereby achieving an accuracy better than 0.5 kcalmol. Benchmark structures and relative energies are computed for the stationary points using the ab initio focal-point analysis (FPA) scheme based on both ECP and Douglas-Kroll approaches providing all-electron relativistic CCSDT(Q)complete-basis-set quality energies. The PESs accurately describe the saddle point of the abstraction reaction and the van der Waals complexes in the entrance and product channels. The SO-corrected PES provides a classical barrier height of 7285(7232 ± 50) cm-1, De values of 867(799 ± 10) and 399(344 ± 10) cm-1 for the complexes CH3-HBr and CH3-BrH, respectively, and reaction endothermicity of 7867(7857 ± 50) cm-1, in excellent agreement with the new, FPA-based benchmark data shown in parentheses. The difference between the Br CH4 asymptotes of the non-SO and SO PESs is 1240 cm-1, in good agreement with the experiment (1228 cm-1). Quasiclassical trajectory calculations based on more than 13 million trajectories for the late-barrier Br CH4(vk= 0, 1) [k= 1, 2, 3, 4] reactions show that the vibrational energy, especially the excitation of the stretching modes, activates the reaction much more efficiently than translational energy, in agreement with the extended Polanyi rules. Angular distributions show dominant backward scattering for the ground-state reaction and forward scattering for the stretching-excited reactions. The reactivity on the non-SO PES is about 3-5 times larger than that on the SO PES in a wide collision energy range of 8000-16 000 cm-1.

AB - Chemically accurate full-dimensional non-spin-orbit and spin-orbit (SO) ground-state potential energy surfaces (PESs) are obtained for the Br + CH 4 → HBr CH3 reaction by fitting 21 574 composite ab initio energy points. The composite method considers electron correlation methods up to CCSD(T), basis sets up to aug-cc-pwCVTZ-PP, correlation of the core electrons, scalar relativistic effects via an effective core potential (ECP), and SO corrections, thereby achieving an accuracy better than 0.5 kcalmol. Benchmark structures and relative energies are computed for the stationary points using the ab initio focal-point analysis (FPA) scheme based on both ECP and Douglas-Kroll approaches providing all-electron relativistic CCSDT(Q)complete-basis-set quality energies. The PESs accurately describe the saddle point of the abstraction reaction and the van der Waals complexes in the entrance and product channels. The SO-corrected PES provides a classical barrier height of 7285(7232 ± 50) cm-1, De values of 867(799 ± 10) and 399(344 ± 10) cm-1 for the complexes CH3-HBr and CH3-BrH, respectively, and reaction endothermicity of 7867(7857 ± 50) cm-1, in excellent agreement with the new, FPA-based benchmark data shown in parentheses. The difference between the Br CH4 asymptotes of the non-SO and SO PESs is 1240 cm-1, in good agreement with the experiment (1228 cm-1). Quasiclassical trajectory calculations based on more than 13 million trajectories for the late-barrier Br CH4(vk= 0, 1) [k= 1, 2, 3, 4] reactions show that the vibrational energy, especially the excitation of the stretching modes, activates the reaction much more efficiently than translational energy, in agreement with the extended Polanyi rules. Angular distributions show dominant backward scattering for the ground-state reaction and forward scattering for the stretching-excited reactions. The reactivity on the non-SO PES is about 3-5 times larger than that on the SO PES in a wide collision energy range of 8000-16 000 cm-1.

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

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

U2 - 10.1063/1.4797467

DO - 10.1063/1.4797467

M3 - Article

VL - 138

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 13

M1 - 134301

ER -