Quantum vibrational analysis of hydrated ions using an ab initio potential

We present full-dimensional potential energy surfaces (PESs) for hydrated chloride based on the sum of ab initio (H(2)O)Cl(-), (H(2)O)(2), and (H(2)O)(3) potentials. The PESs are shown to predict minima and corresponding harmonic frequencies accurately on the basis of comparisons with previous and new ab initio calculations for (H(2)O)(2)Cl(-), (H(2)O)(3)Cl(-), and (H(2)O)(4)Cl(-). An estimate of the effect of the 3-body water interaction is made using a simple 3-body water potential that was recently fit to tens of thousands of ab initio 3-body energies. Anharmonic, coupled vibrational calculations are presented for these clusters, using the "local monomer model" for the high frequency intramolecular modes. This model is tested against previous "exact" calculations for (H(2)O)Cl(-). Radial distribution functions at 0 K obtained from quantum zero-point wave functions are also presented for the (H(2)O)(2)Cl(-) and (H(2)O)(3)Cl(-) clusters.     

Automated exploration of stable isomers of H(+)(H(2)O)(n) (n = 5-7) via ab initio calculations: An application of the anharmonic downward distortion following algorithm

Isomers of protonated water clusters H(+)(H(2)O)(n) (n = 5-7) have been explored on ab initio potential energy surfaces by means of the anharmonic downward distortion following algorithm. Totally, 9, 24, and 131 isomers for n = 5, 6, and 7, respectively, were located by the automatic exploration, and all of known important isomers previously reported by conventional geometry optimization approaches have been included in the present results. Moreover, structure transitions depending on n and temperature, which were observed by experimental studies, could be reproduced via thermodynamic simulation on the basis of the superposition approach and the present isomer sets.     

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.     

Hydration profiles of aromatic amino acids: conformations and vibrations of L-phenylalanine-(H2O)n clusters

IR-UV double resonance spectroscopy and ab initio calculations were employed to investigate the structures and vibrations of the aromatic amino acid, L-phenylalanine-(H(2)O)(n) clusters formed in a supersonic free jet. Our results indicate that up to three water molecules are preferentially bound to both the carbonyl oxygen and the carboxyl hydrogen of L-phenylalanine (L-Phe) in a bridged hydrogen-bonded conformation. As the number of water molecules is increased, the bridge becomes longer. Two isomers are found for L-Phe-(H(2)O)(1), and both of them form a cyclic hydrogen-bond between the carboxyl group and the water molecule. In L-Phe-(H(2)O)(2), only one isomer was identified, in which two water molecules form extended cyclic hydrogen bonds with the carboxyl group. In the calculated structure of L-Phe-(H(2)O)(3) the bridge of water molecules becomes larger and exhibits an extended hydrogen-bond to the pi-system. Finally, in isolated L-Phe, the D conformer was found to be the most stable conformer by the experiment and by the ab initio calculation.     

Hydration of the bisulfate ion: atmospheric implications

Using molecular dynamics configurational sampling combined with ab initio energy calculations, we determined the low energy isomers of the bisulfate hydrates. We calculated the CCSD(T) complete basis set (CBS) binding electronic and Gibbs free energies for 53 low energy isomers of HSO(4)(-)(H(2)O)(n=1-6) and derived the thermodynamics of adding waters sequentially to the bisulfate ion and its hydrates. Comparing the HSO(4)(-)/H(2)O system to the neutral H(2)SO(4)/H(2)O cluster, water binds more strongly to the anion than it does to the neutral molecules. The difference in the binding thermodynamics of HSO(4)(-)/H(2)O and H(2)SO(4)/H(2)O systems decreases with increasing number of waters. The thermodynamics for the formation of HSO(4)(-)(H(2)O)(n=1-5) is favorable at 298.15 K, and that of HSO(4)(-)(H(2)O)(n=1-6) is favorable for T < 273.15 K. The HSO(4)(-) ion is almost always hydrated at temperatures and relative humidity values encountered in the troposphere. Because the bisulfate ion binds more strongly to sulfuric acid than it does to water, it is expected to play a role in ion-induced nucleation by forming a strong complex with sulfuric acid and water, thus facilitating the formation of a critical nucleus.     

First principles study on the solvation and structure of C2O4(2-)(H2O)n, n = 6-12

