Ionic dimers in He droplets: interaction potentials for Li2(+)-He,Na2(+)-He, and K2(+)-He and stability of the smaller clusters

We present post Hartree-Fock calculations of the potential energy surfaces (PESs) for the ground electronic states of the three alkali dimer ions Li(2) (+),Na(2) (+), and K(2) (+) interacting with neutral helium. The calculations were carried out for the frozen molecular equilibrium geometries and for an extensive range of the remaining two Jacobi coordinates, R and theta, for which a total of about 1000 points is generated for each surface. The corresponding raw data were then fitted numerically to produce analytic expressions for the three PESs, which were in turn employed to evaluate the bound states of the three trimers for their J=0 configurations: The final spatial features of such bound states are also discussed in detail. The possible behavior of additional systems with more helium atoms surrounding the ionic dopants is gleaned from further calculations on the structural stability of aggregates with up to six He atoms. The validity of a sum-of-potential approximation to yield realistic total energies of the smaller cluster is briefly discussed vis-a-vis the results from many-body calculations.     

Solvation of Na+, K+, and their dimers in helium

Helium atoms bind strongly to alkali cations which, when embedded in liquid helium, form so-called snowballs. Calculations suggest that helium atoms in the first solvation layer of these snowballs form rigid structures and that their number (n) is well defined, especially for the lighter alkalis. However, experiments have so far failed to accurately determine values of n. We present high-resolution mass spectra of Na(+)He(n), K(+)He(n), Na(2)(+)He(n) and K(2)(+)He(n), formed by electron ionization of doped helium droplets; the data allow for a critical comparison with several theoretical studies. For sodium and potassium monomers the spectra indicate that the value of n is slightly smaller than calculated. Na(2)(+)He(n) displays two distinct anomalies at n=2 and n=6, in agreement with theory; dissociation energies derived from experiment closely track theoretical values. K(2)(+)He(n) distributions are fairly featureless, which also agrees with predictions.     

Fragmentation of ionized doped helium nanodroplets: theoretical evidence for a dopant ejection mechanism

We report a theoretical study of the effect induced by a helium nanodroplet environment on the fragmentation dynamics of a dopant. The dopant is an ionized neon cluster Ne(n) (+) (n=4-6) surrounded by a helium nanodroplet composed of 100 atoms. A newly designed mixed quantum/classical approach is used to take into account both the large helium cluster zero-point energy due to the light mass of the helium atoms and all the nonadiabatic couplings between the Ne(n) (+) potential-energy surfaces. The results reveal that the intermediate ionic dopant can be ejected from the droplet, possibly with some helium atoms still attached, thereby reducing the cooling power of the droplet. Energy relaxation by helium atom evaporation and dissociation, the other mechanism which has been used in most interpretations of doped helium cluster dynamics, also exhibits new features. The kinetic energy distribution of the neutral monomer fragments can be fitted to the sum of two Boltzmann distributions, one with a low kinetic energy and the other with a higher kinetic energy. This indicates that cooling by helium atom evaporation is more efficient than was believed so far, as suggested by recent experiments. The results also reveal the predominance of Ne(2) (+) and He(q)Ne(2) (+) fragments and the absence of bare Ne(+) fragments, in agreement with available experimental data (obtained for larger helium nanodroplets). Moreover, the abundance in fragments with a trimeric neon core is found to increase with the increase in dopant size. Most of the fragmentation is achieved within 10 ps and the only subsequent dynamical process is the relaxation of hot intermediate He(q)Ne(2) (+) species to Ne(2) (+) by helium atom evaporation. The dependence of the ionic fragment distribution on the parent ion electronic state reached by ionization is also investigated. It reveals that He(q)Ne(+) fragments are produced only from the highest electronic state, whereas He(q)Ne(2) (+) fragments originate from all the electronic states. Surprisingly, the highest electronic states also lead to fragments that still contain the original ionic dopant species. A mechanism is conjectured to explain this fragmentation inhibition.     

Vibrational and rotational cooling of NO+ in collisions with He

