Ab initio studies on (H2O)14 - clusters: existence of surface and interior-bound extra electrons

A recent paper by Turi et al. [Science 309, 914 (2005)] suggests that the anionic water clusters smaller than (H2O)(45) (-) (at a low temperature) will only have surface-bound extra electrons and no internally bound electrons. Accordingly, (H2O)(14) (-) cluster isomers should only have surface-bound extra electrons. The ab initio results presented here, however, suggest that the (H2O)(14) (-) cluster isomers can have two distinct types of isomers with almost the same energy. The one type of isomer (type 1) has all the non-H-bonding H atoms (NHB H) directed outward and surface-bound extra electron while the other type (type 2) has a number of NHB H atoms directed toward cavity and has an interior-bound electron, and thus, contradicts the earlier quantum simulation results of Turi et al.     

Theoretical studies on photoelectron and IR spectral properties of Br2.-(H2O)n clusters

We report vertical detachment energy (VDE) and IR spectra of Br2.-.(H2O)n clusters (n=1-8) based on first principles electronic structure calculations. Cluster structures and IR spectra are calculated at Becke's half-and-half hybrid exchange-correlation functional (BHHLYP) with a triple split valence basis function, 6-311++G(d,p). VDE for the hydrated clusters is calculated based on second order Moller-Plesset perturbation (MP2) theory with the same set of basis function. On full geometry optimization, it is observed that conformers having interwater hydrogen bonding among solvent water molecules are more stable than the structures having double or single hydrogen bonded structures between the anionic solute, Br2.-, and solvent water molecules. Moreover, a conformer having cyclic interwater hydrogen bonded network is predicted to be more stable for each size hydrated cluster. It is also noticed that up to four solvent H2O units can reside around the solute in a cyclic interwater hydrogen bonded network. The excess electron in these hydrated clusters is localized over the solute atoms. Weighted average VDE is calculated for each size (n) cluster based on statistical population of the conformers at 150 K. A linear relationship is obtained for VDE versus (n+3)(-1/3) and bulk VDE of Br2.- aqueous solution is calculated as 10.01 eV at MP2 level of theory. BHHLYP density functional is seen to make a systematic overestimation in VDE values by approximately 0.5 eV compared to MP2 data in all the hydrated clusters. It is observed that hydration increases VDE of bromine dimer anion system by approximately 6.4 eV. Calculated IR spectra show that the formation of Br2.--water clusters induces large shifts from the normal O-H stretching bands of isolated water keeping bending modes rather insensitive. Hydrated clusters, Br2.-.(H2O)n, show characteristic sharp features of O-H stretching bands of water in the small size clusters.     

Selecting high stability water clusters by examining locations of free OH bonds: Application in the study of clusters

Computational studies on water clusters can be quite challenging, especially when an irregular cage with non-equivalent oxygen sites are considered which may yield a large number of starting geometries that differ in relative positions of non-H-bonding H (NHB H, free OH) atoms. A detailed study on water octamers suggests that the fewest number of NHB H atoms on neighboring oxygen sites yields the most stable neutral isomer followed by those with increasing number of NHB H atoms on adjacent sites. The least stable cluster has all the NHB H atoms around a ring. By considering a regular cage structure and minimum number of NHB H atoms on adjacent sites, one can readily identify a limited number of starting geometries that are optimized to highly stable isomers. This method has been verified in the identification of the most stable isomer of (H₂O)₈ cubic cage and (H₂O)₂₀ dodecahedral cage. The same method has been applied in the study of cluster isomers.     

Calculation of electron detachment energies for water cluster anions: an appraisal of electronic structure methods, with application to (H2O)20- AND (H2O)24-

