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.     

Quantum solvation dynamics of HCN in a helium-4 droplet

Ultracold nanodroplets of helium-4, containing several thousands of He atoms, offer considerable promise as microscopic cryogenic chambers. Potential applications include the creation of tailor-made chemical or biomolecular complexes and studies of superfluidity in nanoscale systems. Recent experiments have succeeded in interrogating droplets of quantum solvent which consist of as few as 1-20 helium-4 atoms and which contain a single solute molecule. This allows the transition from a floppy, but essentially molecular, complex to a dissolved molecule to be followed and, surprisingly, the transition is found to occur quite rapidly, in some cases for as few as N = 7-20 solvent atoms. For example, in experiments on helium-4 droplets seeded with CO molecules [Tang and McKellar, J. Chem. Phys. 119, 754 (2003)], two series of transitions are observed which correlate with the a-type (Delta K = 0) and b-type (Delta K = +/-1) lines of the binary complex, CO-He (K is the quantum number associated with the projection of the total angular momentum onto the vector connecting the atom and the molecular center of mass). The a-type series, which evolves from the end-over-end rotational motion of the CO-He binary complex, saturates to the nanodroplet limit for as few as 10-15 helium-4 atoms, i.e., the effective moment of inertia of the molecule converges to its asymptotic (solvated) value quite rapidly. In contrast, the b-type series, which evolves from the free-molecule rotational mode, disappears altogether for N approximately 7 atoms. Similar behavior is observed in recent computational studies of HCN(4He)N droplets [Paolini et al., J. Chem. Phys. 123, 114306 (2005)]. In this article the quantum solvation of HCN in small helium-4 droplets is studied using a new fixed-node diffusion Monte Carlo (DMC) procedure. In this approach a Born-Oppenheimer-type separation of radial and angular motions is introduced as a means of computing nodal surfaces of the many-body wave functions which are required in the fixed-node DMC method. Excited rotational energies are calculated for HCN(4He)N droplets with N = 1-20: the adiabatic node approach also allows concrete physical mechanisms to be proposed for the predicted disappearance of the b-type series as well as the rapid convergence of the a-type series to the nanodroplet limit with increasing N. The behavior of the a-type series is traced directly to the mechanics of angular momentum coupling-and decoupling-between identical bosons and the molecular rotor. For very small values of N there exists significant angular momentum coupling between the molecule and the helium atoms: at N approximately 10 solvation appears to be complete as evidenced by significant decoupling of the molecule and solvent angular momenta. The vanishing of the b-type series is predicted to be a result of increasing He-He repulsion as the number of solvent atoms increases.     

Rotational structure of small 4He clusters seeded with HF, HCl, and HBr molecules

Diffusion Monte Carlo calculations are performed for ground and excited rotational states of HX(4He)N, complexes with N相似文献   

Carbon nanotubes in benzene: internal and external solvation

The structure and dynamics of benzene inside and outside of single-walled carbon nanotubes (SWNTs) in the (n,n) armchair configuration are studied via molecular dynamics computer simulations. Irrespective of the nanotube diameter, benzene molecules form cylindrical solvation shell structures on the outside of the nanotubes. Their molecular planes near the SWNTs in the first external solvation shell are oriented parallel to the nanotube surface, forming a π-stacked structure between the two. By contrast, the benzene distributions in the interior of the SWNTs are found to vary markedly with the nanotube diameter. In the case of the (7,7) and (8,8) nanotubes, internal benzene forms a single-file distribution, either in a vertex-to-vertex (n = 7) or face-to-face (n = 8) orientation between two neighboring molecules. Inside a slightly wider (9,9) nanotube channel, however, a cylindrical single-shell distribution of benzene arises. A secondary solvation structure, which begins to appear inside (10,10), develops into a full structure separate from the first internal solvation shell in (12,12). The ring orientation of internal benzene is generally parallel to the nanotube wall for n = 9-12, while it becomes either slanted with respect to (n = 7), or perpendicular to (n = 8), the nanotube axis. The confinement inside the small nanotube pores exerts a strong influence on the dynamics of benzene. Both translational and rotational dynamics inside SWNTs are slower and more anisotropic than in liquid benzene. It is also found that reorientational dynamics of internal benzene deviate dramatically from the rotational diffusion regime and change substantially with the nanotube diameter.     