A quantum mechanical investigation of the vibrational and rotational deactivation of NO(+) in collisions with He atoms in the cold and ultracold regime is presented. Ab initio potential energy calculations are carried out at BCCD(T) level and a new global 3D potential energy surface (PES) is obtained by fitting ab initio points within the reproducing kernel Hilbert space method. As a first test of this PES the bound state energies of the (3)He-NO(+) and (4)He-NO(+) complexes are calculated and compared to previous rigid rotor calculations. The efficiency of the vibrational and the rotational cooling of this molecular ion using a buffer gas of helium is then investigated by performing close coupling scattering calculations for collision energy ranging from 10(-6) to 2000 cm(-1). The calculations are performed for the two isotopes (3)He and (4)He and the results are compared to the available experimental data.     

Correlation holes for the helium dimer

We have investigated the radial electron pair probability distributions (REPPDs) of the helium dimer within the Piris natural orbital functional (PNOF) theory. The analytical formulas to evaluate intracule densities, Fermi, Coulomb, and total correlation holes using our reconstruction functional PNOF-2 [J. Chem. Phys. 126, 214103 (2007)] are derived. The L?wdin's Coulomb holes from PNOF-2 and full configuration interaction calculations are analyzed showing a very similar behavior. New definitions of the Coulomb and Fermi holes based on the cumulant expansion of the two-particle reduced density matrix are presented. The holes are defined in terms of the exact one-particle reduced density matrix and the two-particle cumulant without any reference to the Hartree-Fock state. Through these definitions, we analyze separately the contribution of each component to the total REPPD at several values of the internuclear distance. A straight connection between the Coulomb hole and dispersion interactions is observed.     

Pair potential for helium from symmetry-adapted perturbation theory calculations and from supermolecular data

Symmetry-adapted perturbation theory (SAPT) was applied to the helium dimer for interatomic separations R from 3 to 12 bohrs. The first-order interaction energy and the bulk of the second-order contribution were obtained using Gaussian geminal basis sets and are converged to about 0.1 mK near the minimum and for larger R. The remaining second-order contributions available in the SAPT suite of codes were computed using very large orbital basis sets, up to septuple-zeta quality, augmented by diffuse and midbond functions. The accuracy reached at this level was better than 1 mK in the same region. All the remaining components of the interaction energy were computed using the full configuration interaction method in bases up to sextuple-zeta quality. The latter components, although contributing only 1% near the minimum, have the largest uncertainty of about 10 mK in this region. The total interaction energy at R=5.6 bohrs is -11.000+/-0.011 K. For R< or =6.5 bohrs, the supermolecular (SM) interaction energies computed by us recently turned out to be slightly more accurate. Therefore, we have combined the SM results for R< or =6.5 bohrs with the SAPT results from 7.0 to 12 bohrs to fit analytic functions for the potential and for its error bars. The potential fit uses the best available van der Waals constants C(6) through C(16), including C(11), C(13), and C(15), and is believed to be the best current representation of the Born-Oppenheimer (BO) potential for helium. Using these fits, we found that the BO potential for the helium dimer exhibits the well depth D(e)=11.006+/-0.004 K, the equilibrium distance R(e)=5.608+/-0.012 bohrs, and supports one bound state for (4)He(2) with the dissociation energy D(0)=1.73+/-0.04 mK, and the average interatomic separation R=45.6+/-0.5 A.     

Cs atoms on helium nanodroplets and the immersion of Cs+ into the nanodroplet

