Hg + Br → Hg Br recombination and collision-induced dissociation dynamics

A global potential energy surface has been constructed for the system HgBr+Ar→Hg+Br+Ar to determine temperature dependent rate constants for the collision-induced dissociation (CID) and recombination of Hg and Br atoms. The surface was decomposed using a many-body expansion. Accurate two-body potent...

Full description

Bibliographic Details
Published in:The Journal of Chemical Physics
Main Authors: Shepler, Benjamin C., Balabanov, Nikolai B., Peterson, Kirk A.
Format: Article in Journal/Newspaper
Language:English
Published: AIP Publishing 2007
Subjects:
Online Access:http://dx.doi.org/10.1063/1.2777142
https://pubs.aip.org/aip/jcp/article-pdf/doi/10.1063/1.2777142/15404802/164304_1_online.pdf
Description
Summary:A global potential energy surface has been constructed for the system HgBr+Ar→Hg+Br+Ar to determine temperature dependent rate constants for the collision-induced dissociation (CID) and recombination of Hg and Br atoms. The surface was decomposed using a many-body expansion. Accurate two-body potentials for HgBr, HgAr, and ArBr were calculated using coupled cluster theory with single and double excitations and a perturbative treatment of triple excitations [CCSD(T)], as well as the multireference averaged coupled pair functional method. Correlation consistent basis sets were used to extrapolate to the complete basis set limit and corrections were included to account for scalar and spin-orbit relativistic effects, core-valence correlation, and the Lamb shift. The three-body potential was computed with the CCSD(T) method and triple-zeta quality basis sets. Quasiclassical trajectories using the final analytical potential surface were directly carried out on the CID of HgBr by Ar for a large sampling of initial rotational, vibrational, and collision energies. The recombination rate of Hg and Br atoms is a likely first step in mercury depletion events that have been observed in the Arctic troposphere during polar sunrise. The effective second order rate constant for this process was determined in this work from the calculated CID rate as a function of temperature using the principle of detailed balance, which resulted in k(T)=1.2×10−12cm3molecule−1s−1 at 260K and 1bar pressure.