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.     

Density functional theory based molecular-dynamics study of aqueous fluoride solvation

We use density functional theory based molecular-dynamics simulations to study the aqueous solvation of the fluoride anion. Our studies are focused on the first solvation shell and have resulted in detailed information on its structural and dynamical properties. The fluoride ion leads to the formation of a rigid solvation shell, qualitatively consistent with simulation and experimental studies, classifying fluoride as a "structure making" particle. However, quantitatively we find the solvation shell to be less structured and more mobile than predicted from empirical force-field simulation. The influence on the intramolecular electronical and structural properties of water is minimal, as observed for other halogens. We propose two distinct mechanisms for the exchange of bulk and first solvation shell water molecules.     

Structural and dynamic properties of water within the solvation layer around various conformations of the glycine-based polypeptide

Several conformations of the solvated glycine-based polypeptides were investigated using molecular dynamics simulations. Some properties of water in the neighboring space around these molecules were investigated. It was found that water forms a well-defined layer-the first solvation shell-around the peptide molecule, and thickness of this layer is independent of the peptide structure and is equal to approximately 0.28 nm. Within this layer, water molecules show marked orientations relative to a peptide surface. Using the two-particle contribution to entropy as a measure of structural ordering of water, we found that the first solvation shell contributes 95% or more to the total water ordering around the peptide molecule. In investigating the dynamic properties of water, diffusion coefficients and lifetime of the hydrogen bond, clear differences between solvation layer and the bulk water were observed. It was found that the translational diffusion coefficient, D(T), decreases by 30% or more compared to bulk water; also, the lifetime of the water-water hydrogen bond clearly increases. The rotational diffusion coefficient, however, decreases only slightly, no more than approximately 10%. These differences correspond to the slightly higher energy of the hydrogen bond, and to its slightly distorted geometry. Analyzing the translational dynamics of water in the vicinity of the peptide molecule, it was deduced that the structure of the first solvation shell becomes more rigid than the structure of the bulk water. Investigation of a "pure hydrophobic" form of the polypeptide shows that the structure and the properties of water within the solvation shell are predominantly determined by the hydrophobic effect. The specific interactions between water molecules and various charge groups of the peptide molecule modifies this effect only slightly.     

Be2+在水、甲醇和乙醇中结构性质的从头算分子动力学模拟

采用Car-Parrinello分子动力学(CPMD)方法分别研究了Be²⁺在水、甲醇和乙醇中的溶剂结构性质, 并对Be²⁺的第一溶剂壳结构的实验及理论结果进行了比较. 所得第一溶剂壳结构与已报道的实验和理论结果较为一致. 对径向分布函数、配位数以及角度分布等进行了详细的分析. 结果表明: 在水、甲醇和乙醇中, Be²⁺第一溶剂壳为稳定理想的四面体结构. 在本文的模拟时间尺度内,没有观察到第一溶剂壳中的分子与第二溶剂壳中的分子进行交换, 进一步证明Be²⁺第一溶剂壳为稳定的四配位结构. 根据计算得到的空间分布函数, Be²⁺在溶剂分子的等高面上主要集中分布在溶剂分子接受氢键的方向. 根据氧原子在Be²⁺周围的分布, 壳层分子主要集中分布在Be²⁺周围的四个区域, 进一步证实溶剂壳为四面体对称.     

A first principles molecular dynamics study of lithium atom solvation in binary liquid mixture of water and ammonia: structural, electronic, and dynamical properties

