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

Collisional excitation of HNC by He found to be stronger than for structural isomer HCN in experiments at the low temperatures of interstellar space

HCN and its unstable isomer HNC are widely observed throughout the interstellar medium, with the HNC/HCN abundance ratio correlating strongly with temperature. In very cold environments HNC can even appear more abundant than HCN. Here we use a chirped pulse Fourier transform spectrometer to measure the pressure broadening of HCN and HNC, simultaneously formed in situ by laser photolysis and cooled to low temperatures in uniform supersonic flows of helium. Despite the apparent similarity of these systems, we find the HNC-He cross section to be more than twice as big as the HCN-He cross section at 10 K, confirming earlier quantum calculations. Our experimental results are supported by high-level scattering calculations and are also expected to apply with para-H₂, demonstrating that HCN and HNC have different collisional excitation properties that strongly influence the derived interstellar abundances.

Rotational (de)-excitation of C5 by collision with He at low temperature

An appropriate estimation of the abundance of the observed C₅ radical in the interstellar medium requires accurate radiative and collisional rate coefficients. We present the first two-dimensional potential energy surface (2D-PES) for the ground electronic state of the C₅(X¹Σ⁺)-He(X¹S) van der Waals system, obtained using an explicitly correlated coupled-cluster method with single, double, and perturbative triple excitations (RCCSD(T)-F12). This PES is subsequently used in quantum close-coupling (CC) scattering calculations. Collisional excitation cross-sections of the rotational levels of C₅ by He were calculated for energies up to 1500 cm^-1 using the standard (CC) method. The thermal dependence of the corresponding rate coefficients is given for the low and moderate temperature T ≤ 300 K regime of interstellar molecular clouds. This is the first study on the collisional rate coefficients for this system and may have important implications for the astrophysical detection of C₅(X¹Σ⁺) and modeling of carbon-rich media.

Fitting potential energy and induced dipole surfaces of the van der Waals complex CH4-N2 using non-product quadrature grids

We present an extensive study of the five-dimensional potential energy and induced dipole surfaces of the CH₄-N₂ complex assuming rigid-rotor approximation. Within the supermolecular approach, ab initio calculations of the interaction energies and dipoles were carried out at the CCSD(T)-F12 and CCSD(T) levels of theory using the correlation-consistent aug-cc-pVTZ basis set, respectively. Both potential energy and induced dipole surfaces inherit the symmetry of the molecular system and transform under the A₁⁺ and A₂⁺ irreducible representations of the molecular symmetry group G₄₈, respectively. One can take advantage of the symmetry when fitting the surfaces; first, when constructing angular basis functions and second, when selecting the grid points. The approach to the construction of scalar and vectorial basis functions exploiting the eigenfunction method [Q. Chen, J. Ping and F. Wang, Group Representation Theory for Physicists, World Scientific, 2nd edn, 2002] is developed. We explore the use of Sobolev-type quadrature grids as building blocks of robust quadrature rules adapted to the symmetry of the molecular system. Temperature variations of the cross second virial coefficient and first classical spectral moments of the rototranslational collision-induced band were derived. A reasonable agreement between calculated values and experimental data was found attesting to the high quality of constructed surfaces.

Interaction of the simple carbene c-C3H2with H2: potential energy surface and low-energy scattering

Cyclopropenylidene, c-C₃H₂, is a simple hydrocarbon, ubiquitous in astrophysical gases, and possessing a permanent electric dipole moment.

Potential energy surface of the CO2-N2 van der Waals complex

Four-dimensional potential energy surface (4D-PES) of the atmospherically relevant CO2-N2 van der Waals complex is generated using the explicitly correlated coupled cluster with single, double, and perturbative triple excitation (CCSD(T)-F12) method in conjunction with the augmented correlation consistent triple zeta (aug-cc-pVTZ) basis set. This 4D-PES is mapped along the intermonomer coordinates. An analytic fit of this 4D-PES is performed. Our extensive computations confirm that the most stable form corresponds to a T-shape structure where the nitrogen molecule points towards the carbon atom of CO2. In addition, we located a second isomer and two transition states in the ground state PES of CO2-N2. All of them lay below the CO2 + N2 dissociation limit. This 4D-PES is flat and strongly anisotropic along the intermonomer coordinates. This results in the possibility of the occurrence of large amplitude motions within the complex, such as the inversion of N2, as suggested in the recent spectroscopic experiments. Finally, we show that the experimentally established deviations from the C2v structure at equilibrium for the most stable isomer are due to the zero-point out-of-plane vibration correction.

Rotationally adiabatic pair interactions of para- and ortho-hydrogen with the halogen molecules F2, Cl2, and Br2

The present work is concerned with the weak interactions between hydrogen and halogen molecules, i.e., the interactions of pairs H2-X2 with X = F, Cl, Br, which are dominated by dispersion and quadrupole-quadrupole forces. The global minimum of the four-dimensional (4D) coupled cluster with singles and doubles and perturbative triples (CCSD(T)) pair potentials is always a T shaped structure where H2 acts as the hat of the T, with well depths (De) of 1.3, 2.4, and 3.1 kJ/mol for F2, Cl2, and Br2, respectively. MP2/AVQZ results, in reasonable agreement with CCSD(T) results extrapolated to the basis set limit, are used for detailed scans of the potentials. Due to the large difference in the rotational constants of the monomers, in the adiabatic approximation, one can solve the rotational Schrödinger equation for H2 in the potential of the X2 molecule. This yields effective two-dimensional rotationally adiabatic potential energy surfaces where pH2 and oH2 are point-like particles. These potentials for the H2-X2 complexes have global and local minima for effective linear and T-shaped complexes, respectively, which are separated by 0.4-1.0 kJ/mol, where oH2 binds stronger than pH2 to X2, due to higher alignment to minima structures of the 4D-pair potential. Further, we provide fits of an analytical function to the rotationally adiabatic potentials.