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 fluctuation of molecules in quantum clusters. II. Molecular rotation and superfluidity in OCS-doped helium-4 clusters

In this paper, quantum fluctuations of a carbonyl sulfide molecule in helium-4 clusters are studied as a function of cluster size N in a small-to-large size regime (2or=20, which is larger than the experimental nanodroplet value. Superfluid response of the doped cluster is found to show remarkable anisotropy especially for N相似文献   

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.     

Commensurate-incommensurate transition of 4He adsorbed on a single C60 molecule

Path-integral Monte Carlo calculations have been performed to study (4)He adsorption on a single C(60) molecule. Helium corrugations on the fullerene molecular surface are incorporated with the (4)He-C(60) interaction described by the sum of all (4)He-C interatomic pair potentials. Radial density distributions show a layer-by-layer growth of (4)He with the first adlayer being located at a distance of ~6.3 ? from the center of the C(60) molecule. The monolayer shows different quantum states as the number of (4)He adatoms N varies. For N = 32, we find a commensurate solid, with each of the 32 adsorption sites on the molecular surface being occupied by a single (4)He atom. Various domain-wall structures are observed as more (4)He atoms are added and the first layer crystallizes into an incommensurate solid when it is completely filled. This commensurate-incommensurate transition of the helium monolayer is found to be accompanied by re-entrant superfluid response at a low temperature of 0.31 K with the superfluidity being totally quenched at N = 32, 44, and 48. Finally, the different quantum states observed in the helium monolayer around C(60) are compared with phase diagrams proposed for the corresponding layer on a graphite surface.     

碘乙醇在Ni(100)表面的化学热反应

利用热脱附(TPD)实验和X射线光电子能谱(XPS)研究了碘乙醇在Ni(100)表面的吸附和热反应过程. 实验结果表明, 碘乙醇在100 K温度下以两种分子的形式吸附在Ni(100)的表面, 既有碘原子端的吸附也有碘原子端和羟基端同时吸附在表面. 热分解反应发生在140 K, 伴有少量的乙烯和水产生. 碘乙醇在150 K经过C—I键断裂形成&#8722;O(H)CH2CH2&#8722;和羟乙基两种中间产物. 在160 K温度下&#8722;O(H)CH2CH2&#8722;脱去氢形成&#8722;OCH2CH2&#8722;氧金属环. 中间产物经过进一步分解氧化反应分别在210和250 K产生乙醛, 一部分乙醛从表面脱出, 而其余的则分解成氢气、水和CO.     

Treatment of dilute clusters of methanol and water by ab initio quantum mechanical calculations

Large molecular clusters can be considered as intermediate states between gas and condensed phases, and information about them can help us understand condensed phases. In this paper, ab initio quantum mechanical methods have been used to examine clusters formed of methanol and water molecules. The main goal was to obtain information about the intermolecular interactions and the structure of methanol/water clusters at the molecular level. The large clusters (CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10)) containing one molecule of one component (methanol or water) and many (12, 10) molecules of the other component were considered. M?ller-Plesset perturbation theory (MP2) was used in the calculations. Several representative cluster geometries were optimized, and nearest-neighbor interaction energies were calculated for the geometries obtained in the first step. The results of the calculations were compared to the available experimental information regarding the liquid methanol/water mixtures and to the molecular dynamics and Monte Carlo simulations, and good agreement was found. For the CH(4)O...(H(2)O)(12) cluster, it was shown that the molecules of water can be subdivided into two classes: (i) H bonded to the central methanol molecule and (ii) not H bonded to the central methanol molecule. As expected, these two classes exhibited striking energy differences. Although they are located almost the same distance from the carbon atom of the central methanol molecule, they possess very different intermolecular interaction energies with the central molecule. The H bonding constitutes a dominant factor in the hydration of methanol in dilute aqueous solutions. For the H(2)O...(CH(4)O)(10) cluster, it was shown that the central molecule of water has almost three H bonds with the methanol molecules; this result differs from those in the literature that concluded that the average number of H bonds between a central water molecule and methanol molecules in dilute solutions of water in methanol is about two, with the water molecules being incorporated into the chains of methanol. In contrast, the present predictions revealed that the central water molecule is not incorporated into a chain of methanol molecules, but it can be the center of several (2-3) chains of methanol molecules. The molecules of methanol, which are not H bonded to the central water molecule, have characteristics similar to those of the methane molecules around a central water molecule in the H(2)O...(CH(4))(10) cluster. The ab initio quantum mechanical methods employed in this paper have provided detailed information about the H bonds in the clusters investigated. In particular, they provided full information about two types of H bonds between water and methanol molecules (in which the water or the methanol molecule is the proton donor), including information about their energies and lengths. The average numbers of the two types of H bonds in the CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10) clusters have been calculated. Such information could hardly be obtained with the simulation methods.     

