Structures, Energetics, and Spectra of Cl-(H2O)n Clusters, n = 1−6: Ab Initio Study

Evaluation of Ar tagging toward the vibrational spectra and zero point energy of X-HOH, X-DOH, and X-HOD, for X = F, Cl, Br

In this study, we theoretically evaluated the effect of argon tagging toward the binding energy and vibrational spectra of water halide anion complexes Ar.X-HOH, Ar.X-HOD, and Ar.X-DOH (X = F, Cl, Br). The ionic hydrogen bond (IHB) OH stretching mode was calculated to have a strong peak in the vibrational spectra, and coupling to intermolecular modes as well as bending modes was observed as combination bands and Fermi resonances. We found that the argon tagging affected the IHB OH stretching peak position in Ar.F-H2O, but not in Ar.Cl-H2O and Ar.Br-H2O. Furthermore, D-binding is favored for Cl and Br based on zero point energies, but for F our calculated zero point energies did not show a preference between H- and D-binding. We show that the competition of the energy lowering in the zero point energy of the anharmonic IHB OH (OD) stretching mode versus the intermolecular out-of-plane IHB OH (OD) wagging mode is important for determining the preference between H- and D-binding for these monohydrated halide clusters. We also found that for X-HOD the HOD bending fundamental peak is blue shifted compared to bare HOD, but is redshifted for F-DOH.

Hydrogen bond network structures of protonated short-chain alcohol clusters

Protonated alcohol clusters enable extraction of the physical essence of the nature of hydrogen bond networks.

A Correction Method for Attenuated Total Reflection-Far Ultraviolet Spectra Via the Use of Charge Transfer to Solvent Band Intensities of Iodide in the Ultraviolet Region

Attenuated total reflection (ATR) spectra, which are often used in IR analysis, can be transformed into extinction and refraction spectra by Kramers-Kronig transformation (KKT) with Fresnel equations. However, it is often difficult to obtain correct optical indices due to the inherent instrumental functions. This paper proposes a simple practical method for correction of KKT with two parameters, which include all the effects of the instrumental function. In order to obtain the parameters of the instrumental function, absorption ratios of charge transfer to solvent (CTTS) transitions of aqueous iodide ions observed at 195 nm and 230 nm were used as a standard. The absorption indices calculated from the ATR spectra with the parameters correspond reasonably well to those given by the transmittance spectra not only in the UV region but also in the far-ultraviolet (FUV, 120-200 nm) region. By applying the corrected KKT to the ATR-FUV spectra of aqueous potassium halide solutions in the range of 0-2 M, correct features of the absorption spectra of KCl and KBr, whose CTTS bands are thought to be observed in FUV region, were confirmed. It is possible to use the parameters representing the instrument function as long as the instrument is not changed.

Organic matrix-induced formation of a discrete cyclic [Cl2(H2O)2]2− cluster

Reaction of 6-NH₂-Py-2-NHC(S)NHP(S)(OiPr)₂ with DMFA·2HCl leads to the co-crystallization product of 2,6-di(4H-1,2,4-triazol-4-yl)pyridine and 2,6-diaminopyridinium chloride monohydrate (1) and N-(diisopropoxythiophosphoryl)-N′-(6-(4H-1,2,4-triazol-4-yl))pyridin-2-ylthiocarbamide (2).

Theoretical calculation of the vibrational state dependent photodetachment spectra of X-H2O, X = F, Cl, Br

Vibrational spectra of X^-H₂O (X = F, Cl, Br) were simulated using full dimensional vibrational calculations using quantum chemistry potential energy surfaces. Furthermore, utilizing the reflection approximation, we simulated the photodetachment spectra obtained from different vibrational excited states. From these spectra, we can observe changes in the hydrogen bond interaction between the anion and the neutral XH₂O system. Notably, for F^-H₂O, the excitation of the ionic hydrogen bonded (IHB) OH stretching vibration generates a large tail on the low energy side of the photodetachment spectra compared to the detachment from the zero-point vibration state. This shows that the IHB OH stretching vibration of F^-H₂O causes charge delocalization from F^- to the oxygen atom in H₂O, and that the photodetachment from FHOH^- occurs at lower energies.

Hydrogen-Bonded Networks in Hydride Water Clusters, F-(H2O)n and Cl-(H2O)n: Cubic Form of F-(H2O)7 and Cl-(H2O)7

The anion-water bonds and hydrogen bonds between water molecules in X(-)(H(2)O)(n) (X = F and Cl, n = 3-7) clusters are analyzed by evaluating the charge-transfer (CT) and dispersion terms for every pair of ions and molecules with the perturbation theory based on the locally projected molecular orbitals. In particular, the relative stabilities and the bond strengths in all 11 distinct cubic X(-)(H(2)O)(7) isomers are analyzed by classifying the ligand water (L) with the numbers of the donating (n) and accepting (m) OHs as LD(n)A(m). The number of LD(0)A(2) waters determines the relative stability. It is demonstrated that the strengths of the anion-ligand bonds are strongly influenced by two other hydrogen bonds of the water molecules adjacent to the ligand. When the model theory of Mulliken's charge-transfer interaction is applied to the anion-ligand and water-water hydrogen bonds, the dependence of the bond strengths on the chains of the hydrogen bonds is explained.

