Intermolecular potential energy surfaces and bound states in F–HF

Ab Initio Treatment of the Chemical Reaction Precursor Complex Br(2P)−HCN. 1. Adiabatic and Diabatic Potential Surfaces

The three adiabatic potential surfaces of the Br(2P)-HCN complex that correlate to the 2P ground state of the Br atom were calculated ab initio. With the aid of a geometry-dependent diabatic mixing angle, also calculated ab initio, these adiabatic potential surfaces were transformed into a set of four diabatic potential surfaces required to define the full 3 x 3 matrix of diabatic potentials. Each of these diabatic potential surfaces was expanded in terms of the appropriate spherical harmonics in the atom-linear molecule Jacobi angle theta. The dependence of the expansion coefficients on the distance R between Br and the HCN center of mass and on the CH bond length was fit to an analytic form. For HCN in its equilibrium geometry, the global minimum with De = 800.4 cm(-1) and Re = 6.908a0 corresponds to a linear Br-NCH geometry, with an electronic ground state of Sigma symmetry. A local minimum with De = 415.1 cm-1, Re = 8.730a0, and a twofold degenerate Pi ground state is found for the linear Br-HCN geometry. The binding energy, De, depends strongly on the CH bond length for the Br-HCN complex and much less strongly for the Br-NCH complex, with a longer CH bond giving stronger binding for both complexes. Spin-orbit coupling was included and diabatic states were constructed that correlate to the ground 2P3/2 and excited 2P1/2 spin-orbit states of the Br atom. For the ground spin-orbit state with electronic angular momentum j = (3/2) the minimum in the potential for projection quantum number omega = +/-(3/2) coincides with the local minimum for linear Br-HCN of the spin-free case. The minimum in the potential for projection quantum number omega = +/-(1/2) occurs for linear Br-NCH but is considerably less deep than the global minimum of the spin-free case. According to the lowest spin-orbit coupling included adiabatic potential the two linear isomers, Br-NCH and Br-HCN, are about equally stable. In the subsequent paper, we use these potentials in calculations of the rovibronic states of the Br-HCN complex.

Ab initio treatment of the chemical reaction precursor complex Cl(2P)-HF. 1. Three-dimensional diabatic potential energy surfaces

The three adiabatic potential surfaces of the Cl(2P)-HF complex that correlate with the 2P ground state of the Cl atom were calculated with the ab initio RCCSD(T) method (partially spin-restricted coupled cluster theory including single and double excitations and perturbative correction for the triples). With the aid of a geometry-dependent diabatic mixing angle, calculated by the complete active space self-consistent field (CASSCF) and multireference configuration-interaction (MRCI) methods, these adiabatic potential surfaces were converted to a set of four distinct diabatic potential surfaces required to define the full 3 x 3 matrix of diabatic potentials. Each of these diabatic potential surfaces was expanded in terms of the appropriate spherical harmonics in the angle theta between the HF bond axis r and the Cl-HF intermolecular axis R. The dependence of the expansion coefficients on the Cl-HF distance R and the HF bond length r(HF) was fit to an analytic form. The strongest binding occurs for the hydrogen-bonded linear Cl-HF geometry, with D(e) = 676.5 cm(-1) and R(e) = 6.217 a0 when r(HF) = r(e) = 1.7328 a0. This binding energy D(e) depends strongly on r(HF), with larger r(HF) causing stronger binding. An important contribution to the binding energy is provided by the interaction between the quadrupole moment of the Cl(2P) atom and the dipole of HF. In agreement with this electrostatic picture, the ground state of linear Cl-HF is a 2-fold degenerate electronic Pi state. For the linear Cl-FH geometry the states are in opposite order, i.e., the Sigma state is lower in energy than the Pi state. The following paper in this issue describes full three-dimensional computations of the bound states of the Cl-HF complex, based on the ab initio diabatic potentials of this paper.

Ab Initio Treatment of the Chemical Reaction Precursor Complex Br(2P)−HCN. 2. Bound-State Calculations and Infrared Spectra

