Hydration shell exchange dynamics during ion transfer across the liquid/liquid interface

We examine using molecular dynamics simulations the rate and mechanism of water molecules exchange around the Li(+) and Na(+) ions during ion transfer across the interface between water and nitrobenzene. As the ions are transferred from the water to the organic phase, they keep their first hydration shell and an incomplete second shell. The rate of water exchange between the first shell and the rest of the interfacial water molecule decreases during the transfer, which is consistent with an increase in the barrier along the ion-water potential of mean force. While in bulk water the exchange of water molecules around the Li(+) follows an associative (A) or associative interchange (I(a)) type mechanism, the fraction of exchange events of type A increases at the interface. In contrast, while in bulk water the exchange of water molecules around the six coordinated Na(+) hydrated species mainly follows a dissociative mechanism, the situation at the interface involves an equilibrium interchange between the four- and five-coordinated hydrated ion. Simulation of the reversed process, in which the hydrated Li(+) ion is transferred to the aqueous phase, shows the same general behavior as a function of location from the interface.     

Characteristics of CO3(2-)-water hydrogen bonds in aqueous solution: insights from HF/MM and B3LYP/MM MD simulations

Two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, namely HF/MM and B3LYP/MM, have been performed to investigate the local hydration structure and dynamics of carbonate (CO(3)(2-)) in dilute aqueous solution. With respect to the QM/MM scheme, the QM region, which contains the CO(3)(2-) and its surrounding water molecules, was treated at HF and B3LYP levels of accuracy, respectively, using the DZV+ basis set, while the rest of the system is described by classical MM potentials. For both the HF/MM and B3LYP/MM simulations, it is observed that the hydrogen bonds between CO(3)(2-) oxygens and their nearest-neighbor waters are relatively strong, i.e., compared to water-water hydrogen bonds in the bulk, and that the first shell of each CO(3)(2-) oxygen atom somewhat overlaps with the others, which allows migration of water molecules among the coordinating sites to exist. In addition, it is observed that first-shell waters are either "loosely" or "tightly" bound to the respective CO(3)(2-) oxygen atoms, leading to large fluctuations in the number of first-shell waters, ranging from 1 to 6 (HF/MM) and 2 to 7 (B3LYP/MM), with the prevalent value of 3. Upon comparing the HF and B3LYP methods in describing this hydrated ion, the latter is found to overestimate the hydrogen-bond strength in the CO(3)(2-)-water complexes, resulting in a slightly more compact hydration structure at each of the CO(3)(2-) oxygens.     

Structure and dynamics of the Zr(4+) ion in water

The four-times positively charged zirconium ion in aqueous solution was simulated, using an ab initio quantum mechanical charge field molecular dynamics approach. As no hydrolysis reaction occurred during the simulation time of 10 ps, the target of this study was the evaluation of the structure and dynamics of the monomeric hydrated zirconium(iv) ion. The ion forms three hydration shells. In the first hydration shell the ion is 8-fold coordinated with a maximum probability of the Zr-O distance at 2.25 ?. While no exchanges occurred between the first and second shell, the mean residence time of the water molecules in the second shell is 5.5 ps. A geometry of the first hydration shell in-between a bi-capped trigonal prism and a square antiprism was found and a Zr-O force constant of 188 N m(-1) was evaluated.     

Quantum mechanical charge field molecular dynamics simulation of the TiO2+ ion in aqueous solution

Structural and dynamical properties of the TiO(2+) ion in aqueous solution have been investigated by using the new ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) formalism, which does not require any other potential functions except those for solvent-solvent interactions. Both first and second hydration shell have been treated at Hartree-Fock (HF) quantum mechanical level. A Ti-O bond distance of 1.5 A was observed for the [Ti=O](2+) ion. The first hydration shell of the ion shows a varying coordination number ranging from 5 to 7, five being the dominant one and representing one axial and four equatorial water molecules directly coordinated to Ti, which are located at 2.3 A and 2.1 A, respectively. The flexibility in the coordination number reflects the fast exchange processes, which occur only at the oxo atom, where water ligands are weakly bound through hydrogen bonds. Considering the first shell hydration, the composition of the TiO(2+) hydrate can be characterized as [(H(2)O)(0.7)(H(2)O)(4) (eq)(H(2)O)(ax)](2+). The second shell consists in average of 12 water molecules located at a mean distance of 4.4 A. Several other structural parameters such as radial and angular distribution functions and coordination number distributions were analyzed to fully characterize the hydration structure of the TiO(2+) ion in aqueous solution. For the dynamics of the TiO(2+) ion, different sets of dynamical parameters such as Ti=O, Ti-O(eq), and Ti-O(ax) stretching frequencies and ligands' mean residence times were evaluated. During the simulation time of 15 ps, 3 water exchange processes in the first shell were observed at the oxo atom, corresponding to a mean residence time of 3.6 ps. The ligands' mean residence time for the second shell was determined as 3.5 ps.     