We present benchmark calculations of vertical electron detachment energies (VDEs) for various conformers of (H2O)n-, using both wave function and density functional methods, in sequences of increasingly diffuse Gaussian basis sets. For small clusters (n < or = 6), a systematic examination of VDE convergence reveals that it is possible to converge this quantity to within approximately 0.01 eV of the complete-basis limit, using a highly diffuse but otherwise economical Pople-style basis set of double-zeta quality, with 28 atom-centered basis functions per water molecule. Floating-center basis functions can be useful but are not required to obtain accurate VDEs. Second-order M?ller-Plesset perturbation (MP2) theory suffices to obtain VDEs that are within 0.05 eV of the results from both experiment and coupled-cluster theory, and which always err toward underbinding the extra electron. In contrast to these consistent predictions, VDEs calculated using density functional theory (DFT) vary widely, according to the fraction of Hartree-Fock exchange in a given functional. Common functionals such as BLYP and B3LYP overestimate the VDE by 0.2-0.5 eV, whereas a variant of Becke's "half and half" functional is much closer to coupled-cluster predictions. Exploratory calculations for (H2O)20- and (H2O)24- cast considerable doubt on earlier calculations that were used to assign the photoelectron spectra of these species to particular cluster isomers.     

Mass spectrometry and photoelectron spectroscopy of tetracene cluster anions, (tetracene)n- (n = 1-100): evidence for the highly localized nature of polarization in a cluster analogue of oligoacene crystals

Photoelectron spectroscopy of tetracene cluster anions, (tetracene)n- (n = 1-100), reveals the coexistence of two types of isomers, designated as isomers I and II-1 (n = 10-50) or isomers I and II-2 (n > 60), in a wide size range. The vertical detachment energies (VDEs) of isomer I increase persistently due to polarization and structural relaxation effects, where a monomeric anion core is encompassed with geometrically reorganized neutral molecules. Conversely, a characteristic ion distribution in the mass spectrum of (tetracene)n-ensues from the two-dimensional (2D) herringbone-type ordering of isomer II-1, whose VDEs remain constant at 1.80 eV for n >/= 14. Also, isomer II-2, presumably adopting multilayered structural motifs, exhibits invariable VDEs of 2.0 eV, a manifestation of significant charge screening effects in these isomers. The invariable nature of the VDEs of isomers II-1 and II-2 unambiguously demonstrates a largely localized nature of polarization induced by the excess charge residing in microscopic crystal-like environments. Surprisingly, only 14 tetracene molecules within a 2D herringbone-type layer including an excess charge can provide the charge stabilization energy corresponding to approximately 80% of that of the crystal, and the rest of the energy is provided by polarization of neutral molecules in adjacent layers.     

Photoelectron spectroscopy of the [glycine x (H2O)(1,2)]- clusters: sequential hydration shifts and observation of isomers

The electron binding energies of the small hydrated amino acid anions, [glycine x (H2O)(1,2)]-, are determined using photoelectron spectroscopy. The vertical electron detachment energies (VDEs) are found to increase by approximately 0.12 eV with each additional water molecule such that the higher electron binding isomer of the dihydrate is rather robust, with a VDE value of 0.33 eV. A weak binding isomer of the dihydrate is also recovered, however, with a VDE value (0.14 eV) lower than that of the monohydrate. Unlike the situation in the smaller (n < or = 13) water cluster anions, the [Gly x (H2O)(n > or = 6)]- clusters are observed to photodissociate via water monomer evaporation upon photoexcitation in the O-H stretching region. We discuss this observation in the context of the mechanism responsible for the previously observed [S. Xu, M. Nilles, and K. H. Bowen, Jr., J. Chem. Phys. 119, 10696 (2003)] sudden onset in the cluster formation at [Gly x (H2O)5]-.     

Comprehensive photoelectron spectroscopic study of anionic clusters of anthracene and its alkyl derivatives: electronic structures bridging molecules to bulk

