Single particle and pair dynamics in water-formic acid mixtures containing ionic and neutral solutes: nonideality in dynamical properties

A series of molecular dynamics simulations of water-formic acid mixtures containing either an ionic solute or a neutral hydrophobic solute has been performed to study the extent of nonideality in the dynamics of these solutes for varying composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to show nonideal behavior with variation of composition, and similar nonideality is also observed for the diffusion and orientational relaxation of solvent molecules in these mixtures. The diffusion coefficient of a neutral hydrophobic solute, however, decreases monotonically with increase in water concentration. We have also investigated some of the pair dynamical properties such as water-water and water-formic acid hydrogen bond relaxation and residence dynamics of water molecules in water and formic acid hydration shells. The lifetimes of water-water hydrogen bonds are found to be longer than those between formic acid carbonyl oxygen-water hydrogen bonds, whereas the lifetimes of formic acid hydroxyl hydrogen-water hydrogen bonds are longer than those of water-water hydrogen bonds. In general, the hydrogen bond lifetimes for both water-water and water-formic acid hydrogen bonds are found to decrease with increase in water concentration. Residence times of water molecules also show the same trend with increase in formic acid concentration. Interestingly, these pair dynamical properties show a monotonic dependence on composition without any maximum or minimum and behave almost ideally with respect to changes in the composition of the mixtures. The present calculations are performed with fixed-charge nonpolarizable models of the solvent and solute molecules without taking into account many-body polarization effects in an explicit manner.     

Hydrogen bond and residence dynamics of ion-water and water-water pairs in supercritical aqueous ionic solutions: dependence on ion size and density

We have carried out a series of molecular dynamics simulations to investigate the hydrogen bond and residence dynamics of X(-)-water (X=F, Cl, and I) and pairs in aqueous solutions at a temperature of 673 K. The calculations are done at six different water densities ranging from 1.0 to 0.15 g cm(-3). The hydrogen bonds are defined by using a set of configurational criteria with respect to the anion(oxygen)-oxygen and anion(oxygen)-hydrogen distances and the anion(oxygen)-oxygen-hydrogen angle for an anion(water)-water pair. The F(-)-water hydrogen bonds are found to have a longer lifetime than all other hydrogen bonds considered in the present study. The lifetime of Cl(-)-water hydrogen bonds is shorter than that of F(-)-water hydrogen bonds but longer than the lifetime of water-water hydrogen bonds. The lifetimes of I(-)-water and water-water hydrogen bonds are found to be very similar. Generally, the lifetimes of both anion-water and water-water hydrogen bonds are found to be significantly shorter than those found under ambient conditions. In addition to hydrogen bond lifetimes, we have also calculated the residence times and the orientational relaxation times of water molecules in ion(water) hydration shells and have discussed the correlations of these dynamical quantities with the observed dynamics of anion(water)-water hydrogen bonds as functions of the ion size and density of the supercritical solutions.     

Binding strength of sodium ions in cellulose for different water contents

The interaction strength of sodium ions (Na(+)) with cellulose is investigated from first principles for varying degrees of water content. We find that the interaction of water molecules and Na(+) can be studied independently at the various OH groups in cellulose which we categorize as two different types. In the absence of water, Na(+) forms strong ionic bonds with the OH groups of cellulose. When water molecules are anchored to the OH groups via hydrogen bonds, Na(+) can eventually no longer bind to the OH groups, but will instead interact with the oxygen atoms of the water molecules. Due to the rather weak attachment of the latter to the OH groups, Na(+) becomes effectively more mobile in the fully hydrated cellulose framework. The present study thus represents a significant step toward a first-principles understanding of the experimentally observed dependence of ionic conductivity on the level of hydration in cellulose network.     

Nonideality in diffusion of ionic and hydrophobic solutes and pair dynamics in water-acetone mixtures of varying composition

