Infrared spectra of seeded hydrogen clusters: (paraH2)N - OCS, (orthoH2)N - OCS, and (HD)N - OCS, N = 2 - 7

Infrared spectra of hydrogen-carbonyl sulfide clusters containing paraH2, orthoH2, or HD have been studied in the 2060 cm(-1) region of the C-O stretching vibration. The clusters were formed in pulsed supersonic jet expansions and probed using a tunable infrared diode laser spectrometer. Simple symmetric rotor type spectra were observed and assigned for clusters containing up to N = 7 hydrogen molecules. There was no resolved K structure, and Q-branch features were present for orthoH2 and HD but absent for paraH2. These characteristics can be rationalized in terms of near symmetric rotor structures, very low effective rotational temperatures (0.15 to 0.6 K), and nuclear spin statistics. The observed vibrational shifts were compared with those from recent observations on the same clusters embedded in helium nanodroplets. The observed rotational constants for the paraH2 clusters are in good agreement with a recent quantum Monte Carlo simulation. Some mixed clusters were also observed, such as HD-HD-He-OCS and paraH2 - orthoH2 - OCS.     

Infrared spectra of helium clusters seeded with nitrous oxide, 4HeN-N2O, with N=1-80

High resolution infrared spectra of HeN-N2O clusters are studied in the 2200 cm(-1) region of the N2O nu1 fundamental band. The clusters are produced in a pulsed supersonic jet expansion from a cooled nozzle source and probed using a tunable diode laser operating in a rapid-scan mode. Three isotopic forms are used (14N14N16O, 15N14N16O, and 15N15N16O) in order to support the spectral analyses. For clusters up to N approximately 24, the individual spectra are resolved, assigned, and analyzed together with complementary microwave data. Assignments for larger clusters are uncertain due to overlapping transitions, but an approximate analysis is still possible for N approximately 25-80. Compared to helium clusters containing the related CO2 or OCS molecules, the rotational dynamics of HeN-N2O clusters show similarities but also important differences. In particular, HeN-N2O has more irregular behavior in the range of N=6-17, indicating that conventional molecular structure plays a greater role. In general terms, these differences can be attributed to a greater degree of angular anisotropy in the He-N2O intermolecular potential.     

High resolution infrared spectra of a carbon dioxide molecule solvated with helium atoms

Infrared spectra of He(N)-CO(2) clusters with N up to about 20 have been studied in the region of the CO(2) nu(3) fundamental band ( approximately 2350 cm(-1)) using a tunable diode laser spectrometer and pulsed supersonic jet source with cooled (>-150 degrees C) pinhole or slit nozzles and high backing pressures (<40 atm). Compared to previous studies of He(N)-OCS and -N(2)O clusters, the higher symmetry of CO(2) results in simpler spectra but less information content. Discrete rotation-vibration transitions have been assigned for N=3-17, and their analysis yields the variation of the vibrational band origin and B rotational constant over this size range. The band origin variation is similar to He(N)-OCS, with an initial blueshift up to N=5, followed by a monotonic redshift, consistent with a model where the first five He atoms fill a ring around the equator of the molecule, forcing subsequent He atom density to locate closer to the ends. The B value initially drops as expected for a normal molecule, reaching a minimum for N=5. Its subsequent rise for N=6 to 11 can be interpreted as the transition from a normal (though floppy) molecule to a quantum solvation regime, where the CO(2) molecule starts to rotate separately from the He atoms. For N>13, the B value becomes approximately constant with a value about 17% larger than that measured in much larger helium nanodroplets.     

Vibrational shifts of OCS in mixed clusters of parahydrogen and helium