Interaction of 2,4,6-tris(2-pyrimidyl)-1,3,5-triazine (TPymT) with CoX2(X = Cl, Br) in water: trapping of new self-assembled water–chloride/bromide clusters in a [Co(bpca)2]+host (bpca = bis(2-pyrimidylcarbonyl)amidate anion)

New self-assembled water–chloride/bromide clusters are trapped as counterions in a cobalt-based cationic host.

Temperature dependence of ion mobility signals of halogenated compounds

Ion mobility spectrometry (IMS) as handheld and transportable sensor technique permits the sensitive detection of halogenated compounds with importance in environmental monitoring and process control. The negative ion mobility spectra mostly show one product ion peak which can be attributed to (H(2)O)(n)X(-) ions due to dissociative electron attachments. For minimizing memory effects and contaminations, modern ion mobility spectrometers work at elevated temperatures. In this paper, we investigated the influence of temperature on peak position, resolution and relative abundance of ions formed from halogenated substances. Elevated temperatures affect the peak position in different way. For fluorine- and chlorine-containing product ions, changes in hydration and clustering have a considerable influence on peak position, while these processes are of minor importance for bromine- and iodine-containing product ions. In these cases, the drift time differences mainly result from differences in drift behavior due to differences in gas density, the mean free path of ions and different collision rates. The drift time shift with elevated temperatures provides an enhanced peak-to-peak resolution. Improved separation efficiency can therefore be established with increased temperatures for negative product ions of halogenated compounds. Furthermore, an enhanced sensitivity was found for all compounds with increasing temperatures. However, independent on the temperature, the order of sensitivity is mainly determined by the bonding state of halogen atoms within the molecules.

The calculated infrared spectrum of Cl-H2O using a new full dimensional ab initio potential surface and dipole moment surface

We report a full dimensional, ab initio-based global potential energy surface (PES) and dipole moment surface for Cl-H2O. Both surfaces are symmetric with respect to interchange of the H atoms. The PES is a fit to thousands of electronic energies calculated using the coupled-cluster method [CCSD(T)] with a moderately large basis (aug-cc-pVTZ). Vibrational energies and wave functions are accurately obtained using MULTIMODE. The wave function and dipole moment surface are used to calculate and analyze the pure infrared spectrum at 0 K which is compared with experiment. Vibrational energies and the infrared spectra for DOD and HOD/DOH are also presented.

Absolute hydration free energy scale for alkali and halide ions established from simulations with a polarizable force field

A polarizable potential function for the hydration of alkali and halide ions is developed on the basis of the recent SWM4-DP water model [Lamoureux, G.; MacKerell, A. D., Jr.; Roux, B. J. Chem. Phys. 2003, 119, 5185]. Induced polarization is incorporated using classical Drude oscillators that are treated as auxiliary dynamical degrees of freedom. The ions are represented as polarizable Lennard-Jones centers, whose parameters are optimized to reproduce the binding energies of gas-phase monohydrates and the hydration free energies in the bulk liquid. Systematic exploration of the parameters shows that the monohydrate binding energies can be consistent with a unique hydration free energy scale if the computed hydration free energies incorporate the contribution from the air/water interfacial electrostatic potential (-540 mV for SWM4-DP). The final model, which can satisfyingly reproduce both gas and bulk-phase properties, corresponds to an absolute scale in which the intrinsic hydration free energy of the proton is -247 kcal/mol.

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.

Molecular Simulation Analysis and X-ray Absorption Measurement of Ca2+, K+and Cl-Ions in Solution

This paper presents recent advances in the use of molecular simulations and extended X-ray absorption fine structure (EXAFS) spectroscopy, which enable us to understand solvated ions in solution. We report and discuss the EXAFS spectra and related properties governing solvation processes of different ions in water and methanol. Molecular dynamics (MD) trajectories are coupled to electron scattering simulations to generate the MD-EXAFS spectra, which are found to be in very good agreement with the corresponding experimental measurements. From these simulated spectra, the ion-oxygen distances for the first hydration shell are in agreement with experiment within 0.05-0.1 A. The ionic species studied range from monovalent to divalent, positive and negative: K+, Ca2+, and Cl-. This work demonstrates that the combination of MD-EXAFS and the corresponding experimental measurement provides a powerful tool in the analysis of the solvation structure of aqueous ionic solutions. We also investigate the value of electronic structure analysis of small aqueous clusters as a benchmark to the empirical potentials. In a novel computational approach, we determine the Debye-Waller factors for Ca2+, K+, and Cl- in water by combining the harmonic analysis of data obtained from electronic structure calculations on finite ion-water clusters, providing excellent agreement with the experimental values, and discuss how they compare with results from a harmonic classical statistical mechanical analysis of an empirical potential.

Experimental infrared spectra of Cl−(ROH) (R = H, CH3, CH3CH2) complexes in the gas-phase

