The Ar–C2H2 intermolecular potential from high resolution spectroscopy and ab initio theory: A case for multicenter interactions

Unusual Dynamics and Vibrational Fingerprints of van der Waals Dimers Formed by Linear Molecules and Rare-Gas Atoms

Detailed structural, dynamical, and vibrational analyses have been performed for systems composed of linear triatomic molecules solvated by a single rare-gas atom, He, Ne, or Ar. Among the chromophores of these van der Waals (vdW) dimers, there are four neutral molecules (CO2, CS2, N2O, and OCS) and six molecular cations (HHe2+, HNe2+, HAr2+, HHeNe+, HHeAr+, and HNeAr+), both of apolar and polar nature. Following the exploration of bonding preferences, high-level four-dimensional (4D) potential energy surfaces (PESs) have been developed for 24 vdW dimers, keeping the two intramonomer bond lengths fixed. For these 24 complexes, over 1500 bound vibrational states have been obtained via quasi-variational nuclear-motion computations, employing exact kinetic-energy operators together with the accurate 4D PESs and their 2D/3D cuts. The reduced-dimensional (2D to 4D) dimer models have been compared with full-dimensional (6D) ones in the cases of the neutral CO2·Ar and charged HHe2+·He dimers, corroborating the high accuracy of the 2D to 4D vibrational energies. The reduced-dimensional models suggest that (a) while the equilibrium structures are T-shaped and planar, the effective ground-state structures are nonplanar, (b) certain bound states belong to collinear molecular structures, even when they are not minima, (c) the vdW vibrations are heavily mixed and many states have amplitudes corresponding to both the T-shaped and collinear structures, (d) there are a few dimers, for which even some of the vdW fundamentals lie above the first dissociation limit, and (e) the vdW vibrations are almost fully decoupled from the intramonomer bending motion.

Reaction dynamics within a cluster environment

This perspective article reviews experimental and theoretical works where rare gas clusters and helium nanodroplets are used as a nanoreactor to investigate chemical dynamics in a solvent environment. A historical perspective is presented first followed by specific considerations on the mobility of reactants within these reaction media. The dynamical response of pure clusters and nanodroplets to photoexcitation is shortly reviewed before examining the role of the cluster (or nanodroplet) degrees of freedom in the photodynamics of the guest atoms and molecules.

Symmetry breaking of the bending mode of CO2 in the presence of Ar

The weak infrared spectrum of CO2-Ar corresponding to the (0111) ← (0110) hot band of CO2 is detected in the region of the carbon dioxide ν3 fundamental vibration (≈2340 cm-1), using a tunable OPO laser source to probe a pulsed supersonic slit jet expansion. While this method was previously thought to cool clusters to the lowest rotational states of the ground vibrational state, here we show that under suitable jet expansion conditions, sufficient population remains in the first excited bending mode of CO2 (1-2%) to enable observation of vibrationally hot CO2-Ar, and thus to investigate the symmetry breaking of the intramolecular bending mode of CO2 in the presence of Ar. The bending mode of the CO2 monomer splits into an in-plane and an out-of-plane mode, strongly linked by a Coriolis interaction. Analysis of the spectrum yields a direct measurement of the in-plane/out-of-plane splitting measured to be 0.8770 cm-1. Calculations were carried out to determine if key features of our results, i.e., the sign and magnitude of the shift in the energy for the two intramolecular bending modes, are consistent with a quantum chemical potential energy surface. This aspect of intramolecular interactions has received little previous experimental and theoretical consideration. Therefore, we provide an additional avenue by which to study the intramolecular dynamics of this simplest dimer in its bending modes. Similar results should be possible for other weakly-bound complexes.

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.

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.

Non-covalent C-Cl…π interaction in acetylene-carbon tetrachloride adducts: Matrix isolation infrared and ab initio computational studies