The evolution of the electronic structure of molecular aggregates is investigated using anion photoelectron (PE) spectroscopy for anionic clusters of anthracene (Ac) and its alkyl derivatives: 1-methylanthracene (1MA), 2-methylanthracene (2MA), 9-methylanthracene (9MA), 9,10-dimethylanthracene (DMA), and 2-tert-butylanthracene (2TBA). For their monomer anions (n=1), electron affinities are confined to the range from 0.47 to 0.59 eV and are well reproduced by density functional theory calculations, showing the isoelectronic character of these molecules. For cluster anions (n=2-100) of Ac and 2MA, two types of isomers I and II coexist over a wide size range: isomers I and II-1 (4< or =n<30) or isomers I and II-2 (n> or = approximately 40 for Ac and n> or = approximately 55 for 2MA). However, for the other alkyl-substituted Ac cluster anions (i.e., 1MA, 9MA, DMA, and 2TBA), only isomer I is exclusively formed, and neither isomer II-1 nor II-2 is observed. The vertical detachment energies (VDEs) of isomer I in all the anionic clusters depend almost linearly on n(-1/3). In contrast, the VDEs of isomers II-1 (n> or =14) and II-2 (n=40-100), appeared only in Ac and 2MA cluster anions, remain constant with n and are approximately 0.5 eV lower than those of isomer I. The PE spectra revealed the characteristics of each isomer: isomer I possesses a monomeric anion core that is gradually embedded into the interior of the cluster with increasing n. On the other hand, isomers II-1 and II-2 possess a multimeric (perhaps tetrameric) anion core, but they differ in the number of layers from which they are made up; monolayer (isomer II-1) and multilayers (isomer II-2) of a two-dimensionally ordered, finite herringbone-type structure, in which electron attachment produces only little geometrical rearrangement. Moreover, the agreement of the constant VDEs of isomer II-2 with the bulk data demonstrates the largely localized nature of the electronic polarization around the excess charge in a crystal-like environment, where about 50 molecules provide a charge stabilization energy comparable to the bulk.     

Electronic structure and reactivity of a biradical cluster: Sc3O6(-)

The Sc(3)O(6)(-) cluster anions were produced by laser ablation and studied by reaction with n-butane in a fast flow reactor and by photoelectron spectroscopy. The reactivity experiments indicated that one Sc(3)O(6)(-) cluster can activate two n-butane molecules consecutively with rate constants on the order of 10(-10) cm(3) molecule(-1) s(-1) under near room-temperature conditions, suggesting that the even-electron system Sc(3)O(6)(-) has a highly reactive electronic structure. The photoelectron spectroscopy determined a high vertical detachment energy (VDE) of 5.63 ± 0.08 eV for the Sc(3)O(6)(-) cluster. Density functional computations indicated that the lowest energy isomer of Sc(3)O(6)(-) is an oxygen-centered biradical with a high VDE and is highly reactive toward n-butane, which is in good agreement with the experiments. The Sc(3)O(6)(-) cluster may serve as an ideal model system to provide insight into the real-life chemistry involved with the coupled O(-)˙···O(-)˙ dimers over the surfaces of metal oxide catalysts.     

The structures and electronic states of zinc-water clusters Zn(n)(H2O)(m) (n = 1-32 and m = 1-3)

Ab initio and Density Functional Theory (DFT) calculations have been carried out for zinc-water clusters Zn(n)-(H2O)(m) (n = 1-32 and m = 1-3, where n and m are the numbers of zinc atoms and water molecules, respectively) to elucidate the structure and electronic states of the clusters and the interaction of zinc cluster with water molecules. The binding energies of H2O to zinc clusters were small at n = 2-3 (2.3-4.2 kcal mol(-1)), whereas the energy increased significantly in n = 4 (9.0 kcal mol(-1)). Also, the binding nature of H2O was changed at n = 4. The cluster size dependency of the binding energy of H2O accorded well with that of the natural population of electrons in the 4p orbital of the zinc atom. In the larger clusters (n > 20), it was found that the zinc atoms in surface regions of the zinc cluster have a positive charge, whereas those in the interior region have a negative charge with the large electron population in the 4p orbital. The interaction of H2O with the zinc clusters were discussed on the basis of the theoretical results.     

Clusters of hydrated methane sulfonic acid CH3SO3H.(H2O)n (n = 1-5): a theoretical study

Ab initio and density functional methods have been used to examine the structures and energetics of the hydrated clusters of methane sulfonic acid (MSA), CH3SO3H.(H2O)n (n = 1-5). For small clusters with one or two water molecules, the most stable clusters have strong cyclic hydrogen bonds between the proton of OH group in MSA and the water molecules. With three or more water molecules, the proton transfer from MSA to water becomes possible, forming ion-pair structures between CH3SO3- and H3O+ moieties. For MSA.(H2O)3, the energy difference between the most stable ion pair and neutral structures are less than 1 kJ/mol, thus coexistence of neutral and ion-pair isomers are expected. For larger clusters with four and five water molecules, the ion-pair isomers are more stable (>10 kJ/mol) than the neutral ones; thus, proton transfer takes place. The ion-pair clusters can have direct hydrogen bond between CH3SO3- and H3O+ or indirect one through water molecule. For MSA.(H2O)5, the energy difference between ion pairs with direct and indirect hydrogen bonds are less than 1 kJ/mol; namely, the charge separation and acid ionization is energetically possible. The calculated IR spectra of stable isomers of MSA.(H2O)n clusters clearly demonstrate the significant red shift of OH stretching of MSA and hydrogen-bonded OH stretching of water molecules as the size of cluster increases.     

