Theoretical investigation of weakly-bound complexes of O(3P) with H2

Imaging Resonance Effects in C + H2 Collisions Using a Zeeman Decelerator

An intriguing phenomenon in molecular collisions is the occurrence of scattering resonances, which originate from bound and quasi-bound states supported by the interaction potential at low collision energies. The resonance effects in the scattering behavior are extraordinarily sensitive to the interaction potential, and their observation provides one of the most stringent tests for theoretical models. We present high-resolution measurements of state-resolved angular scattering distributions for inelastic collisions between Zeeman-decelerated C(3P1) atoms and para-H2 molecules at collision energies ranging from 77 cm-1 down to 0.5 cm-1. Rapid variations in the angular distributions were observed, which can be attributed to the consecutive reduction of contributing partial waves and effects of scattering resonances. The measurements showed excellent agreement with distributions predicted by ab initio quantum scattering calculations. However, discrepancies were found at specific collision energies, which most likely originate from an incorrectly predicted quasi-bound state. These observations provide exciting prospects for further high-precision and low-energy investigations of scattering processes that involve paramagnetic species.

Gaseous transport properties of the ground and excited Cr, Co and Ni cations in He: Ab initio study of electronic state chromatography

The electronic state chromatography (ESC) effect allows the differentiation of ions in their ground and metastable states by their gaseous mobilities in the limit of low electrostatic fields. It is investigated here by means of accurate transport calculations with ab initio ion-atom potentials for the Cr, Co and Ni cations in He buffer gas near room temperature. The values for the open-shell ions in degenerate states are shown to be well approximated by using the single isotropic interaction potential. Minimalistic implementation of the multireference configuration interaction (MRCI) method is enough to describe the zero-field transport properties of metastable ions in the 3dm-14s configuration, such as Cr+(a6D), Co+(a5F) and Ni+(4F), due to their weak and almost isotropic interaction with He atom and the low sensitivity of the measured mobilities to the potential well region. By contrast, interactions involving the ions in the ground 3dm states, such as Cr+(a6S), Co+(a3F) and Ni+(2D), are strong and anisotropic; the MRCI potentials poorly describe their transport coefficients. Even the coupled cluster with singles, doubles and non-iterative triples [CCSD(T)] approach taking into account vectorial spin-orbit coupling may not be accurate enough, as shown here for Ni+(2D). The sensitivity of ion mobility and the ESC effect to interaction potentials, similarities in ion-He interactions of the studied ions in distinct configurations, accuracy and possible improvements of the ab initio schemes, and control of the ESC effect by macroscopic parameters are discussed. Extensive sets of improved interaction potentials and transport data are generated.

Extended Mulliken-Hush Method with Applications to the Theoretical Study of Electron Transfer

A novel adiabatic-to-diabatic (ATD) transformation strategy, namely, the extended Mulliken-Hush (XMH) method, is proposed to evaluate diabatic properties including electronic couplings, potential energy surfaces, and their crossings. The XMH method is developed by adopting our recently proposed ATD transformation formula of a general vectorial physical observable, in which a useful ATD transformation is further determined by using an auxiliary dipole between localized frontier orbitals as a simple approximation of the diabatic transition dipole. The XMH method is simple and practical that provides a flexible way to construct diabatic states. To some extent, it can be regarded as an extension of the generalized Mulliken-Hush (GMH) method since the latter takes a stronger approximation, in which the diabatic transition dipole is assumed to be vanishing. Test calculations on the HeH₂⁺ system show that the electronic couplings predicted by the XMH method are closer to the ones calculated by the valence bond block-diagonalization approach than the GMH ones since the XMH method takes into account both the magnitude and direction of the diabatic transition dipole, which is consistent with the properties of this molecule. In the study of electron transfer in the two kinds of donor-bridge-acceptor systems, the XMH method maintains the simplicity of the GMH method and gives reasonable results even when the latter fails, wherein the diabatic transition dipole is nearly perpendicular to the difference of the initial and final adiabatic dipoles. More importantly, the XMH method can be easily combined with high-level electronic structure methods, in which the properties of the ground and excited states may be more accurately calculated, and hence, one may expect that further development of the XMH method would result in a general computational model for studying electron transfer reactions.

