A comprehensive experimental and theoretical study of H2−CO spectra

Theory cracks old data: Rovibrational energy levels of orthoH2-CO derived from experiment

Measurements of rovibrational spectra of clusters provide physical insight only if spectral lines can be assigned to pairs of quantum states, and further insight is obtained if one can deduce the quantitative energy-level pattern. Both steps can be so difficult that some measured spectra remain unassigned, one example is orthoH2-CO. To extend the scope of spectroscopic insights, we propose to use theoretical information in interpretation of spectra. We first performed high accuracy, full-dimensional calculations of the orthoH2-CO spectrum, at the highest practically achievable levels of electronic structure theory and quantum nuclear dynamics. Then, an iterative, theory-guided method developed here allowed us to fully interpret the spectrum of orthoH2-CO, extending the range of van der Waals clusters for which spectroscopy can provide physical insights.

Dispersion Energy from the Time-Independent Coupled-Cluster Polarization Propagator

We present a new method of calculation of the dispersion energy in the second-order symmetry-adapted perturbation theory. Using the Longuet-Higgins integral and time-independent coupled-cluster response theory, one shows that the general expression for the dispersion energy can be written in terms of cluster amplitudes and the excitation operators σ, which can be obtained by solving a linear equation. We introduced an approximate scheme dubbed CCPP2(T) for the dispersion energy accurate to the second order of intramonomer correlation, which includes certain classes to be summed to infinity. Assessment of the accuracy of the CCPP2(T) dispersion energy against the FCI dispersion for He₂ demonstrates its high accuracy. For more complex systems, CCPP2(T) matches the accuracy of the best methods introduced for calculations of dispersion so far. The method can be extended to higher-order levels of excitations, providing a systematically improvable theory of dispersion interaction.

Stereodynamic control of cold rotationally inelastic CO + HD collisions

Quantum control of molecular collision dynamics is an exciting emerging area of cold collisions. Co-expansion of collision partners in a supersonic molecular beam combined with precise control of their quantum states and alignment/orientation using Stark-induced Adiabatic Raman Passage allows exquisite stereodynamic control of the collision outcome. This approach has recently been demonstrated for rotational quenching of HD in collisions with H₂, D₂, and He and D₂ by He. Here we illustrate this approach for HD(v = 0, j = 2) + CO(v = 0, j = 0) → HD(v' = 0, j') + CO(v' = 0, j') collisions through full-dimensional quantum scattering calculations at collision energies near 1 K. It is shown that the collision dynamics at energies between 0.01-1 K are controlled by an interplay of L = 1 and L = 2 partial wave resonances depending on the final rotational levels of the two molecules. Polarized cross sections resolved into magnetic sub-levels of the initial and final rotational quantum numbers of the two molecules also reveal a significant stereodynamic effect in the cold energy regime. Overall, the stereodynamic effect is controlled by both geometric and dynamical factors, with parity conservation playing an important role in modulating these contributions depending on the particular final state.

Kinetic Isotope Effect in Low-Energy Collisions between Hydrogen Isotopologues and Metastable Helium Atoms: Theoretical Calculations Including the Vibrational Excitation of the Molecule

We present very accurate theoretical results of Penning ionization rate coefficients of the excited metastable helium atoms (⁴He(2³S) and ³He(2³S)) colliding with the hydrogen isotopologues (H₂, HD, D₂) in the ground and first excited rotational and vibrational states at subkelvin regime. The calculations are performed using the current best ab initio interaction energy surface, which takes into account the nonrigidity effects of the molecule. The results confirm a recently observed substantial quantum kinetic isotope effect (Nat. Chem. 2014, 6, 332-335) and reveal that the change of the rotational or vibrational state of the molecule can strongly enhance or suppress the reaction. Moreover, we demonstrate the mechanism of the appearance and disappearance of resonances in Penning ionization. The additional model computations, with the morphed interaction energy surface and mass, give better insight into the behavior of the resonances and thereby the reaction dynamics under study. Our theoretical findings are compared with all available measurements, and comprehensive data for prospective experiments are provided.

Probing the location of the unpaired electron in spin-orbit changing collisions of NO with Ar

Understanding the molecular forces that drive a reaction or scattering process lies at the heart of molecular dynamics. Here, we present a combined experimental and theoretical study of the spin-orbit changing scattering dynamics of oriented NO molecules with Ar atoms. Using our crossed molecular beam apparatus, we have recorded velocity-map ion images and extracted differential and integral cross sections of the scattering process in the side-on geometry. We observe an overall preference for collisions close to the N atom in the spin-orbit changing manifold, which is a direct consequence of the location of the unpaired electron on the potential energy surface. In addition, a prominent forward scattered feature is observed for intermediate, even rotational transitions when the atom approaches the molecule from the O-end. The appearance of this peak originates from an attractive well on the A' potential energy surface, which efficiently directs high impact parameter trajectories towards the region of high unpaired electron density near the N-end of the molecule. The ability to orient molecules prior to collision, both experimentally and theoretically, allows us to sample different regions of the potential energy surface(s) and unveil the associated collision pathways.