The preferential solvation of solutes in mixed solvent systems is an interesting phenomenon that plays important roles in solubility and kinetics. In the present study, solvation of a lithium atom in aqueous ammonia solution has been investigated from first principles molecular dynamics simulations. Solvation of alkali metal atoms, like lithium, in aqueous and ammonia media is particularly interesting because the alkali metal atoms release their valence electrons in these media so as to produce solvated electrons and metal counterions. In the present work, first principles simulations are performed employing the Car-Parrinello molecular dynamics method. Spontaneous ionization of the Li atom is found to occur in the mixed solvent system. From the radial distribution functions, it is found that the Li(+) ion is preferentially solvated by water and the coordination number is mostly four in its first solvation shell and exchange of water molecules between the first and second solvation shells is essentially negligible in the time scale of our simulations. The Li(+) ion and the unbound electron are well separated and screened by the polar solvent molecules. Also the unbound electron is primarily captured by the hydrogens of water molecules. The diffusion rates of Li(+) ion and water molecules in its first solvation shell are found to be rather slow. In the bulk phase, the diffusion of water is found to be slower than that of ammonia molecules because of strong ammonia-water hydrogen bonds that participate in solvating ammonia molecules in the mixture. The ratio of first and second rank orientational correlation functions deviate from 3, which suggests a deviation from the ideal Debye-type orientational diffusion. It is found that the hydrogen bond lifetimes of ammonia-ammonia pairs is very short. However, ammonia-water H-bonds are found to be quite strong when ammonia acts as an acceptor and these hydrogen bonds are found to live longer than even water-water hydrogen bonds.     

Probing the transition from hydrophilic to hydrophobic solvation with atomic scale resolution

Picosecond and femtosecond X-ray absorption spectroscopy is used to probe the changes of the solvent shell structure upon electron abstraction of aqueous iodide using an ultrashort laser pulse. The transient L(1,3) edge EXAFS at 50 ps time delay points to the formation of an expanded water cavity around the iodine atom, in good agreement with classical and quantum mechanical/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. These also show that while the hydrogen atoms pointed toward iodide, they predominantly point toward the bulk solvent in the case of iodine, suggesting a hydrophobic behavior. This is further confirmed by quantum chemical (QC) calculations of I(-)/I(0)(H(2)O)(n=1-4) clusters. The L(1) edge sub-picosecond spectra point to the existence of a transient species that is not present at 50 ps. The QC calculations and the QM/MM MD simulations identify this transient species as an I(0)(OH(2)) complex inside the cavity. The simulations show that upon electron abstraction most of the water molecules move away from iodine, while one comes closer to form the complex that lives for 3-4 ps. This time is governed by the reorganization of the main solvation shell, basically the time it takes for the water molecules to reform an H-bond network. Only then is the interaction with the solvation shell strong enough to pull the water molecule of the complex toward the bulk solvent. Overall, much of the behavior at early times is determined by the reorientational dynamics of water molecules and the formation of a complete network of hydrogen bonded molecules in the first solvation shell.     

Solvation structure of hydroxyl radical by Car-Parrinello molecular dynamics

Car-Parrinello molecular dynamics simulations of a hydroxyl radical in liquid water have been performed. Structural and dynamical properties of the solvated structure have been studied in details. The partial atom-atom radial distribution functions for the hydrated hydroxyl do not show drastic differences with the radial distribution functions for liquid water. The OH is found to be a more active hydrogen bond donor and acceptor than the water molecule, but the accepted hydrogen bonds are much weaker than for the hydroxide OH- ion. The first solvation shell of the OH is less structured than the water's one and contains a considerable fraction of water molecules that are not hydrogen bonded to the hydroxyl. Part of them are found to come closer to the solvated radical than the hydrogen bonded molecules do. The lifetime of the hydrogen bonds accepted by the hydroxyl is found to be shorter than the hydrogen bond lifetime in water. A hydrogen transfer between a water molecule and the OH radical has been observed, though it is a much rarer event than a proton transfer between water and an OH- ion. The velocity autocorrelation power spectrum of the hydroxyl hydrogen shows the properties both of the OH radical in clusters and of the OH- ion in liquid.     

The coordination of uranyl in water: a combined quantum chemical and molecular simulation study