We have performed a series of molecular dynamics simulations of water-acetone mixtures containing either an ionic solute or a neutral hydrophobic solute to study the extent of nonideality in the dynamics of these solutes with variation of composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to change nonmonotonically with the composition of the mixtures showing strong nonideality of their dynamics. Also, the extent of nonideality in the diffusion of these charged solutes is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and acetone molecules in these mixtures which show a somewhat similar changes in the solvation characteristics of charged and dipolar solutes with changes of composition of water-acetone mixtures. The diffusion of the hydrophobic solute, however, shows a monotonic increase with increase of acetone concentration showing its different solvation characteristics as compared to the charged and dipolar solutes. The links between the nonideality in diffusion and solvation structures are further confirmed through calculations of the relevant solute-solvent and solvent-solvent radial distribution functions for both ionic and hydrophobic solutes. We have also calculated various pair dynamical properties such as the relaxation of water-water and acetone-water hydrogen bonds and residence dynamics of water molecules in water and acetone hydration shells. The lifetimes of both water-water and acetone-water hydrogen bonds and also the residence times of water molecules are found to increase steadily with increase in acetone concentration. No maximum or minimum was found in the composition dependence of these pair dynamical quantities. The lifetimes of water-water hydrogen bonds are always found to be longer than that of acetone-water hydrogen bonds in these mixtures. The residence times of water molecules are also found to follow a similar trend.     

Adsorption of water molecules on selected charged sodium-chloride clusters

The adsorption of water molecules (H(2)O) on sodium chloride cluster cations and anions was studied at 298 K over a mass range of 100-1200 amu using a custom-built laser desorption ionization reactor and mass spectrometer. Under the conditions used, the cations Na(3)Cl(2)(+) and Na(4)Cl(3)(+) bind up to three water molecules, whereas the larger cations, Na(5)Cl(4)(+) to Na(19)Cl(18)(+), formed hydrates with one or two only. The overall trend is a decrease in hydration with increasing cluster size, with an abrupt drop occurring at the closed-shell Na(14)Cl(13)(+). As compared to the cluster cations, the cluster anions showed almost no adsorption. Among smaller clusters, a weak adsorption of one water molecule was observed for the cluster anions Na(6)Cl(7)(-) and Na(7)Cl(8)(-). In the higher mass region, a substantial adsorption of one water molecule was observed for Na(14)Cl(15)(-). Density functional theory (DFT) computations were carried out for the adsorption of one molecule of H(2)O on the cations Na(n)Cl(n-1)(+), for n = 2-8, and the anions Na(n)Cl(n+1)(-), for n = 1-7. For each ion, the structure of the hydrate, the hydration energy, and the standard-state enthalpy, entropy, and Gibbs energy of hydration at 298 K were computed. In addition, it was useful to compute the distortion energy, defined as the electronic energy lost due to weakening of the Na-Cl bonds upon adsorption of H(2)O. The results show that strong adsorption of a H(2)O molecule occurs for the linear cations only at an end Na ion and for the nonlinear cations only at a corner Na ion bonded to two Cl ions. An unexpected result of the theoretical investigation for the anions is that certain low-energy isomers of Na(6)Cl(7)(-) and Na(7)Cl(8)(-) bind H(2)O strongly enough to produce the observed weak adsorption. The possible implications of these results for the initial hydration of extended NaCl surfaces are discussed.     

Hydrogen bond lifetime dynamics at the interface of a surfactant monolayer

The dynamics of water near the polar headgroups of surfactants in a monolayer adsorbed at the air/water interface is likely to play a decisive role in determining the physical behavior of such organized assemblies. We have carried out an atomistic molecular dynamics (MD) simulation of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature with a surface coverage of that at the critical micelle concentration (78 Angstrom(2)/molecule). Detailed analyses of the lifetime dynamics of surfactant-water (SW) and water-water (WW) hydrogen bonds at the interface have been carried out. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the bound states at the interface and quantification of the dynamic equilibrium between bound and quasi-free water molecules, in terms of time-dependent relaxation rates. It is observed that the water molecules present in the first hydration layer form strong hydrogen bonds with the surfactant headgroups and hence have longer lifetimes. Importantly, it is found that the overall relaxation of the SW hydrogen bonds is faster for those water molecules which form two hydrogen bonds with the surfactant headgroups than those forming one such hydrogen bond. Equally interestingly, it is further noticed that water molecules beyond the first hydration layer form weaker hydrogen bonds than pure bulk water.     

Dynamics of halide ion-water hydrogen bonds in aqueous solutions: dependence on ion size and temperature

