Coupled potential-energy surfaces and quantum reactive scattering for the Cl(2P) + HCl → ClH + Cl(2P) reaction

C. S. Maierle, George C Schatz, Mark S. Gordon, Patrick McCabe, J. N L Connor

Research output: Contribution to journalArticle

62 Citations (Scopus)

Abstract

We construct global multiple-coupled potential-energy surfaces for the three-dimensional Cl(2P) + HCl → ClH + Cl(2P) reaction which we use to generate results from quantum reactive scattering calculations. The three potential surfaces (12A′, 22A′ and 12A″) are derived from multi-configuration self-consistent field (MCSCF) electronic structure calculations for geometries where all atoms are strongly interacting. The MCSCF energies are used to construct diabatic potentials which are fitted with rotated Morse cubic spline functions. The diabatic surfaces are then combined with empirical long-range potentials due to Dubernet and Hutson, which are also used to determine empirical coupling surfaces. Spin-orbit interaction is added perturbatively to the resulting surfaces. Saddle-point properties on these fitted surfaces are scaled to reproduce experimental thermal rate coefficient data based on the quantum scattering calculations. The quantum scattering calculations are based on a coupled-channel hyperspherical coordinate (CCH) method. The basis expansion used in the scattering calculation contains diatomic vibration and rotation, and the atomic fine-structure states, including full angular momentum coupling of rotation with electronic spin and orbital angular momentum. All terms in the combined rovibronic Hamiltonian are evaluated explicitly, including spin-orbit coupling, Coriolis and electrostatic interactions. State-selected and total cumulative reaction probabilities and thermal rate coefficients are calculated using a J-shifting approximation, which starts from converged CCH results for total angular momentum quantum number J = 1/2. The calculated rate coefficients indicate that spin-state change collisions account for a few per cent of the thermal rate coefficient at 300 K. Transition-state ClHCl resonances are lower in energy and broader than those on earlier (single) surfaces, and there is a more extensive progression in the symmetric stretch mode.

Original languageEnglish
Pages (from-to)709-720
Number of pages12
JournalJournal of the Chemical Society - Faraday Transactions
Volume93
Issue number5
Publication statusPublished - Mar 7 1997

Fingerprint

Potential energy surfaces
potential energy
Scattering
Angular momentum
scattering
angular momentum
Orbits
coefficients
Hamiltonians
self consistent fields
Coulomb interactions
spline functions
Splines
Electronic structure
saddle points
configurations
spin-orbit interactions
progressions
quantum numbers
Atoms

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

Coupled potential-energy surfaces and quantum reactive scattering for the Cl(2P) + HCl → ClH + Cl(2P) reaction. / Maierle, C. S.; Schatz, George C; Gordon, Mark S.; McCabe, Patrick; Connor, J. N L.

In: Journal of the Chemical Society - Faraday Transactions, Vol. 93, No. 5, 07.03.1997, p. 709-720.

Research output: Contribution to journalArticle

@article{8c17e32284a346638039f9b5547f0909,
title = "Coupled potential-energy surfaces and quantum reactive scattering for the Cl(2P) + HCl → ClH + Cl(2P) reaction",
abstract = "We construct global multiple-coupled potential-energy surfaces for the three-dimensional Cl(2P) + HCl → ClH + Cl(2P) reaction which we use to generate results from quantum reactive scattering calculations. The three potential surfaces (12A′, 22A′ and 12A″) are derived from multi-configuration self-consistent field (MCSCF) electronic structure calculations for geometries where all atoms are strongly interacting. The MCSCF energies are used to construct diabatic potentials which are fitted with rotated Morse cubic spline functions. The diabatic surfaces are then combined with empirical long-range potentials due to Dubernet and Hutson, which are also used to determine empirical coupling surfaces. Spin-orbit interaction is added perturbatively to the resulting surfaces. Saddle-point properties on these fitted surfaces are scaled to reproduce experimental thermal rate coefficient data based on the quantum scattering calculations. The quantum scattering calculations are based on a coupled-channel hyperspherical coordinate (CCH) method. The basis expansion used in the scattering calculation contains diatomic vibration and rotation, and the atomic fine-structure states, including full angular momentum coupling of rotation with electronic spin and orbital angular momentum. All terms in the combined rovibronic Hamiltonian are evaluated explicitly, including spin-orbit coupling, Coriolis and electrostatic interactions. State-selected and total cumulative reaction probabilities and thermal rate coefficients are calculated using a J-shifting approximation, which starts from converged CCH results for total angular momentum quantum number J = 1/2. The calculated rate coefficients indicate that spin-state change collisions account for a few per cent of the thermal rate coefficient at 300 K. Transition-state ClHCl resonances are lower in energy and broader than those on earlier (single) surfaces, and there is a more extensive progression in the symmetric stretch mode.",
author = "Maierle, {C. S.} and Schatz, {George C} and Gordon, {Mark S.} and Patrick McCabe and Connor, {J. N L}",
year = "1997",
month = "3",
day = "7",
language = "English",
volume = "93",
pages = "709--720",
journal = "Physical Chemistry Chemical Physics",
issn = "1463-9076",
publisher = "Royal Society of Chemistry",
number = "5",

}