Non-covalent halogen-bonding interactions between π cloud of acetylene (C2H2) and chlorine atom of carbon tetrachloride (CCl4) have been investigated using matrix isolation infrared spectroscopy and quantum chemical computations. The structure and the energies of the 1:1 C2H2-CCl4 adducts were computed at the B3LYP, MP2 and M05-2X levels of theory using 6-311++G(d,p) basis set. The computations indicated two minima for the 1:1 C2H2-CCl4 adducts; with the C-Cl…π adduct being the global minimum, where π cloud of C2H2 is the electron donor. The second minimum corresponded to a C-H…Cl adduct, in which C2H2 is the proton donor. The interaction energies for the adducts A and B were found to be nearly identical. Experimentally, both C-Cl…π and C-H…Cl adducts were generated in Ar and N2 matrixes and characterized using infrared spectroscopy. This is the first report on halogen bonded adduct, stabilized through C-Cl…π interaction being identified at low temperatures using matrix isolation infrared spectroscopy. Atoms in Molecules (AIM) and Natural Bond Orbital (NBO) analyses were performed to support the experimental results. The structures of 2:1 ((C2H2)2-CCl4) and 1:2 (C2H2-(CCl4)2) multimers and their identification in the low temperature matrixes were also discussed.

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.

Ab initio characterization of the Mg-HF van der Waals complex

The equilibrium structure and the three-dimensional potential energy surface of the Mg-HF van der Waals complex in its ground electronic state have been determined from accurate ab initio calculations using the coupled-cluster method, CCSD(T), in conjunction with the basis sets of triple- through quintuple-zeta quality. The core-electron correlation, high-order valence-electron correlation, and scalar relativistic effects were investigated. The Mg-HF complex was confirmed to be linear at equilibrium, with a vibrationless dissociation energy (into Mg and HF) D(e) of 280 cm(-1). The vibration-rotation energy levels of two isotopologues, (24)Mg-HF and (24)Mg-DF, were predicted using the variational method. The predicted spectroscopic constants can be useful in a further analysis of high-resolution vibration-rotation spectra of the Mg-HF complex.

High-resolution infrared spectroscopy of Mg-HF and Mg-(HF)2 solvated in helium nanodroplets

High-resolution infrared (IR) spectroscopy is used to investigate the Mg-HF and Mg-(HF)(2) van der Waals complexes. Both complexes are formed and probed within helium nanodroplets. Rotationally resolved zero-field and Stark spectra are assigned to a linear binary complex composed of a Mg atom bound to the hydrogen end of the HF molecule. Although high level ab initio calculations predict a fluorine bonded complex, none of the observed IR bands can be assigned to this complex. The collocation method is employed to determine the bound states on the two-dimensional intermolecular Mg-HF potential energy surface. The ground and first excited state wave functions for this potential surface have zero amplitude in the well corresponding to the fluorine bonded complex, consistent with experiment. The two HF stretching bands of the Mg-(HF)(2) complex are observed and assigned using a combination of the spectral symmetry, ab initio calculations, pick-up cell pressure dependencies, and dipole moment measurements. Comparisons with the helium solvated HF dimer show large changes to the HF stretching frequencies upon the addition of a single Mg atom to the hydrogen side of (HF)(2).

Effect of matrix on IR frequencies of acetylene and acetylene-methanol complex: Infrared matrix isolation andab initiostudy

Effect of nitrogen and argon matrices on the C-H asymmetric stretching and bending infrared frequencies of the acetylene molecule, C(2)H(2), has been studied by matrix isolation experiments as well as by calculations at MP2 level of theory. The complexes of C(2)H(2) in nitrogen and argon matrices, viz., C(2)H(2)(N(2))(m) (with m=2-8) and C(2)H(2)(Ar)(n) (with n=2-10) are theoretically explored. The computed acetylenic C-H asymmetric stretch in C(2)H(2)-nitrogen complexes shows a redshift of 3.0 to 11.9 cm(-1) compared with the frequencies of the free acetylene molecule, and a corresponding blueshift of 7.4 to 26.2 cm(-1) when C(2)H(2) is complexed with argon atoms. The trends in the computed shifts are in good agreement with the experiments. The molecular electrostatic potential minimum of C(2)H(2) becomes more negative when complexed with nitrogen than on complexation with argon. This observation implies a greater basic character for C(2)H(2) in the nitrogen matrix, favoring the formation of H-pi(C(2)H(2)-MeOH) complex as compared to that in the Ar matrix. Experimentally the preferential formation of H-pi(C(2)H(2)-MeOH) complex in the N(2) matrix has indeed been observed.

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.