We have carried out a series of molecular dynamics simulations to investigate the dynamics of X(-)-water (X = F, Cl, Br, and I) and water-water hydrogen bonds in aqueous alkali halide solutions at room temperature and also of Cl(-)-water and water-water hydrogen bonds at seven different temperatures ranging from 238 to 318 K. The hydrogen bonds are defined by using a set of configurational criteria with respect to the anion(oxygen)-oxygen and anion(oxygen)-hydrogen distances and the anion(oxygen)-oxygen-hydrogen angle for an anion(water)-water pair. The results of the hydrogen bond dynamics are obtained for two different cutoff values for the angular criterion. In both cases, similar dynamical behavior of the hydrogen bonds is found with respect to their dependence on ion size and temperature. The fluoride ion-water hydrogen bonds are found to break at a much slower rate than water-water hydrogen bonds, while the lifetimes of chloride and bromide ion-water hydrogen bonds are found to be shorter than those of fluoride ion-water ones but still longer than water-water hydrogen bonds. The short-time dynamics of iodide ion-water hydrogen bonds is found to be slightly faster, while its long-time dynamics is found to be slightly slower than the corresponding water-water hydrogen bond dynamics. Correlations of the observed dynamics of anion(water)-water hydrogen bonds with those of rotational and translational diffusion and residence times of water molecules in ion(water) hydration shells are also discussed. With variation of temperature, the lifetimes of both Cl(-)-water and water-water hydrogen bonds are found to show Arrhenius behavior with a slightly higher activation energy for the Cl(-)-water hydrogen bonds.     

Ion permeation dynamics in carbon nanotubes

Molecular dynamics simulations are carried out to investigate the permeation of ions and water in a membrane consisting of single wall carbon nanotubes possessing no surface charges connecting two reservoirs. Our simulations reveal that there are changes in the first hydration shell of the ions upon confinement in tubes of 0.82 or 0.90 nm effective internal diameter. Although the first minimum in the g(r) is barely changed in the nanotube compared to in the bulk solution, the hydration number of Na(+) ion is reduced by 1.0 (from 4.5 in bulk to 3.5 in the 0.90 nm tube) and the hydration number is reduced further in the 0.82 nm tube. The changes in the hydration shell of Cl(-) ion are negligible, within statistical errors. The water molecules of the first hydration shell of both ions exchange less frequently inside the tube than in the bulk solution. We compare ion trajectories for ions in the same tube under identical reservoir conditions but with different numbers of ions in the tubes. This permits investigation of changes in structure and dynamics which arise from multiple ion occupancy in a carbon nanotube possessing no surface charges. We also investigated the effects of tube flexibility. Ions enter the tubes so as to form a train of ion pairs. We find that the radial distribution profiles of Na(+) ions broaden significantly systematically with increasing number of ion pairs in the tube. The radial distribution profiles of Cl(-) ions change only slightly with increasing number of ions in the tube. Trajectories reveal that Na(+) ions do not pass each other in 0.90 nm tubes, while Cl(-) ions pass each other, as do ions of opposite charge. An ion entering the tube causes the like-charged ions preceding it in the tube to be displaced along the tube axis and positive or negative ions will exit the tube only when one or two other ions of the same charge are present in the tube. Thus, the permeation mechanism involves multiple ions and Coulomb repulsion among the ions plays an essential role.     

Atomistic nature of NaCl nucleation at the solid-liquid interface

The early stage of heterogeneous nucleation of NaCl from supersaturated NaCl aqueous solution at the water-NaCl (001) interface has been investigated by molecular dynamics simulations. The critical size of the nuclei for spontaneous growth was found to be as small as two atoms (a Na(+)-Cl(-) ion pair) at high supersaturation. Due to the presence of a relatively stable water network and the effect of the hydration force at the interface, the stable nuclei formed on the NaCl (001) are found to contain more Na(+) ions than Cl(-) ions. The different deposition characteristics of the Na(+) and Cl(-) solutes lead to a positively charged substrate and thus may introduce another driving force for nucleation besides the level of solution supersaturation. The role of water was further confirmed by comparison with NaCl epitaxy growth in the vacuum.     

Ion condensation behavior and dynamics of water molecules surrounding the sodium poly(methacrylic acid) chain in water: a molecular dynamics study

