Accurate intermolecular ground-state potential-energy surfaces of the HCCH–He, Ne, and Ar van der Waals complexes

Calculation of the Transport and Relaxation Properties of the Ar···HCl van der Waals Complex Using a New Potential Energy Surface: Comparison of the Classical and Full Quantum Mechanical Kinetic Theory Results with Molecular Dynamics Simulations

The intermolecular potential energy surface (PES) of the Ar···HCl complex was calculated at the RCCSD(T)/aug-cc-pvQz-BF level of theory. The obtained potential was expanded in terms of Legendre polynomials and fitted to a mathematical model. The fitting results are highly correlated with the ab initio PES data with SD = 5.9 × 10^-3 cm^-1 and average absolute deviation (AAD) = 4.0 × 10^-6 cm^-1. The interaction second virial coefficients (B₁₂) in the temperature range of 190-480 K were calculated by considering classical and first quantum corrections and compared with the available experimental data. A reasonable agreement with the experimental and calculated B₁₂ was obtained. The PES was also used to obtain the rovibrational energy levels, and the spectroscopic rovibrational constants were obtained. It was found that the D₀ values differ ∼2.25 cm^-1 from the experimental values of the ground rovibrational state. Furthermore, the obtained potential was used to calculate the transport and relaxation properties using full quantum close-coupling (CC) formalism and the classical kinetic theory methods based on the Mason-Monchik approximation (MMA). It was found that the deviation between MMA and CC calculations is increased with increasing the temperature due to the higher influence of the rotational degrees of freedom on the transport properties. Also, the contribution of the inelastic (off-diagonal) transitions for diffusion coefficient is higher than the viscosity. Furthermore, the classical molecular dynamics simulations were performed using LJ(12,6) and Vashishta models, to calculate the interaction diffusion and viscosity coefficients, and compared with the results of the full quantum CC calculations. The obtained results confirm that the Vashishta model is better fitted to the ab initio potentials and is more accurate than LJ(12,6) in calculation of the diffusion coefficients.

Large amplitude motion within acetylene-rare gas complexes hosted in helium droplets

Near-infrared spectroscopy of the C2H2-Ar, Kr complexes was performed in the spectral region overlapping the ν3/ν2 + ν4 + ν5 Fermi-type resonance of C2H2. The experiment was conducted along the HElium NanoDroplet Isolation (HENDI) technique in order to study the coupling dynamics between a floppy molecular system (C2H2-Ar and C2H2-Kr) and a mesoscopic quantum liquid (the droplet). Calculations were performed using a spectral element based close-coupling program and state-of-the-art 2-dimensional potential energy surfaces to determine the bound states of the C2H2-Ar and C2H2-Kr complexes and simulate the observed spectra. This furnished a quantitative basis to unravel how the superfluid and non-superfluid components of the droplet affect the rotation and the deformation dynamics of the hosted complex.

Ab initio study of the CO–N2 complex: a new highly accurate intermolecular potential energy surface and rovibrational spectrum

We present a highly accurate ab initio intermolecular potential-energy surface and rovibrational spectrum for the CO–N₂ complex.

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.

Large amplitude motion of the acetylene molecule within acetylene-neon complexes hosted in helium droplets

Superfluid helium droplets provide an ideal environment for spectroscopic studies with rotational resolution. Nevertheless, the molecular rotation is hindered because the embedded molecules are surrounded by a non-superfluid component. The present work explores the dynamical role of this component in the hindered rotation of C2H2 within the C2H2-Ne complex. A HENDI experiment was built and near-infrared spectroscopy of C2H2-Ne and C2H2 was performed in the spectral region overlapping the ν3/ν2 + ν4 + ν5 Fermi-type resonance of C2H2. The comparison between measured and simulated spectra helped to address the above issue.

The infrared spectrum of the Ne-C2D2 complex

Infrared spectra of Ne-C2D2 are observed in the region of the ν3 fundamental band (asymmetric C-D stretch, ≈2440 cm(-1)) using a tunable optical parametric oscillator to probe a pulsed supersonic slit jet expansion from a cooled nozzle. Like helium-acetylene, this system lies close to the free rotor limit, making analysis tricky because stronger transitions tend to pile up close to monomer (C2D2) rotation-vibration transitions. Assignments are aided by predicted rotational energies calculated from a published ab initio intermolecular potential energy surface. The analysis extends up to the j = 3←2 band, where j labels C2D2 rotation within the dimer, and is much more complete than the limited infrared assignments previously reported for Ne-C2H2 and Ne-C2HD. Two previous microwave transitions within the j = 1 state of Ne-C2D2 are reassigned. Coriolis model fits to the theoretical levels and to the spectrum are compared. Since the variations observed as a function of C2D2 vibrational excitation are comparable to those noted between theory and experiment, it is evident that more detailed testing of theory will require vibrational averaging over the acetylene intramolecular modes.