On the possible "supersolid" character of parahydrogen clusters

We present results of a theoretical study of structural and superfluid properties of parahydrogen (p-H(2)) clusters comprising 25, 26, and 27 molecules at low temperature. The microscopic model utilized here is based on the Silvera-Goldman pair potential. Numerical results are obtained by means of quantum Monte Carlo simulations, making use of the continuous-space worm algorithm. The clusters are superfluid in the low temperature limit, but display markedly different physical behaviors. For N = 25 and 27, superfluidity at low temperature arises as clusters melt, that is, become progressively liquid-like as a result of quantum effects. On the other hand, for N = 26, the cluster remains rigid and solid-like. We argue that the cluster (p-H(2))(26) can be regarded as a mesoscopic "supersolid". This physical picture is supported by results of simulations in which a single p-H(2) molecule in the cluster is isotopically substituted.     

Clusters containing open-shell molecules: minimum-energy structures and low-lying isomers of ArnCH (X 2 pi), n = 1 to 15

The size evolution of the equilibrium structures of open-shell ArnCH (X 2 pi) Van der Waals clusters is investigated for n = 1 to 15. We describe a method for combining pair potentials for Ar-CH and Ar-Ar interactions to obtain potential energy surfaces for ArnCH clusters. For each cluster size considered, the global and a few energetically close local minima are calculated using simulated annealing followed by a direct minimization scheme. Ar2CH is found to have an unusually stable planar structure, which persists as a motif in larger ArnCH clusters and has a strong effect on their optimal geometries. The lowest-energy isomers of ArnCH with n = 3 to 11 have all Ar atoms in a shell around CH. The only exception is Ar4CH, where the fully solvated isomer is 3 cm-1 higher in energy than the optimal isomer with CH bound to the surface of the Ar4 tetrahedron. For n = 7 to 11, the minimum-energy structure of ArnCH derives from the global minimum of the Arn + 1 cluster, by replacing the Ar atom at the bottom of the pentagonal bipyramid with CH. The lowest-energy structure of Ar12CH is that of the optimal icosahedral Ar13 cluster, with CH replacing one of the Ar atoms on the cluster surface. This structure supports the proposition based on the spectroscopic data, that for ArnCH clusters with about 10 to 50 Ar atoms CH resides on the surface of Arn.     

Ion-molecule reactions in 4He droplets: flying nano-cryo-reactors

Ion-molecule reactions are studied inside large (approximately equal to 10(4) atoms) very cold (0.37 K) superfluid (4)He droplets by mass spectrometric detection of the product ions. He+ ions initially formed inside the droplets by electron impact ionization undergo charge transfer with either embedded D(2), N(2), or CH(4). For D(2) this charge transfer process was studied in detail by varying the pickup pressure. For either N(2) or CH(4) the reagent ions were formed by this charge transfer and the reaction pathways of the secondary reactions N(2) (+)+D(2), CH(4) (+)+D(2), and CH(3) (+)+D(2) each with an additionally embedded D(2) molecule were also determined from the pickup pressure dependencies. In several cases, notably He.N(2) (+) and CH(3)D(2) (+) reaction intermediates are observed. The analysis is facilitated by the tendency for molecular ion products to appear without (or with only very few) attached He atoms whereas the atomic ion products usually appear in the mass spectra with several attached He atoms, e.g., He(m).D+ ions with up to m=19.     

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 and ab initio calculations for the Cl--(CH4)n (n = 1-10) anion clusters