Molecular dynamics study of the solvation of calcium carbonate in water

We performed molecular dynamics simulations of diluted solutions of calcium carbonate in water. To this end, we combined and tested previous polarizable models. The carbonate anion forms long-living hydrogen bonds with water and shows an amphiphilic character, in which the water molecules are expelled in a region close to its C(3) symmetry axis. The calcium cation forms a strongly bound ion pair with the carbonate. The first hydration shell around the CaCO(3) pair is found to be very similar to the location of the water molecules surrounding CaCO(3) in ikaite, the hydrated mineral.     

Characterization of structure and dynamics of an aqueous scandium(iii) ion by an extended ab initio QM/MM molecular dynamics simulation

Hydration structure and dynamics of an aqueous Sc(iii) solution were characterized by means of an extended ab initio quantum mechanical/molecular dynamical (QM/MM) molecular dynamics simulation at Hartree-Fock level. A monocapped trigonal prismatic structure composed of seven water molecules surrounding scandium(iii) ion was proposed by the QM/MM simulation including the quantum mechanical effects for the first and second hydration shells. The mean Sc(iii)-O bond length of 2.14 ? was identified for six prism water molecules with one capping water located at around 2.26 ?, reproducing well the X-ray diffraction data. The Sc(iii)-O stretching frequency of 432 cm(-1) corresponding to a force constant of 130 N m(-1), evaluated from the enlarged QM/MM simulation, is in good agreement with the experimentally determined value of 430 cm(-1) (128 N m(-1)). Various water exchange processes in the second hydration shell of the hydrated Sc(iii) ion predict a mean ligand residence time of 7.3 ps.     

Hydration of the Calcium Dication: Direct Evidence for Second Shell Formation from Infrared Spectroscopy

Infrared laser action spectroscopy in a Fourier‐transform ion cyclotron resonance mass spectrometer is used in conjunction with ab initio calculations to investigate doubly charged, hydrated clusters of calcium formed by electrospray ionization. Six water molecules coordinate directly to the calcium dication, whereas the seventh water molecule is incorporated into a second solvation shell. Spectral features indicate the presence of multiple structures of Ca(H₂O)₇²⁺ in which outer‐shell water molecules accept either one (single acceptor) or two (double acceptor) hydrogen bonds from inner‐shell water molecules. Double‐acceptor water molecules are predominately observed in the second solvent shells of clusters containing eight or nine water molecules. Increased hydration results in spectroscopic signatures consistent with additional second‐shell water molecules, particularly the appearance of inner‐shell water molecules that donate two hydrogen bonds (double donor) to the second solvent shell. This is the first reported use of infrared spectroscopy to investigate shell structure of a hydrated multiply charged cation in the gas phase and illustrates the effectiveness of this method to probe the structures of hydrated ions.     

Structure and dynamics of Au+ ion in aqueous solution: ab initio QM/MM MD simulations

Structure and dynamics of hydrated Au(+) have been investigated by means of molecular dynamics simulations based on ab initio quantum mechanical molecular mechanical forces at Hartree-Fock level for the treatment of the first hydration shell. The outer region of the system was described using a newly constructed classical three-body corrected potential. The structure was evaluated in terms of radial and angular distribution functions and coordination number distributions. Water exchange processes between coordination shells and bulk indicate a very labile structure of the first hydration shell whose average coordination number of 4.7 is a mixture of 3-, 4-, 5-, 6-, and 7-coordinated species. Fast water exchange reactions between first and second hydration shell occur, and the second hydration shell is exceptionally large. Therefore, the mean residence time of water molecules in the first hydration shell (5.6 ps/7.5 ps for t*= 0.5 ps/2.0 ps) is shorter than that in the second shell (9.4 ps/21.2 ps for t*= 0.5 ps/2.0 ps), leading to a quite specific picture of a "structure-breaking" effect.     