Explicitly correlated ab initio potential energy surface and predicted rovibrational spectra for H2O-N2 and D2O-N2 complexes

An ab initio intermolecular potential energy surface (PES) for the van der Waals complex of H₂O-N₂ that explicitly incorporates the intramolecular Q₂ bending normal mode of the H₂O monomer is presented. The electronic structure computations have been carried out at the explicitly correlated coupled cluster theory [CCSD(T)-F12] with an augmented correlation-consistent triple zeta basis set and an additional bond function. Analytic five-dimensional intermolecular PESs for ν₂(H₂O) = 0 and 1 are obtained by fitting to the multi-dimensional Morse/long-range potential function form. These fits to 40 890 points have the root-mean-square (rms) discrepancy of 0.88 cm^-1 for interaction energies less than 2000.0 cm^-1. The resulting vibrationally averaged PESs provide good representations of the experimental microwave and infrared data: for microwave transitions of H₂O-N₂, the rms discrepancy is only 0.0003 cm^-1, and for infrared transitions of the A₁ symmetry of the H₂O(ν₂ = 1 ← 0)-N₂, the rms discrepancy is 0.001 cm^-1. The calculated infrared band origin shifts associated with the ν₂ bending vibration of water are 2.210 cm^-1 and 1.323 cm^-1 for H₂O-N₂ and D₂O-N₂, respectively, in good agreement with the experimental values of 2.254 cm^-1 and 1.266 cm^-1. The benchmark tests and comparisons of the predicted spectral properties are carried out between CCSD(T)-F12a and CCSD(T)-F12b approaches.

SAPT codes for calculations of intermolecular interaction energies

Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies. The set of SAPT codes that is described here, the current version named SAPT2020, includes virtually all variants of SAPT developed so far, among them two-body SAPT based on perturbative, coupled cluster, and density functional theory descriptions of monomers, three-body SAPT, and two-body SAPT for some classes of open-shell monomers. The properties of systems governed by noncovalent interactions can be predicted only if potential energy surfaces (force fields) are available. SAPT is the preferred approach for generating such surfaces since it is seamlessly connected to the asymptotic expansion of interaction energy. SAPT2020 includes codes for automatic development of such surfaces, enabling generation of complete dimer surfaces with a rigid monomer approximation for dimers containing about one hundred atoms. These codes can also be used to obtain surfaces including internal degrees of freedom of monomers.

Evidence of Nonrigidity Effects in the Description of Low-Energy Anisotropic Molecular Collisions of Hydrogen Molecules with Excited Metastable Helium Atoms

Cold collisions serve as a sensitive probe of the interaction potential. In the recent study of Klein et al. (Nature Phys.2017, 13, 35–38), the one-parameter scaling of the interaction potential was necessary to obtain agreement between theoretical and observed patterns of the orbiting resonances for excited metastable helium atoms colliding with hydrogen molecules. Here, we show that the effect of nonrigidity of the H₂ molecule on the resonant structure, absent in the previous study, is critical to predict the correct positions of the resonances in that case. We have complemented the theoretical description of the interaction potential and revised reaction rate coefficients by proper inclusion of the flexibility of the molecule. The calculated reaction rate coefficients are in remarkable agreement with the experimental data without empirical adjustment of the interaction potential. We have shown that even state-of-the-art calculations of the interaction energy cannot ensure agreement with the experiment if such an important physical effect as flexibility of the interacting molecule is neglected. Our findings about the significance of the nonrigidity effects can be especially crucial in cold chemistry, where the quantum nature of molecules is pronounced.

Automatic Generation of Flexible-Monomer Intermolecular Potential Energy Surfaces