The infrared spectrum of the He-C2D2 complex

Spectra of the helium-acetylene complex are elusive because this weakly bound system lies close to the free rotor limit. Previously, limited assignments of He-C2D2 transitions in the R(0) region of the ν3 fundamental band (≈2440 cm(-1)) were published. Here, new He-C2D2 infrared spectra of this band are obtained using a tunable optical parametric oscillator laser source to probe a pulsed supersonic slit jet expansion from a cooled nozzle, and the analysis is extended to the weaker and more difficult P(1) and R(1) regions. A term value approach is used to obtain a consistent set of "experimental" energy levels. These are compared directly with calculations using two recently reported ab initio intermolecular potential energy surfaces, which exhibit small but significant differences. Rovibrational energies for the He-C2H2 complex are also calculated using both surfaces. A Coriolis model, useful for predicting spectral intensities, is used to interpret the energy level patterns, and a comparison with the isoelectronic complex He-CO is made.

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.

Theoretical prediction of the linear isomers for rare gas-carbon disulfide complexes: He-CS₂, Ne-CS₂, and Ar-CS₂

Theoretical studies of the potential energy surfaces (PESs) and bound states are performed for rare gas-carbon disulfide complexes, He-CS2, Ne-CS2, and Ar-CS2. Three two-dimensional intermolecular PESs are constructed from ab initio data points which are calculated at the CCSD(T) level with aug-cc-pVTZ basis set supplemented with bond functions. We find that the three PESs have very similar features and each PES can be characterized by a global T-shaped minimum, two equivalent local linear minima, and the saddle points between them. The T-shaped isomer is energetically more stable than the linear isomer for each complex. The linear isomers, which have not been observed in experiment so far, are predicted from our PESs and further identified by bound state calculations. Moreover, we assign several intermolecular vibrational states for both the T-shaped and linear isomers of the three complexes via the analysis of wavefunctions. The corresponding vibrational frequencies are calculated from the bound state energies for these assigned states. These frequencies could be helpful for further experimental studies, especially for the linear isomers. We also calculate the rovibrational transition frequencies for the three T-shaped isomers and the pure rotational transition frequencies for the linear isomers, respectively. The accuracy of the PESs is validated by the good agreement between theoretical and experimental results for the rovibrational transition frequencies and spectroscopic parameters.

Intermolecular interaction potentials for the Ar–C2H2, Kr–C2H2, and Xe–C2H2 weakly bound complexes: Information from molecular beam scattering, pressure broadening coefficients, and rovibrational spectroscopy

Integral cross sections and pressure broadening coefficients have been measured for the acetylene-krypton complex, by molecular beam scattering and by high resolution IR spectroscopy, respectively. A new potential energy surface (PES) is proposed to describe structure and dynamical properties of this prototypical weakly bound complex. The PES has been parametrized exploiting a novel atom-bond pairwise additive scheme and has been fitted to the experimental data. A similar PES has been obtained for the acetylene-xenon system by a proper scaling of the interaction parameters of the krypton case, based on empirical considerations. These PESs together with that recently proposed by the same authors [J. Phys. Chem. 109, 8471 (2005)] for the acetylene-argon case have been employed for close coupling calculations of the pressure broadening cross sections and for a characterization of the rovibrational structure of the complexes.

HXeCCH in Ar and Kr matrices

HXeCCH molecule is prepared in Ar and Kr matrices and characterized by IR absorption spectroscopy. The experiments show that HXeCCH can be made in another host than the polarizable Xe environment. The H-Xe stretching absorption of HXeCCH in Ar and Kr is blueshifted from the value measured in solid Xe. The maximum blueshifts are +44.9 and +32.3 cm(-1) in Ar and Kr, respectively, indicating stabilization of the H-Xe bond. HXeCCH has a doublet H-Xe stretching absorption measured in Xe, Kr, and Ar matrices with a splitting of 5.7, 13, and 14 cm(-1), respectively. Ab initio calculations for the 1:1 HXeCCHcdots, three dots, centeredNg complexes (Ng = Ar, Kr, or Xe) are used to analyze the interaction of the hosts with the embedded molecule. These calculations support the matrix-site model where the band splitting observed experimentally is caused by specific interactions of the HXeCCH molecule with noble-gas atoms in certain local morphologies. However, the 1:1 complexation is unable to explain the observed blueshifts of the H-Xe stretching band in Ar and Kr matrices compared to a Xe matrix. More sophisticated computational approach is needed to account in detail the effects of solid environment.