The structures and energies of hydrated oxalate clusters, C2O4(2-)(H2O)n, n = 6-12, are obtained by density functional theory (DFT) calculations and compared to SO4(2-)(H2O)n. Although the evolution of the cluster structure with size is similar to that of SO4(2-)(H2O)n, there are a number of important and distinctive futures in C2O4(2-)(H2O)n, including the separation of the two charges due to the C-C bond in C2O4(2-), the lower symmetry around C2O4(2-), and the torsion along the C-C bond, that affect both the structure and the solvation energy. The solvation dynamics for the isomers of C2O4(2-)(H2O)12 are also examined by DFT based ab initio molecular dynamics.     

Infrared predissociation spectroscopy of M+ (C6H6)(1-4)(H2O)(1-2)Ar(0-1) cluster ions, M = Li, Na

Infrared predissociation (IRPD) spectra of Li(+)(C(6)H(6))(1-4)(H(2)O)(1-2)Ar(0-1) and Na(+)(C(6)H(6))(2-4)(H(2)O)(1-2)Ar(1) are presented along with ab initio calculations. The results indicate that the global minimum energy structure for Li(+)(C(6)H(6))(2)(H(2)O)(2) has each water forming a π-hydrogen bond with the same benzene molecule. This bonding motif is preserved in Li(+)(C(6)H(6))(3-4)(H(2)O)(2)Ar(0-1) with the additional benzene ligands binding to the available free OH groups. Argon tagging allows high-energy Li(+)(C(6)H(6))(2-4)(H(2)O)(2)Ar isomers containing water-water hydrogen bonds to be trapped and detected. The monohydrated, Li(+) containing clusters contain benzene-water interactions with varying strength as indicated by shifts in OH stretching frequencies. The IRPD spectra of M(+)(C(6)H(6))(1-4)(H(2)O)(1-2)Ar are very different for lithium-bearing versus sodium-bearing cluster ions emphasizing the important role of ion size in determining the most favorable balance of competing noncovalent interactions.     

Characterization of hydrated Na+(phenol) and K+(phenol) complexes using infrared spectroscopy

Hydrated alkali metal ion-phenol complexes were studied to model these species in aqueous solution for M=Na and K. IR predissociation spectroscopy in the O-H stretch region was used to analyze the structures of M+(Phenol)(H2O)n cluster ions, for n = 1-4. The onset of hydrogen bonding was observed to occur at n=4. Ab initio calculations were used to qualitatively explore the types of hydrogen-bonded structures of the M+(Phenol)(H2O)4 isomers. By combining the ab initio calculations and IR spectra, several different structures were identified for each metal ion. In contrast to benzene, detailed in a previous study of Na+(Benzene)n(H2O)m [J. Chem. Phys. 110, 8429 (1999)], phenol is able to bind directly to Na+ even in the presence of four waters. This is likely the result of the sigma-type interaction between the phenol oxygen and the ion. With K+, the dominant isomers are those in which the phenol O-H group is involved in a hydrogen bond with the water molecules, while with Na+, the dominant isomers are those in which the phenol O-H group is free and the water molecules are hydrogen-bonded to each other. Spectra and ab initio calculations for the M+(Phenol)Ar cluster ions for M=Na and K are reported to characterize the free phenol O-H stretch in the M+(Phenol) complex. While pi-type configurations were observed for binary M+(Phenol) complexes, sigma-type configurations appear to dominate the hydrated cluster ions.     

Cu3C4-: a new sandwich molecule with two revolving C2(2-) units

A combined photoelectron spectroscopy (PES) and ab initio study was carried out on a novel copper carbide cluster in the gas phase: Cu(3)C(4)(-). It was generated in a laser vaporization cluster source and appeared to exhibit enhanced stability among the Cu(3)C(n)(-) series. Its PES spectra were obtained at several photon energies, showing numerous well-resolved bands. Extensive ab initio calculations were performed on Cu(3)C(4)(-), and two isomers were identified: a C(2) structure ((1)A) with a Cu(3)(3+) triangular group sandwiched by two C(2)(2-) units and a linear CuCCCuCCCu structure (D(infinity)(h), (1)Sigma(g)(+)). A comparison of ab initio PES spectra with experimental data showed that the sandwich Cu(3)C(4)(-) cluster was solely responsible for the observed spectra and the linear isomer was not present, suggesting that the C(2) structure is the global minimum in accordance with CCSD(T)/6-311+G predictions. Interestingly, a relatively low barrier (0.4-0.6 kcal/mol) was found for the internal rotation of the C(2)(2-) units in the sandwich Cu(3)C(4)(-). To test different levels of theory in describing the Cu(m)C(n)(-) systems and lay foundations for the validity of the theoretical methods, extensive calculations at a variety of levels were also carried out on a simpler copper carbide species CuC(2)(-), where two isomers were found to be close in energy: a linear one (C(infinity)(v), (1)Sigma(+)) and a triangular one (C(2)(v), (1)A(1)). The calculated electronic transitions for CuC(2)(-) were also compared with the PES data, in which both isomers were present.     

