Electronic structure, statistical mechanical simulations, and EXAFS spectroscopy of aqueous potassium

We investigate the solvation structure of aqueous potassium ions, using a combination of electronic structure calculations, statistical mechanical simulations with a derived polarizable empirical potential and experimental measurement of the extended X-ray absorption fine structure (EXAFS) spectra. The potassium K-edge (at 3,608 eV) EXAFS spectra were acquired on the bending magnet of sector 20 at the Advanced Photon Source, at ambient conditions and for the concentrations of 1 and 4 m KCl. We focus on the coordination distances and the degree of disorder of the first hydration shell as determined by electronic structure calculations, molecular dynamics simulations and experimental measurement. Finally, we characterize the changes of the structure in the first hydration shell with increasing temperature as predicted by molecular simulation     

Hydrogen and higher shell contributions in Zn2+, Ni2+, and Co2+ aqueous solutions: an X-ray absorption fine structure and molecular dynamics study

A detailed investigation of the hydration structure of Zn2+, Ni2+, and Co2+ in water solutions has been carried out combining X-ray absorption fine structure (EXAFS) spectroscopy and Molecular Dynamics (MD) simulations. The first quantitative analysis of EXAFS from hydrogen atoms in 3d transition metal ions in aqueous solutions has been carried out and the ion-hydrogen interactions have been found to provide a detectable contribution to the EXAFS spectra. An accurate determination of the structural parameters associated with the first hydration shell has been performed and compared with previous experimental results. No evidence of significant contributions from the second hydration shell to the EXAFS signal has been found for these solutions, while the inclusion of the hydrogen signal has been found to be important in performing a quantitative analysis of the experimental data. The high-frequency contribution present in the EXAFS spectra has been found to be due to multiple scattering (MS) effects inside the ion-oxygen first coordination shell. MD has been used to generate three-body distribution functions from which a reliable analysis of the MS contributions to the EXAFS spectra of these systems has been carried out.     

Development and validation of an integrated computational approach for the study of ionic species in solution by means of effective two-body potentials. The case of Zn2+, Ni2+, and Co2+ in aqueous solutions

In this paper we have developed an effective computational procedure for the structural and dynamical investigation of ions in aqueous solutions. Quantum mechanical potential energy surfaces for the interaction of a transition metal ion with a water molecule have been calculated taking into account the effect of bulk solvent by the polarizable continuum model (PCM). The effective ion-water interactions have been fitted by suitable analytical potentials, and have been utilized in molecular dynamics (MD) simulations to obtain structural and dynamical properties of the ionic aqueous solutions. This procedure has been successfully applied to the Co2+-H2O open-shell system and, for the first time, Co-oxygen and Co-hydrogen pair potential functions have been determined and employed in MD simulations. The reliability of the whole procedure has been assessed by applying it also to the Zn2+ and Ni2+ aqueous solutions, and the structural and dynamical properties of the three systems have been calculated by means of MD simulations and have been found to be in very good agreement with experimental results. The structural parameters of the first solvation shells issuing from the MD simulations provide an effective complement to extended X-ray absorption fine structure (EXAFS) experiments.     

Hydration of the calcium ion. An EXAFS, large-angle x-ray scattering, and molecular dynamics simulation study.

The structure of the hydrated calcium(II) ion in aqueous solution has been studied by means of extended X-ray absorption fine structure spectroscopy (EXAFS), large-angle X-ray scattering (LAXS), and molecular dynamics (MD) methods. The EXAFS data displayed a broad and asymmetric distribution of the Ca-O bond distances with the centroid at 2.46(2) A. LAXS studies on four aqueous calcium halide solutions (1.5-2 mol dm(-)(3)) gave a mean Ca-O bond distance of 2.46(1) A. This is consistent with a hydration number of 8 determined from correlations between mean distances and coordination numbers from crystal structures. The LAXS studies showed a second coordination sphere with a mean Ca.O(II) distance of 4.58(5) A, and for the hydrated halide ions the distances Cl.O 3.25(1) A, Br.O 3.36(1) A, and I.O 3.61(1) A were obtained. Molecular dynamics simulations of CaCl(2)(aq) were performed using three different Ca(2+)-OH(2) pair potentials. The potential from the GROMOS program gave results in agreement with experiments, i.e., a coordination number of 8 and an average Ca-O distance of 2.46 A, and was used for further comparisons. Theoretical EXAFS oscillations were computed for individual MD snapshots and showed very large variations, though the simulated average spectrum from 2000 snapshots gave satisfactory agreement with the experimental EXAFS spectra. The effect of thermal motions of the coordinated atoms is inherent in the MD simulation method. Thermal disorder parameters evaluated from simulated spatial atom distribution functions of the oxygen atoms coordinated to the calcium ion were in close agreement with those from the current LAXS and EXAFS analyses. The combined results are consistent with a root-mean-square displacement from the mean Ca-O distance of 0.09(2) A in aqueous solution at 300 K.     

