Infrared photodissociation spectroscopy of Na(NH3)n clusters: probing the solvent coordination

The first mass-selective vibrational spectra have been recorded for Na(NH3)n clusters. Infrared spectra have been obtained for n = 3-8 in the N-H stretching region. The spectroscopic work has been supported by ab initio calculations carried out at both the DFT(B3LYP) and MP2 levels, using a 6-311++G(d,p) basis set. The calculations reveal that the lowest energy isomer for n or= 7 is indicative of molecules entering a second solvation shell, i.e., the inner solvation shell around the sodium atom can accommodate a maximum of six NH3 molecules.     

Structures and infrared spectra of fluoride-hydrogen sulfide clusters from ab initio calculations: F(-)-(H2S)n, n = 1-5

Clusters formed between a fluoride anion and several hydrogen sulfide molecules have been investigated via ab initio calculations at the MP2 level of theory, using Dunning's augmented correlation consistent basis sets. Optimised geometries, vibrational frequencies, and enthalpy changes for the ligand association reactions are presented for clusters with up to five H2S ligands interacting with a F- anion. The minimum energy structure for the 1:1 F(-)-H2S complex features proton transfer from the H2S to the F- anion, forming a planar C(s) symmetry FH...SH- structure. For the F(-)-(H2S)2 cluster, the FH...SH- core remains and is solvated by a perturbed H2S ligand. For the larger F(-)-(H2S)(3-5) clusters, in addition to the FH...SH(-)-(H2S)n cluster forms, other minima featuring a 'solvated F-' anion are predicted. Calculated infrared spectra for the minima of each cluster size are presented to aid in assigning spectra from future experimental studies.     

Ab initio molecular dynamics studies of ionic dissolution and precipitation of sodium chloride and silver chloride in water clusters, NaCl(H2O)n and AgCl(H2O)n, n = 6, 10, and 14

An ab initio molecular dynamics method was used to compare the ionic dissolution of soluble sodium chloride (NaCl) in water clusters with the highly insoluble silver chloride (AgCl). The investigations focused on the solvation structures, dynamics, and energetics of the contact ion pair (CIP) and of the solvent-separated ion pair (SSIP) in NaCl(H(2)O)(n) and AgCl(H(2)O)(n) with cluster sizes of n = 6, 10 and 14. We found that the minimum cluster size required to stabilize the SSIP configuration in NaCl(H(2)O)(n) is temperature-dependent. For n = 6, both configurations are present as two distinct local minima on the free-energy profile at 100 K, whereas SSIP is unstable at 300 K. Both configurations, separated by a low barrier (<10 kJ mol(-1)), are identifiable on the free energy profiles of NaCl(H(2)O)(n) for n = 10 and 14 at 300 K, with the Na(+)/Cl(-) pairs being internally solvated in the water cluster and the SSIP configuration being slightly higher in energy (<5 kJ mol(-1)). In agreement with the low bulk solubility of AgCl, no SSIP minimum is observed on the free-energy profiles of finite AgCl(H(2)O)(n) clusters. The AgCl interaction is more covalent in nature, and is less affected by the water solvent. Unlike NaCl, AgCl is mainly solvated on the surface in finite water clusters, and ionic dissolution requires a significant reorganization of the solvent structure.     

Application of the valence bond mixing configuration diagrams to hypervalency in trihalide anions: a challenge to the Rundle-Pimentel model