All-atom molecular dynamics simulations are used to study the condensation behavior of monovalent (Na(+)) and multivalent (Ca(2+)) salt counterions associated with the co-ions (Cl(-)) surrounding the charged poly(methacrylic acid) (PMAA) chain in water. The study is extended to the influences on chain conformation, local arrangement, and dynamics of water in the highly diluted aqueous solutions. We find that even when the salt ions are monovalent, they attract more than one charged monomer and act as a bridging agent within the chain, as the multivalent salt ions. In principle, the salt ions bridge between not only the "non-adjacent" but also the "adjacent" charged monomers, leading to a more coil-like and a locally stretched conformation, respectively. With an increase in the salt concentration, the amount of coiled-type condensed ions increase and reach a maximum when the chain conformation becomes the most collapsed; whereas, the stretched-type shows an opposite trend. Our results show that the attractive interactions through the condensed salt ions between the non-adjacent monomers are responsible for the conformational collapse. When the salt concentration increases high enough, a significant increase for the stretched-type condensed ions makes an expansion effect on the chain. These stretched-type salt ions, followed by the adsorption of the co-ions and water molecules, tend to form a multilayer organization outside surrounding the PMAA chain. Thus, the expansion degree of the chain conformation is greatly limited. When only the monovalent Na(+) ions are present in the solutions, water molecules are primarily adsorbed into either the condensed Na(+) ions or the COO(-) groups. These adsorbed water molecules form hydrogen bonds with each other and enhance the local bridging behavior associated with the Na(+) condensation on the resultant chain conformation. With an increase in the amount of multivalent Ca(2+) salt ions, more water molecules are bonded directly with the condensed Ca(2+) ions. In this case, only the condensed Ca(2+) ions provide a strong bridging effect within the polymer chain. We observe a significant shift towards a higher frequency of the oxygen vibration spectrum and only a slight shift towards a higher frequency of the hydrogen spectrum for the water molecules associated with the ion condensation.     

Properties of a water layer on hydrophilic and hydrophobic self-assembled monolayer surfaces: A molecular dynamics study

The microscopic behaviors of a water layer on different hydrophilic and hydrophobic surfaces of well ordered self-assembled monolayers (SAMs) are studied by molecular dynamics simulations. The SAMs consist of 18-carbon alkyl chains bound to a silicon(111) substrate, and the characteristic of its surface is tuned from hydrophobic to hydrophilic by using different terminal functional groups ( CH 3 , COOH). In the simulation, the properties of water membranes adjacent to the surfaces of SAMs were reported by comparing pure water in mobility, structure, and orientational ordering of water molecules. The results suggest that the mobility of water molecules adjacent to hydrophilic surface becomes weaker and the molecules have a better ordering. The distribution of hydrogen bonds indicates that the number of water-water hydrogen bonds per water molecule tends to be lower. However, the mobility of water molecules and distribution of hydrogen bonds of a water membrane in hydropho- bic system are nearly the same as those in pure water system. In addition, hydrogen bonds are mainly formed between the hydroxyl of the COOH group and water molecules in a hydrophilic system, which is helpful in understanding the structure of interfacial water.     

Ab initio molecular dynamics study of glycine intramolecular proton transfer in water

We use ab initio molecular-dynamics simulations to quantify structural and thermodynamic properties of a model proton transfer reaction that converts a neutral glycine molecule, stable in the gas phase, to the zwitterion that predominates in aqueous solution. We compute the potential of mean force associated with the direct intramolecular proton transfer event in glycine. Structural analyses show that the average hydration number (N(w)) of glycine is not constant along the reaction coordinate, but rather progresses from N(w) = 5 in the neutral molecule to N(w) = 8 for the zwitterion. We report the free-energy difference between the neutral and charged glycine molecules, and the free-energy barrier to proton transfer. Finally, we identify the approximations inherent in our method and estimate the corresponding corrections to our reported thermodynamic predictions.     

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.     

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.     

Ions in water: the microscopic structure of concentrated hydroxide solutions

Neutron-diffraction data on aqueous solutions of hydroxides, at solute concentrations ranging from 1 solute per 12 water molecules to 1 solute per 3 water molecules, are analyzed by means of a Monte Carlo simulation (empirical potential structure refinement), in order to determine the hydration shell of the OH- in the presence of the smaller alkali metal ions. It is demonstrated that the symmetry argument between H+ and OH- cannot be used, at least in the liquid phase at such high concentrations, for determining the hydroxide hydration shell. Water molecules in the hydration shell of K+ orient their dipole moment at about 45 degrees from the K+-water oxygen director, instead of radially as in the case of the Li+ and Na+ hydration shells. The K+-water oxygen radial distribution function shows a shallower first minimum compared to the other cation-water oxygen functions. The influence of the solutes on the water-water radial distribution functions is shown to have an effect on the water structure equivalent to an increase in the pressure of the water, depending on both ion concentration and ionic radius. The changes of the water structure in the presence of charged solutes and the differences among the hydration shells of the different cations are used to present a qualitative explanation of the observed cation mobility.     

Monolayer of aerosol-OT surfactants adsorbed at the air/water interface: an atomistic computer simulation study