We present a detailed theoretical study of the solvation structure and solvent induced vibrational shifts for an OCS molecule embedded in pure parahydrogen clusters and in mixed parahydrogen/helium clusters. The use of two recent OCS-(parahydrogen) and OCS-helium ab initio potential energy surfaces having explicit dependence on the asymmetric stretch of the OCS molecule allows calculation of the frequency shift of the OCS nu(3) vibration as a function of the cluster size and composition. We present results for clusters containing up to a full first solvation shell of parahydrogen (N=17 molecules), and up to M=128-N helium atoms. Due to the greater interaction strength of parahydrogen than helium with OCS, in the mixed clusters the parahydrogen molecules always displace He atoms in the first solvation shell around OCS and form multiple axial rings as in the pure parahydrogen clusters. In the pure clusters, the chemical potential of parahydrogen shows several magic numbers (N=8,11,14) that reflect an enhanced stability of axial rings containing one less molecule than required for complete filling at N=17. Only the N=14 magic number survives in the mixed clusters, as a result of different filling orders of the rings and greater delocalization of both components. The OCS vibration shows a redshift in both pure and mixed clusters, with N-dependent values that are in good agreement with the available experimental data. The dependence of the frequency shift on the cluster size and its composition is analyzed in terms of the parahydrogen and helium density distributions around the OCS molecule as a function of N and M. The frequency shift is found to be strongly dependent on the detailed distribution of the parahydrogen molecules in the pure parahydrogen clusters, and to be larger but show a smoother dependence on N in the presence of additional helium, consistent with the more delocalized nature of the mixed clusters.     

Computational spectroscopy of helium-solvated molecules: effective inertia, from small He clusters toward the nanodroplet regime

Accurate computer simulations of the rotational dynamics of linear molecules solvated in He clusters indicate that the large-size (nanodroplet) regime is attained quickly for light rotors (HCN) and slowly for heavy ones (OCS, N2O, and CO2), thus challenging previously reported results. Those results spurred the view that the different behavior of light rotors with respect to heavy ones-including a smaller reduction of inertia upon solvation of the former-would result from the lack of adiabatic following of the He density upon molecular rotation. We have performed computer experiments in which the rotational dynamics of OCS and HCN molecules was simulated using a fictitious inertia appropriate to the other molecule. These experiments indicate that the approach to the nanodroplet regime, as well as the reduction of the molecular inertia upon solvation, is determined by the anistropy of the potential, more than by the molecular weight. Our findings are in agreement with recent infrared and/or microwave experimental data which, however, are not yet totally conclusive by themselves.     

OCS in para-hydrogen clusters: rotational dynamics and superfluidity

We present a detailed analysis of the rotational excitations of the linear OCS molecule solvated by a variable number of para-hydrogen molecules (9 < or = N < or = 17). The effective rotational constant extracted from the fit of the rotational energy levels decreases up to N = 13, indicating near-rigid coupling between OCS rotations and para-hydrogen motion. Departure from rigidity is instead seen for larger clusters with 14 < or = N < or = 17. Path-integral Monte Carlo calculations show that the N dependence of the effective rotational constant can be explained in terms of a partial superfluid response of para-hydrogen to rotations about an axis perpendicular to the OCS axis. Complete para-hydrogen superfluid response to rotations about the OCS axis is found for N > or = 10.     

Spectroscopy of OCS-hydrogen clusters in He droplets

Clusters of para-hydrogen (pH2) and ortho-deuterium (oD2) have been assembled around an OCS chromophore molecule inside He droplets in a molecular beam and studied via IR diode laser depletion spectroscopy (nu approximately 2060 cm-1). The superfluid 4He droplets provide a gentle host ensuring a constant low temperature of either T = 0.38 K for 4He droplets or T = 0.15 K for both the pure 3He and mixed 4He-3He droplets. The spectra show well resolved rotational structure of the vibrational bands for each attached hydrogen molecule in the range n = 1-8. With only one (n = 1) attached pH2, HD or an oD2 molecule the best fit rotational constants were used to determine the structure of the complex, which was found to be in surprisingly good agreement with quantum chemical calculations for the free complex. With n = 5 and 6 the Q-branch disappears for the pH2 clusters but not for the oD2 clusters which is consistent with a donut model. The moments of inertia of the pH2 and the oD2 complexes are explained by a new model in which each of the 18 attached helium atoms in a shell surrounding the OCS molecule are assigned a mass of 0.55, while each attached H2 and D2 molecule has an effective mass of about 10 and 12 u, respectively.     

Infrared spectra of CO2-doped hydrogen clusters, (H2)N-CO2