The X(3)(-) hypercoordinated anions (H, F, Cl, Br, I) are studied by means of the breathing-orbital valence bond ab initio method. The valence bond wave functions describe the different X(3)(-) complexes in terms of only six valence bond structures and yield energies relative to the two exit channels, X(2) + X(-) and X(2)(-) + X(*), in very good agreement with reference CCSD(T) calculations. Although H(3)(-) is unstable and dissociates to H(2) + H(-), all the trihalogen anions are stable intermediates, Br(3)(-) and I(3)(-) being more stable than F(3)(-) and Cl(3)(-). As a challenge to the traditional Rundle-Pimentel model, the different energies of the hypercoordinated species relative to the normal-valent dissociation products X(2) + X(-) are interpreted in terms of valence bond configuration mixing diagrams and found to correlate with a single parameter of the X(2) molecule, its singlet-triplet energy gap. Examination of the six-structure wave functions show that H(3)(-), Cl(3)(-), Br(3)(-), and I(3)(-) share the same bonding picture and can be mainly described in terms of the interplay of two Lewis structures. On the other hand, F(3)(-) is bonded in a different way and possesses a significant three-electron bonding character that is responsible for the dissociation of this complex to F(2)(-) + F(*), instead of the more stable products F(2) + F(-). This counterintuitive preference for the thermodynamically disfavored exit channel is found to be an experimental manifestation of the large charge-shift resonance energy that generally characterizes fluorine-containing bonds.     

VBEFP: a valence bond approach that incorporates effective fragment potential method

An ab initio explicit solvation valence bond (VB) method, called VBEFP, is presented. The VBEFP method is one type of QM/MM approach in which the QM part of system is treated within the ab initio valence bond scheme and the solvent water molecules are accounted by the effective fragment potential (EFP) method, which is a polarized force field approach developed by Gordon et al. (J. Chem. Phys. 1996, 105, 1968). This hybrid method enables one to take the first-solvation shell and heterogeneous solvation effects into account explicitly with VB wave function. Therefore, the nature of chemical bonding and the mechanism of chemical reactions with explicit solvent environments can be explored at the ab inito VB level. In this paper, the hydrated metal-ligand complexes [M(2+)L](H(2)O)(n) (M(2+): Mg(2+), Zn(2+); L: NH(3), CH(2)O) are studied by the VBEFP method. Resonance energy and bond order are computed, and the influence of the solvent coordination and hydrogen bonding to the metal-ligand bonding are explored in the paper.     

Structure and dissociation characteristics of metal chloride anion clusters containing redox-active metal ions studied by laser desorption and electrospray ionization mass spectrometry and ab initio calculations

Copper chloride anion clusters with both copper oxidation states can be made by laser desorption of CuCl(2) crystals. We have used this method to study the dissociation characteristics of such cluster ions. The stability and the structure of the observed complexes were probed by ab initio calculations. These calculations show that many of these complexes are bridged structures. Thus, for the Cu(2)Cl(4) dimer anion, formally [ClCu-Cl-CuCl(2)](-) , with putative mixed copper oxidation states, the two copper ions become equivalent through bridging. Such bridging does not occur when redox inactive metal ions are present as in [ClCu-Cl-CaCl(2)](-) . By observing the dissociation characteristics of a variety of metal chloride cluster anions produced from binary mixtures, the following Cl(-) affinity order is obtained: FeCl(3) > CuCl > CaCl(2) > FeCl(2) > AgCl ≈ CuCl(2) ≈ ZnCl(2) > LiCl. Ab initio calculations on the Cl(-) affinity of selected chlorides confirm this order as do Cl(-) affinity estimates from the experimentally known vertical electron detachment energies of the superhalogens CaCl(3)(-) and LiCl(2)(-) . An equimolar mixture of CuCl(2) and FeCl(3) produces an intense cluster ion, which, from (65)Cu labeling experiments, is best described as FeCl(4)(-)···Cu(+)···(-)Cl(4) Fe, a Cu(+) bound superhalogen FeCl(4)(-) dimer. The Cu(+) ion can be replaced by the redox inactive alkali cations and by Ag(+) but these metal ion bound FeCl(4)(-) dimers show an entirely different fragmentation behavior which is attributed to the absence of bridging. Electrospray ionization (ESI) of CuCl(2) produces an extended series of (CuCl(2))(n) Cl(-) anions (n = 1-11) and so in ESI very limited reduction of Cu(2+) takes place. The (CuCl(2))(n) Cl(-) anions show an abundant dissociation via loss of neutral Cu(2)Cl(4) which according to our ab initio calculations is 9 kcal/mol more stable than two CuCl(2).     

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties.     

Ab initio molecular dynamics study of a mixture of HF(aq) and HCl(aq)

