Otoniel Denis-Alpizar

The interaction of He with vibrating HCN: potential energy surface, bound states and rotationally inelastic cross sections

A four-dimensional potential energy surface representing the interaction between He and hydrogen cyanide (HCN) subjected to bending vibrational motion is presented. Ab initio calculations were carried out at the coupled-cluster level with single and double excitations and a perturbative treatment of triple excitations, using a quadruple-zeta basis set and mid-bond functions. The global minimum is found in the linear He-HCN configuration with the H atom pointing towards helium at the intermolecular separation of 7.94 a0. The corresponding well depth is 30.35 cm(-1). First, the quality of the new potential has been tested by performing two comparisons with previous theoretical and experimental works. (i) The rovibrational energy levels of the He-HCN complex for a rigid linear configuration of the HCN molecule have been calculated. The dissociation energy is 8.99 cm(-1), which is slightly smaller than the semi-empirical value of 9.42 cm(-1). The transitions frequencies are found to be in good agreement with the experimental data. (ii) We performed close coupling calculations of the rotational de-excitation of rigid linear HCN in collision with He and observed a close similarity with the theoretical data published in a recent study. Second, the effects of the vibrational bending of HCN have been investigated, both for the bound levels of the He-HCN system and for the rotationally inelastic cross sections. This was performed with an approximate method using the average of the interaction potential over the vibrational bending wavefunction. While this improves slightly the comparison of calculated transitions frequencies with experiment, the cross sections remain very close to those obtained with rigid linear HCN.

Ro-vibrational relaxation of HCN in collisions with He: Rigid bender treatment of the bending-rotation interaction

We present a new theoretical method to treat atom-rigid bender inelastic collisions at the Close Coupling (RB-CC) level in the space fixed frame. The coupling between rotation and bending is treated exactly within the rigid bender approximation and we obtain the cross section for the rotational transition between levels belonging to different bending levels. The results of this approach are compared with those obtained when using the rigid bender averaged approximation (RBAA) introduced in our previous work dedicated to this system. We discuss the validity of this approximation and of the previous studies based on rigid linear HCN. We find that l-type transitions cross sections have to be calculated at the RB-CC level for the He-HCN collision while pure rotational transitions cross sections may be calculated accurately at the RBAA level.

On the use of explicitly correlated treatment methods for the generation of accurate polyatomic –He/H2 interaction potential energy surfaces: The case of C3–He complex and generalization

Through the study of the C3(X1Σg (+)) (1)Σg (+)) + He((1)S) astrophysical relevant system using standard (CCSD(T)) and explicitly correlated (CCSD(T)-F12) coupled cluster approaches, we show that the CCSD(T)-F12/aug-cc-pVTZ level represents a good compromise between accuracy and low computational cost for the generation of multi-dimensional potential energy surfaces (PESs) over both intra- and inter-monomer degrees of freedom. Indeed, the CCSD(T)-F12/aug-cc-pVTZ 2D-PES for linear C3 and the CCSD(T)-F12/aug-cc-pVTZ 4D-PES for bent C3 configurations gently approach those mapped at the CCSD(T)/aug-cc-pVXZ (X = T,Q) + bond functions level, whereas a strong reduction of computational effort is observed. After exact dynamical computations, the pattern of the rovibrational levels of the intermediate C3-He complex and the rotational and rovibrational (de-) excitation of C3 by He derived using both sets of PESs agree quite well. Since C3 shows a floppy character, the interaction PES is defined in four dimensions to obtain realistic collisional parameters. The C-C-C bending mode, which fundamental lies at 63 cm(-1) and can be excited at very low temperatures is explicitly considered as independent coordinate. Our work suggests hence that CCSD(T)-F12/aug-cc-pVTZ methodology is the key method for the generation of accurate polyatomic - He/H2 multi-dimensional PESs.

A new ab initio potential energy surface for the collisional excitation of HCN by para- and ortho-H2.

