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.     

Formation of giant solvation shells around fulleride anions in liquid ammonia

Here, we measure the solvation structure of fulleride C605- anions in potassium ammonia solution using neutron diffraction. We find a very strong solvation structure consisting of two shells of ammonia densely packed around the anion. The system's structure is driven by the propensity of ammonia molecules to direct one of their hydrogen bonds to the center of the anion while retaining axial hydrogen bonding within the shells. This permits high concentrations of solvent separated fulleride anions.     

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

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

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     

Electron tunneling in lithium-ammonia solutions probed by frequency-dependent electron spin relaxation studies

Electron transfer or quantum tunneling dynamics for excess or solvated electrons in dilute lithium-ammonia solutions have been studied by pulse electron paramagnetic resonance (EPR) spectroscopy at both X- (9.7 GHz) and W-band (94 GHz) frequencies. The electron spin-lattice (T(1)) and spin-spin (T(2)) relaxation data indicate an extremely fast transfer or quantum tunneling rate of the solvated electron in these solutions which serves to modulate the hyperfine (Fermi-contact) interaction with nitrogen nuclei in the solvation shells of ammonia molecules surrounding the localized, solvated electron. The donor and acceptor states of the solvated electron in these solutions are the initial and final electron solvation sites found before, and after, the transfer or tunneling process. To interpret and model our electron spin relaxation data from the two observation EPR frequencies requires a consideration of a multiexponential correlation function. The electron transfer or tunneling process that we monitor through the correlation time of the nitrogen Fermi-contact interaction has a time scale of (1-10) × 10(-12) s over a temperature range 230-290 K in our most dilute solution of lithium in ammonia. Two types of electron-solvent interaction mechanisms are proposed to account for our experimental findings. The dominant electron spin relaxation mechanism results from an electron tunneling process characterized by a variable donor-acceptor distance or range (consistent with such a rapidly fluctuating liquid structure) in which the solvent shell that ultimately accepts the transferring electron is formed from random, thermal fluctuations of the liquid structure in, and around, a natural hole or Bjerrum-like defect vacancy in the liquid. Following transfer and capture of the tunneling electron, further solvent-cage relaxation with a time scale of ～10(-13) s results in a minor contribution to the electron spin relaxation times. This investigation illustrates the great potential of multifrequency EPR measurements to interrogate the microscopic nature and dynamics of ultrafast electron transfer or quantum-tunneling processes in liquids. Our results also impact on the universal issue of the role of a host solvent (or host matrix, e.g. a semiconductor) in mediating long-range electron transfer processes and we discuss the implications of our results with a range of other materials and systems exhibiting the phenomenon of electron transfer.     

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.     

Molecular orbital theory of the hydrogen bond. 26. The hydration of uracil

Ab initio SCF calculations with the STO -3G basis set have been performed to determine the structure and stability of a 6:1 water:uracil heptamer in which water molecules are hydrogen bonded to uracil at each of the six hydrogen-bonding sites in the uracil molecular plane. The structure of the heptamer describes a stable arrangement of these six water molecules, which are the primary solvent molecules in the first solvation shell, and is suggestive of the arrangement of secondary solvent molecules in that shell in the nonpolar region of the uracil molecular plane. The stabilization energy of the heptamer is 49.6 kcal/mol, or 8.3 kcal/mol per water molecule. The hydrogen bonds between uracil and water are the primary factor in the stabilization of the complex, although water–water interactions and nonadditivity effects are also significant.     

The structure of calcium-ammonia solutions by neutron diffraction