In an effort to elucidate their structures, mass-selected Cl--(CH4)n (n = 1-10) clusters are probed using infrared spectroscopy in the CH stretch region (2800-3100 cm(-1)). Accompanying ab initio calculations at the MP2/6-311++G(2df,2p) level for the n = 1-3 clusters suggest that methane molecules prefer to attach to the chloride anion by single linear H-bonds and sit adjacent to one another. These conclusions are supported by the agreement between experimental and calculated vibrational band frequencies and intensities. Infrared spectra in the CH stretch region for Cl--(CH4)n clusters containing up to ten CH4 ligands are remarkably simple, each being dominated by a single narrow peak associated with stretching motion of hydrogen-bonded CH groups. The observations are consistent with cluster structures in which at least ten equivalent methane molecules can be accommodated in the first solvation shell about a chloride anion.     

Theoretical infrared spectrum of the ethanol hexamer

The hydrogen bond network of ethanol clusters is among the most complex hydrogen bond networks of molecular clusters. One of the reasons of its complexity arises from the number of possible ethanol monomers (there are three isoenergetic isomers of the ethanol monomer). This leads to difficulties in the exploration of potential energy surfaces (PESs) of ethanol clusters. In this work, we have explored the PES of the ethanol hexamer at the MP2/aug-cc-pVDZ level of theory. We have provided structures and their relative stability at 0 K and for temperatures ranging from 20 to 400 K in the gas phase. These structures are used to compute the theoretical infrared (IR) spectrum of the ethanol hexamer at the MP2/aug-cc-pVDZ level of theory. As a result, 98 different structures have been investigated, and six isomers are reported to be the most isoenergetically stable structures of the ethanol hexamer. These isomers are folded cyclic structures in which the stability is enhanced by the implication of CH⋯O interactions. Our investigations show that the PES of the ethanol hexamer is very flat, yielding several isoenergetic structures. Furthermore, we have noted that several isomers contribute to the population of the ethanol hexamer at high temperatures. As far as the IR spectroscopic study is concerned, we have found that the IR spectra of the most stable structures are in good agreement with the experiment. Considering this agreement, these structures are used to assign the experimental peaks in the CH-stretching region. We concluded that the stability of the structures of the ethanol hexamer is related both to OH⋯O hydrogen bonds and CH⋯O interactions. Overall, we have found that the IR spectrum of the ethanol hexamer, calculated from the contribution of all the possible stable structures weighted by their probability, excellently reproduce the experimental spectrum of the ethanol hexamer.     

Ab initio studies of alkyl radical reactions: Combination and disproportionation reactions of CH3 with C2H5, and the decomposition of chemically activated C3H8

This paper reports the first quantitative ab initio prediction of the disproportionation/combination ratio of alkyl+alkyl reactions using CH3+C2H5 as an example. The reaction has been investigated by the modified Gaussian-2 method with variational transition state or Rice-Ramsperger-Kassel-Marcus calculations for several channels producing (1) CH4+CH2CH2, (2) C3H8, (3) CH4CH3CH, (4) H2+CH3CHCH2, (5) H2+CH3CCH3, and (6) C2H6+CH2 by H-abstraction and association/decomposition mechanisms through singlet and triplet potential energy paths. Significantly, the disproportionation reaction (1) producing CH4+C2H4 was found to occur primarily by the lowest energy path via a loose hydrogen-bonding singlet molecular complex, H3CHC2H4, with a 3.5 kcal/mol binding energy and a small decomposition barrier (1.9 kcal/mol), instead of a direct H-abstraction process. Bimolecular reaction rate constants for the formation of the above products have been calculated in the temperature range 300-3000 K. At 1 atm, formation of C3H8 is dominant below 1200 K. Over 1200 K, the disproportionation reaction becomes competitive. The sum of products (3)-(6) accounts for less than 0.3% below 1500 K and it reaches around 1%-4% above 2000 K. The predicted rate constant for the disproportionation reaction with multiple reflections above the complex well, k1=5.04 x T(0.41) exp(429/T) at 200-600 K and k1=1.96 x 10(-20) T(2.45) exp(1470/T) cm3 molecule(-1) s(-1) at 600-3000 K, agrees closely with experimental values. Similarly, the predicted high-pressure rate constants for the combination reaction forming C3H8 and its reverse dissociation reaction in the temperature range 300-3000 K, k2(infinity)=2.41 x 10(-10) T(-0.34) exp(259/T) cm3 molecule(-1) s(-1) and k(-2)(infinity)=8.89 x 10(22) T(-1.67)exp(-46 037/T) s(-1), respectively, are also in good agreement with available experimental data.     

Temperature and structural changes of water clusters in vacuum due to evaporation