We report the non-desorption of cesium (Cs) atoms on the surface of helium nanodroplets (He(N)) in their 6(2)P(1/2) ((2)Π(1/2)) state upon photo-excitation as well as the immersion of Cs(+) into the He(N) upon photo-ionization via the 6(2)P(1/2) ((2)Π(1/2)) state. Cesium atoms on the surface of helium nanodroplets are excited with a laser to the 6(2)P states. We compare laser-induced fluorescence (LIF) spectra with a desorption-sensitive method (Langmuir-Taylor detection) for different excitation energies. Dispersed fluorescence spectra show a broadening of the emission spectrum only when Cs-He(N) is excited with photon energies close to the atomic D(1)-line, which implies an attractive character of the excited state system (Cs?-He(N)) potential energy curve. The experimental data are compared with a calculation of the potential energy curves of the Cs atom as a function of its distance R from the center of the He(N) in a pseudo-diatomic model. Calculated Franck-Condon factors for emission from the 6(2)P(1/2) ((2)Π(1/2)) to the 6(2)S(1/2) ((2)Σ(1/2)) state help to explain the experimental data. The stability of the Cs?-He(N) system allows to form Cs(+) snowballs in the He(N), where we use the non-desorbing 6(2)P(1/2) ((2)Π(1/2)) state as a springboard for ionization in a two-step ionization scheme. Subsequent immersion of positively charged Cs ions is observed in time-of-flight mass spectra, where masses up to several thousand amu were monitored. Only ionization via the 6(2)P(1/2) ((2)Π(1/2)) state gives rise to a very high yield of immersed Cs(+) in contrast to an ionization scheme via the 6(2)P(3/2) ((2)Π(3/2)) state. When resonant two-photon ionization is applied to cesium dimers on He droplets, Cs(2) (+)-He(N) aggregates are observed in time-of-flight mass spectra.     

Dissociation energy of the water dimer from quantum Monte Carlo calculations

We report a study of the electronic dissociation energy of the water dimer using quantum Monte Carlo techniques. We have performed variational quantum Monte Carlo and diffusion quantum Monte Carlo (DMC) calculations of the electronic ground state of the water monomer and dimer using all-electron and pseudopotential approaches. We have used Slater-Jastrow trial wave functions with B3LYP type single-particle orbitals, into which we have incorporated backflow correlations. When backflow correlations are introduced, the total energy of the water monomer decreases by about 4-5 mhartree, yielding a DMC energy of -76.428 30(5) hartree, which is only 10 mhartree above the experimental value. In our pseudopotential DMC calculations, we have compared the total energies of the water monomer and dimer obtained using the locality approximation with those from the variational scheme recently proposed by Casula [Phys. Rev. B 74, 161102(R) (2006)]. The time step errors in the Casula scheme are larger, and the extrapolation of the energy to zero time step always lies above the result obtained with the locality approximation. However, the errors cancel when energy differences are taken, yielding electronic dissociation energies within error bars of each other. The dissociation energies obtained in our various all-electron and pseudopotential calculations range between 5.03(7) and 5.47(9) kcalmol and are in good agreement with experiment. Our calculations give monomer dipole moments which range between 1.897(2) and 1.909(4) D and dimer dipole moments which range between 2.628(6) and 2.672(5) D.     

Hydrogen bonds in 1,4-dioxane/ammonia binary clusters

With synchrotron radiation, we have studied the photoionization and dissociation of 1,4-dioxane/ammonia clusters in a supersonic expansion. The observed major product ions are the 1,4-dioxane cation M(+) and protonated cluster ions M(NH(3))(n)H(+) (where M=1,4-dioxane), and the intensities of the unprotonated cluster ions M(NH(3))(n) (+) are much lower. Fully optimized geometries and energies of the neutral cluster M(NH(3))(2) and related cluster ions have been obtained using the ab initio molecular orbital method and density functional theory. The potential energy surface of the excited state of M(NH(3))(2) (+) was also calculated. With these results, the mechanisms of different photoionization-dissociation channels have been suggested. The most probable channel is electron ejection from the highest occupied molecular orbital, followed by the dissociation into M(+) and (NH(3))(2). For another main channel, after removing an electron from the second highest occupied molecular orbital, the intracluster proton transfer process takes place to form the stable unprotonated cluster ion M(NH(3))H(+)-NH(2), which usually leads to the dissociated protonated cluster ion M(NH(3))H(+) and a radical NH(2).     

Ionization of the ammonia dimer: Proton transfer in the ionic state