The effect of nonadiabaticity on the C+ + HF reaction

The chemistry of fluorine in the interstellar medium is particularly simple, with only a few key species and important reactions. Of the latter, the rate of the reaction of C⁺ ions with HF is not well established but is one of the key reactions that sets the relative abundance of HF and the CF⁺ ion, the two fluorine-bearing species that have been observed in interstellar clouds. The C⁺ + HF → CF⁺ + H reaction proceeds through a deeply bound HCF⁺ well. In this work, statistical methods, namely, the statistical adiabatic channel method originally developed by Quack and Troe and the quantum statistical method of Manolopoulos and co-workers, are applied to compute the total cross section as a function of energy for this reaction. This reaction proceeds on the ground 1² A' potential energy surface (PES), and there are also two non-reactive PES's, 1² A″ and 2² A', correlating with the C⁺(² P _1/2,3/2) + HF reactants. Two sets of scattering calculations were carried out, namely, a single-surface calculation on the 1² A' PES and the one in which all three PES's and the spin-orbit splitting of C⁺ are included in the description of the entrance channel. In the latter, reactivity of the spin-orbit excited ² P _3/2 level can be computed, and not just assumed to be zero, as in the single-state adiabatic approximation.

Probing Nonadiabatic Effects in Low-Energy C(3 P j) + H2 Collisions

Nonadiabatic effects are of fundamental interest in collision dynamics. In particular, inelastic collisions between open-shell atoms and molecules, such as the collisional excitation of C(³ P _j) by H₂, are governed by nonadiabatic and spin-orbit couplings that are the sole responsible of collisional energy transfer. Here, we study collisions between carbon in its ground state C(³ P _j=0) and molecular hydrogen (H₂) at low collision energies that result in spin-orbit excitation to C(³ P _j=1) and C(³ P _j=2). State-to-state integral cross sections are obtained experimentally from crossed-beam experiments with a source of almost pure beam of C(³ P _j=0) and theoretically from highly accurate quantum calculations. We observe very good agreement between experimental and theoretical data that demonstrates our ability to model nonadiabatic dynamics. New rate coefficients at temperatures relevant to astrochemical modeling are also provided. They should lead to an increase of the abundance of atomic C(³ P) derived from the observations of interstellar clouds and a decrease of the efficiency of the cooling of the interstellar gas due to carbon atoms.

State-to-State Quantum Dynamics of Reactions O((3)P) + HD (v = 0-1, j = 0) → OH+D and OD+H: Reaction Mechanism and Vibrational Excitation

Time-dependent quantum wave packet dynamics calculations have been performed in order to characterize the dynamics and mechanism of O((3)P) + HD (v = 0-1, j = 0) → OH+D and OD+H reactive collisions using the adiabatic potential energy surface by Rogers et al. [J. Phys. Chem. A 2000, 104, 2308] Special attention has been paid to the calculations and discussion of the state resolved integral and differential cross sections and the product state distributions. In addition, the intramolecular isotopic branching ratio has been determined. The results revealed that the OD + H is the favored product channel and the product OH has the same quantum number v as the reactant HD. For low collision energy, the product angular distributions concentrate in the backward region being consistent with a rebounding mechanism. In the case of higher collision energy, the stripping collisions with larger impact parameters tend to produce sideways and forward scatterings, especially for the HD vibrationally excited state. The cross section and intramolecular isotopic branching ratio are in agreement with the previous theoretical results. A cartoon depiction collision model is built and works well for our calculation results.

Quantum reactive scattering of O(3P)+H2 at collision energies up to 4.4 eV

We report the results of quantum scattering calculations for the O((3)P)+H2 reaction for a range of collision energies from 0.4 to 4.4 eV, important for astrophysical and atmospheric processes. The total and state-to-state reactive cross sections are calculated using a fully quantum time-independent coupled-channel approach on recent potential energy surfaces of (3)A' and (3)A″ symmetry. A larger basis set than in the previous studies was used to ensure single-surface convergence at higher energies. Our results agree well with the published data at lower energies and indicate the breakdown of reduced dimensionality approach at collision energies higher than 1.5 eV. Differential cross sections and momentum transfer cross sections are also reported.