The microscopic structures of calcium-ammonia solutions have been established by using neutron diffraction. Total structure factors measured at 230 K reveal immediately the evolution of an uncommonly intense diffraction prepeak in the metallic solutions. As concentration is increased from 4 mole percent metal to 10 mole percent metal (i.e., saturation), this feature intensifies and shifts from 0.6 to 0.9 A(-1). It is therefore evidence of well developed intermediate-range ordering among the solvated cations, and is a microstructural signature of the observed strong phase separation of metallic (concentrated) and nonmetallic (dilute) solutions. The technique of isotopic labelling of *N by 15N was then used in conjunction with difference analysis to focus on the solvent structure in metallic solutions at 4 and 10 mole percent metal. These nitrogen-centered functions are analyzed in conjunction with classical Monte Carlo computer simulation techniques, to provide us with detailed insight into the calcium solvation and the extent of hydrogen bonding. We find that calcium is solvated by approximately 6-7 ammonia molecules, with a Ca-N distance of around 2.45 A. There is evidence of hydrogen bonding among the solvent molecules, even in the saturated 10 mole percent metal solution.     

Structure formation features of water and concentrated aqueous lithium halide solutions at low temperatures from the data of integral equation method

The structure of water and the influence of halide ions on the structure formation of concentrated LiX : H₂O (1 : 5; X = Cl, Br, I) solutions at low temperatures were studied by the method of integral equations. Based on the results obtained, supercooling of pure water is expected to significantly enhance the tetrahedral ordering of its molecules, strengthen hydrogen bonding in the system, and decrease the number of the nearest-neighbor water molecules. The effects for the solutions on lowering the temperature include a partial restoration of the tetrahedral network of H-bonds of the solvent molecules, insignificant increase in the number of the nearest-neighbor water molecules, enhancement of the coordination ability of Li⁺ cation, strengthening of hydrogen bonding between anions and water molecules in the first hydration shell, increase in the number of solvent-separated ion pairs, and weakening of the temperature effect on these structural parameters in the following order of solutions: LiCl > LiBr > LiI. The probability of contact ion pair formation in the systems studied should appreciably decrease. The temperature should to a greater extent influence the associative ability of larger anions.     

Local order in aqueous NaCl solutions and pure water: X-ray scattering and molecular dynamics simulations study

The microstructures of pure water and aqueous NaCl solutions over a wide range of salt concentrations (0-4 m) under ambient conditions are characterized by X-ray scattering and molecular dynamics (MD) simulations. MD simulations are performed with the rigid SPC water model as a solvent, while the ions are treated as charged Lennard-Jones particles. Simulated data show that the first peaks in the O...O and O...H pair correlation functions clearly decrease in height with increasing salt concentration. Simultaneously, the location of the second O...O peak, the signature of the so-called tetrahedral structure of water, gradually disappears. Consequently, the degree of hydrogen bonding in liquid water decreases when compared to pure fluid. MD results also show that the hydration number around the cation decreases as the salt concentration increases, which is most likely because some water molecules in the first hydration shell are occasionally substituted by chlorine. In addition, the fraction of contact ion pairs increases and that of solvent-separated ion pairs decreases. Experimental data are analyzed to deduce the structure factors and the pair correlation functions of each system. X-ray results clearly show a perturbation of the association structure of the solvent and highlight the appearance of new interactions between ions and water. A model of intermolecular arrangement via MD results is then proposed to describe the local order in each system, as deduced from X-ray scattering data.     

Monte Carlo simulations of the solution structure of simple alcohols in water-acetonitrile mixtures

Monte Carlo simulations have been performed to explore the solution structure of ethyl, isopropyl, isobutyl, and tertiary butyl alcohols in pure water, pure acetonitrile, and different mixtures of the two solvents. The explicit solvent studies in NpT ensembles at T = 298 K illustrate that the solute "discriminates" the solvent's components and that the composition of the first solvation shell differs from that of the bulk solution. Since the polarizable continuum dielectric method (PCM) does not presently model the solvation of molecules with both polar and apolar sites in mixed protic solvents, we suggest a direction for further program development wherein a continuum dielectric method would accept more than one solvent and the solute sites would be solvated by user-defined solvent components. The prevailing solvation model will be determined upon the lowest free energy calculated for a particular solvation pattern of the solute having a specific conformational/tautomeric state. Characterization of equilibrium hydrogen-bond formation becomes a complicated problem that depends on the chemical properties of the solute and its conformation, as well as upon the varying nature of the first solvation shell. For example, while the number of hydrogen bonds to secondary and tertiary alcohol solutes are nearly constant in pure water and in water-acetonitrile mixtures with at least 50% water content, the number of hydrogen bonds to primary alcohols gradually decreases for most of their conformations when acetonitrile content is increased. Nonetheless, the calculations indicate that O-H...O(water) hydrogen bonds are still possible in a small fraction of the arrangements for the solution models with water content of 30% or less. The isopentene solute does not form any observable hydrogen bonds, despite having an electron-rich, double-bond site.     

