A hierarchical family of three-dimensional potential energy surfaces for He-CO

Zero- and high-pressure mechanisms in the complex forming reactions of OH with methanol and formaldehyde at low temperatures

A recent Ring Polymer Molecular Dynamics study of the reactions of OH with methanol and formaldehyde, at zero pressure and below 100 K, has shown the formation of long lived complexes, with long lifetimes, longer than 100 ns for the lower temperatures studied, 20-100 K (del Mazo-Sevillano et al., 2019). These long lifetimes support the existence of multi collision events with the He buffer-gas atoms under experimental conditions, as suggested by several transition state theory studies of these reactions. In this work we study these secondary collisions, as a dynamical approach to study pressure effects on these reactions. For this purpose, the potential energy surfaces of He with H2CO, OH, H2O and HCO are calculated at highly accurate ab initio level. The stability of some of the complexes is studied using Path Integral Molecular dynamics techniques, determining that OH-H2CO complexes can be formed up to 100 K or higher temperatures, while the weaker He-H2CO complexes dissociate at approximately 50 K. The predicted IR intensity spectra shows new features which could help the identification of the OH-H2CO complex. Finally, the He-H2CO + OH and OH-H2CO + He collisions are studied using quassi-classical trajectories, finding that the cross section to produce HCO + H2O products increases with decreasing collision energy, and that it is ten times higher in the He-H2CO + OH case.

Quasiparticle Approach to Molecules Interacting with Quantum Solvents

Understanding the behavior of molecules interacting with superfluid helium represents a formidable challenge and, in general, requires approaches relying on large-scale numerical simulations. Here, we demonstrate that experimental data collected over the last 20 years provide evidence that molecules immersed in superfluid helium form recently predicted angulon quasiparticles [Phys. Rev. Lett. 114, 203001 (2015)PRLTAO0031-900710.1103/PhysRevLett.114.203001]. Most important, casting the many-body problem in terms of angulons amounts to a drastic simplification and yields effective molecular moments of inertia as straightforward analytic solutions of a simple microscopic Hamiltonian. The outcome of the angulon theory is in good agreement with experiment for a broad range of molecular impurities, from heavy to medium-mass to light species. These results pave the way to understanding molecular rotation in liquid and crystalline phases in terms of the angulon quasiparticle.

van der Waals potential energy surfaces from the exchange-hole dipole moment dispersion model

The potential energy surfaces (PESs) of 28 simple van der Waals complexes, each consisting of a rare-gas (Rg) atom interacting with a linear molecule, are calculated using the exchange-hole dipole moment (XDM) dispersion model in conjunction with three base density functionals (HFPBE, PW86PBE, and a commensurate hybrid functional). Results are compared with literature coupled-cluster reference data. The quality of the computed PESs is assessed based on the positions of the global minima and the corresponding binding energies. Only the hybrid functional is found to provide generally reliable PESs. Dispersion-corrected HFPBE strongly underestimates the equilibrium intermolecular separations and predicts different global minima than the reference PESs for Rg–HCl, Rg–HBr, and two of the Rg–HCN complexes. Analysis of the binding-energy errors reveals that the performance of HFPBE degrades as the size of the Rg atoms increase down the group, while the performance of PW86PBE is significantly worse for strongly-polar molecules. PW86PBE, and to a lesser extent the hybrid, strongly overbind Kr–HF due to charge-transfer error. Despite this, the XDM-corrected hybrid functional displays the best overall error statistics and provides binding energies to within ca. 10 cm–1 of the coupled-cluster reference data at a greatly reduced computational cost.

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.

Inelastic Scattering of CO with He: Polarization Dependent Differential State-to-State Cross Sections