Infrared multiple photon dissociation spectra for the chloride ion solvated by either water, methanol or ethanol have been recorded using an FTICR spectrometer coupled to a free-electron laser, and are presented here along with assignments to the observed bands. The assignments made to the Cl(-)/H(2)O, Cl(-)(CH(3)OH), and Cl(-)(CH(3)CH(2)OH) spectra are based on comparison with the neutral H(2)O, CH(3)OH, and CH(3)CH(2)OH spectra, respectively. This work confirms that a band observed around 1400 cm(-1) in the Cl(-)(H(2)O) spectrum is not due to the Ar tag in Ar predissociation spectra. The carrier of this band is, most likely, the first overtone of the OHCl bend. Based on the position of the overtone in the IRMPD spectrum, 1375 cm(-1), the fundamental must occur very close to 700 cm(-1) and observation of this band should aid theoretical treatments of the spectrum of this complex. B3LYP/6-311++G(2df,2pd) calculations are shown to reproduce the IRMPD spectra of all three solvated chloride species. They also predict that attaching one or two Ar atoms to the Cl(-)(H(2)O) complex results in a shift of no more than a few wavenumbers in the fundamental bands for the bare complex, in agreement with previous experiment. For both alcohol-Cl(-) complexes, the S(N)2 "backside attack" isomers are not observed and Cl(-) is predicted theoretically, and confirmed experimentally, to be bound to the hydroxyl hydrogen. For Cl(-)(CH(3)CH(2)OH), the trans and gauche conformers are similar in energy, with the gauche conformer predicted to be thermodynamically favoured. The experimental infrared spectrum agrees well with that predicted for the gauche conformer but a mixture of gauche and anti conformers cannot be ruled out based on the experimental spectra nor on the computed thermochemistry.

Microsolvation of F- in Water

Microsolvation of F- in water is studied by ionization and double ionization spectra of (H2O)1-3F- calculated by ab initio methods. It is shown that the presence of the fluorine electrons introduces many-body properties in the spectra which cannot be reproduced by the presence of a negative point charge. The increase of the solvation shell increases the complexity in particular of the double ionization spectra. Ionization and double ionization energies slowly increase with continued solvation, and many-body effects in the inner valence spectra become more prominent.

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.

A first principles study on the solvation and structure of SO[sub 4][sup 2−](H[sub 2]O)[sub n], n=6–12

The structures of hydrated sulfate clusters, SO(4) (2-)(H(2)O)(n) with n=6-12, are obtained by density functional theory calculations. For SO(4) (2-)(H(2)O)(12), two structures with symmetric distribution of H(2)O molecules around the sulfate group are favored in energy. The structures for the smaller clusters, SO(4) (2-)(H(2)O)(n) with n=6-11, are obtained by taking away one H(2)O molecule successively from the two symmetric SO(4) (2-)(H(2)O)(12) isomers. The hydrogen bonding between the sulfate O atoms and H(2)O molecules are strong. So are the hydrogen bonds among H(2)O molecules, which are facilitated by the structure of the polyatomic sulfate group. The solvation energy is quite large (often exceeding 15 kcal/mol). The patterns for structural and energy changes as the cluster size increases are very different from the well studied hydrated halide ions, although the competition between solute-solvent and solvent-solvent interactions is again an important factor. Ab initio molecular dynamics simulations also show "crowding" effects in the first solvation of SO(4) (2-)(H(2)O)(12) at raised temperature.

Ion Solvation Thermodynamics from Simulation with a Polarizable Force Field

Thermodynamic measurements of the solvation of salts and electrolytes are relatively straightforward, but it is not possible to separate total solvation free energies into distinct cation and anion contributions without reference to an additional extrathermodynamic assumption. The present work attempts to resolve this difficulty using molecular dynamics simulations with the AMOEBA polarizable force field and perturbation techniques to directly compute absolute solvation free energies for potassium, sodium, and chloride ions in liquid water and formamide. Corresponding calculations are also performed with two widely used nonpolarizable force fields. The simulations with the polarizable force field accurately reproduce in vacuo quantum mechanical results, experimental ion-cluster solvation enthalpies, and experimental solvation free energies for whole salts, while the other force fields do not. The results indicate that calculations with a polarizable force field can capture the thermodynamics of ion solvation and that the solvation free energies of the individual ions differ by several kilocalories from commonly cited values.

An environmental pseudopotential approach to molecular interactions: Implementation in MOLPRO

We present the implementation into the MOLPRO package of a model for the interaction of a central system with its surrounding environment. The properties of a target system enclosed by a noncovalently bound environment or solvent are modeled as those of a system embedded into the effective pseudopotential arising from the exact electrostatic Coulomb potential and the approximated exchange-repulsion potential. For the latter we use the charge-density overlap model, which relates the exchange-repulsion interaction energy between two species with the overlap of their ground-state electron charge densities. The solutions of the modified Hartree-Fock equations for the target system are obtained self-consistently. This way the exchange-induction effects arising from the converged electron-charge density of the embedded system are implicitly included. Inclusion of the correlation effects is provided by the use of post-Hartree-Fock and density-functional techniques available in the MOLPRO package. The computational and conceptual advantages provided by this approach are shown in the calculation of the dipole polarizabilities of halide and chalcogenide anions in different environments.