Benchmark structures and binding energies of small water clusters with anharmonicity corrections

For (H(2)O)(n) where n = 1-10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.     

Extra electron in (H2O)24- cluster isomers: a theoretical study

The isomers of (H(2)O)(24) (-) tetrakaidecahedral cluster are studied by applying the Becke-3-parameter density functional theory and Lee-Yang-Parr correlation functional (B3LYP) and 6-311++G** basis set. Three isomers are selected on the basis of stabilization energy values. The vertical electron dissociation energies (VDE) of these isomers are 1.353, 0.404, and 0.258 eV, respectively. The experimental VDE value of 1.31 eV [J. Chem. Phys. 92, 3980 (1990)] for this cluster size is in excellent agreement with that calculated for isomer 1, suggesting the dominance of this isomer in the experiment. Four water molecules in this isomer share most of the -1 charge. These four water molecules have non-H-bonding H (NHB H) atoms turned toward the cavity, and the inward turned H atoms exhibit a significant lowering of O-H stretch frequency compared to that of a monomer. Isomers 2 and 3 have all 12 NHB H atoms projected outward and have the -1 charge distributed among 7-8 water molecules on the cluster surface.     

On the structure of the matrix isolated water trimer

Infrared spectra of partially deuterated water trimers have been investigated. It is found that HDO(H(2)O)(2) has a single, bound OD stretching fundamental, (HDO)(2)H(2)O two bound OD stretches. (HDO)(3) has a single, bound OD stretch and (H(2)O)(3) has a pair of bound OH stretches. Ab initio and discrete Fourier transform (DFT) calculations predict that the water trimer has C(1) symmetry with six different, isoenergetic minima. These calculations consequently give three numerically different OD stretches for HDO(H(2)O)(2), six for (HDO)(2)H(2)O, three for (HDO)(3), and three bound OH stretches for (H(2)O)(3). The connection between the observations and the pseudorotation of the trimer is discussed with the help of Wales' pseudorotation model. It is found that pseudorotation is sufficiently fast to average the effective symmetry of the A(3) trimer to C(3h) and to eliminate the difference between the different ab initio minima for A(2)B. The only exception is (H(2)O)(3) where the splitting between the different bound OH stretches is largest. Here a doublet is observed due to incomplete averaging. DFT calculations indicate that the D-bonded form of HDO(H(2)O)(2) is between 50 and 60 cm(-1) more stable than the H-bonded form. The energy difference is determined by differences in zero point vibration energy of intermolecular librations of the two forms. Attempts to measure the energy difference indicate that the energy difference is larger, of the order of 100 cm(-1).     

Unraveling the mechanisms of O2 activation by size-selected gold clusters: transition from superoxo to peroxo chemisorption

The activation of dioxygen is a key step in CO oxidation catalyzed by gold nanoparticles. It is known that small gold cluster anions with even-numbered atoms can molecularly chemisorb O(2) via one-electron transfer from Au(n)(-) to O(2), whereas clusters with odd-numbered atoms are inert toward O(2). Here we report spectroscopic evidence of two modes of O(2) activation by the small even-sized Au(n)(-) clusters: superoxo and peroxo chemisorption. Photoelectron spectroscopy of O(2)Au(8)(-) revealed two distinct isomers, which can be converted from one to the other depending on the reaction time. Ab initio calculations show that there are two close-lying molecular O(2)-chemisorbed isomers for O(2)Au(8)(-): the lower energy isomer involves a peroxo-type binding of O(2) onto Au(8)(-), while the superoxo chemisorption is a slightly higher energy isomer. The computed detachment transitions of the superoxo and peroxo species are in good agreement with the experimental observation. The current work shows that there is a superoxo to peroxo chemisorption transition of O(2) on gold clusters at Au(8)(-): O(2)Au(n)(-) (n = 2, 4, 6) involves superoxo binding and n = 10, 12, 14, 18 involves peroxo binding, whereas the superoxo binding re-emerges at n = 20 due to the high symmetry tetrahedral structure of Au(20), which has a very low electron affinity. Hence, the two-dimensional (2D) Au(8)(-) is the smallest anionic gold nanoparticle that prefers peroxo binding with O(2). At Au(12)(-), although both 2D and 3D isomers coexist in the cluster beam, the 3D isomer prefers the peroxo binding with O(2).     