Ultrafast Vibrational Dynamics and Local Interactions of Hydrated DNA

Biochemical processes occur mainly in aqueous environments, where interactions with water molecules play a key role for both the structure and function of biomolecules. Deoxyribonucleic acid (DNA), the basic carrier of genetic information, is characterized by an equilibrium double helix structure which is held together by intermolecular hydrogen bonds between base pairs and hydrated by an environment of water molecules with fluctuating hydrogen bonds. Basic vibrational motions of hydrated DNA and the fastest changes in the DNA–water interactions and hydration geometries occur in less than 1 ps. These processes can be accessed by mapping the vibrational dynamics of DNA and water in a time‐resolved way by nonlinear ultrafast vibrational spectroscopy. Recent studies provide a detailed understanding of DNA vibrations and their dynamics, and give insight into nonequilibrium properties and structures of hydrated DNA.     

Structure and dynamics of phosphate ion in aqueous solution: an ab initio QMCF MD study

A simulation of phosphate in aqueous solution was carried out employing the new QMCF MD approach which offers the possibility to investigate composite systems with the accuracy of a QMMM method but without the time consuming creation of solute-solvent potential functions. The data of the simulations give a clear picture of the hydration shells of the phosphate anion. The first shell consists of 13 water molecules and each oxygen of the phosphate forms in average three hydrogens bonds to different solvent molecules. Several structural parameters such as radial distribution functions and coordination number distributions allow to fully characterize the embedding of the highly charged phosphate ion in the solvent water. The dynamics of the hydration structure of phosphate are described by mean residence times of the solvent molecules in the first hydration shell and the water exchange rate.     

Ab initio molecular dynamics study of formate ion hydration

We perform ab initio molecular dynamics simulations of the aqueous formate ion. The mean number of water molecules in the first solvation shell, or the hydration number, of each formate oxygen is found to be consistent with recent experiments. Our ab initio pair correlation functions, however, differ significantly from many classical force field results and hybrid quantum mechanics/molecular mechanics predictions. They yield roughly one less hydrogen bond between each formate oxygen and water than force field or hybrid methods predict. Both the BLYP and PW91 exchange correlation functionals give qualitatively similar results. The time dependence of the hydration numbers are examined, and Wannier function techniques are used to analyze electronic configurations along the molecular dynamics trajectory.     

Tl(I)‐the strongest structure‐breaking metal ion in water? A quantum mechanical/molecular mechanical simulation study

Structural and dynamical properties of the Tl(I) ion in dilute aqueous solution have been investigated by ab initio quantum mechanics in combination with molecular mechanics. The first shell plus a part of the second shell were treated by quantum mechanics at Hartree-Fock level, the rest of the system was described by an ab initio constructed potential. The radial distribution functions indicate two different bond lengths (2.79 and 3.16 A) in the first hydration shell, in good agreement with large-angle X-ray scattering and extended X-ray absorption fine structure spectroscopy results. The average first shell coordination number was found as 5.9, and several other structural parameters such as coordination number distributions, angular distribution functions, and tilt- and theta-angle distributions were evaluated. The ion-ligand vibration spectrum and reorientational times were obtained via velocity auto correlation functions. The Tl-O stretching force constant is very weak with 5.0 N m(-1). During the simulation, numerous water exchange processes took place between first and second hydration shell and between second shell and bulk. The mean ligand residence times for the first and second shell were determined as 1.3 and 1.5 ps, respectively, indicating Tl(I) to be a typical "structure-breaker". The calculated hydration energy of -84 +/- 16 kcal mol(-1) agrees well with the experimental value of -81 kcal mol(-1). All data obtained for structure and dynamics of hydrated Tl(I) characterize this ion as a very special case among all monovalent metal ions, being the most potent "structure-breaker", but at the same time forming a distinct second hydration shell and thus having a far-reaching influence on the solvent structure.     

Molecular dynamics simulation of glycine zwitterion in aqueous solution