Clusters of para-H(2) and/or ortho-H(2) containing a single carbon dioxide molecule are studied by high resolution infrared spectroscopy in the 2300 cm(-1) region of the CO(2) ν(3) fundamental band. The (H(2))(N)-CO(2) clusters are formed in a pulsed supersonic jet expansion from a cooled nozzle and probed using a rapid scan tunable diode laser. Simple symmetric rotor type spectra are observed with little or no resolved K-structure, and prominent Q-branch features for ortho-H(2) but not para-H(2). Observed rotational constants and vibrational shifts are reported for ortho-H(2) up to N = 7 and para-H(2) up to N = 15, with the N > 7 assignments only made possible with the help of theoretical simulations. The para-H(2) cluster with N = 12 shows clear evidence for superfluid effects, in good agreement with theory. The presence of larger clusters with N > 15 is evident in the spectra, but specific assignments are not possible. Mixed para- + ortho-H(2) cluster transitions are well predicted by linear interpolation between corresponding pure cluster line positions.     

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.     

Infrared spectra of seeded hydrogen clusters: (para-H2)N-N2O and (ortho-H2)N-N2O, N = 2-13

High-resolution infrared spectra of clusters containing para-H2 and/or ortho-H2 and a single nitrous oxide molecule are studied in the 2225-cm(-1) region of the upsilon1 fundamental band of N2O. The clusters are formed in pulsed supersonic jet expansions from a cooled nozzle and probed using a tunable infrared diode laser spectrometer. The simple symmetric rotor-type spectra generally show no resolved K structure, with prominent Q-branch features for ortho-H2 but not para-H2 clusters. The observed vibrational shifts and rotational constants are reported. There is no obvious indication of superfluid effects for para-H2 clusters up to N=13. Sharp transitions due to even larger clusters are observed, but no definite assignments are possible. Mixed (para-H2)N-(ortho-H2)M-N2O cluster line positions can be well predicted by linear interpolation between the corresponding transitions of the pure clusters.     

Infrared spectra of CO2-doped 4He clusters, 4HeN-CO2, with N=1-60

High resolution spectra of (4)He(N)-CO(2) clusters are studied in the region of the CO(2) nu(3) fundamental band (approximately 2300 cm(-1)). The clusters are produced in a pulsed supersonic jet expansion from a cooled nozzle source and probed by direct absorption using a tunable diode laser operating in a rapid-scan mode. Four carbon dioxide isotopes ((16)O(12)C(16)O, (16)O(13)C(16)O, (18)O(13)C(18)O, and (16)O(13)C(18)O) are used to support the analysis, and because additional rotational transitions are allowed for the asymmetric one ((16)O(13)C(18)O). Resolved R(0) (J=1<--0) rotation-vibration transitions are observed for clusters up to N=60. A detailed rotational analysis is possible up to N approximately 20 and, with some assumptions, to N approximately 37 and beyond. The derived rotational constants (B values) vary smoothly with N and show evidence for broad oscillations similar to those already reported for He(N)-OCS and He(N)-N(2)O. Possible indications of a disruption are observed in the J=2 levels of larger clusters (N>22) which could be caused by interactions with a "dark" helium cluster modes.     

Infrared absorption of gaseous ClCS detected with time-resolved Fourier-transform spectroscopy

A transient infrared absorption spectrum of gaseous ClCS was detected with a step-scan Fourier-transform spectrometer coupled with a multipass absorption cell. ClCS was produced upon irradiating a flowing mixture of Cl2CS and N2 or CO2 with a KrF excimer laser at 248 nm. A transient band in the region of 1160-1220 cm-1, which diminished on prolonged reaction, is assigned to the C-S stretching (nu1) mode of ClCS. Calculations with density-functional theory (B3P86 and B3LYP/aug-cc-pVTZ) predict the geometry, vibrational wave numbers, and rotational parameters of ClCS. The rotational contour of the spectrum of ClCS simulated based on predicted rotational parameters agrees satisfactorily with experimental observation; from spectral simulation, the band origin is determined to be at 1194.4 cm-1. Reaction kinetics involving ClCS, CS, and CS2 are discussed.     

Isotope effects in the infrared spectra of OCS-He complexes and clusters

Infrared spectra of the OCS-He van der Waals complex and of OCS-He(N) clusters have been studied in the region of the OCS nu1 fundamental band using a tunable diode laser to probe a pulsed supersonic slit jet. For the complex, the spectrum of the normal isotope, 16O12C32S-4He, has been considerably extended and the 34S- and 13C-substituted forms have been recorded for the first time. The data could be analyzed satisfactorily using a conventional asymmetric rotor Hamiltonian with sextic centrifugal distortion terms. For the clusters, the 34S- and 13C-substituted forms have been observed and assigned for N = 2-7, including some transitions with higher J values than previously reported for the normal isotope, e.g., R5. The observed vibrational shifts, relative to the free OCS molecule, were very similar to those of the normal isotope, and most of the difference could be explained by simple scaling. These results constitute a subtle and precise probe of intermolecular forces and dynamical effects in a system which is of current interest for cluster studies.     