Probing isomer interconversion in anionic water clusters using an Ar-mediated pump-probe approach: Combining vibrational predissociation and velocity-map photoelectron imaging spectroscopies

We present the first results from an experiment designed to explore barriers for interconversion between isomers of cluster anions using an Ar-cluster mediated pump-probe technique. In this approach, anions are generated with many Ar atoms attached, and one of the isomers present is selectively excited by tuning an infrared laser to one of the isomer's characteristic vibrational resonances. The excited cluster is then cooled by evaporation of Ar atoms, and the isomer distribution in the lighter daughter ions is measured after secondary mass selection by recording their photoelectron spectra using velocity-map imaging. We apply the method to the water hexamer anion, (H(2)O)(6) (-), which is known to occur in two isomeric forms with different electron-binding energies. We find that conversion of the high-binding (type I) form to the low-binding (type II) isomer is not efficiently driven in (H(2)O)(6) (-) with excitation energies in the 0.4 eV range even though it is possible to create both isomers in abundance in the ion source. This observation is discussed in the context of the competition between isomerization and electron autodetachment, which depends on the relative positions of the neutral and ionic potential surfaces along the isomerization pathway. Application of the method to the more complex heptamer ion, however, does reveal that interconversion is available among the highest binding isomer classes (I and I(')).     

MRCI investigation of different isomers of Ni2O2H2(+)

The geometrical and electronic structures of different isomers of Ni(2)O(2)H(2)(+) are investigated by multireference configuration interaction (MRCI) calculations using natural atomic orbital basis sets. The lowest-lying isomer, Ni(2)(OH)(2)(+), has a rhombic shape with two OH groups bridging the Ni atoms. The next isomer in energetic order with a relative energy of 0.29 eV consists of a linear NiONi(OH(2))(+) chain. Other structures with a rhombic shape, (NiH)(2)O(2)(+), with H bound to the Ni atoms have considerably higher energies, above 4 eV. Especially the low-lying isomers are characterised by a large number of low-lying electronic terms. The product Ni(2)O(2)H(2)(+) of the reaction of Ni(2)O(2)(+) with small alkanes is likely to have the rhombic Ni(2)(OH)(2)(+) structure. The reaction energy of the reaction Ni(2)O(2)(+) + H(2)→ Ni(2)(OH)(2)(+) is estimated to be about -3.5 eV.     

Observation of Au2H- impurity in pure gold clusters and implications for the anomalous Au-Au distances in gold nanowires

Au(2)H(-) was recognized and confirmed as a minor contamination to typical photoelectron spectra of Au(2) (-), produced by laser vaporization of a pure Au target using an ultrahigh purity helium carrier gas. The hydrogen source was shown to be from trace H impurities present in the bulk gold target. Carefully designed experiments using H(2)- and D(2)-seeded helium carrier gas were used to study the electronic structure of Au(2)H(-) and Au(2)D(-) using photoelectron spectroscopy and density functional calculations. Well-resolved photoelectron spectra with vibrational resolution were obtained for Au(2)H(-) and Au(2)D(-). Two isomers were observed both experimentally and theoretically. The ground state of Au(2)H(-) turned out to be linear with a terminal H atom [Au-Au-H](-) ((1)A(1),C(infinity v)), whereas a linear [Au-H-Au](-) ((1)A(1),D(infinity h)) structure with a bridging H atom was found to be a minor isomer 0.6 eV higher in energy. Calculated electron detachment energies for both isomers agree well with the experimental spectra, confirming their existence in the cluster beam. The observation and confirmation of H impurity in pure gold clusters and the 3.44 A Au-Au distance in the [Au-H-Au](-) isomer presented in the current work provide indirect experimental evidence that the anomalous 3.6 A Au-Au distances observed in gold nanowires is due to an "invisible" hydrogen impurity atom.     