An atomistic molecular dynamics (MD) simulation has been carried out to investigate the structural and dynamical properties of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature and at a surface coverage corresponding to that at its critical micelle concentration (78 A(2)/molecule). The estimated thickness of the adsorbed layer is in good agreement with neutron reflection data. The study shows that the surfactants exhibit diffusive motion in the plane of the interface. It is observed that the surfactant monolayer has a strong influence in restricting both the translational and reorientational motions of the water molecules close to the interface. A drastic difference in the dipolar reorientational motion of water molecules in the aqueous layer is observed with a small variation of the distance from the surfactant headgroups. It has been observed that the water molecules in the first hydration layer (region 1) form strong hydrogen bonds with surfactant headgoups. This results in the slower structural relaxation of water-water hydrogen bonds in the first hydration layer compared to that in the pure bulk water. Most interestingly, we notice that the water molecules present in the layer immediately after the first hydration layer form weaker hydrogen bonds and thus relax faster than even pure bulk water.     

Sensitivity of hydrogen bond lifetime dynamics to the presence of ethanol at the interface of a phospholipid bilayer

Atomistic molecular dynamics simulations of a fully hydrated liquid crystalline lamellar phase of a dimyrystoylphosphatidylcholine lipid bilayer containing ethanol at 1:1 composition as well as of the pure lamellar phase of the bilayer have been performed. Detailed analyses have been carried out to investigate the effects of ethanol, if any, on the lifetime dynamics of lipid-water and water-water hydrogen bonds in the hydration layer of the lipid headgroups. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed the identification of the bound states at the bilayer interface and the quantification of the dynamic equilibrium between the bound and the free water molecules, in terms of time-dependent relaxation rates. The calculations show that the overall relaxation of phosphate-water hydrogen bonds is faster in the presence of ethanol. Studies of the residence time and the number fluctuation of the hydration layer water molecules reveal that the presence of ethanol molecules decreases the rigidity of the lipid hydration layer.     

Study of intermolecular interactions in aqueous solutions by millimeter spectroscopy

The absorption of millimeter electromagnetic radiation (v=1.4, 1.71, and 5 cm⁻¹) by aqueous solutions of glycine (pH 6.1–6.2) in the concentration range of 0.5–2.5 mol L⁻¹ was measured. It was found that the absorbing ability of the water present in the solutions, is higher than that of pure water. This phenomenon is explained by the presence of a center of negative hydration in the structure of the glycine zwitterion, which results in an increase in the rotational mobility of water molecules immobilized in the hydrate shell of the glycine zwitterion. For Part 5, see Ref. 1. Deceased. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 7, pp. 1305–1307, July, 1997.     

Density functional study of ion hydration for the alkali metal ions (Li+, Na+, K+) and the halide ions (F-, Br-, Cl-)

We performed first principles density functional calculations to study the effect of monovalent ions M+ (M = Li,Na,K) and A- (A = F,Cl,Br) in water with the aim of characterizing the local molecular properties of hydration. For this reason, several ion-water clusters, up to five or six water molecules were considered; such structures were optimized, and the Wannier analysis was then applied to determine the average molecular dipole moment of water. We found that with an increasing number of water molecules, the molecular polarization is determined by the water-water interaction rather than the water-ion interaction, as one would intuitively expect. These results are consistent with those obtained in previous density functional calculations and with other results obtained by employing classical polarizable water models. The main message of this work is that as one increases the number of molecules the average dipole moment of all water molecules and the ones in the first shell tends to the same value as the average of a similar sized cluster of pure water. This supports the use of nonpolarizable classical models of water in classical atomistic simulations.     

Hydration of alkali ions from first principles molecular dynamics revisited

Structural and dynamical properties of the hydration of Li(+), Na(+), and K(+) in liquid water at ambient conditions were studied by first principles molecular dynamics. Our simulations successfully captured the different hydration behavior shown by the three alkali ions as observed in experiments. The present analyses of the dependence of the self-diffusion coefficient and rotational correlation time of water on the ion concentration suggest that Li(+) (K(+)) is certainly categorized as a structure maker (breaker), whereas Na(+) acts as a weak structure breaker. An analysis of the relevant electronic structures, based on maximally localized Wannier functions, revealed that the dipole moment of H(2)O molecules in the first solvation shell of Na(+) and K(+) decreases by about 0.1 D compared to that in the bulk, due to a contraction of the oxygen lone pair orbital pointing toward the metal ion.