A joint theoretical and experimental study of state-to-state rotationally inelastic polarization dependent differential cross sections (PDDCSs) for CO (v = 0, j = 0, 1, 2) molecules colliding with helium is reported for collision energies of 513 and 840 cm(-1). In a crossed molecular beam experiment, velocity map imaging (VMI) with state-selective detection by (2 + 1) and (1 + 1') resonance enhanced multiphoton ionization (REMPI) is used to probe rotational excitation of CO due to scattering. By taking account of the known fractions of the j = 0, 1, and 2 states of CO in the rotationally cold molecular beam (Trot ≈ 3 K), close-coupling theory based on high-quality ab initio potential energy surfaces for the CO-He interaction is used to simulate the differential cross sections for the mixed initial states. With polarization-sensitive 1 + 1' REMPI detection and a direct analysis procedure described by Suits et al. ( J. Phys, Chem. A 2015 , 119 , 5925 ), alignment moments are extracted from the images and the latter are compared with images simulated by theory using the calculated DCS and alignment moments. In general, good agreement of theory with the experimental results is found, indicating the reliability of the experiment in reproducing state-to-state differential and polarization-dependent differential cross sections.

Quantum dynamical resonances in low-energy CO(j = 0) + He inelastic collisions

In molecular collisions, long-lived complexes may be formed that correspond to quasi-bound states in the van der Waals potential and give rise to peaks in the collision energy-dependent cross-sections. They are known as 'resonances' and their experimental detection remains difficult because their signatures are extremely challenging to resolve. Here, we show a complete characterization of quantum-dynamical resonances occurring in CO-He inelastic collisions with rotational CO(j = 0->1) excitation. Crossed-beam scattering experiments were performed at collision energies as low as 4 cm(-1), equivalent to a temperature of 4 K. Resonance structures in the measured cross-sections were identified by comparison with quantum-mechanical calculations. The excellent agreement found confirms that the potential energy surfaces describing the CO-He van der Waals interaction are perfectly suitable for calculating state-to-state (de)excitation rate coefficients at the very low temperatures needed in chemical modelling of the interstellar medium. We also computed these rate coefficients.

Formation of coherent rotational wavepackets in small molecule-helium clusters using impulsive alignment

We show that rotational line spectra of molecular clusters with near zero permanent dipole moments can be observed using impulsive alignment. Aligned rotational wavepackets were generated by non-resonant interaction with intense femtosecond laser pump pulses and then probed using Coulomb explosion by a second, time-delayed femtosecond laser pulse. By means of a Fourier transform a rich spectrum of rotational eigenstates was derived. For the smallest cluster, C(2)H(2)-He, we were able to establish essentially all rotational eigenstates up to the dissociation threshold on the basis of theoretical level predictions. The C(2)H(2)-He complex is found to exhibit distinct features of large amplitude motion and very early onset of free internal rotor energy level structure.

Low-temperature rate coefficients for vibrational relaxation of (3)Σ(u)(+)Rb2 molecules by (3)He and (4)He atoms

We present quantum-scattering calculations of (4)He and (3)He colliding with (87)Rb2. For both helium isotopes, the elastic and inelastic rate coefficients are strongly influenced by the J = 1 partial wave. For the lighter isotope, a strong resonance feature of the J = 1 partial wave is responsible for an extremely efficient vibrational relaxation process. We also perform bound-state calculations of the Rb2He complex for even Rb permutation symmetry and nonzero total angular momentum. The global Rb2He (3)Σu(+) potential-energy surface used supports four bound states for (4)He and a single one for (3)He. We propose an analysis of the (87)Rb2(4)He spectrum separating the contributions of Rb2 rotation and helium motion.

Fixed node diffusion Monte Carlo using a genetic algorithm: a study of the CO-(4)He(N) complex, N = 1…10

The diffusion Monte Carlo (DMC) method is a widely used algorithm for computing both ground and excited states of many-particle systems; for states without nodes the algorithm is numerically exact. In the presence of nodes approximations must be introduced, for example, the fixed-node approximation. Recently we have developed a genetic algorithm (GA) based approach which allows the computation of nodal surfaces on-the-fly [Ramilowski and Farrelly, Phys. Chem. Chem. Phys., 2010, 12, 12450]. Here GA-DMC is applied to the computation of rovibrational states of CO-(4)He(N) complexes with N≤ 10. These complexes have been the subject of recent high resolution microwave and millimeter-wave studies which traced the onset of microscopic superfluidity in a doped (4)He droplet, one atom at a time, up to N = 10 [Surin et al., Phys. Rev. Lett., 2008, 101, 233401; Raston et al., Phys. Chem. Chem. Phys., 2010, 12, 8260]. The frequencies of the a-type (microwave) series, which correlate with end-over-end rotation in the CO-(4)He dimer, decrease from N = 1 to 3 and then smoothly increase. This signifies the transition from a molecular complex to a quantum solvated system. The frequencies of the b-type (millimeter-wave) series, which evolves from free rotation of the rigid CO molecule, initially increase from N = 0 to N∼ 6 before starting to decrease with increasing N. An interesting feature of the b-type series, originally observed in the high resolution infra-red (IR) experiments of Tang and McKellar [J. Chem. Phys., 2003, 119, 754] is that, for N = 7, two lines are observed. The GA-DMC algorithm is found to be in good agreement with experimental results and possibly detects the small (∼0.7 cm(-1)) splitting in the b-series line at N = 7. Advantages and disadvantages of GA-DMC are discussed.

A new ab initio potential energy surface for the NeCO complex with the vibrational coordinate dependence

A new high quality three-dimensional potential energy surface for the Ne-CO van der Waals complex is developed using the CCSD(T) method and avqz∕avqz+33221 basis set. The ab initio calculation is performed in a total of 1365 configurations with supermolecule method. There is a single global minimum located in a nearly T-shaped geometry. The global minimum energy is -49.4090 cm(-1) at R(e)=6.40a(0) and θ(e)=82.5(∘) for V(00). Using the three-dimensional potential energy surface, we have calculated bound rovibrational energy levels up to J = 10 including the Coriolis coupling terms. Compared with the experimental transition frequencies, the theoretical results are in good agreement with the experimental results.

Theoretical and experimental study of the rovibrational spectrum of He(2)-CO

We report calculated microwave and infrared ro-vibrational transitions of the van der Waals complex He(2)-CO. The calculations were done using a product basis and a Lanczos eigensolver, together with He-CO and He-He potential functions taken from the literature. The results are found to be in good agreement with previously reported experimental results, and they enable the experimental assignments to be clarified, augmented, and (in one case) corrected. Unlike some other van der Waals complexes with two He atoms such as He(2)-N(2)O, it is not possible to associate a set of energy levels with the "torsional" motion of the two He atoms on a ring encircling the dopant (in this case CO). Although we assume that the dipole moment is along the CO axis we find nonetheless that many transitions have appreciable intensity due to ro-vibrational coupling.

An exchange-Coulomb model potential energy surface for the Ne-CO interaction. I. Calculation of Ne-CO van der Waals spectra

Exchange-Coulomb model potential energy surfaces have been developed for the Ne-CO interaction. The initial model is a three-dimensional potential energy surface based upon computed Heitler-London interaction energies and literature results for the long-range induction and dispersion energies, all as functions of interspecies distance, the orientation of CO relative to the interspecies axis, and the bond length of the CO molecule. Both a rigid-rotor model potential energy surface, obtained by setting the CO bond length equal to its experimental spectroscopic equilibrium value, and a vibrationally averaged model potential energy surface, obtained by averaging the stretching dependence over the ground vibrational motion of the CO molecule, have been constructed from the full data set. Adjustable parameters in each model potential energy surface have been determined through fitting a selected subset of pure rotational transition frequencies calculated for the (20)Ne-(12)C(12)O isotopolog to precisely known experimental values. Both potential energy surfaces provide calculated results for a wide range of available experimental microwave, millimeter-wave, and midinfrared Ne-CO transition frequencies that are generally far superior to those obtained using the best current literature potential energy surfaces. The vibrationally averaged CO ground state potential energy surface, employed together with a potential energy surface obtained from it by replacing the ground vibrational state average of the CO stretching dependence of the potential energy surface by an average over the first excited CO vibrational state, has been found to be particularly useful for computing and/or interpreting mid-IR transition frequencies in the Ne-CO dimer.

Rovibrational structure of the Xe-CO complex based on a new three-dimensional ab initio potential

The first three-dimensional interaction potential energy surface of the Xe-CO complex is developed using the single and double excitation coupled cluster theory with noniterative treatment of triple excitations. Mixed basis sets, aug-cc-pVQZ for the C and O atoms and aug-cc-pVQZ-PP for the Xe atom, including an additional (3s3p2d2f1g) set of midbond functions are used. The calculated single point energies at five fixed r(co) values are fitted to an analytic two-dimensional potential model, and further the five model potentials are used to construct the three-dimensional potential energy surface by interpolating along (r-r(e)). Dynamical calculations with the vibrationally averaged potentials are performed to determine the energy levels and the frequencies of various rovibrational transitions. Our results agree well with the experiment. For example, the IR transitions of 508 lines are precisely reproduced with only a total rms error of 0.105 cm(-1).

Computational Spectroscopy of Carbon Monoxide Isotopomers in Helium Clusters

The rotational excitation spectrum, including the vibrational shift of the rotational band, of several CO isotopomers solvated in He clusters has been calculated. Reptation quantum Monte Carlo simulations are used in conjunction with an accurate He-CO potential energy surface, which quantitatively describes the rovibrational spectrum of the binary complex. Our simulations, when compared with number-selective infrared spectra taken for different isotopomers, help discriminate among the alternative assignments proposed for cluster sizes around 15 He atoms. The origin of the vibrational band has a red shift that is nearly linear with the cluster size within the first solvation shell and is almost constant up to the largest cluster studied, well beyond completion of the second solvation shell. A blue upturn at even larger sizes would be needed to attain the nanodroplet limit, as recently estimated from the isotopic dependence of the measured R(0) transitions.

Aqueous Microsolvation of Mercury Halide Species†

The effects of aqueous solvation on the thermochemistry of reactions between mercury and small halogen molecules has been investigated by the microsolvation approach using ab initio and density functional theory (DFT) calculations. The structures, vibrational frequencies, and binding energies of 1, 2, and 3 water molecules with mercury-halide (HgBr2, HgBrCl, HgCl2, HgBr, and HgCl) and related mercury and halogen species (Br2, BrCl, Cl2, Cl, Hg, and Br) have been computed with second order Møller-Plesset perturbation theory (MP2) and the B3LYP density functional method. Accurate incremental water binding energies have been obtained at the complete basis set (CBS) limit using sequences of correlation consistent basis sets, including higher order correlation effects estimated from coupled cluster calculations. The resulting energetics were used to calculate the influence of water molecules on the thermochemistry of a number of reactions between mercury and small halogen-containing molecules. In general, the presence of water favors the formation of oxidized mercury halide species.

Helium clusters seeded with CO molecules: New results for HeN–C1813O and the approach to the nanodroplet limit

Infrared spectra of helium clusters seeded with doubly substituted carbon monoxide molecules, 13C18O, have been studied in order to complement recent helium nanocluster results and to determine whether additional isotopic data would help to separate vibrational and rotational contributions to the observed transitions. The experiments were made by direct infrared absorption in pulsed supersonic jet expansions using a tunable diode laser probe in the region of the fundamental band (approximately 2045 cm-1 for 13C18O). Even with data on the R0 transitions from four CO isotopomers, it was found that a clear and consistent separation of vibration and rotation could not be achieved for HeN-CO clusters in the size range N approximately 10-20. Isotope shifts observed for clusters with 13C18O (relative to 12C16O) were found to be close to the sums of the shifts previously determined for 13C16O and 12C18O. The new measurements generally supported previous assignments of cluster size, but some modifications for the range N=14-16 are suggested here. New measurements for HeN-12C16O under conditions favoring larger clusters (high backing pressure and low jet temperature) showed that individual transitions could be resolved even at N approximately 50. For larger clusters, a partly resolved "lump" of transitions was observed to approach the nanodroplet limit.