Rovibronic energy levels and properties of the Br(2P)-HCN complex were obtained from three-dimensional calculations, with HCN kept linear and the CN bond frozen. All diabatic states that correlate to the 2P3/2 and 2P1/2 states of the Br atom were included and spin-orbit coupling was taken into account. The 3 x 3 matrix of diabatic potential surfaces was taken from the preceding paper (paper 1). In agreement with experiment, we found two linear isomers, Br-NCH and Br-HCN. The calculated binding energies are very similar: D0 = 352.4 cm(-1) and D0 = 349.1 cm(-1), respectively. We established, also in agreement with experiment, that the ground electronic state of Br-NCH has |Omega| = (1/2) and that Br-HCN has a ground state with |Omega| = (3/2), where the quantum number, Omega, is the projection of the total angular momentum, J, of the complex on the intermolecular axis R. This picture can be understood as being caused by the electrostatic interaction between the quadrupole of the Br(2P) atom and the dipole of HCN, combined with the very strong spin-orbit coupling in Br. We predicted the frequencies of the van der Waals modes of both isomers and found a direct Renner-Teller splitting of the bend mode in Br-HCN and a smaller, indirect, splitting in Br-NCH. The red shift of the CH stretch frequency in the complex, relative to free HCN, was calculated to be 1.98 cm(-1) for Br-NCH and 23.11 cm(-1) for Br-HCN, in good agreement with the values measured in helium nanodroplets. Finally, with the use of the same potential surfaces, we modeled the Cl(2P)-HCN complex and found that the experimentally observed linear Cl-NCH isomer is considerably more stable than the (not observed) Cl-HCN isomer. This was explained mainly as an effect of the substantially smaller spin-orbit coupling in Cl, relative to Br.

A high-resolution infrared spectroscopic investigation of the halogen atom–HCN entrance channel complexes solvated in superfluid helium droplets

Rotationally resolved infrared spectra are reported for the X-HCN (X = Cl, Br, I) binary complexes solvated in helium nanodroplets. These results are directly compared with those obtained previously for the corresponding X-HF complexes [J. M. Merritt, J. Küpper and R. E. Miller, Phys. Chem. Chem. Phys., 2005, 7, 67]. For bromine and iodine atoms complexed with HCN, two linear structures are observed and assigned to the (2)Sigma(1/2) and (2)Pi(3/2) ground electronic states of the nitrogen and hydrogen bound geometries, respectively. Experiments for HCN + chlorine atoms give rise to only a single band which is attributed to the nitrogen bound isomer. That the hydrogen bound isomer is not stabilized is rationalized in terms of a lowering of the isomerization barrier by spin-orbit coupling. Theoretical calculations with and without spin-orbit coupling have also been performed and are compared with our experimental results. The possibility of stabilizing high-energy structures containing multiple radicals is discussed, motivated by preliminary spectroscopic evidence for the di-radical Br-HCCCN-Br complex. Spectra for the corresponding molecular halogen HCN-X(2) complexes are also presented.

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.

Potential energy surface for H2O(3A″) from accurateab initiodata with inclusion of long-range interactions

A new potential energy surface for the O((3)P) + H(2) system in the lowest (3)A(") state is built using ab initio data calculated by Rogers et al. [J. Phys. Chem. A 104, 2308 (2000)] and the double many-body expansion formalism. It incorporates a semiempirical model of long-range interactions, which should play an important role at low collision energies. Preliminary quasiclassical trajectory results at 12.6 kcal/mol collision energy, show that the deeper van der Waals region described in this new surface translates into a four times higher cross section than that of Rogers' (3)A(") surface. To confirm this hypothesis, a second surface was calibrated. The two surfaces are fitted with rmsd<0.5 kcal/mol and differ mainly on the depth of the van der Waals region. That difference in the van der Waals region corresponds to a 22% lower cross section of the less deep surface, which is still three times higher than the equivalent results from Rogers' (3)A(") surface. This study reflects the importance of a correct description of van der Waals forces on potential energy surfaces.