The coordination environment of uranyl in water has been studied using a combined quantum mechanical and molecular dynamics approach. Multiconfigurational wave function calculations have been performed to generate pair potentials between uranyl and water. The quantum chemically determined energies have been used to fit parameters in a polarizable force field with an added charge transfer term. Molecular dynamics simulations have been performed for the uranyl ion and up to 400 water molecules. The results show a uranyl ion with five water molecules coordinated in the equatorial plane. The U-O(H(2)O) distance is 2.40 A, which is close to the experimental estimates. A second coordination shell starts at about 4.7 A from the uranium atom. No hydrogen bonding is found between the uranyl oxygens and water. Exchange of waters between the first and second solvation shell is found to occur through a path intermediate between association and interchange. This is the first fully ab initio determination of the solvation of the uranyl ion in water.     

Quantum mechanical/molecular mechanical simulations of the Tl(III) ion in water

Classical molecular dynamics (MD) and combined quantum mechanical/molecular mechanical (QM/MM) MD simulations have been performed to investigate the structural and dynamical properties of the Tl(III) ion in water. A six-coordinate hydration structure with a maximum probability of the Tl-O distance at 2.21 A was observed, which is in good agreement with X-ray data. The librational and vibrational spectra of water molecules in the first hydration shell are blue-shifted compared with those of pure liquid water, and the Tl-O stretching force constant was evaluated as 148 Nm(-1). Both structural and dynamical properties show a distortion of the first solvation shell structure. The second shell ligands' mean residence time was determined as 12.8 ps. The Tl(III) ion can be classified as "structure forming" ion; the calculated hydration energy of -986 +/- 9 kcal mol agrees well with the experimental value of -986 kcal mol.     

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.     

Ions in water: the microscopic structure of concentrated NaOH solutions

A neutron diffraction experiment with isotopic H/D substitution on four concentrated NaOH/H(2)O solutions is presented. The full set of partial structure factors is extracted, by combining the diffraction data with a Monte Carlo simulation. These allow to investigate both the changes of the water structure in the presence of ions and their solvation shells. It is found that the interaction with the solute affects the tetrahedral network of hydrogen bonded water molecules in a manner similar to the application of high pressure to pure water. The solvation shell of the OH(-) ions has an almost concentration independent structure, although with concentration dependent coordination numbers. The hydrogen site coordinates a water molecule through a weak bond, while the oxygen site forms strong hydrogen bonds with a number of molecules that is on the average very close to four at the higher water concentrations and decreases to about three at the lowest one. The competition between hydrogen bond interaction and Coulomb forces in determining the orientation of water molecules within the cation solvation shell is visible in the behavior of the g(NaHw)(r) function     

A density functional theory based study of the microscopic structure and dynamics of aqueous HCl solutions

The aqueous solvation of hydrochloric acid is studied using density functional theory based molecular dynamics simulations at two concentrations. The large simulation boxes that we use allow us to investigate larger-scale structures such as the water-bridged chloride ion network. We find a strong concentration dependence for almost all structural and dynamical properties. Excess protons are mostly present both as Eigen and Zundel structures, either as a direct hydronium-chloride contact-ion pair or a solvent-separated ion pair. Increasing the concentration has a detrimental effect on the natural hydrogen bonded network of water molecules. This effect is visible in our studies as a decrease in the persistence time of the solvation shells around the chloride ions. Also the number of proton hops, determined by a new and well defined identification procedure, suffers from the breakdown of the natural hydrogen bond network.     

Effects of ion concentration on the hydrogen bonded structure of water in the vicinity of ions in aqueous NaCl solutions

Molecular dynamics simulations of dilute and concentrated aqueous NaCl solutions are carried out to investigate the changes of the hydrogen bonded structures in the vicinity of ions for different ion concentrations. An analysis of the hydrogen bond population in the first and second solvation shells of the ions and in the bulk water is done. Although essentially no effect of ions on the hydrogen bonding is observed beyond the first solvation shell of the ions for the dilute solutions, for the concentrated solutions a noticeable change in the average number of water-water hydrogen bonds is observed in the second solvation shells of the ions and even beyond. However, the changes in the average number of hydrogen bonds are found to be relatively less when both water-water and ion-water hydrogen bonds are counted. Thus, the changes in the total number of hydrogen bonds per water are not very dramatic beyond the first solvation shell even for concentrated solutions.     

