Charge transfer to solvent (CTTS) energies of small X−(H2O)n=1–4 (X=F, Cl, Br, I) clusters: Ab initio study

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.

Ab initiostudy of hydrated potassium halides KX(H2O)1–6 (X=F,Cl,Br,I)

The ionic dissociation of salts was examined with a theoretical study of KX (X=F,Cl,Br,I) hydrated by up to six water molecules KX(H2O)n (n=1-6). Calculations were done using the density functional theory and second order Møller-Plesset (MP2) perturbational theory. To provide more conclusive results, single point energy calculations using the coupled cluster theory with single, double, and perturbative triple excitations were performed on the MP2 optimized geometries. The dissociation feature of the salts was examined in terms of K-X bond lengths and K-X stretch frequencies. In general, the successive incorporation of water molecules to the cluster lengthens the K-X distance, and consequently the corresponding frequency decreases. Near 0 K, the KX salt ion pairs can be partly separated by more than five water molecules. The pentahydrated KX salt is partly dissociated, though these partly dissociated structures are almost isoenergetic to the undissociated ones for KFKCl. For the hexahydrated complexes, KF is undissociated, KClKBr is partly dissociated, and KI is dissociated (though this dissociated structure is nearly isoenergetic to a partly dissociated one). On the other hand, at room temperature, the penta- and hexahydrated undissociated structures which have less hydrogen bonds are likely to be more stable than the partly dissociated ones because of the entropy effect. Therefore, the dissociation at room temperature could take place for higher clusters than the hexahydrated ones.

Structures, energetics, and spectra of aqua-cesium (I) complexes: An ab initio and experimental study

The design of cesium-selective ionophores must include the nature of cesium-water interactions. The authors have carried out extensive ab initio and density functional theory calculations of hydrated cesium cations to obtain reasonably accurate energetics, thermodynamic quantities, and IR spectra. An extensive search was made to find the most stable structures. Since water...water interactions are important in the aqua-Cs+ clusters, the authors investigated the vibrational frequency shifts as a function of the number of water molecules and the frequency characteristics with and without the presence of outer-shell water molecules. The predicted vibrational frequencies were then compared with the infrared photodissociation spectra of argon-tagged hydrated cesium cluster ions. This comparison allowed the identification of specific hydrogen-bonding structures present in the experimental spectra.

Structure, energy, and IR spectra of I2∙−.nH2O clusters (n=1–8): A theoretical study

The authors report theoretical results on structure, bonding, energy, and infrared spectra of iodine dimer radical anion hydrated clusters, I(2) (-).nH(2)O (n=1-8), based on a systematic study following density functional theory. Several initial guess structures are considered for each size cluster to locate minimum energy conformers with a Gaussian 6-311++G(d,p) split valence basis function (triple split valence 6-311 basis set is applied for iodine). It is observed that three different types of hydrogen bonded structures, namely, symmetrical double hydrogen bonding, single hydrogen bonding, and interwater hydrogen bonding structures, are possible in these hydrated clusters. But conformers having interwater hydrogen bonding arrangements are more stable compared to those of double or single hydrogen bonded structures. It is also noticed that up to four solvent H(2)O units can reside around the solute in interwater hydrogen bonding network. At the maximum six H(2)O units are independently linked to the dimer anion having four double hydrogen bonding and two single hydrogen bonding, suggesting the hydration number of I(2) (-) to be 6. However, conformers having H(2)O units independently linked to the iodine dimer anion are not the most stable structures. In all these hydrated clusters, the odd electron is found to be localized over two I atoms and the two atoms are bound by a three-electron hemi bond. The solvation, interaction, and vertical detachment energies are calculated for all I(2) (-).nH(2)O clusters. Energy of interaction and vertical detachment energy profiles show stepwise saturation, indicating geometrical shell closing in the hydrated clusters, but solvation energy profile fails to show such behavior. A linear correlation is observed between the calculated energy of interaction and vertical detachment energy. It is observed that formation of I(2) (-)-water cluster induces significant shifts from the normal O-H stretching modes of isolated H(2)O. However, bending mode of H(2)O remains insensitive to the successive addition of solvent H(2)O units. Weighted average energy profiles and IR spectra are reported for all the hydrated clusters based on the statistical population of individual conformers at room temperature.

De novo design approach based on nanorecognition toward development of functional molecules/materials and nanosensors/nanodevices