Theoretical study of negatively charged Fe(-)-(H2O)(n ≤ 6) clusters

Interactions of a singly negatively charged iron atom with water molecules, Fe(-)-(H(2)O)(n≤6), in the gas phase were studied by means of density functional theory. All-electron calculations were performed using the B3LYP functional and the 6-311++G(2d,2p) basis set for the Fe, O, and H atoms. In the lowest total energy states of Fe(-)-(H(2)O)(n), the metal-hydrogen bonding is stronger than the metal-oxygen one, producing low-symmetry structures because the water molecules are directly attached to the metal by basically one of their hydrogen atoms, whereas the other ones are involved in a network of hydrogen bonds, which together with the Fe(δ-)-H(δ+) bonding accounts for the nascent hydration of the Fe(-) anion. For Fe(-)-(H(2)O)(3≤n), three-, four-, five-, and six-membered rings of water molecules are bonded to the metal, which is located at the surface of the cluster in such a way as to reduce the repulsion with the oxygen atoms. Nevertheless, internal isomers appear also, lying less than 3 or 5 kcal/mol for n = 2-3 or n = 4-6. These results are in contrast with those of classical TM(+)-(H(2)O)(n) complexes, where the direct TM(+)-O bonding usually produces high symmetry structures with the metal defining the center of the complex. They show also that the Fe(-) anions, as the TM(+) ions, have great capability for the adsorption of water molecules, forming Fe(-)-(H(2)O)(n) structures stabilized by Fe(δ-)-H(δ+) and H-bond interactions.     

Theoretical characterization of four distinct isomer types in hydrated-electron clusters, and proposed assignments for photoelectron spectra of water cluster anions

Water cluster anions, (H(2)O)(N)(-), are examined using mixed quantum/classical molecular dynamics based on a one-electron pseudopotential model that incorporates many-body polarization and predicts vertical electron detachment energies (VDEs) with an accuracy of ~0.1 eV. By varying the initial conditions under which the clusters are formed, we are able to identify four distinct isomer types that exhibit different size-dependent VDEs. On the basis of a strong correlation between the electron's radius of gyration and its optical absorption maximum, and extrapolating to the bulk limit (N → ∞), our analysis supports the assignment of the "isomer Ib" data series, observed in photoelectron spectra of very cold clusters, as arising from cavity-bound (H(2)O)(N)(-) cluster isomers. The "isomer I" data reported in warmer experiments are assigned to surface-bound isomers in smaller clusters, transitioning to partially embedded isomers in larger clusters. The partially embedded isomers are characterized by a partially formed solvent cavity at the cluster surface, and they are spectroscopically quite similar to internalized cavity isomers. These assignments are consistent with various experimental data, and our theoretical characterization of these isomers sheds new light on a long-standing assignment problem.     

Vertical detachment energies of anionic thymidine: Microhydration effects

