This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design.

Precision vibrational spectroscopy of the molecular hydrogen ions is of significant interest for determining fundamental constants, for searching for new forces, and for testing quantum electrodynamics calculations. Future experiments can profit from the ability of preparing molecular hydrogen ions at ultralow kinetic energy and in preselected internal states, with respect to vibration, rotation, and spin degrees of freedom. For the homonuclear ions (H 2 + , D 2 + ), direct laser cooling of the rotational degree of freedom is not feasible. We show by quantum calculations that rotational cooling by cold He buffer gas is an effective approach. For this purpose we have computed the energy-dependent cross sections for rotationally elastic and inelastic collisions, h 2 + (v = 0, N) + He → h 2 + (v = 0, N ) + He (where h = H, D), using ab initio coupled-channel calculations. We find that rotational cooling to the lowest rotational state is possible within tens of seconds under experimentally realistic conditions. We furthermore describe possible protocols for the preparation of a single quantum state, where also the spin state is well defined.

We present quantum calculations for the inelastic collisions between H molecules, in rotationally excited internal states, and He atoms. This work is motivated by the possibility of experiments in which the molecular ions are stored and translationally cooled in an ion trap and a He buffer gas is added for deactivation of the internal rotational population, in particular at low (cryogenic) translational temperatures. We carry out an accurate representation of the forces at play from an ab initio description of the relevant potential energy surface, with the molecular ion in its ground vibrational state, and obtain the cross sections for state-changing rotationally inelastic collisions by solving the coupled channel quantum scattering equations. The presence of hyperfine and fine structure effects in both ortho- and para-H molecules is investigated and compared to the results where such a contribution is disregarded. An analysis of possible propensity rules that may predict the relative probabilities of inelastic events involving rotational state-changing is also carried out, together with the corresponding elastic cross sections from several initial rotational states. Temperature-dependent rotationally inelastic rates are then computed and discussed in terms of relative state-changing collisional efficiency under trap conditions. The results provide the essential input data for modeling different aspects of the experimental setups which can finally produce internally cold molecular ions interacting with a buffer gas.

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.

No abstract

A resolution-of-the-identity (RI) approximation for two-electron integrals over Gaussian basis functions with a complex-scaled exponent is presented. Such functions are used in non-Hermitian quantum mechanics to represent electronic resonances by L2 integrable wave functions with complex energies. We have implemented this new RI approximation for second-order Møller-Plesset perturbation (MP2) theory as well as for the Coulomb and exchange contributions in Hartree-Fock (HF) theory. We discuss the differences to the standard RI approximation of Hermitian quantum mechanics and demonstrate the utility of the non-Hermitian RI-MP2 and RI-HF methods by computations of the orientation-dependent ionization rates of CO, C6H6, and C10H8 in static electric fields. Our results illustrate that RI-MP2 correctly describes correlation effects in molecular electronic resonances while the computational cost is low enough to allow for investigations of medium-sized molecules.

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.

The collisional cooling of the internal rotational states of the nonlinear anion NH2- (1A1), occurring at the low temperature of a cold ion trap under helium buffer gas cooling, is examined via quantum dynamics calculations and ion decay rate measurements. The calculations employ a novel ab initio potential energy surface that describes the interaction anisotropy and range of action between the molecular anions and the neutral He atoms. The state changing integral cross sections are employed to obtain the state-to-state rate coefficients, separately for the ortho- and the para-NH2- ions. These rates are in turn used to compute the state population evolution in the trap for both species, once photodetachment by a laser is initiated in the trap. The present work shows results for the combined losses of both species after the photodetachment laser is switched on and analyzes the differences of loss kinetics between the two hyperfine isomers.