TY - JOUR

T1 - Coupled potential-energy surfaces and quantum reactive scattering for the Cl(2P) + HCl → ClH + Cl(2P) reaction

AU - Maierle, C. S.

AU - Schatz, George C

AU - Gordon, Mark S.

AU - McCabe, Patrick

AU - Connor, J. N L

PY - 1997/3/7

Y1 - 1997/3/7

N2 - We construct global multiple-coupled potential-energy surfaces for the three-dimensional Cl(2P) + HCl → ClH + Cl(2P) reaction which we use to generate results from quantum reactive scattering calculations. The three potential surfaces (12A′, 22A′ and 12A″) are derived from multi-configuration self-consistent field (MCSCF) electronic structure calculations for geometries where all atoms are strongly interacting. The MCSCF energies are used to construct diabatic potentials which are fitted with rotated Morse cubic spline functions. The diabatic surfaces are then combined with empirical long-range potentials due to Dubernet and Hutson, which are also used to determine empirical coupling surfaces. Spin-orbit interaction is added perturbatively to the resulting surfaces. Saddle-point properties on these fitted surfaces are scaled to reproduce experimental thermal rate coefficient data based on the quantum scattering calculations. The quantum scattering calculations are based on a coupled-channel hyperspherical coordinate (CCH) method. The basis expansion used in the scattering calculation contains diatomic vibration and rotation, and the atomic fine-structure states, including full angular momentum coupling of rotation with electronic spin and orbital angular momentum. All terms in the combined rovibronic Hamiltonian are evaluated explicitly, including spin-orbit coupling, Coriolis and electrostatic interactions. State-selected and total cumulative reaction probabilities and thermal rate coefficients are calculated using a J-shifting approximation, which starts from converged CCH results for total angular momentum quantum number J = 1/2. The calculated rate coefficients indicate that spin-state change collisions account for a few per cent of the thermal rate coefficient at 300 K. Transition-state ClHCl resonances are lower in energy and broader than those on earlier (single) surfaces, and there is a more extensive progression in the symmetric stretch mode.

AB - We construct global multiple-coupled potential-energy surfaces for the three-dimensional Cl(2P) + HCl → ClH + Cl(2P) reaction which we use to generate results from quantum reactive scattering calculations. The three potential surfaces (12A′, 22A′ and 12A″) are derived from multi-configuration self-consistent field (MCSCF) electronic structure calculations for geometries where all atoms are strongly interacting. The MCSCF energies are used to construct diabatic potentials which are fitted with rotated Morse cubic spline functions. The diabatic surfaces are then combined with empirical long-range potentials due to Dubernet and Hutson, which are also used to determine empirical coupling surfaces. Spin-orbit interaction is added perturbatively to the resulting surfaces. Saddle-point properties on these fitted surfaces are scaled to reproduce experimental thermal rate coefficient data based on the quantum scattering calculations. The quantum scattering calculations are based on a coupled-channel hyperspherical coordinate (CCH) method. The basis expansion used in the scattering calculation contains diatomic vibration and rotation, and the atomic fine-structure states, including full angular momentum coupling of rotation with electronic spin and orbital angular momentum. All terms in the combined rovibronic Hamiltonian are evaluated explicitly, including spin-orbit coupling, Coriolis and electrostatic interactions. State-selected and total cumulative reaction probabilities and thermal rate coefficients are calculated using a J-shifting approximation, which starts from converged CCH results for total angular momentum quantum number J = 1/2. The calculated rate coefficients indicate that spin-state change collisions account for a few per cent of the thermal rate coefficient at 300 K. Transition-state ClHCl resonances are lower in energy and broader than those on earlier (single) surfaces, and there is a more extensive progression in the symmetric stretch mode.

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

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

M3 - Article

VL - 93

SP - 709

EP - 720

JO - Physical Chemistry Chemical Physics

JF - Physical Chemistry Chemical Physics

SN - 1463-9076

IS - 5

ER -