The adiabatic (I_a) and the vertical (I_v) ionization energy of the ammonia dimer as well as the appearance energies of NH⁺₄ and (NH₃)⁺₂ have been obtained from ab initio potential energy surface calculations. The comparison with the experimental result leads to the following conclusions: (1) the real I_a of (NH₃)₂ is not observed in the experiment due to the unfavorable Franck—Condon factor of the equilibrium dimer cation which has the structure of H₃NH⁺ …NH₂. (2) The value of I_a is estimated to be about 8.6 eV, almost 1 eV lower than the AE of (NH₃)⁺₂. (3) Consequently, the two kinds of dissociation energy of the dimer cation into NH₃ + NH⁺₃ and NH⁺₄ + NH₂, respectively, are about 1 eV larger than the “experimental” estimate from the apparent AE values. The present interpretation solves the question of the large discrepancy in the dissociation energies of the ammonia dimer cation.     

Collision Energetics in a Tandem Time-of-Flight (TOF/TOF) Mass Spectrometer with a Curved-Field Reflectron

Collisions of fullerene ions (C(60) (+)) with helium and neon were carried out over a range of laboratory energies (3-20 keV) on a unique tandem time-of-flight (TOF/TOF) mass spectrometer equipped with a curved-field reflectron (CFR). The CFR enables focusing of product ions over a wide kinetic energy range. Thus, ions extracted from a laser desorption/ionization (LDI) source are not decelerated prior to collision, and collision energies in the laboratory frame are determined by the source extraction voltages. Comparison of product ion mass spectra obtained following collisions with inert gases show a time (and apparent mass) shift for product ions relative to those observed in spectra obtained by metastable dissociation (unimolecular decay), consistent with impulse collision models, in which interactions of helium with fullerene in the high energy range are primarily with a single carbon atom. In addition, within a narrow range of kinetic energies an additional peak corresponding to the capture of helium is observed for fragment ions C(50) (+), C(52) (+), C(54) (+), C(56) (+) and C(58) (+).     

Dissociative recombination of the weakly bound NO-dimer cation: cross sections and three-body dynamics

Dissociative recombination (DR) of the dimer ion (NO)(2) (+) has been studied at the heavy-ion storage ring CRYRING at the Manne Siegbahn Laboratory, Stockholm. The experiments were aimed at determining details on the strongly enhanced thermal rate coefficient for the dimer, interpreting the dissociation dynamics of the dimer ion, and studying the degree of similarity to the behavior in the monomer. The DR rate reveals that the very large efficiency of the dimer rate with respect to the monomer is limited to electron energies below 0.2 eV. The fragmentation products reveal that the breakup into the three-body channel NO+O+N dominates with a probability of 0.69+/-0.02. The second most important channel yields NO+NO fragments with a probability of 0.23+/-0.03. Furthermore, the dominant three-body breakup yields electronic and vibrational ground-state products, NO(upsilon=0)+N((4)S)+O((3)P), in about 45% of the cases. The internal product-state distribution of the NO fragment shows a similarity with the product-state distribution as predicted by the Franck-Condon overlap between a NO moiety of the dimer ion and a free NO. The dissociation dynamics seem to be independent of the NO internal energy. Finally, the dissociation dynamics reveal a correlation between the kinetic energy of the NO fragment and the degree of conservation of linear momentum between the O and N product atoms. The observations support a mechanism in which the recoil takes place along one of the NO bonds in the dimer.     

Communication: rigorous calculation of dissociation energies (D0) of the water trimer, (H2O)3 and (D2O)3