Relative stability of planar versus double-ring tubular isomers of neutral and anionic boron cluster B20 and B20-

High-level ab initio molecular-orbital methods have been employed to determine the relative stability among four neutral and anionic B20 isomers, particularly the double-ring tubular isomer versus three low-lying planar isomers. Calculations with the fourth-order Moller-Plessset perturbation theory [MP4(SDQ)] and Dunning's correlation consistent polarized valence triple zeta basis set as well as with the coupled-cluster method including single, double, and noniteratively perturbative triple excitations and the 6-311G(d) basis set show that the double-ring tubular isomer is appreciably lower in energy than the three planar isomers and is thus likely the global minimum of neutral B20 cluster. In contrast, calculations with the MP4(SDQ) level of theory and 6-311+G(d) basis set show that the double-ring anion isomer is appreciably higher in energy than two of the three planar isomers. In addition, the temperature effects on the relative stability of both 10B20- and 11B20- anion isomers are examined using the density-functional theory. It is found that the three planar anion isomers become increasingly more stable than the double-ring isomer with increasing the temperature. These results are consistent with the previous conclusion based on a joint experimental/simulated anion photoelectron spectroscopy study [B. Kiran et al., Proc. Natl. Acad. Sci. U.S.A. 102, 961 (2005)], that is, the double-ring anion isomer is notably absent from the experimental spectra. The high stability of the double-ring neutral isomer of B20 can be attributed in part to the strong aromaticity as characterized by its large negative nucleus-independent chemical shift. The high-level ab initio calculations suggest that the planar-to-tubular structural transition starts at B20 for neutral clusters but should occur beyond the size of B20- for the anion clusters.     

Ab initio study of the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2) and related C(4)H(10)(2+) isomers.

Ab initio calculations at the MP4(SDTQ)/6-311G//MP2/6-31G level were performed to study the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2), and related C(4)H(10)(2+) isomers. Two doubly hydrogen bridged diborane type trans 1 and cis 2 isomers were located as minima. The trans isomer was found to be more favorable than cis isomer by only 0.6 kcal/mol. Several other minima for C(4)H(10)(2+) were also located. However, the global energy minimum corresponds to C-H (C(4) position) protonated 2-butyl cation 10. Structure 10 was computed to be substantially more stable than 1 by 31.7 kcal/mol. The structure 10 was found to be lower in energy than 2-butyl cation 13 by 34.4 kcal/mol.     

Decomposition of protonated formic acid: One transition state—Two product channels

The unimolecular chemistry of protonated formic acid, [HCOOH]H(+), has been investigated by analyzing the fragmentation of metastable ions (MI) during flight in a sector mass spectrometer, and by proton transfer to formic acid in a Fourier-transform ion cyclotron resonance (FT-ICR) mass spectrometer. High level ab initio calculations have been used to model the relevant parts of the potential energy surface (PES). In addition, ab initio direct dynamics calculations have been conducted, tracing out 60 different reaction trajectories. The only stable isomer in the mass spectrometric experiments is HC(OH)(2)(+), which is the precursor to both observed ionic products, HCO(+) and H(3)O(+), via the same saddle point of the potential energy surface. The detailed motion of the dissociating molecule during passage of the post-transition state region of the PES therefore determines which product ion is formed. After passing the TS a transient HC(O)OH(2)(+) molecule is first formed. High total energy increases the probability that the nascent water molecule will have sufficient speed to escape the HCO(+) moiety. Otherwise, typically at low energies, the two units recombine, upon which intra-complex proton transfer is very likely. Eventually, this will give the more stable H(3)O(+).     

Theoretical study of the solvation of fluorine and chlorine anions by water

The solvation of fluoride and chloride anions (F(-) and Cl(-), respectively) by water has been studied using effective fragment potentials (EFPs) for the water molecules and ab initio quantum mechanics for the anions. In particular, the number of water molecules required to fully surround each anion has been investigated. Monte Carlo calculations have been used in an attempt to find the solvated system X(-)(H(2)O)(n) (X = F, Cl) with the lowest energy for each value of n. It is predicted that 18 water molecules are required to form a complete solvation shell around a Cl(-) anion, where "complete solvation" is interpreted as an ion that is completely surrounded by solvent molecules. Although fewer water molecules may fully solvate the Cl(-) anion, such structures are higher in energy than partially solvated molecules, up to n > or = 18. Calculations on the F(-) anion suggest that 15 water molecules are required for a complete solvation shell. The EFP predictions are in good agreement with the relative energies predicted by ab initio energy calculations at the EFP geometries.     