A method is developed for automatic generation of nonreactive intermolecular two-body potential energy surfaces (PESs) including intramonomer degrees of freedom. This method, called flex-autoPES, is an extension of the autoPES method developed earlier, which assumes rigid monomers. In both cases, the whole PES development proceeds without any human intervention. The functional form used is a sum of products of site-site functions (both atomic and off-atomic sites can be used). The leading terms with sites involving different monomers are of physically motivated form. The long-range part of a PES is computed from monomer properties without using any dimer information. The close-range part is fitted to dimer interaction energies computed using electronic structure methods. Virtually any method can be used in such calculations, but the use of symmetry-adapted perturbation theory provides a seamless connection to the long-range part of the PES. The performance of the flex-autoPES code was tested by developing a full-dimensional PES for the water dimer and PESs including only some soft intramonomer degrees of freedom for the ethylene glycol dimer and for the ethylene glycol-water dimer. In the case of the water dimer, the root-mean-square error (RMSE) of the PES from the data points with negative total energies is 0.03 kcal/mol, and we expect this PES to be more accurate than any previously published PES of this type. For the ethylene glycol dimer and the ethylene glycol-water dimers, the analogous RMSEs are 0.25 and 0.1 kcal/mol, respectively.

Full quantum calculation of the rovibrational states and intensities for a symmetric top-linear molecule dimer: Hamiltonian, basis set, and matrix elements

The rovibrational energy levels and intensities of the CH₃F-H₂ dimer have been obtained using our recent global intermolecular potential energy surface [X.-L. Zhang et al., J. Chem. Phys. 148, 124302 (2018)]. The Hamiltonian, basis set, and matrix elements are derived and given for a symmetric top-linear molecule complex. This approach to the generation of energy levels and wavefunctions can readily be utilized for studying the rovibrational spectra of other van der Waals complexes composed of a symmetric top molecule and a linear molecule, and may readily be extended to other complexes of nonlinear molecules and linear molecules. To confirm our method, the rovibrational levels of the H₂O-H₂ dimer have been computed and shown to be in good agreement with experiment and with previous theoretical results. The rovibrational Schrödinger equation has been solved using a Lanczos algorithm together with an uncoupled product basis set. As expected, dimers containing ortho-H₂ are more strongly bound than dimers containing para-H₂. Energies and wavefunctions of the discrete rovibrational levels of CH₃F-paraH₂ complexes obtained from the direct vibrationally averaged 5-dimensional potentials are in good agreement with the results of the reduced 3-dimensional adiabatic-hindered-rotor (AHR) approximation. Accurate calculations of the transition line strengths for the orthoCH₃F-paraH₂ complex are also carried out, and are consistent with results obtained using the AHR approximation. The microwave spectrum associated with the orthoCH₃F-orthoH₂ dimer has been predicted for the first time.

Double resonance rotational spectroscopy of He-HCO. Phys Chem Chem Phys 2019;21:3440-3445. [PMID: 30191208 DOI: 10.1039/c8cp04532a]

The ground state of He-HCO+ is investigated using a recently developed double resonance technique, consisting of a rotational transition followed by a vibrational transition into a dissociative state. In order to derive precise predictions for the rotational states, the high resolution infrared predissociation spectroscopy of the v1 C-H stretching mode is revisited. Eleven pure rotational transitions are measured via the double resonance method. A least squares fit of these transitions to a standard linear rotor Hamiltonian reveals that the semirigid rotor model cannot fully describe the loosely bound He-HCO+ complex. The novel double resonance technique is compared with other action spectroscopic schemes, and some potential future applications are presented.

Molecular dimers of methane clathrates: ab initio potential energy surfaces and variational vibrational states

Motivated by the energetic and environmental relevance of methane clathrates, highly accurate ab initio potential energy surfaces (PESs) have been developed for the three possible dimers of the methane and water molecules: (H₂O)₂, CH₄·H₂O, and (CH₄)₂.

A general variational approach for computing rovibrational resonances of polyatomic molecules. Application to the weakly bound H2He+ and H2⋅CO systems

The quasi-variational quantum chemical protocol and code GENIUSH [E. Mátyus et al., J. Chem. Phys. 130, 134112 (2009) and C. Fábri et al., J. Chem. Phys. 134, 074105 (2011)] has been augmented with the complex absorbing potential (CAP) technique, yielding a method for the determination of rovibrational resonance states. Due to the effective implementation of the CAP technique within GENIUSH, the GENIUSH-CAP code is a powerful tool for the study of important dynamical features of arbitrary-sized molecular systems with arbitrary composition above their first dissociation limit. The GENIUSH-CAP code has been tested and validated on the H₂He⁺ cation: the computed resonance energies and lifetimes are compared to those obtained with a previously developed triatomic rovibrational resonance-computing code, D²FOPI-CCS [T. Szidarovszky and A. G. Császár Mol. Phys. 111, 2131 (2013)], utilizing the complex coordinate scaling method. A unique feature of the GENIUSH-CAP protocol is that it allows the simple implementation of reduced-dimensional dynamical models. To prove this, resonance energies and lifetimes of the H₂⋅CO van der Waals complex have been computed utilizing a four-dimensional model (freezing the two monomer stretches), and a related potential energy surface, of the complex.