For the design of functional molecules and nanodevices, it is very useful to utilize nanorecognition (which is governed mainly by interaction forces such as hydrogen bonding, ionic interaction, π-H/π-π interactions, and metallic interactions) and nanodynamics (involving capture, transport, and release of electrons, photons, or protons). The manifestation of these interaction forces has led us to the design and realization of diverse ionophores/receptors, organic nanotubes, nanowires, molecular mechanical devices, molecular switches, enzyme mimetics, protein folding/unfolding, etc. In this review, we begin with a brief discussion of the interaction forces, followed by some of our representative applications. We discuss ionophores with chemo-sensing capability for biologically important cations and anions and explain how the understanding of hydrogen bonding and π-interactions has led to the design of self-assembled nanotubes from calix[4]hydroquinone (CHQ). The binding study of neutral and cationic transition metals with the redox system of hydroquinone (HQ) and quinone (Q) predicts what kind of nanostructures would form. Finally, we look into the conformational changes between stacked and edge-to-face conformers in π-benzoquinone-benzene complexes controlled by alternating electrochemical potential. The resulting flapping motion illustrates a promising pathway toward the design of mobile nanomechanical devices.

Interaction induced shifts in O–H stretching frequency of water in halide-ion water clusters: A microscopic approach with a bond descriptor

Based on conceptual density functional theory, two new descriptors, the normalized atom-condensed Fukui functions (NFFs) and bond deformation kernel (BDK), are proposed to describe the O-H stretching frequency shifts in water halide-ion clusters by including local polarizations in the system. NFF, along with BDK, can be successfully used to describe interactions between different atoms, especially in cases where polarizations are quite important.

Hydration Phenomena of Sodium and Potassium Hydroxides by Water Molecules

The hydrated structures, dissociation energies, thermodynamic quantities, infrared spectra, and electronic properties of alkali-metal hydroxides (MOH, M = Na and K) hydrated by up to six water molecules [MOH(H(2)O)(n=1-6)], are investigated by using the density functional theory and Møller-Plesset second-order perturbation theory. Further accurate analysis based on the coupled cluster theory with singles, doubles, and perturbative triples excitations is more consistent with the MP2 results. NaOH shows a peculiar trend in dissociation: it begins to form a partially dissociated structure for n = 3, and it dissociates for n = 4 and 6, whereas it is undissociated for n = 5. However, for n = 5, the dissociated structure is nearly isoenergetic to the undissociated structure. For KOH, it begins to show partial dissociation for n = 5, and complete dissociation for n = 6.

Microhydration shell structure in Cl2∙−∙nH2O clusters: A theoretical study

We present the results of a detailed study on structure and electronic properties of hydrated cluster Cl2*-.nH2O (n = 1-7) based on a nonlocal density functional, namely, Becke's [J. Chem. Phys. 98, 1372 (1993)] half and half hybrid exchange-correlation functional with a split valence 6-311++G(d,p) basis function. Geometry optimizations for all the clusters are carried out with various possible initial guess structures without any symmetry restriction. Several minimum energy structures (conformers) are predicted with a small difference in total energy. There is a competition between the binding of solvent H2O units with Cl2*- dimer radical anion directly through ion-molecule interaction and forming interwater hydrogen-bonding network in Cl2*-.nH2O (n > or = 2) hydrated cluster. Structure having interwater H-bonded network is more stable over the structure where H2O units are connected to the solute dimer radical anion Cl2*- rather independently either by single or double H bonding in a particular size (n) of hydrated cluster Cl2*-.nH2O. At the maximum four solvent H2O units reside in interwater H-bonding network present in these hydrated clusters. It is observed that up to six H2O units are independently linked to the anion having four double H bondings and two single H bondings suggesting the primary hydration number of Cl2*- to be 6. In all these clusters, the odd electron is found to be mostly localized over the two Cl atoms and these two atoms are bound by a three-electron hemibond. Calculated interaction (between solute and different water clusters) and vertical detachment energy profiles show saturation at n = 6 in the hydrated cluster Cl2*-.nH2O (n = 1-7). However, calculated solvation energy increases with the increase in number of solvent H2O molecules in the cluster. Interaction energy varies linearly with vertical detachment energy for the hydrated clusters Cl2*-.nH2O (n < or = 6). Calculation of the vibration frequencies show that the formation of Cl2*(-)-water clusters induces significant shifts from the normal stretching modes of isolated water. A clear difference in the pattern of IR spectra is observed in the O-H stretching region of water from hexa- to heptahydrated cluster.