We have studied a mixture of HF and HCl molecules in water using Car-Parrinello ab initio molecular dynamics (CPMD). We have done simulations with 1 HF and 3 HCl molecules, 3 HF and 4 HCl, 6 HF and 8 HCl (6/8 simulation), and 14 HF molecules. All simulations consist of 32 molecules, and they were 10-96 ps long. The HF dissociation probability was around 30%, and HCl's was more than 90%. The solvation of the HF molecule was much better than the solvation of HCl. The solvation environment of F, both the F- ion and the F in HF, did not depend much on the acids concentration, whereas the Cl coordination numbers were rather sensitive to the concentration. In the 6/8 simulation, all XH-Y (X, Y = F, Cl) type molecules were observed and the FH-F was the most probable. In general, the molecular structures in mixed aqueous acid systems were similar to the pure HF(aq) and HCl(aq) systems.     

Chloride complexation by uranyl in a room temperature ionic liquid. A computational study

The stepwise addition of 1 to 4 Cl(-) anions to the uranyl cation has been studied via potential of mean force (PMF) calculations in the [BMI][Tf 2N] ionic liquid based on the 1-butyl-3-methylimidazolium cation (BMI(+)) and the bis(trifluoromethylsulfonyl)imide anion (Tf2N(-)). According to these calculations, the four Cl(-) complexation reactions are favored and UO2Cl4(2-) is the most stable chloride complex in [BMI][Tf2N]. The solvation of the different chloro-complexes is found to evolve from purely anionic (ca. 5 Tf2N(-) ions around UO2(2+)) to purely cationic (ca. 8.5 BMI (+) cations around UO2Cl4(2-)), with onion-type alternation of solvent shells. We next compare the solvation of the UO2Cl4(2-) complex to its reduced analogue UO2Cl4(3-) in the [BMI][Tf2N] and [MeBu3N][Tf2N] liquids that possess the same anion, but differ by their cation (imidazolium BMI(+) versus ammonium MeBu3N(+)). The overall solvation structure of both complexes is found to be similar in both liquids with a first solvation shell formed exclusively of solvent cations (about 9 BMI(+) cations or 7 MeBu3N(+) cations). However, a given complex is better solvated by the [BMI][Tf2N] liquid, due to hydrogen bonding interactions between Cl(-) ligands and imidazolium-ring C-H protons. According to free energy calculations, the gain in solvation energy upon reduction of UO2Cl4(2-) to UO2Cl4(3-) is found to be larger in [BMI][Tf2N] than in [MeBu3N][Tf2N], which is fully consistent with recent experimental results (Inorg. Chem. 2006, 45, 10419).     

The isobutylene-isobutane alkylation process in liquid HF revisited

Details on the mechanism of HF catalyzed isobutylene-isobutane alkylation were investigated. On the basis of available experimental data and high-level quantum chemical calculations, a detailed reaction mechanism is proposed taking into account solvation effects of the medium. On the basis of our computational results, we explain why the density of the liquid media and stirring rates are the most important parameters to achieve maximum yield of alkylate, in agreement with experimental findings. The ab initio Car-Parrinello molecular dynamics calculations show that isobutylene is irreversibly protonated in the liquid HF medium at higher densities, leading to the ion pair formation, which is shown to be a minimum on the potential energy surface after optimization using periodic boundary conditions. The HF medium solvates preferentially the fluoride anion, which is found as solvated [FHF](-) or solvated F(-.)(HF)(3). On the other hand, the tert-butyl cation is weakly solvated, where the closest HF molecules appear at a distance of about 2.9 Angstrom with the fluorine termination of an HF chain.     

Interaction of Co(m)O(-) (m = 1-3) with water: anion photoelectron spectroscopy and density functional calculations

We investigated the reactions between cobalt-oxides and water molecules using photoelectron spectroscopy and density functional calculations. It has been confirmed by both experimental observation and theoretical calculations that dihydroxide anions, Co(m)(OH)(2)(-) (m = 1-3), were formed when Co(m)O(-) clusters interact with the first water molecule. Addition of more water molecules produced solvated dihydroxide anions, Co(m)(OH)(2)(H(2)O)(n)(-) (m = 1-3). Hydrated dihydroxide anions, Co(m)(OH)(2)(H(2)O)(n)(-), are more stable than their corresponding hydrated metal-oxide anions, Co(m)O(H(2)O)(n+1)(-).     

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.     