Protonation and complex formation of 5-sulfosalicylate in NaCl,CaCl2 and MgCl2 aqueous media. Speciation in synthetic seawater

The binding capacity of 5-sulfosalicylic acid (ssa) towards cationic macro-components of natural waters has been investigated in different ionic media (NaCl, MgCl2 and CaCl2 aqueous solutions) and in the ionic strength range 0 < or = I < or = 1 mol dm-3. In order to contribute to the speciation of this multi-sites ligand, measurements have been carried out also in a synthetic seawater (SSWE) containing the major components of seawater (Na+, K+, Ca2+, Mg2+, Cl- and SO4(2-)). Measurements have been performed by potentiometry ([H+]-glass electrode), at t = 25 degrees C. A critical analysis on the experimental and literature data is also given.     

Solvation of uranyl(II), europium(III) and europium(II) cations in "basic" room-temperature ionic liquids: a theoretical study

We report a molecular dynamics study of the solvation of UO2(2+), Eu3+ and Eu2+ ions in two "basic" (Lewis acidity) room-temperature ionic liquids (IL) composed of the 1-ethyl-3-methylimidazolium cation (EMI+) and a mixture of AlCl4- and Cl- anions, in which the Cl-/AlCl4- ratio is about 1 and 3, respectively. The study reveals the importance of the [UO2Cl4]2- species, which spontaneously form during most simulations, and that the first solvation shell of europium is filled with Cl- and AlCl4- ions embedded in a cationic EMI+ shell. The stability of the [UO2Cl4]2- and [Eu(III)Cl6]3- complexes is supported by quantum mechanical calculations, according to which the uranyl and europium cations intrinsically prefer Cl- to the AlCl4- ion. In the gas phase, however, [Eu(III)Cl6]3- and [Eu(II)Cl6]4- complexes are predicted to be metastable and to lose two to three Cl- ions. This contrasts with the results of simulations of complexes in ILs, in which the "solvation" of the europium complexes increases with the number of coordinated chlorides, leading to an equilibrium between different chloro species. The behavior of the hydrated [Eu(OH2)8]3+ complex is considered in the basic liquids; the complex exchanges H2O molecules with Cl- ions to form mixed [EuCl3(OH2)4] and [EuCl4(OH2)3]- complexes. The results of the simulations allow us to better understand the microscopic nature and solvation of lanthanide and actinide complexes in "basic" ionic liquids.     

Halide anion recognition in water by a hexaprotonated octaaza-cryptand: a molecular dynamics investigation

Based on molecular dynamics simulations, we describe the F- versus Cl- complexation by an hexaprotonated cryptand L6+ in aqueous solution, in order to elucidate their structures, solvation properties and the status of external halide counterions. In water, F- and Cl- simulated inclusive complexes adopt a structure somewhat different from the solid state structure of the F- complex: The anion binding involves two diammonium bridges only, and the accompanying counterions are dissociated from the +5 charged complex. A remarkable result is obtained for the dissociated L6+,3F- ,3Cl- system, where spontaneous complexation of F- (the anion which forms the most stable complex with L6+) takes place during the dynamics. The resulting complex is of facial type; this suggests that the equilibrium involves multiple binding modes and structures in aqueous solution. The question of F-/Cl- binding selectivity is investigated by free energy perturbations simulations which nicely reproduce the spectacular preference for F- over Cl-. Two different methodologies used for the treatment of electrostatics (standard versus Ewald calculations) yield similar conclusions.     

Modeling of ionic equilibria of trace metals (Cu2+, Zn2+, Cd2+) in concentrated aqueous electrolyte solutions at 25 degrees C

This work presents a model of activity coefficients and a database for ionic equilibria of Cu2+, Zn2+, and Cd2+ in (H+, Na+, K+, Mg2+, Ca2+)(OH-, Cl-, NO-3, ClO-4, HSO-4, SO2-(4), HCO-3, CO(2-)3) aqueous media valid up to 6-12 m ionic strength. The activity coefficient of a dissolved species is represented by empirical equation [Formula: see text] , where Agamma is Debye-Hückel constant (1.17 at 25 degrees C), gammai and zi are activity coefficient and charge of a dissolved species i, I is molal ionic strength, bij is model parameter, and mj is molal concentration of dissolved species. The model is applicable to the modeling of ionic equilibria, as well as to simulation of solubility of salts in mixed electrolyte solutions.     

Crystal lattice properties fully determine short-range interaction parameters for alkali and halide ions