DYNAMICS OF SINGLE-MOLECULE ROTATIONS ON SURFACES DEPEND ON SYMMETRY, INTERACTIONS AND MOLECULAR SIZES

Rotating surface-mounted molecules have attracted attention of many research groups as a way to develop new nanoscale devices and materials. However, mechanisms of motion of these rotors at the single-molecule level are still not well understood. Theoretical and experimental studies on thioether molecular rotors on gold surfaces suggest that the size of the molecules, their flexibility and steric repulsions with the surface are important for dynamics of the system. A complex combination of these factors leads to the observation that the rotation speeds have not been hindered by increasing the length of the alkyl chains. However, experiments on diferrocene derivatives indicated that a significant increase in the rotational barriers for longer molecules. We present here a comprehensive theoretical study that combines molecular dynamics simulations and simple models to investigate what factors influence single-molecule rotations on the surfaces. Our results suggest that rotational dynamics is determined by the size and by the symmetry of the molecules and surfaces, and by interactions with surfaces. Our theoretical predictions are in excellent agreement with current experimental observations.     

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.     

Infrared laser spectroscopy of the CH3-HCN radical complex stabilized in helium nanodroplets

The CH3-HCN and CD3-HCN radical complexes have been formed in helium nanodroplets by sequential pickup of a CH3 (CD3) radical and a HCN molecule and have been studied by high-resolution infrared laser spectroscopy. The complexes have a hydrogen-bonded structure with C3v symmetry, as inferred from the analysis of their rotationally resolved nu = 1 <-- 0 H-CN vibrational bands. The A rotational constants of the complexes are found to change significantly upon vibrational excitation of the C-H stretch of HCN within the complex, DeltaA = A'-A" = -0.04 cm(-1) (for CH3-HCN), whereas the B rotational constants are found to be 2.9 times smaller than that predicted by theory. The reduction in B can be attributed to the effects of helium solvation, whereas the large DeltaA is found to be a sensitive probe of the vibrational averaging dynamics of such weakly bound systems. The complex has a permanent electric dipole moment of 3.1 +/- 0.2 D, as measured by Stark spectroscopy. A vibration-vibration resonance is observed to couple the excited C-H stretching vibration of HCN within the complex to the lower-frequency C-H stretches of the methyl radical. Deuteration of the methyl radical was used to detune these levels from resonance, increasing the lifetime of the complex by a factor of 2. Ab initio calculations for the energies and molecular parameters of the stationary points on the CN+CH4 --> HCN+CH3 potential-energy surface are also presented.     

Photodissociation of laboratory oriented molecules: Revealing molecular frame properties of nonaxial recoil

We report the photodissociation of laboratory oriented OCS molecules. A molecular beam of OCS molecules is hexapole state-selected and spatially oriented in the electric field of a velocity map imaging lens. The oriented OCS molecules are dissociated at 230 nm with the linear polarization set at 45 degrees to the orientation direction of the OCS molecules. The CO(nu=0,J) photofragments are quantum state-selectively ionized by the same 230 nm pulse and the angular distribution is measured using the velocity map imaging technique. The observed CO(nu=0,J) images are strongly asymmetric and the degree of asymmetry varies with the CO rotational state J. From the observed asymmetry in the laboratory frame we can directly extract the molecular frame angles between the final photofragment recoil velocity and the permanent dipole moment and the transition dipole moment. The data for CO fragments with high rotational excitation reveal that the dissociation dynamics is highly nonaxial, even though conventional wisdom suggests that the nearly limiting beta parameter results from fast axial recoil dynamics. From our data we can extract the relative contribution of parallel and perpendicular transitions at 230 nm excitation.     

Rovibrational spectra for the HCCCN*HCN and HCN*HCCCN binary complexes in 4He droplets