A theoretical investigation of the interactions between water molecules and ionic liquids

Quantum chemical calculations have been used to investigate the interaction between water molecules and ionic liquids based on the imidazolium cation with the anions [Cl(-)], [Br(-)], [BF(4)(-)], and [PF(6)(-)]. The predicted geometries and interaction energies implied that the water molecules interact with the Cl(-), Br(-), and BF(4)(-0 anions to form X(-)...W (X = Cl or Br, W = H(2)O), 2X-...2W, BF(4)(-)...W, and W...BF(4)(-)...W complexes. The hydrophobic PF(6)(-) anion could not form a stable complex with the water molecules at the density functional theory (DFT) level. Further studies indicate that the cation could also form a strong interaction with the water molecules. The 1-ethyl-3-methylimidazolium cation (Emim(+)) has been used as a model cation to investigate the interaction between a water molecule and a cation. In addition, the interaction between the ion pairs and the water was studied by using 1-ethyl-3-methylimidazolium chloride (Emim x Cl) as a model ionic liquid. The strengths of the interactions in these categories follow the trend anion-W > cation-W > ion pair-W.     

Effect of solvation on the vertical ionization energy of thymine: from microhydration to bulk

The effect of hydration on the vertical ionization energy (VIE) of thymine was characterized using equation-of-motion ionization potential coupled-cluster (EOM-IP-CCSD) and effective fragment potential (EFP) methods. We considered several microsolvated clusters as well as thymine solvated in bulk water. The VIE in bulk water was computed by averaging over solvent-solute configurations obtained from equilibrium molecular dynamics trajectories at 300 K. The effect of microsolvation was analyzed and contrasted against the combined effect of the first solvation shell in bulk water. Microsolvation reduces the ionization energy (IE) by about 0.1 eV per water molecule, while the first solvation shell increases the IE by 0.1 eV. The subsequent solvation lowers the IE, and the bulk value of the solvent-induced shift of thymine's VIE is approximately -0.9 eV. The combined effect of the first solvation shell was explained in terms of specific solute-solvent interactions, which were investigated using model structures. The convergence of IE to the bulk value requires the hydration sphere of approximately 13.5 ? radius. The performance of the EOM-IP-CCSD/EFP scheme was benchmarked against full EOM-IP-CCSD using microhydrated structures. The errors were found to be less than 0.01-0.02 eV. The relative importance of the polarization and higher multipole moments in EFP model was also investigated.     

Ab initio potential and dipole moment surfaces for water. II. Local-monomer calculations of the infrared spectra of water clusters

We employ recent flexible ab initio potential energy and dipole surfaces [Y. Wang, X. Huang, B. C. Shepler, B. J. Braams, and J. M. Bowman, J. Chem. Phys. 134, 094509 (2011)] to the calculation of IR spectra of the intramolecular modes of water clusters. We use a quantum approach that begins with a partitioned normal-mode analysis of perturbed monomers, and then obtains solutions of the corresponding Schro?dinger equations for the fully coupled intramolecular modes of each perturbed monomer. For water clusters, these modes are the two stretches and the bend. This approach is tested against benchmark calculations for the water dimer and trimer and then applied to the water clusters (H(2)O)(n) for n = 6-10 and n = 20. Comparisons of the spectra are made with previous ab initio harmonic and empirical potential calculations and available experiments.     

Hemibonding of hydroxyl radical and halide anion in aqueous solution

Molecular geometries and properties of the possible reaction products between the hydroxyl radical and the halide anions in aqueous solution were investigated. The formation of two-center three-electron bonding (hemibonding) between the hydroxyl radical and halide anions (Cl, Br, I) was examined by density functional theory (DFT) calculation with a range-separated hybrid (RSH) exchange-correlation functional. The long-range corrected hybrid functional (LC-ωPBE), which have given quantitatively satisfactory results for odd electron systems and excited states, was examined by test calculations for dihalogen radical anions (X(2)(-); X = Cl, Br, I) and hydroxyl radical-water clusters. Equilibrium geometries with hemibonding between the hydroxyl radical and halide anions were located by including four hydrogen-bonded water molecules. Excitation energies and oscillator strengths of σ-σ* transitions calculated by the time-dependent DFT method showed good agreement with observed values. Calculated values of the free energy of reaction on the formation of hydroxyl halide radical anion from the hydroxyl radical and halide anion were endothermic for chloride but exothermic for bromide and iodide, which is consistent with experimental values of equilibrium constants.     

Gradients of the polarization energy in the effective fragment potential method

The effective fragment potential (EFP) method is an ab initio based polarizable classical method in which the intermolecular interaction parameters are obtained from preparative ab initio calculations on isolated molecules. The polarization energy in the EFP method is modeled with asymmetric anisotropic dipole polarizability tensors located at the centroids of localized bond and lone pair orbitals of the molecules. Analytic expressions for the translational and rotational gradients (forces and torques) of the EFP polarization energy have been derived and implemented. Periodic boundary conditions (the minimum image convention) and switching functions have also been implemented for the polarization energy, as well as for other EFP interaction terms. With these improvements, molecular dynamics simulations can be performed with the EFP method for various chemical systems.     

Solvation effects on angular distributions in H- (NH3)n and NH2(-) (NH3)n photodetachment: role of solute electronic structure

We report 355 and 532 nm photoelectron imaging results for H(-)(NH(3))(n) and NH(2)(-)(NH(3))(n), n = 0-5. The photoelectron spectra are consistent with the electrostatic picture of a charged solute (H(-) or NH(2)(-)) solvated by n ammonia molecules. For a given number of solvent molecules, the NH(2)(-) core anion is stabilized more strongly than H(-), yet the photoelectron angular distributions for solvated H(-) deviate more strongly from the unsolvated limit than those for solvated NH(2)(-). Hence, we conclude that solvation effects on photoelectron angular distributions are dependent on the electronic structure of the anion, i.e., the type of the initial orbital of the photodetached electron, rather than merely the strength of solvation interactions. We also find evidence of photofragmentation and autodetachment of NH(2)(-)(NH(3))(2-5), as well as autodetachment of H(-)(NH(3))(5), upon 532 nm excitation of these species.     

Interaction of NaCl with solid water

The interaction of NaCl with solid water, deposited on tungsten at 80 K, was investigated with metastable impact electron spectroscopy (MIES) and ultraviolet photoelectron spectroscopy (UPS) (He I). We have studied the ionization of Cl(3p) and the 1b(1), 3a(1), and 1b(2) bands of molecular water. The results are supplemented by first-principles density functional theory (DFT) calculations of the electronic structure of solvated Cl(-) ions. We have prepared NaCl/water interfaces at 80 K, NaCl layers on thin films of solid water, and H(2)O ad-layers on thin NaCl films; they were annealed between 80 and 300 K. At 80 K, closed layers of NaCl on H(2)O, and vice versa, are obtained; no interpenetration of the two components H(2)O and NaCl was observed. However, ionic dissociation of NaCl takes place when H(2)O and NaCl are in direct contact. Above 115 K solvation of the ionic species Cl(-) becomes significant. Our results are compatible with a transition of Cl(-) species from an interface site (Cl in direct contact with the NaCl lattice) to an energetically favored configuration, where Cl species are solvated. The DFT calculations show that Cl(-) species, surrounded by their solvation shell, are nevertheless by some extent accessed by MIES because the Cl(3p)-charge cloud extends through the solvation shell. Water desorption is noticeable around 145 K, but is not complete before 170 K, about 15 K higher than for pure solid water. Above 150 K the NaCl-induced modification of the water network gives rise to gas phase like structures in the water spectra. In particular, the 3a(1) emission turns into a well-defined peak. This suggests that under these conditions water molecules interact mainly with Cl(-) rather than among themselves. Above 170 K only Cl is detected on the surface and desorbs around 450 K.