Hydration and Dissociation of Hydrogen Fluoric Acid (HF)

The hydration and dissociation phenomena of HF(H(2)O)(n)() (n < or = 10) clusters have been studied by using both the density functional theory with the 6-311++G[sp] basis set and the Møller-Plesset second-order perturbation theory with the aug-cc-pVDZ+(2s2p/2s) basis set. The structures for n > or = 8 are first reported here. The dissociated form of the hydrogen-fluoric acid in HF(H(2)O)(n) clusters is found to be less stable at 0 K than the undissociated form until n = 10. HF may not be dissociated at 0 K solely by water molecules because the HF H bond is stronger than the OH H bond, against the expectation that the dissociated HF(H(2)O)(n) would be more stable than the undissociated one in the presence of a number of water molecules. The dissociation would be possible for only a fraction of a number of hydrated HF clusters by the Boltzmann distribution at finite temperatures. This is in sharp contrast to other hydrogen halide acids (HCl, HBr, HI) showing the dissociation phenomena at 0 K for n > or = 4. The IR spectra of dissociated and undissociated structures of HF(H(2)O)(n) are compared. The structures and binding energies of HF(H(2)O)(n) are found to be similar to those of (H(2)O)(n+1). It is interesting that HF(H(2)O)(n=5,6,10) are slightly less stable compared with other sizes of clusters, just like the fact that (H(2)O)(n=6,7,11) are slightly less stable. The present study would be useful for the experimental/spectroscopic investigation of not only the dissociation phenomena of HF but also the similarity of the HF-water clusters to the water clusters.

Dynamical Calculations of Charge-Transfer-to-Solvent Excited States of Small I-(CH3CN)n Clusters

Relaxation dynamics of photoexcited charge-transfer-to-solvent (CTTS) states for the I(-)(CH(3)CN)(n) (n = 2 and 3) clusters has been theoretically studied using electronic structure methods. First, we have calculated several lowest singlet and triplet potential energy surfaces using the multireference configuration interaction method. It was found that the character of the singlet CTTS excited-state potential surfaces is very similar to that of the triplet CTTS states. Due to a small singlet-triplet splitting, the lowest triplet potential energy surface was used as a good model to understand the dynamics of the photoexcited singlet CTTS states. We have carried out direct molecular dynamics simulations on the lowest triplet surface at the B3LYP level. When an I(-) anion is exteriorly solvated by CH(3)CN molecules, we found that the (CH(3)CN)(n)(-) anion cluster is effectively produced. In addition, when the I(-) anion is placed in the interior in I(-)(CH(3)CN)(n) clusters, photoexcitation gives an acetonitrile monomer anion plus neutral monomers. However, if the initial geometric configuration is distorted from the minimum structure, we also found that the (CH(3)CN)(2)(-) anion cluster, where an excess electron is internally trapped, is formed via I(-)(CH(3)CN)(2) + hnu --> I + (CH(3)CN)(2)(-) process.

Photofragment coincidence imaging of small I−(H2O)n clusters excited to the charge-transfer-to-solvent state

The photodissociation dynamics of small I-(H2O)n(n=2-5) clusters excited to their charge-transfer-to-solvent (CTTS) states have been studied using photofragment coincidence imaging. Upon excitation to the CTTS state, two photodissociation channels were observed. The major channel (approximately 90%) is a two-body process forming neutral I+(H2O)n photofragments, and the minor channel is a three-body process forming I+(H2O)n-1+H2O fragments. Both processes display translational energy [P(ET)] distributions peaking at ET=0 with little available energy partitioned into translation. Clusters excited to the detachment continuum rather than to the CTTS state display the same two channels with similar P(ET) distributions. The observation of similar P(ET) distributions from the two sets of experiments suggests that in the CTTS experiments, I atom loss occurs after autodetachment of the excited [I(H2O)n-]* cluster or, less probably, that the presence of the excess electron has little effect on the departing I atom.

Ab initio study of hydrated sodium halides NaX(H2O)1–6 (X=F, Cl, Br, and I)