Accurate models of alkali and halide ions in aqueous solution are necessary for computer simulations of a broad variety of systems. Previous efforts to develop ion force fields have generally focused on reproducing experimental measurements of aqueous solution properties such as hydration free energies and ion-water distribution functions. This dependency limits transferability of the resulting parameters because of the variety and known limitations of water models. We present a solvent-independent approach to calibrating ion parameters based exclusively on crystal lattice properties. Our procedure relies on minimization of lattice sums to calculate lattice energies and interionic distances instead of equilibrium ensemble simulations of dense fluids. The gain in computational efficiency enables simultaneous optimization of all parameters for Li+, Na+, K+, Rb+, Cs+, F-, Cl-, Br-, and I- subject to constraints that enforce consistency with periodic table trends. We demonstrate the method by presenting lattice-derived parameters for the primitive model and the Lennard-Jones model with Lorentz-Berthelot mixing rules. The resulting parameters successfully reproduce the lattice properties used to derive them and are free from the influence of any water model. To assess the transferability of the Lennard-Jones parameters to aqueous systems, we used them to estimate hydration free energies and found that the results were in quantitative agreement with experimentally measured values. These lattice-derived parameters are applicable in simulations where coupling of ion parameters to a particular solvent model is undesirable. The simplicity and low computational demands of the calibration procedure make it suitable for parametrization of crystallizable ions in a variety of force fields.     

Structure and Raman spectrum of clavulanic acid in aqueous solution

The calculation of the vibrational Raman spectrum of enzyme-bound beta-lactamase inhibitors may be of help to understand the mechanisms responsible for bacterial drug resistance. Here, we present a study of the solvation structure and the vibrational properties of clavulanate, an important beta-lactamase inhibitor, in aqueous solution as obtained from full quantum and hybrid empirical/quantum molecular dynamics simulations at ambient conditions. The analysis of the vibrational density of states indicates that hybrid empirical/quantum mechanical simulations are able to properly describe the vibrational levels of clavulanate in solution. In addition, we propose a computationally efficient protocol to calculate the vibrational Raman effect for large solute molecules in water, which is able to faithfully reproduce the experimentally recorded clavulanate Raman spectrum and discloses the possibility to employ hybrid simulations to assign the experimental Raman spectra of inhibitors bound to beta-lactamases.     

Interaction of water with the model ionic liquid [bmim][BF4]: molecular dynamics simulations and comparison with NMR data

We report on molecular dynamics simulations of the ionic liquid [bmim][BF 4] and its mixtures with water, from zero up to 0.5 mol fraction of water. All of the simulations are carried out with two published force fields. The results are compared with each other and with published as well as new NMR data on the same mixtures, whenever possible. We perform extensive analyses of structural quantities, such as pair correlation functions, nearest-neighbor analysis and size distribution of the water clusters formed at higher concentrations. We show that the water clusters are formed almost exclusively by linear chains of hydrogen-bonded molecules. There is a nanoscale structuring of the mixtures but no macroscopic phase separation among the components, in agreement with experiment. Roughly, we identify two solvation regimes. At low water content, the ions are selectively coordinated by individual water molecules, but their ionic network is largely unperturbed. At high water content, the ionic network is somewhat disrupted or swollen in a nonspecific way by the water clusters.     

The effects of charge transfer on the aqueous solvation of ions

Ab initio-based charge partitioning of ionic systems results in ions with non-integer charges. This charge-transfer (CT) effect alters both short- and long-range interactions. Until recently, the effects of CT have been mostly neglected in molecular dynamics (MD) simulations. The method presented in this paper for including charge transfer between ions and water is consistent with ab initio charge partitioning and does not add significant time to the simulation. The ions of sodium, potassium, and chloride are parameterized to reproduce dimer properties and aqueous structures. The average charges of the ions from MD simulations (0.900, 0.919, and -0.775 for Na(+), K(+), and Cl(-), respectively) are consistent with quantum calculations. The hydration free energies calculated for these ions are in agreement with experimental estimates, which shows that the interactions are described accurately. The ions also have diffusion constants in good agreement with experiment. Inclusion of CT results in interesting properties for the waters in the first solvation shell of the ions. For all ions studied, the first shell waters acquire a partial negative charge, due to the difference between water-water and water-ion charge-transfer amounts. CT also reduces asymmetry in the solvation shell of the chloride anion, which could have important consequences for the behavior of chloride near the air-water interface.     

Effective interaction potentials for alkali and alkaline earth metal ions in SPC/E water and prediction of mean ion activity coefficients