Conformational characterization of lanthanide(III)-DOTA complexes by ab initio investigation in vacuo and in aqueous solution

The conformational behavior of four [Ln(DOTA)(H(2)O)](-) systems (Ln = La, Gd, Ho, and Lu) has been characterized by means of ab initio calculations performed in vacuo and in aqueous solution, the latter by using the polarizable continuum model (PCM). Calculated molecular geometries and conformational energies of the [Ln(DOTA)(H(2)O)](-) systems, and interconversion mechanisms, barriers, and (13)C NMR spectra of the [Lu(DOTA)](-) complex are compared with experimental values. For each system, geometry optimizations, performed in vacuo and in solution at the HF/3-21G level and using a 46+4f(n) core electron effective core potential (ECP) for lanthanides, provide two minima corresponding to a square antiprismatic (A) and an inverted antiprismatic (IA) coordination geometry. All the systems are nonacoordinated, with the exception of the IA isomer of the Lu complex that, from in solution calculations, is octacoordinated, in agreement with experimental data. On comparing the in vacuo relative free energies calculated at different theory levels it can be seen that the nonacoordinated species dominates at the beginning of the lanthanide series while the octacoordinated one does so at the end. Furthermore, on passing along the series the IA isomer becomes less and less favored with respect to A and for the Lu complex a stabilization of the IAisomer is observed in solution (but not in vacuo), in agreement with the experimental data. Investigation of the A<-->IA isomerization process in the [Lu(DOTA)](-) system provides two different interconversion mechanisms: a single-step process, involving the simultaneous rotation of the acetate arms, and a multistep path, involving the inversion of the cyclen cycle configuration. While in vacuo the energy barrier for the acetate arm rotation is higher than that involved in the ring inversion (23.1 and 13.1 kcal mol(-)(1) at the B3LYP/6-311G level, respectively), in solution the two mechanisms present comparable barriers (14.7 and 13.5 kcal mol(-)(1)), in fairly good agreement with the experimental values. The NMR shielding constants for the two isomers of the [Lu(DOTA)](-) complex have been calculated by means of the ab initio GIAO and CSGT methods, and using a 46-core-electron ECP for Lu. The calculated (13)C NMR chemical shifts are in close agreement with the experimental values (rms 3.3 ppm, at the HF/6-311G level) and confirm the structural assignment of the two isomers based on experimental NMR spectra in solution. The results demonstrate that our computational approach is able to predict several physicochemical properties of lanthanide complexes, allowing a better characterization of this class of compounds for their application as contrast agents in medical magnetic resonance imaging (MRI).     

Micro-hydration of the MgNO3+ cation in the gas phase.

Coordination complexes of the magnesium nitrate cation with water [MgNO(3)(H(2)O)(n)](+) up to n=7 are investigated by experiment and theory. The fragmentation patterns of [MgNO(3)(H(2)O)(n)](+) clusters generated via electrospray ionization indicate a considerable change in stability between n=3 and 4. Further, ion-molecule reactions of mass-selected [MgNO(3)(H(2)O)(n)](+) cations with D(2)O reveal the occurrence of consecutive replacement of water ligands by heavy water, and in this respect the complexes with n=4 and 5 are somewhat more reactive than their smaller homologs with n=1-3 as well as the larger clusters with n=6 and 7. For the latter two ions, the theory suggests the existence of isomers, such as complexes with monodentate nitrato ligands as well as solvent-separated ion pairs with a common solvation shell. The reactions observed and the ion thermochemistry are discussed in the context of ab initio calculations, which also reveal the structures of the various hydrated cation complexes.     

An ab Initio Study of Tetrathiosquaric Acid

Inthepastfewdecades,sincetheoxocarbonanionsC.O.'-wererecognizedasmembersofanewclassofstabilizedcarbocyclicnonbenzenoidaromaticcompounds"',theknowledgeoftheirchemicalandphysicalpropertieshasbeendevelopingrapidly'Theiruniqueelectronicstructures,highdegreeofsymmetry,andestheticallybeauti-fulgeometrygeneratedafreshimpetustostudytheeffectofreplacingtheoxygensbyvariousfunctionalgroups,suchassulfur',andselenium[11.Analogoustotheoxocar-bondianions,thethioxocarbondianionscanbecharacterizedbythegeneralf…