Microsolvation of the sodium and iodide ions and their ion pair in acetonitrile clusters: a theoretical study

The structural and thermodynamic properties of Na+(CH3CN)n, I-(CH3CN)n, and NaI(CH3CN)n clusters have been investigated by means of room-temperature Monte Carlo simulations with model potentials developed to reproduce the properties of small clusters predicted by quantum chemistry. Ions are found to adopt an interior solvation shell structure, with a first solvation shell containing approximately 6 and approximately 8 acetonitrile molecules for large Na+(CH3CN)n and I-(CH3CN)n clusters, respectively. Structural features of Na+(CH3CN)n are found to be similar to those of Na+(H2O)n clusters, but those of I-(CH3CN)n contrast with those of I-(H2O)n, for which "surface" solvation structures were observed. The potential of mean force calculations demonstrates that the NaI ion pair is thermodynamically stable with respect to ground-state ionic dissociation in acetonitrile clusters. The properties of NaI(CH3CN)n clusters exhibit some similarities with NaI(H2O)n clusters, with the existence of contact ion pair and solvent-separated ion pair structures, but, in contrast to water clusters, both types of ion pairs adopt a well-defined interior ionic solvation shell structure in acetonitrile clusters. Whereas contact ion pair species are thermodynamically favored in small clusters, solvent-separated ion pairs tend to become thermodynamically more stable above a cluster size of approximately 26. Hence, ground-state charge separation appears to occur at larger cluster sizes for acetonitrile clusters than for water clusters. We propose that the lack of a large Na+(CH3CN)n product signal in NaI(CH3CN)n multiphoton ionization experiments could arise from extensive stabilization of the ground ionic state by the solvent and possible inhibition of the photoexcitation mechanism, which may be less pronounced for NaI(H2O)n clusters because of surface solvation structures. Alternatively, increased solvent evaporation resulting from larger excess energies upon photoexcitation or major solvent reorganization on the ionized state could account for the observed solvent-selectivity in NaI cluster multiphoton ionization.     

The elimination of a hydrogen atom in Na(H2O)n

By a systematic examination on Na(H2O)n, with n = 4-7, 9, 10, and 15, we demonstrate that a hydrogen loss reaction can be initiated by a single sodium atom with water molecules. This reaction is similar to the well-known size-dependent intracluster hydrogen loss in Mg+(H2O)n, which is isoelectronic to Na(H2O)n. However, with one less charge on Na(H2O)n than that on Mg+(H2O)n, the hydrogen loss for Na(H2O)n is characterized by a higher barrier and a more flexible solvation shell around the metal ion, although the reaction should be accessible, as the lowest barrier is around 8 kcal/mol. Interestingly, the hydroxide ion OH- produced in the process is stabilized by the solvation of H2O molecules and the formation of an ion pair Na+(H2O)4(H2O)n-l-4[OH-(H2O)l]. The activation barrier is reduced as the unpaired electron in Na(H2O)n moves to higher solvation shells with increasing cluster size, and the reaction is not switched off for larger clusters. This is in sharp contrast to the reaction for Mg+(H2O)n, in which the OH- ion is stabilized by direct coordination with Mg2+ and the reaction is switched off for n > 17, as the unpaired electron moved to higher solvation shells. Such a contrast illustrates the important link between microsolvation environment and chemical reactivity in solvation clusters.     

Strength of hydrogen bonds of water depends on local environment