We have studied the dissociation phenomena of sodium halides by water molecules. The structures, binding energies, electronic properties, and IR spectroscopic features have been investigated by using the density-functional theory, second-order Moller-Plesset perturbation theory, and coupled clusters theory with single, double, and perturbative triplet excitations. In the case that the sodium halides are hydrated by three water molecules, the most stable structures show the partial (or half) dissociation feature. The dissociated structures are first found for NaX(H2O)(n=5) for X=BrI, though these structures are slightly higher in energy than the global minimum-energy structure. In the case of hexahydrated sodium halides the global minimum-energy structures (which are different from the structures reported in any previous work) are found to be dissociated (X=F/I) or partially/half dissociated (X=Cl/Br), while other nearly isoenergetic structures are undissociated, and the dissociated cubical structures are higher in energy than the corresponding global minimum-energy structure.

A comparative ab initio study of Br2•− and Br2 water clusters

The work presents ab initio results on structure and electronic properties of Br2*-.nH2O(n=1-10) and Br2.nH2O(n=1-8) hydrated clusters to study the effects of an excess electron on the microhydration of the halide dimer. A nonlocal density functional, namely, Becke's half-and-half hybrid exchange-correlation functional is found to perform well on the present systems with a split valence 6-31++G(d,p) basis function. Geometry optimizations for all the clusters are carried out with several initial guess structures and without imposing any symmetry restriction. Br2*-.nH2O clusters prefer to have symmetrical double hydrogen-bonding structures. Results on Br2.nH2O(n>or=2) cluster show that the O atom of one H2O is oriented towards one Br atom and the H atom of another H2O is directed to other Br atom making Br2 to exist as Br+-Br- entity in the cluster. The binding and solvation energies are calculated for the Br2*-.nH2O and Br2.nH2O clusters. Calculations of the vibrational frequencies show that the formation of Br2*- and Br2 water clusters induces significant shifts from the normal modes of isolated water. Excited-state calculations are carried out on Br2*-.nH2O clusters following configuration interaction with single electron excitation procedure and UV-VIS absorption profiles are simulated. There is an excellent agreement between the present theoretical UV-VIS spectra of Br2*-.10H2O cluster and the reported transient optical spectra for Br2*- in aqueous solution.

Dissolution Nature of Cesium Fluoride by Water Molecules

The structures, stabilities, thermodynamic quantities, dissociation energies, infrared spectra, and electronic properties of CsF hydrated by water molecules are investigated by using density functional theory, Møller-Plesset second-order perturbation theory (MP2), coupled cluster theory with singles, doubles, and perturbative triples excitations (CCSD(T)), and ab initio molecular dynamic (AIMD) simulations. It is revealed that at 0 K three water molecules (as a global minimum structure) begin to half-dissociate the Cs-F, and six water molecules (though not a global minimum energy structure) can dissociate it. By the combination of the accurate CCSD(T) conformational energies for Cs(H2O)6 at 0 K with the AIMD thermal energy contribution, it reveals that the half-dissociated structure is the most stable at 0 K, but this structure (which is still the most stable) changes to the dissociated structure above 50 K. The spectra of CsF(H2O)(1-6) from MP2 calculations and the power spectra of CsF(H2O)6 from 50 and 100 K AIMD simulations are also reported.

Ab Initio Molecular Dynamics Simulations of an Excited State of X-(H2O)3 (X = Cl, I) Complex

Upon excitation of Cl(-)(H(2)O)(3) and I(-)(H(2)O)(3) clusters, the electron transfers from the anionic precursor to the solvent, and then the excess electron is stabilized by polar solvent molecules. This process has been investigated using ab initio molecular dynamics (AIMD) simulations of excited states of Cl(-)(H(2)O)(3) and I(-)(H(2)O)(3) clusters. The AIMD simulation results of Cl(-)(H(2)O)(3) and I(-)(H(2)O)(3) are compared, and they are found to be similar. Because the role of the halogen atom in the photoexcitation mechanism is controversial, we also carried out AIMD simulations for the ground-state bare excess electron -- water trimer [e(-)(H(2)O)(3)] at 300 K, the results of which are similar to those for the excited state of X(-)(H(2)O)(3) with zero kinetic energy at the initial excitation. This indicates that the rearrangement of the complex is closely related to that of e(-)(H(2)O)(3), whereas the role of the halide anion is not as important.

Dissolution nature of the lithium hydroxide by water molecules

The structures, stabilities, thermodynamic quantities, dissociation energies, infrared spectra, and electronic properties of LiOH hydrated by up to seven water molecules are investigated by using the density-functional theory and the Møller-Plesset second-order perturbation theory (MP2). Further accurate analysis based on the coupled-cluster theory with singles, doubles, and perturbative triples excitations agrees with the MP2 results. The Li-OH stretch mode significantly shifts with the increase of water molecules, and it eventually disappears upon dissociation. It is revealed that seven water molecules are needed for the stable dissociation of LiOH (as a completely dissociated conformation), in contrast to the cases of RbOH and CsOH which require four and three water molecules, respectively.