Excess protons in water-acetone mixtures

Using molecular dynamics experiments, we analyze equilibrium and dynamical characteristics related to the solvation of excess protons in water-acetone mixtures. Our approach is based on the implementation of an extended valence-bond Hamiltonian, which incorporates translocation of the excess charge between neighboring water molecules. Different mixtures have been analyzed, starting from the pure water case down to solutions with a water molar fraction x(w) = 0.25. In all cases, we have verified that the structure of the first solvation shell of the H(3)O(+) moiety remains practically unchanged, compared to the one observed in pure water. This shell is composed by three water molecules acting as hydrogen bond acceptors, with no evidence of hydrogen bond donor-like connectivity. Moreover, the increment in the acetone concentration leads to a gradual stabilization of Eigen-like [H(3)O[middle dot](H(2)O)(3)](+) configurations, in detriment of Zundel-like [H[middle dot](H(2)O)(2)](+) ones. Rates of proton transfer and proton diffusion coefficients have been recorded at various water-acetone relative concentrations. In both cases, we have found a transition region, in the vicinity of x(w) ～ 0.8, where the concentration dependences of the two magnitudes change at a quantitative level. A crude estimate shows that, at this tagged concentration, the volumes "occupied" by the two solvents become comparable. The origins of this transition separating water-rich from acetone-rich realms is rationalized in terms of modifications operated in the nearby, second solvation shell, which in the latter solutions, normally includes at least, one acetone molecule. Our results would suggest that one possible mechanism controlling the proton transfer in acetone-rich solutions is the exchange of one of these tagged acetone molecules, by nearby water ones. This exchange would give rise to Zundel-like structures, exhibiting a symmetric, first solvation shell composed exclusively by water molecules, and would facilitate the transfer between neighboring water molecules along the resonant complex.     

Structural studies on the hydration of L-glutamic acid in solution

A combination of neutron diffraction augmented with isotopic substitution and computer modeling using empirical potential structure refinement has been used to extract detailed structural information for L-glutamic acid dissolved in 2 M NaOH solution. This work shows that the tetrahedral hydrogen bonding network in water is severely disrupted by the addition of glutamic acid and NaOH, with the number of water-water hydrogen bonds being reduced from 1.8 bonds per water molecule in pure water to 1.4 bonds per water molecule in the present solution. In the glutamic acid molecule, each carboxylate oxygen atom forms an average of three hydrogen bonds with the surrounding water solvent with one of these hydrogens being shared between the two oxygen atoms on each carboxylate group, while each amine hydrogen forms a single hydrogen bond with the surrounding water solvent. Additionally, the average conformation of the glutamic acid molecules in these solutions is extracted.     

Ions in water: the microscopic structure of a concentrated HCl solution