A classical molecular dynamics method was used to study the modifications of the solution structure and the properties of glycine zwitterion in aqueous solution due to the increase of glycine zwitterion concentration and the incorporation of Na(+) and Cl(-) ions to the solution. The glycine zwitterion had fundamentally a hydrophilic behavior at infinite dilution, establishing around six hydrogen bonds with the water molecules that surrounded it, which formed a strong hydration layer. Because of the increase of glycine zwitterion concentration, the hydration structure became more compact and the quantity of water molecules bound to this molecule decreased. The Na(+) ion bound to the CO(2) group of glycine, while the Cl(-) ion bound mainly to the NH(3) group of this molecule. The integration of the ions to the hydration layer of the glycine zwitterion produced modifications in the orientational correlation between atoms of glycine zwitterion and water that surrounded them and an increase of the approaches between the glycine zwitterion molecules. The incorporation of ions to the solution also produced changes in the water-water orientational correlation. Decreases of the water-water hydrogen bonds and diffusion coefficient of all molecules were observed when the glycine zwitterion concentration increased and when the ions were incorporated to the solution.     

Unraveling the Role of Water in the Stereoselective Step of Aqueous Proline‐Catalyzed Aldol Reactions

A multiscale computational study was performed with the aim of tracing the source of stereoselectivity and disclosing the role of water in the stereoselective step of propionaldehyde aldol self‐condensation catalyzed by proline amide in water, a reaction that serves as a model for aqueous organocatalytic aldol condensations. Solvent mixing and hydration behavior were assessed by classical molecular dynamics simulations, which show that the reaction between propanal and the corresponding enamine takes place in a fully hydrated environment. First‐principles molecular dynamics simulations were used to study the free‐energy profile of four possible reaction paths, each of which yields a different stereoisomer, and high‐level static first‐principles calculations were employed to characterize the transition states for microsolvated species. The first solvation shell of the oxygen atom of the electrophilic aldehyde at the transition states contains two water molecules, each of which donates one hydrogen bond to the nascent alkoxide and thereby largely stabilizes its excess electron density. The stereoselectivity originates in an extra hydrogen bond donated by the amido group of proline amide in two reaction paths.     

Hydrogen Bonding in Liquid Water and in the Hydration Shell of Salts

A near‐IR spectral study on pure water and aqueous salt solutions is used to investigate stoichiometric concentrations of different types of hydrogen‐bonded water species in liquid water and in water comprising the hydration shell of salts. Analysis of the thermodynamics of hydrogen‐bond formation signifies that hydrogen‐bond making and breaking processes are dominated by enthalpy with non‐negligible heat capacity effects, as revealed by the temperature dependence of standard molar enthalpies of hydrogen‐bond formation and from analysis of the linear enthalpy–entropy compensation effects. A generalized method is proposed for the simultaneous calculation of the spectrum of water in the hydration shell and hydration number of solutes. Resolved spectra of water in the hydration shell of different salts clearly differentiate hydrogen bonding of water in the hydration shell around cations and anions. A comparison of resolved liquid water spectra and resolved hydration‐shell spectra of ions highlights that the ordering of absorption frequencies of different kinds of hydrogen‐bonded water species is also preserved in the bound state with significant changes in band position, band width, and band intensity because of the polarization of water molecules in the vicinity of ions.     

Nature of hydration shells of a polyoxy‐anion with a large cationic centre: The case of iodate ion in water

The structural nature of the solvation shells of an iodate ion, which is known to be a polyoxy‐anion with a large cationic centre, is investigated by means of Born–Oppenheimer molecular dynamics (BOMD) simulations using BLYP and the dispersion corrected BLYP‐D3 functionals. The iodate ion is found to have two distinct solvation regions around the positively charged iodine (iodine solvation shell or ISS) and the negatively charged oxygens (oxygen solvation shell or OSS). We have looked at the spatial, orientational, and hydrogen bond distributions of water in the two solvation regions. It is found that the water orientational profile in the ISS is typical of a cation hydration shell. The hydrogen bonded structure of water in the OSS is found to be very similar to that of the bulk water structure. Thus, the iodate ion essentially behaves like a positively charged iodine ion in water as if there is no anionic part. This explains why the cationic character of the iodate ion was prominently seen in earlier studies. The arrangement of water molecules in the two solvation shells and in the intervening regions around the iodate ion is further resolved by looking at structural cross‐correlations. The electronic properties of the solvation shells are also looked at by calculating the solute–solvent orbital overlap and dipole moments of the solute and solvation shell water. We have also performed BOMD simulations of iodate ion‐water clusters at experimentally relevant conditions. The simulation results are found to be in agreement with experimental results. © 2018 Wiley Periodicals, Inc.     

Hydration of sodium, potassium, and chloride ions in solution and the concept of structure maker/breaker