Using a recent, full-dimensional, ab initio potential energy surface [Y. Wang, X. Huang, B. C. Shepler, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 134, 094509 (2011)] together with rigorous diffusion Monte Carlo calculations of the zero-point energy of the water trimer, we report dissociation energies, D(0), to form one monomer plus the water dimer and three monomers. The calculations make use of essentially exact zero-point energies for the water trimer, dimer, and monomer, and benchmark values of the electronic dissociation energies, D(e), of the water trimer [J. A. Anderson, K. Crager, L. Fedoroff, and G. S. Tschumper, J. Chem. Phys. 121, 11023 (2004)]. The D(0) results are 3855 and 2726 cm(-1) for the 3H(2)O and H(2)O + (H(2)O)(2) dissociation channels, respectively, and 4206 and 2947 cm(-1) for 3D(2)O and D(2)O + (D(2)O)(2) dissociation channels, respectively. The results have estimated uncertainties of 20 and 30 cm(-1) for the monomer plus dimer and three monomer of dissociation channels, respectively.     

Thermochemistry of N-N containing ions: a threshold photoelectron photoion coincidence spectroscopy study of ionized methyl- and tetramethylhydrazine

The unimolecular dissociation reactions of the methylhydrazine (MH) and tetramethylhydrazine (TMH) radical cations have been investigated using tandem mass spectrometry and threshold photoelectron photoion coincidence spectroscopy in the photon energy ranges 9.60-31.95 eV (for the MH ion) and 7.74-29.94 eV (for the TMH ion). Methylhydrazine ions (CH3NHNH2(+*)) have three low-energy dissociation channels: hydrogen atom loss to form CH2NHNH2(+) (m/z 45), loss of a methyl radical to form NHNH2(+) (m/z 31), and loss of methane to form the fragment ion m/z 30, N2H2(+*). Tetramethylhydrazine ions only exhibit two dissociation reactions near threshold: that of methyl radical loss to form (CH3)2NNCH3(+) (m/z 73) and of methane loss to form the fragment ion m/z 72 with the empirical formula C3H8N2(+*). The experimental breakdown curves were modeled with Rice-Ramsperger-Kassel-Marcus theory, and it was found that, particularly for methyl radical loss, variational transition state theory was needed to obtain satisfactory fits to the data. The 0 K enthalpies of formation (delta(f)H0) for all fragment ions (m/z 73, m/z 72, m/z 45, m/z 31, and m/z 30) have been determined from the 0 K activation energies (E0) obtained from the fitting procedure: delta(f)H0[(CH3)2NNCH3(+)] = 833 +/- 5 kJ mol(-1), delta(f)H0 [C3H8N2(+*)] = 1064 +/- 5 kJ mol(-1), delta(f)H0[CH2NHNH2(+)] = 862 +/- 5 kJ mol(-1), delta(f)H0[NHNH2(+)] = 959 +/- 5 kJ mol(-1), and delta(f)H0[N2H2(+*)] = 1155 +/- 5 kJ mol(-1). The breakdown curves have been measured from threshold up to h nu approximately 32 eV for both hydrazine ions. As the photon energy increases, other dissociation products are observed and their appearance energies are reported.     

Calculation of intermolecular interactions in the benzene dimer using coupled-cluster and local electron correlation methods

Potential energy curves for the parallel-displaced, T-shaped and sandwich structures of the benzene dimer are computed with density fitted local second-order M?ller-Plesset perturbation theory (DF-LMP2) as well as with the spin-component scaled (SCS) variant of DF-LMP2. While DF-LMP2 strongly overestimates the dispersion interaction, in common with canonical MP2, the DF-SCS-LMP2 interaction energies are in excellent agreement with the best available literature values along the entire potential energy curves. The DF-SCS-LMP2 dissociation energies for the three structures are also compared with new complete basis set estimates of the interaction energies obtained from accurate coupled cluster (CCSD(T)) and DF-SCS-MP2 calculations. Since LMP2 is essentially free of basis set superposition errors, counterpoise corrections are not required. As a result, DF-SCS-LMP2 is computationally inexpensive and represents an attractive method for the study of larger pi-stacked systems such as truncated sections of DNA.     

Quantum Monte Carlo calculations of the dissociation energies of three-electron hemibonded radical cationic dimers