Vibrational excitation and relaxation of five polyatomic molecules in an electrical discharge

Vibrational excitation and relaxation of five linear polyatomic molecules, OCS, OC3S, HC3N, HC5N, and SiC2S, have been studied by Fourier transform microwave spectroscopy in a supersonic expansion after the application of a low-current dc electric discharge. For each chain, the populations in bending and stretching modes have been characterized as a function of the applied discharge current; for stable OCS and HC3N, vibrational populations were studied as well in the absence of a discharge. With no discharge present the derived vibrational temperatures are slightly below T, the temperature of the gas before the supersonic expansion (i.e., 300 K). In the presence of the discharge, vibrational excitation occurs via inelastic collisions with the electrons and the vibrational temperatures rise as the applied current increases. Global vibrational relaxation is governed by rapid vibration-vibration (VV) energy transfer and slow vibration-translation (VT) energy transfer. The latter process is rate-determining and depends primarily on the wave number of the vibration. Vibrational modes with wave numbers near and below kT/hc (where T = 300 K and kT/hc-210 cm(-1)) are efficiently cooled by VT transfer because a sufficient number of collisions occur in the initial stages of the supersonic expansion. Vibrational modes with wave numbers around 450 cm(-l) appear to be inefficiently cooled in the molecular beam; at these energies VV and VT rates are probably comparable. For high-frequency vibrations, VV energy transfer dominates. For the longer chains OC3S and HC5N, higher-lying modes are generally not detectable and vibrational temperatures of most lower-lying modes were found to be lower than those of OCS and HC3N, suggesting that as the size of the molecules increases, intermode VV transfer becomes more efficient, plausibly due to the higher density of vibrational levels. New high resolution spectroscopic data have been obtained for several vibrationally excited states of OC3S, HC3N, and HC5N. Rotational lines of the 13C and 15N isotopic species of HC5N have been measured, yielding improved rotational and centrifugal distortion constants; 14N nitrogen quadrupole coupling constants for the isotopic species of HC5N with 13C have been determined for the first time.     

Solvation structure and rotational dynamics of LiH in 4He clusters

We present results of path integral Monte Carlo simulations of LiH solvated in superfluid 4He clusters of size up to N = 100. Despite the light mass of LiH and the strongly anisotropic LiH-He potential with a large repulsion at the hydrogen end, LiH is solvated inside the cluster for sufficiently large N. Using path integral correlation function analysis, we have determined the dipole (J = 1) rotational excitations of the cluster and a corresponding effective rotational constant Beff of the solvated LiH. We predict that Beff is greatly reduced with respect to the gas-phase rotational constant B, to a value of only about 6% of B. This exceptionally large reduction of the rotational constant is due to the highly anisotropic 4He solvation structure around LiH. It does not follow the previously established trend of a relatively small B reduction for light molecules, showing the strongest reduction of all molecules in 4He to date. Comparison of the calculated rotational spectra of LiH in helium obeying Bose and Boltzmann statistics, respectively, demonstrates that the Bose statistics of helium is an essential requirement for obtaining well-defined molecule rotational spectra in helium-4.     

High resolution infrared spectroscopy of carbon dioxide clusters up to (CO2)13

Thirteen specific infrared bands in the 2350 cm(-1) region are assigned to carbon dioxide clusters, (CO(2))(N), with N = 6, 7, 9, 10, 11, 12 and 13. The spectra are observed in direct absorption using a tuneable infrared laser to probe a pulsed supersonic jet expansion of a dilute mixture of CO(2) in He carrier gas. Assignments are aided by cluster structure calculations made using two reliable CO(2) intermolecular potential functions. For (CO(2))(6), two highly symmetric isomers are observed, one with S(6) symmetry (probably the more stable form), and the other with S(4) symmetry. (CO(2))(13) is also symmetric (S(6)), but the remaining clusters are asymmetric tops with no symmetry elements. The observed rotational constants tend to be slightly (≈2%) smaller than those from the predicted structures. The bands have increasing vibrational blueshifts with increasing cluster size, similar to those predicted by the resonant dipole-dipole interaction model but significantly larger in magnitude.     