The potential of mean force (PMF) acting between two simple ions surrounded by SPC/E water have been determined by molecular dynamics (MD) simulations using a spherical cavity approach. Such effective ion-ion potentials were obtained for Me-Me, Me-Cl-, and Cl(-)-Cl- pairs, where Me is a Li+, Na+, K+, Mg2+, Ca2+, Sr2+, and Ba2+ cation. The ionic sizes estimated from the effective potentials are not pairwise additive, a feature in the frequently used primitive model for electrolytes. The effective potentials were used in Monte Carlo (MC) simulations with implicit water to calculate mean ion activity coefficients of LiCl, NaCl, KCl, MgCl2, CaCl2, SrCl2, and BaCl2. Predicted activities were compared with experimental ones in the electrolyte concentration range 0.1-1 M. A qualitative agreement for LiCl and a satisfactory agreement for NaCl were found, whereas the predictions for KCl by two K+ models were less coherent. In the case of alkaline earth metal ions, all experimental activities were successfully reproduced at c = 0.1 M. However, at higher concentrations, similar deviations occurred for all divalent cations, suggesting that the dependence of the permittivity on the salt concentration and the polarization deficiency arising from the ordering of water molecules in the ion hydration shells are important in such systems.     

Electrolytes in a nanometer slab-confinement: ion-specific structure and solvation forces

We study the liquid structure and solvation forces of dense monovalent electrolytes (LiCl, NaCl, CsCl, and NaI) in a nanometer slab-confinement by explicit-water molecular dynamics (MD) simulations, implicit-water Monte Carlo (MC) simulations, and modified Poisson-Boltzmann (PB) theories. In order to consistently coarse-grain and to account for specific hydration effects in the implicit methods, realistic ion-ion and ion-surface pair potentials have been derived from infinite-dilution MD simulations. The electrolyte structure calculated from MC simulations is in good agreement with the corresponding MD simulations, thereby validating the coarse-graining approach. The agreement improves if a realistic, MD-derived dielectric constant is employed, which partially corrects for (water-mediated) many-body effects. Further analysis of the ionic structure and solvation pressure demonstrates that nonlocal extensions to PB (NPB) perform well for a wide parameter range when compared to MC simulations, whereas all local extensions mostly fail. A Barker-Henderson mapping of the ions onto a charged, asymmetric, and nonadditive binary hard-sphere mixture shows that the strength of structural correlations is strongly related to the magnitude and sign of the salt-specific nonadditivity. Furthermore, a grand canonical NPB analysis shows that the Donnan effect is dominated by steric correlations, whereas solvation forces and overcharging effects are mainly governed by ion-surface interactions. However, steric corrections to solvation forces are strongly repulsive for high concentrations and low surface charges, while overcharging can also be triggered by steric interactions in strongly correlated systems. Generally, we find that ion-surface and ion-ion correlations are strongly coupled and that coarse-grained methods should include both, the latter nonlocally and nonadditive (as given by our specific ionic diameters), when studying electrolytes in highly inhomogeneous situations.     

First principles and experimental 1H NMR signatures of solvated ions: The case of HCl(aq).

A combined experimental and ab initio study is presented of the 1H NMR chemical shift distribution of aqueous hydrogen chloride solution as a function of acid concentration, based on Car-Parrinello molecular dynamics simulations and fully periodic NMR chemical-shift calculations. The agreement of computed and experimental spectra is very good. From first-principles calculations, we can show that the individual contributions of Eigen and Zundel ions, regular water molecules, and the chlorine solvation shell to the NMR line are very distinct and almost independent of the acid concentration. From the computed instantaneous NMR distributions, it is further possible to characterize the average variation in hydrogen-bond strength of the different complexes.     

Structural characterization of zinc(II) chloride in aqueous solution and in the protic ionic liquid ethyl ammonium nitrate by x-ray absorption spectroscopy

Extended x-ray absorption fine structure (EXAFS) spectroscopy has been used to investigate the species and structures existing in a series of ZnCl(2)-H(2)O-NaCl solutions with different chloride/zinc ratios and in a solution of ZnCl(2) in the protic ionic liquid ethyl ammonium nitrate (EAN). The average coordination numbers and distances of zinc species were determined from the analysis of the EXAFS data. In aqueous solution the number of chloride ions tightly bounded to Zn(2+) is significantly related to the chloride/zinc ratio, and no inner complex formation between Zn(2+) and Cl(-) ions has been detected for low ZnCl(2) concentration (0.1 and 0.2 M). Conversely, in the same concentration range (0.13 M) the ZnCl(2) species do not dissociate in EAN and the Zn(2+) first coordination shell has two chloride ions and is completed by two oxygen atoms of the nitrate anion. The results of this investigation show that notwithstanding the existence of similar characteristics between EAN and water, the solvation properties of the two solvents are markedly different.