van der Waals interactions and dipole polarizabilities of lanthanides: Tm(F2)–He and Yb(S1)–He potentials

Anisotropic dipole polarizabilities of Tm(2F), Tm+2(2F), and Yb(1S) are calculated using the finite-field multireference averaged quadratic coupled cluster (MR-AQCC) (Tm and Tm+2) and RCCSD(T) (Yb) methods with small-core relativistic pseudopotentials ECP28MWB combined with the augmented ANO basis sets. The lanthanide atoms are strongly polarizable with the scalar part originating from the 6s electrons and the tensorial part from the open 4f shells. The adiabatic interaction potentials 2Sigma+, 2Pi, 2Delta, and 2Phi of Tm(2F)-He and Tm+2(2F)-He were examined by the multireference approaches, multireference configuration interaction and MR-AQCC, using the basis sets designed in the polarizability calculations. A closed-shell lanthanide system Yb(1S)-He was included for comparison. The Tm-He 2Sigma+, 2Pi, 2Delta, and 2Phi interaction potentials are very shallow and nearly degenerate (within 0.01 cm(-1)), with the well depths in the range of 2.35-2.36 cm(-1) at R=6.17 A. The basis-set saturated well depths are expected to be larger by ca. 25%, as estimated using the bond-function augmented basis set. The interactions of lanthanide atoms with He are one order of magnitude less anisotropic than those involving first-row transition metal atoms. The suppression of anisotropy is chiefly attributed to the screening effected by the 6s shell. When these electrons are removed as in the di-cation complex Tm+2(2F)-He, the potentials deepen to a thousand wave number range and their anisotropy is enhanced 500-fold.

An ab initio investigation of the O(3P)–H2(1Σ+g) van der Waals well

We report an ab initio study of the van der Waals region of the O(3P)-H2 potential energy surface based on RCCSD(T) calculations with an aug-cc-pVQZ basis supplemented by bond functions. In addition, an open-shell implementation of symmetry-adapted perturbation theory (SAPT) is used to corroborate the RCCSD(T) calculations and to investigate the relative magnitudes of the various contributions to the van der Waals interaction. We also investigate the effect of the spin-orbit coupling on the position and depth of the van der Waals well. We predict the van der Waals minimum to occur in perpendicular geometry, and located at a closer distance than a secondary well in colinear geometry. The potentials obtained in the present study confirm the previous calculations of Alexander [M. H. Alexander, J. Chem. Phys., 1998, 108, 4467], but disagree with the earlier work of Harding and co-workers [Z. Li, V. A. Apkarian and L. B. Harding, J. Chem. Phys., 1997, 106, 942] as well as with recently refitted surfaces of Brandão and coworkers [J. Brandão, C. Mogo and B. C. Silva, J. Chem. Phys., 2004, 121, 8861]. Inclusion of spin-orbit coupling reduces the depth of the van der Waals minimum without causing a change in its position.

Interactions in open-shell clusters: ab initio study of pre-reactive complex O(3P) + HCl

Van der Waals interactions between the ground-state triplet O(3P) atom and the closed-shell HCl molecule are investigated in the pre-reactive region. Three adiabatic (two of A'' symmetry and one of A' symmetry) and four non-relativistic diabatic potential energy surfaces are obtained by combining a restricted open-shell coupled cluster approach with the multireference configuration interaction method. The lower A'' adiabatic potential surface has a single minimum (D(e) = 589 cm(-1)) for a linear O...HCl configuration. The upper A'' potential has a weak (D(e) = 65 cm(-1)) minimum for a linear HCl...O configuration. The A' adiabatic potential has a weak (124 cm(-1)) T-shaped minimum. Adiabatic potentials intersect once in the O...HCl linear configuration and twice in the linear HCl...O geometry. The role of electrostatic interactions in shaping these potentials is discussed. The effects of spin-orbit coupling on this interaction are also investigated assuming a constant value of the SO parameter.