Rotational dynamics of CO solvated in small He clusters: a quantum Monte Carlo study

The rotational dynamics of CO single molecules solvated in small He clusters (CO @ HeN) has been studied using reptation quantum Monte Carlo simulations for cluster sizes up to N = 30. Our results are in good agreement with the rotovibrational features of the infrared spectrum recently determined for this system and provide a deep insight into the relation between the structure of the cluster and its dynamics. Simulations for large N also provide a prediction of the effective moment of inertia of CO in the He nanodroplet regime, which has not been measured so far.     

Rotational excitations of N2O in small helium clusters and the role of Bose permutation symmetry

We present a detailed study of the energetics, structures, and Bose properties of small clusters of (4)He containing a single nitrous oxide (N(2)O) molecule, from N=1 (4)He up to sizes corresponding to completion of the first solvation shell around N(2)O (N=16 (4)He). Ground state properties are calculated using the importance-sampled rigid-body diffusion Monte Carlo method, rotational excited state calculations are made with the projection operator imaginary time spectral evolution method, and Bose permutation exchange and associated superfluid properties are calculated with the finite temperature path integral method. For N< or =5 the helium atoms are seen to form an equatorial ring around the molecular axis, at N=6 helium density starts to occupy the second (local) minimum of the N(2)O-He interaction at the oxygen side of the molecule, and N=9 is the critical size at which there is onset of helium solvation all along the molecular axis. For N> or =8 six (4)He atoms are distributed in a symmetric, quasirigid ring around N(2)O. Path integral calculations show essentially complete superfluid response to rotation about the molecular axis for N> or =5, and a rise of the perpendicular superfluid response from zero to appreciable values for N> or =8. Rotational excited states are computed for three values of the total angular momentum, J=1-3, and the energy levels fitted to obtain effective spectroscopic constants that show excellent agreement with the experimentally observed N dependence of the effective rotational constant B(eff). The non-monotonic behavior of the rotational constant is seen to be due to the onset of long (4)He permutation exchanges and associated perpendicular superfluid response of the clusters for N> or =8. We provide a detailed analysis of the role of the helium solvation structure and superfluid properties in determining the effective rotational constants.     

Rotational spectra of the van der Waals complexes of molecular hydrogen and OCS

The a- and b-type rotational transitions of the weakly bound complexes formed by molecular hydrogen and OCS, para-H2-OCS, ortho-H2-OCS, HD-OCS, para-D2-OCS, and ortho-D2-OCS, have been measured by Fourier transform microwave spectroscopy. All five species have ground rotational states with total rotational angular momentum J=0, regardless of whether the hydrogen rotational angular momentum is j=0 as in para-H2, ortho-D2, and HD or j=1 as in ortho-H2 and para-D2. This indicates quenching of the hydrogen angular momentum for the ortho-H2 and para-D2 species by the anisotropy of the intermolecular potential. The ground states of these complexes are slightly asymmetric prolate tops, with the hydrogen center of mass located on the side of the OCS, giving a planar T-shaped molecular geometry. The hydrogen spatial distribution is spherical in the three j=0 species, while it is bilobal and oriented nearly parallel to the OCS in the ground state of the two j=1 species. The j=1 species show strong Coriolis coupling with unobserved low-lying excited states. The abundance of para-H2-OCS relative to ortho-H2-OCS increases exponentially with decreasing normal H2 component in H2He gas mixtures, making the observation of para-H2-OCS in the presence of the more strongly bound ortho-H2-OCS dependent on using lower concentrations of H2. The determined rotational constants are A=22 401.889(4) MHz, B=5993.774(2) MHz, and C=4602.038(2) MHz for para-H2-OCS; A=22 942.218(6) MHz, B=5675.156(7) MHz, and C=4542.960(7) MHz for ortho-H2-OCS; A=15 970.010(3) MHz, B=5847.595(1) MHz, and C=4177.699(1) MHz for HD-OCS; A=12 829.2875(9) MHz, B=5671.3573(7) MHz, and C=3846.7041(6) MHz for ortho-D2-OCS; and A=13 046.800(3) MHz, B=5454.612(2) MHz, and C=3834.590(2) MHz for para-D2-OCS.