A neutron diffraction experiment with isotopic H/D substitution on a concentrated HCl/H2O solution is presented. The full set of partial structure factors is extracted, by combining the diffraction data with a Monte Carlo simulation. This allows us to investigate both the changes of the water structure in the presence of ions and their solvation shell, overcoming the limitations of standard diffraction experiments. It is found that the interaction with the solutes affects the tetrahedral network of hydrogen bonded water molecules, in a manner similar to the application of an external pressure to pure water, although HCl seems less effective than other solutes, such as NaOH, at the same concentration. Consistent with experimental and theoretical data, the number of water molecules in the solution is not sufficient to completely dissociate the acid molecule. As a consequence, both dissociated H+ and Cl- ions and undissociated HCl molecules coexist in the sample, and this mixture is correctly reproduced in the simulation box. In particular, the hydrated H+ ions, forming a H3O+ complex, participate in three strong and short hydrogen bonds, while a well-defined hydration shell is found around the chlorine ion. These results are not consistent with the findings of early diffraction experiments on the same system and could only be obtained by combining high quality experimental data with a proper computer simulation.     

Solvation effects on the intracluster elimination channels in M+(L)n, where M+= Mg+ and Ca+, L = CH3OH, and NH3, and n = 2-6

The methanol and ammonia solvated Ca (+) or Mg (+) clusters are known to go through intracluster H or CH 3 eliminations which are typically switched on just below n = 6. By first principles calculations at the B3LYP/6-311+G** level, we have identified the transition structures, activation barriers, and energy changes in these reactions for clusters with 2-6 solvent molecules. The activation barrier is crucial to explain the previously reported experimental results. While increasing number of solvent molecules stabilizes a transition structure, the increasing presence of solvent molecules in the first solvation shell makes it difficult for the metal ion to assist the bond breaking through its interaction with the departing H atom or CH 3 group. The balance of these two factors determines whether a particular elimination channel could be switched on.     

Infrared spectroscopy of Li(NH3)n clusters for n=4-7

Infrared spectra of Li(NH3)(n) clusters as a function of size are reported for the first time. Spectra have been recorded in the N-H stretching region for n=4-->7 using a mass-selective photodissociation technique. For the n=4 cluster, three distinct IR absorption bands are seen over a relatively narrow region, whereas the larger clusters yield additional features at higher frequencies. Ab initio calculations have been carried out in support of these experiments for the specific cases of n=4 and 5 for various isomers of these clusters. The bands observed in the spectrum for Li(NH3)(4) can all be attributed to N-H stretching vibrations from solvent molecules in the first solvation shell. The appearance of higher frequency N-H stretching bands for n > or =5 is assigned to the presence of ammonia molecules located in a second solvent shell. These data provide strong support for previous suggestions, based on gas phase photoionization measurements, that the first solvation shell for Li(NH3)(n) is complete at n=4. They are also consistent with neutron diffraction studies of concentrated lithium/liquid ammonia solutions, where Li(NH3)(4) is found to be the basic structural motif.     

Heat capacity effects associated with the hydrophobic hydration and interaction of simple solutes: a detailed structural and energetical analysis based on molecular dynamics simulations

We examine the SPCE [H. J. C. Berendsen et al., J. Chem. Phys. 91, 6269 (1987)] and TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys 112, 8910 (2000)] water models using a temperature series of molecular-dynamics simulations in order to study heat-capacity effects associated with the hydrophobic hydration and interaction of xenon particles. The temperature interval between 275 and 375 K along the 0.1-MPa isobar is studied. For all investigated models and state points we calculate the excess chemical potential for xenon employing the Widom particle insertion technique. The solvation enthalpy and excess heat capacity is obtained from the temperature dependence of the chemical potentials and, alternatively, directly by Ewald summation, as well as a reaction field based method. All three methods provide consistent results. In addition, the reaction field technique allows a separation of the solvation enthalpy into solute/solvent and solvent/solvent parts. We find that the solvent/solvent contribution to the excess heat capacity is dominating, being about one order of magnitude larger than the solute/solvent part. This observation is attributed to the enlarged heat capacity of the water molecules in the hydration shell. A detailed spatial analysis of the heat capacity of the water molecules around a pair of xenon particles at different separations reveals that even more enhanced heat capacity of the water located in the bisector plane between two adjacent xenon atoms is responsible for the maximum of the heat capacity found for the desolvation barrier distance, recently reported by Shimizu and Chan [J. Am. Chem. Soc. 123, 2083 (2001)]. The about 60% enlarged heat capacity of water in the concave part of the joint xenon-xenon hydration shell is the result of a counterplay of strengthened hydrogen bonds and an enhanced breaking of hydrogen bonds with increasing temperature. Differences between the two models with respect to the heat capacity in the xenon-xenon contact state are attributed to the different water model bulk heat capacities, and to the different spatial extension of the structure effect introduced by the hydrophobic particles. Similarities between the different states of water in the joint xenon-xenon hydration shell and the properties of stretched water are discussed.     

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