Rovibrational spectra are measured for the HCCCN*HCN and HCN*HCCCN binary complexes in helium droplets at low temperature. Though no Q-branch is observed in the infrared spectrum of the linear HCN*HCCCN dimer, which is consistent with previous experimental results obtained for other linear molecules, a prominent Q-branch is found in the corresponding infrared spectrum of the HCCCN*HCN complex. This Q-branch, which is reminiscent of the spectrum of a parallel band of a prolate symmetric top, implies that some component of the total angular momentum is parallel to the molecular axis. The appearance of this particular spectroscopic feature is analyzed here in terms of a nonsuperfluid helium density induced by the molecular interactions. Finite temperature path integral Monte Carlo simulations are performed using potential energy surfaces calculated with second-order M?ller-Plesset perturbation theory, to investigate the structural and superfluid properties of both HCCCN*HCN(4He)N and HCN*HCCCN(4He)N clusters with N < or = 200. Explicit calculation of local and global nonsuperfluid densities demonstrates that this difference in the rovibrational spectra of the HCCCN*HCN and HCN*HCCCN binary complexes in helium can be accounted for by local differences in the superfluid response to rotations about the molecular axis, i.e., different parallel nonsuperfluid densities. The parallel and perpendicular nonsuperfluid densities are found to be correlated with the locations and strengths of extrema in the dimer interaction potentials with helium, differences between which derive from the variable extent of polarization of the CN bond in cyanoacetylene and the hydrogen-bonded CH unit in the two isomers. Calculation of the corresponding helium moments of inertia and effective rotational constants of the binary complexes yields overall good agreement with the experimental values.     

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.     

Fragmentation of HCN in optically selected mass spectrometry: nonthermal ion cooling in helium nanodroplets

A technique that combines infrared laser spectroscopy and helium nanodroplet mass spectrometry, which we refer to as optically selected mass spectrometry, is used to study the efficiency of ion cooling in helium. Electron-impact ionization is used to form He(+) ions within the droplets, which go on to transfer their charge to the HCN dopant molecules. Depending upon the droplet size, the newly formed ion either fragments or is cooled by the helium before fragmentation can occur. Comparisons with gas-phase fragmentation data suggest that the cooling provided by the helium is highly nonthermal. An "explosive" model is proposed for the cooling process, given that the initially hot ion is embedded in such a cold solvent.     

Exponential probe rotation in glass-forming liquids

Using time resolved optical depolarization, we have studied the rotational behavior of molecular probes in supercooled liquids near the glass transition temperature T(g). Simultaneously, the dynamics of the liquid immediately surrounding these rigid probes is measured by triplet state solvation experiments. This direct comparison of solute and solvent dynamics is particularly suited for assessing the origin of exponential orientational correlation functions of probe molecules embedded in liquids which exhibit highly nonexponential structural relaxation. Polarization angle dependent Stokes shift correlation functions demonstrate that probe rotation time and solvent response time are locally correlated quantities in the case of smaller probe molecules. Varying the size of both guest and host molecules shows that the size ratio determines the rotational behavior of the probes. The results are indicative of time averaging being at the origin of exponential rotation of probes whose rotational time constant is slower than solvent relaxation by a factor of 20 or more.     

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.     

Spectroscopic studies of OCS-doped 4He clusters with 9-72 helium atoms: observation of broad oscillations in the rotational moment of inertia

High-resolution spectra of HeN-OCS clusters with N up to 39 in the microwave region and up to 72 in the infrared region were observed with apparatus-limited line widths of about 15 kHz and 0.001 cm(-1), respectively. The cold (approximately 0.2 K) clusters were produced in pulsed supersonic jet expansions of very dilute OCS + He mixtures and probed using a microwave Fourier transform spectrometer or a tunable infrared diode laser spectrometer. Consistent analyses of the microwave and infrared data yield band origins for the carbonyl stretching vibration, together with rotational parameters for the ground and excited vibrational states. The rotational constant, B, passes through a minimum at N = 9 and then rises as the He atoms uncouple from the OCS rotational motion as a result of superfluid effects. There are broad unexpected oscillations in B, with maxima at N = 24 and 47 and minima at N = 36 and 62. The change in B upon vibrational excitation, which is negative for the OCS molecule, converges to positive values for N > 15. These results help to bridge the gap between individual molecules and bulk matter with atom-by-atom resolution over a significant range of cluster sizes.     

Dynamics of a supercooled ionic liquid studied by optical and dielectric spectroscopy