In-depth knowledge of water-water potential is important for devising and evaluating simple water models if they are to accurately describe water properties and reflect various solvation phenomena. Water-water potential depends upon inter-molecular distance, relative orientation of water molecules, and also local environment. When placed at a favorable distance in a favorable orientation, water molecules exhibit a particularly strong attractive interaction called hydrogen bond. Although hydrogen bond is very important for its effects on the elements of life, industrial applications, and bulk water properties, there is no scientific consensus on its true nature and origin. Using quantum-mechanical methods, hydrogen bond strength was calculated in different local environments. A simple empirical linear relationship was discovered between maximum hydrogen bond strength and the number of water molecules in the local environment. The local environment effect was shown to be considerable even on the second coordination shell. Additionally, a negative linear correlation was found between maximum hydrogen bond strength and the distance, at which it was observed. These results provide novel insights into the nature of hydrogen bonding.     

Spectral signatures and molecular origin of acid dissociation intermediates

The existence of a broad, mid-infrared absorption ranging from 1000 to 3000 cm(-1) is usually interpreted as a signature for the existence of protonated water networks. Herein, we use cryogenic mixtures of water and hydrogen fluoride (HF) and show experimental and computational evidence that similarly wide absorptions can be generated by a broad distribution of proton-shared and ion pair complexes. In the present case, we demonstrate that the broadening is mainly inhomogeneous, reflecting the fact that the topology of the first solvation shell determines the local degree of ionization and the shared-proton asymmetric stretching frequency within H2O x HF complexes. The extreme sensitivity of the proton transfer potential energy hypersurface to local hydrogen bonding topologies modulates its vibrational frequency from 2800 down to approximately 1300 cm(-1), the latter value being characteristic of solvation geometries that yield similar condensed-phase proton affinities for H2O and fluoride. By linking the local degree of ionization to the solvation pattern, we are able to propose a mechanism of ionization for HF in aqueous solutions and to explain some of their unusual properties at large concentrations. However, an important conclusion of broad scientific interest is our prediction that spectral signatures that are normally attributed to protonated water networks could also reveal the presence of strong hydrogen bonds between un-ionized acids and water molecules, with important consequences to spectroscopic investigations of biologically relevant proton channels and pumps.     

Influence of the ring shape of different cyclodextrins on complexation properties

The complexation properties of various cyclodextrins depend on their ring size and the molecular shape of the guest molecules, but also on the hydrophobicity of the interior. Convenient systems for the characterization of the environment (solvation shell or complexing host molecules) are alkylaminomethyl-naphthols, molecules with an intramolcular hydrogen bond allowing proton transfer, which is extremely sensitive to solvation phenomena.     

Theoretical study of the solvation of HgCl2, HgClOH, Hg(OH)2 and HgCl3(-): a density functional theory cluster approach

The determination of the solvation shell of Hg(II)-containing molecules and especially the interaction between Hg(II) and water molecules is the first requirement to understand the transmembrane passage of Hg into the cell. We report a systematic DFT study by stepwise solvation of HgCl(2) including up to 24 water molecules. In order to include pH and salinity effects, the solvation patterns of HgClOH, Hg(OH)(2) and HgCl(3)(-) were also studied using 24 water molecules. In all cases the hydrogen bond network is crucial to allow orbital-driven interactions between Hg(II) and the water molecules. DFT Born-Oppenheimer molecular dynamics simulations starting from the stable HgCl(2)-(H(2)O)(24) structure revealed that an HgCl(2)-(H(2)O)(3) trigonal bipyramid effective solute appears and then the remaining 21 water molecules build a complete first solvation shell, in the form of a water-clathrate. In the HgCl(2), HgClOH, Hg(OH)(2)-(H(2)O)(24) optimized structures Hg also directly interacts with 3 water molecules from an orbital point of view (three Hg-O donor-acceptor type bonds). All the other interactions are through hydrogen bonding. The cluster-derived solvation energies of HgCl(2), HgClOH and Hg(OH)(2) are estimated to be -34.4, -40.1 and -47.2 kcal mol(-1), respectively.