This paper presents a study on evaporation of pure water clusters. Molecular dynamics simulations between 20 ns and 3 micros of clusters ranging from 125 to 4096 molecules in vacuum were performed. Three different models (SPC, TIP4P, and TIP5P) were used to simulate water, starting at temperatures of 250, 275, and 300 K. We monitored the temperature, the number of hydrogen bonds, the tetrahedral order, the evaporation, the radial distribution functions, and the diffusion coefficients. The three models behave very similarly as far as temperature and evaporation are concerned. Clusters starting at a higher temperature show a higher initial evaporation rate and therefore reach the point where evaporation stop (around 240 K) sooner. The radius of the clusters is decreased by 0.16-0.22 nm after 0.5 micros (larger clusters tend to decrease their radius slightly more), which corresponds to around one evaporated molecule per nm(2). The cluster temperature seems to converge towards 215 K independent of cluster size, when starting at 275 K. We observe only small structural changes, but the clusters modeled by TIP5P show a larger percentage of molecules with low diffusion coefficient as t-->infinity, than those using the two other water models. TIP4P seems to be more structured and more hydrogen bonds are formed than in the other models as the temperature falls. The cooling rates are in good agreement with experimental results, and evaporation rates agree well with a phenomenological expression based on experimental observations.     

Finite temperature path integral Monte Carlo simulations of structural and dynamical properties of Ar(N)-CO(2) clusters

We report finite temperature quantum mechanical simulations of structural and dynamical properties of Ar(N)-CO(2) clusters using a path integral Monte Carlo algorithm. The simulations are based on a newly developed analytical Ar-CO(2) interaction potential obtained by fitting ab initio results to an anisotropic two-dimensional Morse∕Long-range function. The calculated distributions of argon atoms around the CO(2) molecule in Ar(N)-CO(2) clusters with different sizes are consistent to the previous studies of the configurations of the clusters. A first-order perturbation theory is used to quantitatively predict the CO(2) vibrational frequency shift in different clusters. The first-solvation shell is completed at N = 17. Interestingly, our simulations for larger Ar(N)-CO(2) clusters showed several different structures of the argon shell around the doped CO(2) molecule. The observed two distinct peaks (2338.8 and 2344.5 cm(-1)) in the υ(3) band of CO(2) may be due to the different arrangements of argon atoms around the dopant molecule.     

N_2H_4-CH_3OH氢键团簇体系的从头计算

用从头计算法研究了 (N2 H4-CH3OH)氢键团簇体系。分别在HF/6 31G 和HF/6 31G 水平上对它们的中性和离子团簇进行几何全优化 ,得到了 3种中性混合团簇稳定构型和离子混合团簇稳定构型 ,并对其能量和稳定性进行了比较。讨论了 3种不同构型离子团簇可能的解离通道。给出了质子化混合团簇的稳定构型 ,并对其可能的解离通道进行了讨论。文中最后计算出N2 H4,CH3OH ,(N2 H4-CH3OH)团簇的质子亲和能 (PA) ,分别为 :2 0 6.7kcal/mol,1 78.3kcal/mol,2 2 7.5kcal/mol,其中质子亲和能PAcalc[N2 H4]与实验值PAexp[N2 H4]=2 0 4 .8kcal/mol符合得很好。     

Substituent effects in five oxo-centered trinuclear rhodium(III) clusters

We here report the rates of water substitution by methanol-d(4) for four new oxo-centered trinuclear rhodium(III) clusters with different carboxylate-bridging ligands, [Rh3(mu3-O)(mu-O2CR)6(OH2)3]+ (R = CH2CH3, CH2CH2Cl, CH2Cl, and CHCl2), and [Rh3(mu3-O)(mu-O2CCH3)6(OH2)3]+. By varying the R group alkyl chain, water substitution rates were found to span almost 3 orders of magnitude (k(298K) = 1.2 x 10(-2)-2.3 x 10(-5) s(-1)) and reflect the following trend R = CH2CH3 > CH3 > CH2CH2Cl > CH(2)Cl > CHCl2. Activation parameters for substitution point toward a dissociative activation pathway (DeltaH = 99-115 kJ mol(-1); DeltaS = 48-52 J mol(-1) K(-1)), indicating that there is little association with the incoming methanol molecule during the formation of the transition-state complex. Because the mechanism for substitution in all five trimers has a considerable dissociative character, substitution rates are likely very similar to water exchange rates. These data suggest that the kinetic reactivity of the ligated waters is heavily influenced by the inductive ability of the aliphatic substituents, but yet the mechanism of substitution remains virtually unchanged. Structural data are also reported for the four new rhodium(III) trimer salts as well as 103Rh NMR spectra. We find that 103Rh NMR chemical shifts span more than 200 ppm and mirror the same reactivity trend found for the rates of water substitution (103Rh delta (9406-9620 ppm): R = CH2CH3 < CH3 < CH2CH2Cl < CH2Cl < CHCl2). Taken together, these data suggest a means for estimating water exchange rates for other oxo-centered rhodium(III) trimers from chemical shift data alone.     