Density functional theory has been employed to investigate microhydration effects on the vertical detachment energy (VDE) of the thymidine anion by considering the various structures of its monohydrates. Structures were located using a random searching procedure. Among 14 distinct structures of the anionic thymidine monohydrate, the low-energy structures, in general, have the water molecule bound to the thymine base unit. The negative charge developed on the thymine moiety increases the strength of the intermolecular hydrogen bonding between the water and base units. The computed VDE values of the thymidine monohydrate anions are predicted to range from 0.67 to 1.60 eV and the lowest-energy structure has a VDE of 1.32 eV. The VDEs of the monohydrates of the thymidine anion, where the N(1)[Single Bond]H hydrogen of thymine has been replaced by a 2(')-deoxyribose ring, are greater by ～0.30?eV, compared to those of the monohydrates of the thymine anion. The results of the present study are in excellent agreement with the accompanying experimental results of Bowen and co-workers [J. Chem. Phys. 133, 144304 (2010)].     

Addition of water to Al5O4- determined by anion photoelectron spectroscopy and quantum chemical calculations

The anion photoelectron spectra of Al5O4- and Al5O5H2- are presented and interpreted within the context of quantum chemical calculations on these species. Experimentally, the electron affinities of these two molecules are determined to be 3.50(5) eV and 3.10(10) eV for the bare and hydrated cluster, respectively. The spectra show at least three electronic transitions crowded into a 1 eV energy window. Calculations on Al5O4- predict a highly symmetric near-planar structure with a singlet ground state. The neutral structure calculated to be most structurally similar to the ground state structure of the anion is predicted to lie 0.15 eV above the ground state structure of the neutral. The lowest energy neutral isomer does not have significant Franck-Condon overlap with the ground state of the anion. Dissociative addition of water to Al5O4- is energetically favored over physisorption. The ground state structure for the Al5O4- +H(2)O product forms when water adds to the central Al atom in Al5O4- with -H migration to one of the neighboring O atoms. Again, the ground state structures for the anion and neutral are very different, and the PE spectrum represents transitions to a higher-lying neutral structure from the ground state anion structure.     

烟酸二聚体的结构与性质

用密度泛函理论B3LYP方法选取6-311+G(d,p)基组对烟酸-烟酸复合物进行了量子化学计算研究, 通过在相同水平下的频率振动分析发现了势能面上存在7个极小值点, 其最稳定构型1对应一N…H—O型强氢键, 其结合能在消除基组重叠误差后为-48.3 kJ·mol-1. 通过自然键轨道(NBO)分析, 研究了电荷转移及轨道相互作用. 通过自洽反应场(SCRF)理论中的极化连续模型(PCM)在介电常数分别为1.0(真空)、2.247(苯)、10.36(二氯乙烷)、20.7(丙酮)、32.63(甲醇)、78.39(水)的不同溶剂环境下重新优化烟酸复合物势能面上最稳定构型1, 研究了溶剂对烟酸复合物几何构型及结合能的影响. 发现溶剂化作用增大了烟酸复合物分子间的结合能, 导致N…H距离减小. 当溶液介电常数在1.0-32.63范围时, 溶剂效应十分显著, 当介电常数大于32.63后, 溶剂化作用趋于稳定.     

Atomic and Molecular Chemisorption of Oxygen in WO4- Clusters

Density functional theory (DFT) calculations were performed to study the monotungsten-oxide WO4 cluster in the anionic and neutral charge states. The results show the two most stable WO4- isomers have C2v and D2d symmetries and both have the four oxygen atoms attached to the tungsten W monomer atomically. The WO4- species previously suggested with a molecular adsorption of di-oxygen is found to be a metastable isomer of WO4-, whose geometric, vibrational properties and electron affinities are in good agreement with the ultraviolet photoelectron spectroscopy (UPS) experimental results. The reason why this metastable isomer could be observed in the experiment is given by a molecule formation mechanism. The UPS spectrums compare well with the excitation spectrum computed by time-dependent DFT method.     

Identification of two distinct electron binding motifs in the anionic water clusters: a vibrational spectroscopic study of the (H2O)6- isomers

Photoelectron spectroscopy of the water cluster anions, (H2O)n-, has revealed that several isomeric forms are present for most sizes, and here, we use vibrational spectroscopy to address the structure of the (H2O)6- isomer that more weakly binds the extra electron. To overcome the severe line broadening that occurs in the OH stretching region of this isomer caused by fast electron autodetachment, we concentrate on the low-energy bending modes of the perdeutero isotopomer. Sharp spectroscopic signatures are recovered for two isomers using argon predissociation spectroscopy, and the resulting bands are heavily overlapped. To extract their independent contributions to the observed spectra, we exploit the substantial dependence of their relative populations on the number of attached argon atoms in the (D2O)6-.Ar(m) clusters, determined by photoelectron spectroscopy. The vibrational spectra of each isomer can then be isolated by spectral subtraction, which is implemented with a covariance mapping approach. The resulting band patterns establish that the more weakly binding isomer does not display the characteristic electron-binding motif common to the more strongly bound isomer class. Whereas the strongly binding isomer features a single water molecule pointing toward the excess electron cloud with both of its hydrogen atoms, the spectrum of the more weakly binding isomer suggests a structure where the electron is bound by a number of dangling OH groups corresponding to water molecules in acceptor-donor binding sites.