Scattering resonances in bimolecular collisions between NO radicals and H2 challenge the theoretical gold standard

Over the last 25 years, the formalism known as coupled-cluster (CC) theory has emerged as the method of choice for the ab initio calculation of intermolecular interaction potentials. The implementation known as CCSD(T) is often referred to as the gold standard in quantum chemistry. It gives excellent agreement with experimental observations for a variety of energy-transfer processes in molecular collisions, and it is used to calibrate density functional theory. Here, we present measurements of low-energy collisions between NO radicals and H₂ molecules with a resolution that challenges the most sophisticated quantum chemistry calculations at the CCSD(T) level. Using hitherto-unexplored anti-seeding techniques to reduce the collision energy in a crossed-beam inelastic-scattering experiment, a resonance structure near 14 cm^-1 is clearly resolved in the state-to-state integral cross-section, and a unique resonance fingerprint is observed in the corresponding differential cross-section. This resonance structure discriminates between two NO-H₂ potentials calculated at the CCSD(T) level and pushes the required accuracy beyond the gold standard.

All-dimensional H2-CO potential: Validation with fully quantum second virial coefficients

We use a new high-accuracy all-dimensional potential to compute the cross second virial coefficient B₁₂(T) between molecular hydrogen and carbon monoxide. The path-integral method is used to fully account for quantum effects. Values are calculated from 10 K to 2000 K and the uncertainty of the potential is propagated into uncertainties of B₁₂. Our calculated B₁₂(T) are in excellent agreement with most of the limited experimental data available, but cover a much wider range of temperatures and have lower uncertainties. Similar to recently reported findings from scattering calculations, we find that the reduced-dimensionality potential obtained by averaging over the rovibrational motion of the monomers gives results that are a good approximation to those obtained when flexibility is fully taken into account. Also, the four-dimensional approximation with monomers taken at their vibrationally averaged bond lengths works well. This finding is important, since full-dimensional potentials are difficult to develop even for triatomic monomers and are not currently possible to obtain for larger molecules. Likewise, most types of accurate quantum mechanical calculations, e.g., spectral or scattering, are severely limited in the number of dimensions that can be handled.

On the importance of full-dimensionality in low-energy molecular scattering calculations

Scattering of H₂ on CO is of great importance in astrophysics and also is a benchmark system for comparing theory to experiment. We present here a new 6-dimensional potential energy surface for the ground electronic state of H₂-CO with an estimated uncertainty of about 0.6 cm⁻¹ in the global minimum region, several times smaller than achieved earlier. This potential has been used in nearly exact 6-dimensional quantum scattering calculations to compute state-to-state cross-sections measured in low-energy crossed-beam experiments. Excellent agreement between theory and experiment has been achieved in all cases. We also show that the fully 6-dimensional approach is not needed with the current accuracy of experimental data since an equally good agreement with experiment was obtained using only a 4-dimensional treatment, which validates the rigid-rotor approach widely used in scattering calculations. This finding, which disagrees with some literature statements, is important since for larger systems full-dimensional scattering calculations are currently not possible.

Observation of quantum dynamical resonances in near cold inelastic collisions of astrophysical molecules

Quantum resonances in inelastic collisions, predicted by theory and detected at low energies in a crossed-beam experiment, are reviewed.

This mini review summarizes experimental findings of quantum dynamical resonances in inelastic collisions at energies equivalent to temperatures of a few to a few tens of Kelvin, corresponding to physical conditions prevailing in dense molecular clouds of the interstellar medium. Information obtained is thus relevant to collision energy transfer modelling in such media. Crossed-beam scattering experiments performed at Bordeaux university for inelastic collisions of important astrophysical molecules such as CO with H₂ or He and O₂ with H₂ are described. The peaks that show up in the collision energy dependence of the state-to-state integral cross sections for the lowest rotational excitation transitions reveal the quantum nature of such processes. They are ascribed as shape and Feshbach resonances by comparison with the results of close coupling quantum mechanical calculations performed concomitantly on accurate potential energy surfaces.

Ab initio water pair potential with flexible monomers

A potential energy surface for the water dimer with explicit dependence on monomer coordinates is presented. The surface was fitted to a set of previously published interaction energies computed on a grid of over a quarter million points in the 12-dimensional configurational space using symmetry-adapted perturbation theory and coupled-cluster methods. The present fit removes small errors in published fits, and its accuracy is critically evaluated. The minimum and saddle-point structures of the potential surface were found to be very close to predictions from direct ab initio optimizations. The computed second virial coefficients agreed well with experimental values. At low temperatures, the effects of monomer flexibility in the virial coefficients were found to be much smaller than the quantum effects.