We have measured the dynamics of solvation of a triplet state probe, quinoxaline, in the glass-forming ionic liquid dibutylammonium formate near its glass transition temperature Tg=153 K. The Stokes-shift correlation function displays a relaxation time dispersion of considerable magnitude and the optical line width changes systematically along the solvation coordinate. The solvent dynamics in the viscous regime is compared with the rotational behavior of the solute and with the dielectric relaxation of the ionic liquid. Among the different quantities derived from the dielectric experiments, the relaxation of the macroscopic electric field, i.e., the modulus Mt, matches best the solvent response Ct regarding time scale and stretching exponent. Many other properties are reminiscent of the behavior of polar molecular liquids which lack the ionic character.     

Ultrafast solvent response upon a change of the solute size in non-polar supercritical fluids

Non-polar solvation dynamics has been investigated using steady-state absorption and emission spectroscopy of the NO A ²Σ⁺(3sσ) Rydberg state in fluid Ar over a wide range of densities spanning the supercritical regime. Equilibrium molecular dynamics simulations were implemented to derive a new isotropic NO A(3sσ)–Ar pair potential which was further used to investigate the role of local density enhancements on the solvation process by non-equilibrium molecular dynamics simulations. These density inhomogeneities were found to have no influence on the solvation dynamics. Furthermore, the latter was shown to take place in a strongly non-linear regime, especially at low temperatures. This process results from the dramatic change of solute–solvent short range interaction associated with the large solute size change upon excitation to the Rydberg state.     

Solid State Machinery of Multiple Dynamic Elements in a Metal–Organic Framework

Engineering coordinated rotational motion in porous architectures enables the fabrication of molecular machines in solids. A flexible two-fold interpenetrated pillared Metal-Organic Framework precisely organizes fast mobile elements such as bicyclopentane (BCP) (10⁷ Hz regime at 85 K), two distinct pyridyl rotors and E-azo group involved in pedal-like motion. Reciprocal sliding of the two sub-networks, switched by chemical stimuli, modulated the sizes of the channels and finally the overall dynamical machinery. Actually, iodine-vapor adsorption drives a dramatic structural rearrangement, displacing the two distinct subnets in a concerted piston-like motion. Unconventionally, BCP mobility increases, exploring ultra-fast dynamics (10⁷ Hz) at temperatures as low as 44 K, while the pyridyl rotors diverge into a faster and slower dynamical regime by symmetry lowering. Indeed, one pillar ring gained greater rotary freedom as carried by the azo-group in a crank-like motion. A peculiar behavior was stimulated by pressurized CO_2, which regulates BCP dynamics upon incremental site occupation. The rotary dynamics is intrinsically coupled to the framework flexibility as demonstrated by complementary experimental evidence (multinuclear solid-state NMR down to very low temperatures, synchrotron radiation XRD, gas sorption) and computational modelling, which helps elucidate the highly sophisticated rotor-structure interplay.     

Computer simulations of the solvation dynamics of Coumarin 153 in dimethylsulfoxide

We report molecular dynamics (MD) simulations of the solvation dynamics of Coumarin 153 in liquid dimethylsulfoxide using two distinct sets of partial charges for the coumarin probe. The excited state dipole moment of the coumarin and the dynamic Stokes shift in solution depend significantly on the type of charge distributions used. Nevertheless, the overall characteristics of the solvation responses obtained from both sets of charges are very similar and show good agreement with time-dependent Stokes shift experiments. Microscopic details of the solvent reorganization around the probe are discussed in light of the charge transfer upon photoexcitation.     

Path integral methods for rotating molecules in superfluids

We present a path integral Monte Carlo (PIMC) methodology for quantum simulation of molecular rotations in superfluid environments such as helium and para-hydrogen that combines the sampling of rotational degrees of freedom for a molecular impurity with multilevel Metropolis sampling of Bose permutation exchanges for the solvating species. We show how the present methodology can be applied to the evaluation of imaginary time rotational correlation functions of the molecular impurity, from which the effective rotational constants can be extracted. The combined rotation/permutation sampling approach allows for the first time explicit assessment of the effect of Bose permutations on molecular rotation dynamics, and the converse, i.e., the effect of molecular rotations on permutation exchanges and local superfluidity. We present detailed studies showing that the effect of Bose permutations in the solvating environment is more significant for the dynamics of heavy than light molecules in helium, and that Bose permutation exchanges are slightly enhanced locally by molecular rotation. Finally, the examples studied here reveal a size dependence of rotational excitations for molecules possessing a strongly anisotropic interaction with helium in 4HeN clusters between N approximately 20 and N approximately 10(3).     

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.