Neutron diffraction data with hydrogen isotope substitution on aqueous solutions of NaCl and KCl at concentrations ranging from high dilution to near-saturation are analyzed using the Empirical Potential Structure Refinement technique. Information on both the ion hydration shells and the microscopic structure of the solvent is extracted. Apart from obvious effects due to the different radii of the three ions investigated, it is found that water molecules in the hydration shell of K+ are orientationally more disordered than those hydrating a Na+ ion and are inclined to orient their dipole moments tangentially to the hydration sphere. Cl- ions form instead hydrogen-bonded bridges with water molecules and are readily accommodated into the H-bond network of water. The results are used to show that concepts such as structure maker/breaker, largely based on thermodynamic data, are not helpful in understanding how these ions interact with water at the molecular level.     

Hydrated germanium (II): Irregular structural and dynamical properties revealed by a quantum mechanical charge field molecular dynamics study

Structural and dynamical properties of Ge (II) in aqueous solution have been investigated using the novel ab initio quantum mechanical charge field (QMCF) molecular dynamics (MD) formalism. The first and second hydration shells were treated by ab initio quantum mechanics at restricted Hartree–Fock (RHF) level using the cc‐pVDZ‐PP basis set for Ge (II) and Dunning double‐ζ plus polarization basis sets for O and H. Besides ligand exchange processes and mean ligand residence times to observe dynamics, tilt‐ and theta‐angle distributions along with an advanced structural parameter, namely radial and angular distribution functions (RAD) for different regions were also evaluated. The combined radial and angular distribution depicted through surface plot and contour map is presented to provide a detailed insight into the density distribution of water molecules around the Ge²⁺ ion. A strongly distorted hydration structure with two trigonal pyramidal substructures within the first hydration shell is observed, which demonstrates the lone‐pair influence and provides a new basis for the interpretation of the catalytic and pharmacological properties of germanium coordination compounds. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

A study of the hydration of the alkali metal ions in aqueous solution

The hydration of the alkali metal ions in aqueous solution has been studied by large angle X-ray scattering (LAXS) and double difference infrared spectroscopy (DDIR). The structures of the dimethyl sulfoxide solvated alkali metal ions in solution have been determined to support the studies in aqueous solution. The results of the LAXS and DDIR measurements show that the sodium, potassium, rubidium and cesium ions all are weakly hydrated with only a single shell of water molecules. The smaller lithium ion is more strongly hydrated, most probably with a second hydration shell present. The influence of the rubidium and cesium ions on the water structure was found to be very weak, and it was not possible to quantify this effect in a reliable way due to insufficient separation of the O-D stretching bands of partially deuterated water bound to these metal ions and the O-D stretching bands of the bulk water. Aqueous solutions of sodium, potassium and cesium iodide and cesium and lithium hydroxide have been studied by LAXS and M-O bond distances have been determined fairly accurately except for lithium. However, the number of water molecules binding to the alkali metal ions is very difficult to determine from the LAXS measurements as the number of distances and the temperature factor are strongly correlated. A thorough analysis of M-O bond distances in solid alkali metal compounds with ligands binding through oxygen has been made from available structure databases. There is relatively strong correlation between M-O bond distances and coordination numbers also for the alkali metal ions even though the M-O interactions are weak and the number of complexes of potassium, rubidium and cesium with well-defined coordination geometry is very small. The mean M-O bond distance in the hydrated sodium, potassium, rubidium and cesium ions in aqueous solution have been determined to be 2.43(2), 2.81(1), 2.98(1) and 3.07(1) ?, which corresponds to six-, seven-, eight- and eight-coordination. These coordination numbers are supported by the linear relationship of the hydration enthalpies and the M-O bond distances. This correlation indicates that the hydrated lithium ion is four-coordinate in aqueous solution. New ionic radii are proposed for four- and six-coordinate lithium(I), 0.60 and 0.79 ?, respectively, as well as for five- and six-coordinate sodium(I), 1.02 and 1.07 ?, respectively. The ionic radii for six- and seven-coordinate K(+), 1.38 and 1.46 ?, respectively, and eight-coordinate Rb(+) and Cs(+), 1.64 and 1.73 ?, respectively, are confirmed from previous studies. The M-O bond distances in dimethyl sulfoxide solvated sodium, potassium, rubidium and cesium ions in solution are very similar to those observed in aqueous solution.