We report variational and diffusion quantum Monte Carlo (VMC and DMC) calculations of the dissociation energies of the three-electron hemibonded radical cationic dimers of He, NH3, H2O, HF, and Ne. These systems are particularly difficult for standard density-functional methods such as the local-density approximation and the generalized gradient approximation. We have performed both all-electron (AE) and pseudopotential (PP) calculations using Slater-Jastrow wave functions with Hartree-Fock single-particle orbitals. Our results are in good agreement with coupled-cluster CCSD(T) calculations. We have also studied the relative stability of the hemibonded and hydrogen-bonded water radical dimer isomers. Our calculations indicate that the latter isomer is more stable, in agreement with post-Hartree-Fock methods. The excellent agreement between our AE and PP results demonstrates the high quality of the PPs used within our VMC and DMC calculations.     

Actinide sulfides in the gas phase: experimental and theoretical studies of the thermochemistry of AnS (An = Ac, Th, Pa, U, Np, Pu, Am and Cm)

The gas-phase thermochemistry of actinide monosulfides, AnS, was investigated experimentally and theoretically. Fourier transform ion cyclotron resonance mass spectrometry was employed to study the reactivity of An(+) and AnO(+) (An = Th, Pa, U, Np, Pu, Am and Cm) with CS(2) and COS, as well as the reactivity of the produced AnS(+) with oxidants (COS, CO(2), CH(2)O and NO). From these experiments, An(+)-S bond dissociation energies could be bracketed. Density functional theory studies of the energetics of neutral and monocationic AnS (An = Ac, Th, Pa, U, Np, Pu, Am and Cm) provided values for bond dissociation energies and ionization energies; the computed energetics of neutral and monocationic AnO were also obtained for comparison. The theoretical data, together with comparisons with known An(+)-O bond dissociation energies and M(+)-S and M(+)-O dissociation energies for the early transition metals, allowed for the refining of the An(+)-S bond dissociation energy ranges obtained from experiment. Examination of the reactivity of AnS(+) with dienes, coupled to comparisons with reactivities of the AnO(+) analogues, systematic considerations and the theoretical results, allowed for the estimation of the ionization energies of the AnS; the bond dissociation energies of neutral AnS were consequently derived. Estimates for the case of AcS were also made, based on correlations of the data for the other An and the electronic energetics of neutral and ionic An. The nature of the bonding in the elementary molecular actinide chalcogenides (oxides and sulfides) is discussed, based on both the experimental data and the computed electronic structures. DFT calculations of ionization energies for the actinide atoms and the diatomic sulfides and oxides are relatively reliable, but the calculation of bond dissociation energies is not uniformly satisfactory, either with DFT or CCSD(T). A key conclusion from both the experimental and theoretical results is that the 5f electrons do not substantially participate in actinide-sulfur bonding. We emphasize that actinides form strikingly strong bonds with both oxygen and sulfur.     

Characterization of Tyrosine Nitration and Cysteine Nitrosylation Modifications by Metastable Atom-Activation Dissociation Mass Spectrometry

The fragmentation behavior of nitrated and S-nitrosylated peptides were studied using collision induced dissociation (CID) and metastable atom-activated dissociation mass spectrometry (MAD-MS). Various charge states, such as 1+, 2+, 3+, 2–, of modified and unmodified peptides were exposed to a beam of high kinetic energy helium (He) metastable atoms resulting in extensive backbone fragmentation with significant retention of the post-translation modifications (PTMs). Whereas the high electron affinity of the nitrotyrosine moiety quenches radical chemistry and fragmentation in electron capture dissociation (ECD) and electron transfer dissociation (ETD), MAD does produce numerous backbone cleavages in the vicinity of the modification. Fragment ions of nitrosylated cysteine modifications typically exhibit more abundant neutral losses than nitrated tyrosine modifications because of the extremely labile nature of the nitrosylated cysteine residues. However, compared with CID, MAD produced between 66% and 86% more fragment ions, which preserved the labile –NO modification. MAD was also able to differentiate I/L residues in the modified peptides. MAD is able to induce radical ion chemistry even in the presence of strong radical traps and therefore offers unique advantages to ECD, ETD, and CID for determination of PTMs such as nitrated and S-nitrosylated peptides.