Why the hydration energy of Au+ is larger for the second water molecule than the first one: Skewed orbitals overlap

Using molecular-orbital analysis, we have elucidated the quantum-chemical origin of the intriguing phenomena in sequential hydration energies of the gold cation, which is known to be the most conspicuous among all transition metals. The hydration energy of Au+ with the second water molecule is found to be much larger than that with the first water molecule. Owing to the large relativistic effect of gold (i.e., significant lowering of the 6s orbital energy and significant raising of the 5d orbital energy), the highest occupied molecular orbital of the hydrated gold cation has a large portion of the 6s orbital. As the electron density of the 6s orbital populates in a large outer spherical shell far off the gold nucleus, the p orbitals (or sp hybridized lone-pair orbitals) of the water molecules are able to overlap with the outer part of the 6s orbital in the dihydrated gold cation, resulting in the unusual skewed overlap of p-6s-p orbitals (not the atom-to-atom bond overlap). No previous molecular-orbital analysis has reported this peculiar skewed orbitals overlap. Since this skewed orbitals overlap is saturated with two water molecules, this property is responsible for the low coordination number of the gold ion.

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.

Study of the stabilization energies of halide-water clusters: An application of first-principles interaction potentials based on a polarizable and flexible model

The aim of this work is to compute the stabilization energy E(stab)(n) of [X(H(2)O)(n)](-) (X identical with F, Br, and I for n=1-60) clusters from Monte Carlo simulations using first-principles ab initio potentials. Stabilization energy of [X(H(2)O)(n)](-) clusters is defined as the difference between the vertical photodeachment energy of the cluster and the electron affinity of the isolated halide. On one hand, a study about the relation between cluster structure and the E(stab)(n) value, as well as the dependence of the latter with temperature is performed, on the other hand, a test on the reliability of our recently developed first-principles halide ion-water interaction potentials is carried out. Two different approximations were applied: (1) the Koopmans' theorem and (2) calculation of the difference between the interaction energy of [X(H(2)O)(n)](-) and [X(H(2)O)(n)] clusters using the same ab initio interaction potentials. The developed methodology allows for using the same interaction potentials in the case of the ionic and neutral clusters with the proviso that the charge of the halide anion was switched off in the latter. That is, no specific parametrization of the interaction potentials to fit the magnitude under study was done. The good agreement between our predicted E(stab)(n) and experimental data allows us to validate the first-principles interaction potentials developed elsewhere and used in this study, and supports the fact that this magnitude is mainly determined by electrostatic factors, which can be described by our interaction potentials. No relation between the value of E(stab)(n) and the structure of clusters has been found. The diversity of E(stab)(n) values found for different clusters with similar interaction energy indicates the need for statistical information to properly estimate the stabilization energy of the halide anions. The effect of temperature in the prediction of the E(stab)(n) is not significant as long as it was high enough to avoid cluster trapping into local equilibrium configurations which guarantees an appropriate sampling of the configurational space. Parallel tempering method was applied in particular cases to guarantee satisfactory sampling of clusters at low temperature.

Molecular aspects of halide ion hydration: the cluster approach

This review provides a historical context for our understanding of the hydration shell surrounding halide ions and illustrates how the cluster systems can be used, in combination with theory, to elucidate the behavior of water molecules in direct contact with the anion. We discuss how vibrational predissociation spectroscopy, carried out with weakly bound argon atoms, has been employed to deduce the morphology of the small water networks attached to anions in the primary steps of hydration. We emphasize the importance of charge-transfer in the binary interaction, and discuss how this process affects the structures of the larger networks. Finally, we survey how the negatively charged water clusters (H2O)n(-) are providing a molecular-level perspective on how diffuse excess electrons interact with the water networks.

Theory of dipole-bound anions

The problem of the binding of an excess electron to polar molecules and their clusters has long fascinated researchers. Although excess electrons bound to such species tend to be very extended spatially and to have little spatial overlap with the valence electrons of the neutral molecules, inclusion of electron correlation effects is essential for quantitatively describing the electron binding. The major electron correlation contribution may be viewed as a dispersion interaction between the excess electron and the electrons of the molecule or cluster. Recent work using a one-electron Drude model to describe excess electrons interacting with polar molecules is reviewed.