Evidence for partial quenching of orbital angular momentum upon complex formation in the infrared spectrum of OH-acetylene

The entrance channel leading to the addition reaction between the hydroxyl radical and acetylene has been examined by spectroscopic characterization of the asymmetric CH stretching band of the pi-hydrogen bonded OH-acetylene reactant complex. The infrared action spectrum observed at 3278.6 cm(-1) (origin) consists of seven peaks of various intensities and widths, and is very different from those previously reported for closed-shell HF/HCl-acetylene complexes. The unusual spectrum arises from a partial quenching of the OH orbital angular momentum in the complex, which in turn is caused by a significant splitting of the OH monomer orbital degeneracy into (2)A(') and (2)A(") electronic states. The magnitude of the (2)A(')-(2)A(") splitting as well as the A rotational constant for the OH-acetylene complex are determined from the analysis of this b-type infrared band. The most populated OH product rotational state, j(OH)=9/2, is consistent with intramolecular vibrational energy transfer to the nu2 C triple bonded C stretching mode of the departing acetylene fragment. The lifting of the OH orbital degeneracy and partial quenching of its electronic orbital angular momentum indicate that the electronic changes accompanying the evolution of reactants into products have begun to occur in the reactant complex.

Collision Cross Sections, Pressure-Broadening Coefficients and Second Virial Coefficients for the Acetylene-Argon Complex: Experiments and Calculations on a New Potential Energy Surface

Integral cross sections and pressure-broadening coefficients have been measured by molecular beam scattering and by high-resolution infrared spectroscopy, respectively, for the acetylene-argon system. 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. Calculations of the scattering cross sections (both differential and integral), pressure-broadening, and second virial coefficients have been performed using both the present and also the most recent ab initio PES available in the literature. Analysis of the new experimental data indicates that the anisotropy of the interaction in the well region should be larger than that obtained in ab initio calculations. This is also in line with previous spectroscopic results.

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

Accurate ground-state intermolecular potential-energy surfaces are obtained for the HCCH-He, Ne, and Ar van der Waals complexes. The interaction energies are calculated at the coupled cluster singles and doubles including connected triple excitations level and fitted to analytic functions. For the three complexes we start with systematic basis set studies carried out at several intermolecular geometries, and using augmented correlation consistent polarized valence basis sets x-aug-cc-pVXZ (x=-,d; X=D,T,Q,5), also extended with a set of 3s3p2d1f1g midbond functions. The aug-cc-pVQZ-33211 surfaces of HCCH-He, Ne, and Ar complexes are characterized by absolute minima of -24.22, -50.20, and -122.17 cm(-1) at distances R between the rare-gas atom and the HCCH centers of mass of 4.35, 3.95, and 3.99 A, respectively; and at angles between the vector R and the HCCH main symmetry axis of 0 degrees , 43.3 degrees , and 60.6 degrees . The results are compared and considerably improve those previously available.

Quantum vibrational state-dependent potentials for classical many-body simulations

In this paper, we present a method for constructing simple state-dependent many-body potentials for quantum vibrations in a classical bath. The approach is based on an adiabatic separation between high-frequency quantum vibrational modes of the solute and the lower frequency classical motion of the solvent, and on a first-order perturbation theory description of the dependence of the quantum energies on bath configuration. In the simplest realization of the method, the delocalized quantum probability density of the vibrational mode is approximated by a sum of two delta functions, with positions and weights chosen to represent the lowest three moments of the exact distribution. Thus, in the many-body description of the system, each atom describing the quantum vibration is represented by a pair of particles. These quantum particles are held in rigid relative position and interact with the bath via potentials the magnitudes of which are modified by the delta-function weights. The resulting approach allows the classical molecular dynamics of molecules in arbitrary quantum vibrational states to be simulated with a little more effort than a purely classical description. The applicability of the method is illustrated in many-body simulations of the dephasing of vibrational superposition states of I(2) in a cryogenic krypton matrix, yielding results in good agreement with experiment.