Metal core bonding motifs of monodisperse icosahedral Au13 and larger Au monolayer-protected clusters as revealed by X-ray absorption spectroscopy and transmission electron microscopy

The atomic metal core structures of the subnanometer clusters Au13[PPh3]4[S(CH2)11CH3]2Cl2 (1) and Au13[PPh3]4[S(CH2)11CH3]4 (2) were characterized using advanced methods of electron microscopy and X-ray absorption spectroscopy. The number of gold atoms in the cores of these two clusters was determined quantitatively using high-angle annular dark field scanning transmission electron microscopy. Multiple-scattering-path analyses of extended X-ray absorption fine structure (EXAFS) spectra suggest that the Au metal cores of each of these complexes adopt an icosahedral structure with a relaxation of the icosahedral strain. Data from microscopy and spectroscopy studies extended to larger thiolate-protected gold clusters showing a broader distribution in nanoparticle core sizes (183 +/- 116 Au atoms) reveal a bulklike fcc structure. These results further support a model for the monolayer-protected clusters (MPCs) in which the thiolate ligands bond preferentially at 3-fold atomic sites on the nanoparticle surface, establishing an average composition for the MPC of Au180[S(CH2)11CH3]40. Results from EXAFS measurements of a gold(I) dodecanethiolate polymer are presented that offer an alternative explanation for observations in previous reports that were interpreted as indicating Au MPC structures consisting of a Au core, Au2S shell, and thiolate monolayer.     

Production of methyl-oxonium ion and its complexes in the core-excited (HC(O)OCH3)n clusters: H/H+ transfer from the alpha carbonyl

Dissociation of free methyl-formate (MF), HC(O)OCH3, and its clusters (MF)n, (HC(O)OCH3)n, induced by core-level excitation was studied near the oxygen K edge by time-of-flight fragment-mass spectroscopy. Besides the protonated clusters, (MF)nH+ with n < or = 15, we identified the production for another series of (MF)mCH3OH2+ with m < or = 14 as well as methyl-oxonium ion, CH3OH2+, characteristic of hydrogen transfer reactions in the cationic clusters. Here; specifically labeled methyl-formate-d (MFD), DC(O)OCH3 was also used to examine the core-excited dissociation mechanisms. Deuterium-labeled experiments indicated that MFD+ with low internal energies, partially generated after the core excitation, produces CH3OD+ via a site-specific deuterium transfer from the alpha carbonyl in the molecular cation and that CH3OD2+ can be formed via the successive transfer of another deuterium from the neighbor molecule in the clusters. The deuteron (proton) transfer was also found to take place preferentially from the alpha carbonyl of the neighbor molecule for the production of deuteronated (MFD)nD+, (protonated (MF)nH+), clusters. The minimal energy requirement paths were examined for dimer (MF)2+ cation to support the present dissociation mechanisms of core-excited (MF)n clusters using ab initio molecular-orbital calculations.     

气体分子对甲烷水合物稳定性的影响

通过B3LYP方法, 在6-31G(d,p)水平下, 分别优化了结构I型甲烷水合物十二面体和十四面体晶穴结构. 结果表明, CH4分子使晶穴的相互作用能降低, 增强了晶穴的稳定性. 计算了晶穴中甲烷分子C—H键的对称伸缩振动频率, 计算结果与实验值相符合. 研究发现CH4分子影响晶穴中氧原子的电荷分布, 从而增强了氢键的稳定性. 通过分子动力学方法研究水合物晶胞中气体分子的占有率对水合物稳定性的影响, 进一步说明气体分子对水合物晶穴稳定性的重要作用.