We present a new four-dimensional potential energy surface for the collisional excitation of HCN by H2. Ab initio calculations of the HCN-H2 van der Waals complex, considering both molecules as rigid rotors, were carried out at the explicitly correlated coupled cluster with single, double, and perturbative triple excitations [CCSD(T)-F12a] level of theory using an augmented correlation-consistent triple zeta (aVTZ) basis set. The equilibrium structure is linear HCN-H2 with the nitrogen pointing towards H2 at an intermolecular separation of 7.20 a0. The corresponding well depth is -195.20 cm(-1). A secondary minimum of -183.59 cm(-1) was found for a T-shape configuration with the H of HCN pointing to the center of mass of H2. We also determine the rovibrational energy levels of the HCN-para-H2 and HCN-ortho-H2 complexes. The calculated dissociation energies for the para and ortho complexes are 37.79 cm(-1) and 60.26 cm(-1), respectively. The calculated ro-vibrational transitions in the HCN-H2 complex are found to agree by more than 0.5% with the available experimental data, confirming the accuracy of the potential energy surface.

Rotational excitation of HCN by para- and ortho-H2

Rotational excitation of the hydrogen cyanide (HCN) molecule by collisions with para-H2(j = 0, 2) and ortho-H2(j = 1) is investigated at low temperatures using a quantum time independent approach. Both molecules are treated as rigid rotors. The scattering calculations are based on a highly correlated ab initio 4-dimensional (4D) potential energy surface recently published. Rotationally inelastic cross sections among the 13 first rotational levels of HCN were obtained using a pure quantum close coupling approach for total energies up to 1200 cm(-1). The corresponding thermal rate coefficients were computed for temperatures ranging from 5 to 100 K. The HCN rate coefficients are strongly dependent on the rotational level of the H2 molecule. In particular, the rate coefficients for collisions with para-H2(j = 0) are significantly lower than those for collisions with ortho-H2(j = 1) and para-H2(j = 2). Propensity rules in favor of even Δj transitions were found for HCN in collisions with para-H2(j = 0) whereas propensity rules in favor of odd Δj transitions were found for HCN in collisions with H2(j ⩾ 1). The new rate coefficients were compared with previously published HCN-para-H2(j = 0) rate coefficients. Significant differences were found due the inclusion of the H2 rotational structure in the scattering calculations. These new rate coefficients will be crucial to improve the estimation of the HCN abundance in the interstellar medium.

Rovibrational energy transfer in the He-C3 collision: Potential energy surface and bound states

We present a four-dimensional potential energy surface (PES) for the collision of C3 with He. Ab initio calculations were carried out at the coupled-cluster level with single and double excitations and a perturbative treatment of triple excitations, using a quadruple-zeta basis set and mid-bond functions. The global minimum of the potential energy is found to be -26.9 cm(-1) and corresponds to an almost T-shaped structure of the van der Waals complex along with a slightly bent configuration of C3. This PES is used to determine the rovibrational energy levels of the He-C3 complex using the rigid monomer approximation (RMA) and the recently developed atom-rigid bender approach at the Close Coupling level (RB-CC). The calculated dissociation energies are -9.56 cm(-1) and -9.73 cm(-1), respectively at the RMA and RB-CC levels. This is the first theoretical prediction of the bound levels of the He-C3 complex with the bending motion.

Potential energy surface and rovibrational energy levels of the H2-CS van der Waals complex.

Owing to its large dipole, astrophysicists use carbon monosulfide (CS) as a tracer of molecular gas in the interstellar medium, often in regions where H(2) is the most abundant collider. Predictions of the rovibrational energy levels of the weakly bound complex CS-H(2) (not yet observed) and also of rate coefficients for rotational transitions of CS in collision with H(2) should help to interpret the observed spectra. This paper deals with the first goal, i.e., the calculation of the rovibrational energy levels. A new four-dimensional intermolecular potential energy surface for the H(2)-CS complex is presented. Ab initio potential energy calculations were carried out at the coupled-cluster level with single and double excitations and a perturbative treatment of triple excitations, using a quadruple-zeta basis set and midbond functions. The potential energy surface was obtained by an analytic fit of the ab initio data. The equilibrium structure of the H(2)-CS complex is found to be linear with the carbon pointing toward H(2) at the intermolecular separation of 8.6 a(o). The corresponding well depth is -173 cm(-1). The potential was used to calculate the rovibrational energy levels of the para-H(2)-CS and ortho-H(2)-CS complexes. The present work provides the first theoretical predictions of these levels. The calculated dissociation energies are found to be 35.9 cm(-1) and 49.9 cm(-1), respectively, for the para and ortho complexes. The second virial coefficient for the H(2)-CS pair has also been calculated for a large range of temperature. These results could be used to assign future experimental spectra and to check the accuracy of the potential energy surface.