We study the solvation of iodide in water using density functional theory based molecular-dynamics simulations. Detailed analysis of the structural and dynamical properties of the first solvation shell is presented, showing a disruptive influence of the ion on the local water structure. Iodide-water hydrogen bonding is weak, compared to water-water hydrogen bonds. This effective repulsive ion-water interaction leads to the formation of a quite unstructured solvation shell. The dynamics of water molecules surrounding the iodide is relatively fast. The intramolecular structural and electronical properties of water molecules around the ion are not affected.     

Ab initio studies on the mechanism of the size-dependent hydrogen-loss reaction in Mg+(H2O)n

The mechanism of size-dependent intracluster hydrogen loss in the cluster ions Mg(+)(H(2)O)(n), which is switched on around n=6, and off around n=14, was studied by ab initio calculations at the MP2/6-31G* and MP2/6-31G** levels for n=1-6. The reaction proceeds by Mg(+)-assisted breaking of an H-O bond in one of the H(2)O molecules. The reaction barrier is dependent on both the cluster size and the solvation structure. As n increases from 1 to 6, there is a dramatic drop in the reaction barrier, from greater than 70 kcal mol(-1) for n=1 to less than 10 kcal mol(-1) for n=6. In the transition structures, the Mg atom is close to the oxidation state of +2, and H(2)O molecules in the first solvation shell are much more effective in stabilizing the transition structures and lowering the reaction barriers than H(2)O molecules in the other solvation shells. While the reaction barrier for trimer core structures with only three H(2)O molecules in the first shell is greater than 24 kcal mol(-1), even for Mg(+)(H(2)O)(6), it drops considerably for clusters with four-six H(2)O molecules in the first shell. The more highly coordinated complexes have comparable or slightly higher energy than the trimer core structures, and the presence of such high coordination number complexes is the underlying kinetic factor for the switching on of the hydrogen-loss reaction around n=6. For clusters with trimer core structures, the hydrogen loss reaction is much easier when it is preceded by an isomerization step that increases the coordination number around Mg(+). Delocalization of the electron on the singly occupied molecular orbital (SOMO) away from the Mg(+) ion is observed for the hexamer core structure, while at the same time this isomer is the most reactive for the hydrogen-loss reaction, with an energy barrier of only 2.7 kcal mol(-1) at the MP2/6-31G** level.     

Effects of solute structure on local solvation and solvent interactions: results from UV/vis spectroscopy and molecular dynamics simulations

Solvation of heterocyclic amines in CO(2)-expanded methanol (MeOH) has been explored with UV/vis spectroscopy and molecular dynamics (MD) simulations. A synergistic study of experiments and simulations allows exploration of solute and solvent effects on solvation and the molecular interactions that affect absorption. MeOH-nitrogen hydrogen bonds hinder the n-pi* transition; however, CO(2) addition causes a blue shift relative to MeOH because of Lewis acid/base interactions with nitrogen. Effects of solute structure are considered, and very different absorption spectra are obtained as nitrogen positions change. MD simulations provide detailed solvent clustering behavior around the solute molecules and show that the local solvent environment and ultimately the spectra are sensitive to the solute structure. This work demonstrates the importance of atomic-level information in determining the structure-property relationships between solute structure, local salvation, and solvatochromism.