Structural and dynamical properties of the Hg2+ aqua ion: a molecular dynamics study

Molecular dynamics simulations of the Hg2+ ion in aqueous solution have been carried out using an effective two-body potential derived from quantum mechanical calculations. A stable heptacoordinated structure of the Hg2+ first hydration shell has been observed and confirmed by extended X-ray absorption fine structure (EXAFS) experimental data. The structural properties of the Hg2+ hydration shells have been investigated using radial and angular distribution functions, while the dynamical behavior has been discussed in terms of reorientational correlation functions, mean residence times of water molecules in the first and second hydration shells, and self-diffusion coefficients. The effect of water-water interactions on the Hg2+ hydration properties has been evaluated using the SPC/E and TIP5P water models.     

Development of new Cd2+ and Pb2+ Lennard-Jones parameters for liquid simulations

We present new Lennard-Jones parameters for Cd2+ and Pb2+ ion-water interactions and describe a general methodology to obtain these parameters for any ion. Our strategy is based on the adjustment of ion parameters to reproduce simultaneously experimental absolute hydration free energy and structural properties, namely, g(r) and coordination numbers, obtained from X-ray liquid scattering and quantum mechanical/molecular mechanical (QM/MM) calculations. The validation of the obtained parameters is made by the calculation of dynamical properties and comparing them with experimental values and theoretical results from the literature. The transferability of parameters is checked by the calculation of thermodynamic, structural, and dynamical properties cited above with four different water models. The results obtained for Cd2+ and Pb2+ show an overall agreement with reference values. The absolute hydration free energy calculated with the TIP3P, SPC/E, SPC, and TIP4P water models presents, respectively, percent differences of 3.8, 3.0, 4.3, and 7.2% for lead(II) and 9.8, 8.4, 10.2, and 14.1% for cadmium(II) when compared with experimental values. Ion-water mean distance and coordination numbers for the first coordination shell are in good agreement with experimental and QM/MM results for both ions. Cd2+ shows a lesser diffusion coefficient compared to that of Pb2+ despite its smaller atomic radius, indicating a more persistent first coordination shell for the cadmium(II) ion, a result confirmed with calculations of the mean residence time of water molecules in the first coordination shell.     

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.     

Exploring the capabilities of X-ray absorption spectroscopy for determining the structure of electrolyte solutions: computed spectra for Cr(3+) or Rh(3+) in water based on molecular dynamics

Extended X-ray absorption fine structure (EXAFS) spectra of Cr(3+) and Rh(3+) in aqueous solution are analyzed and compared with computed spectra derived from structural results obtained by molecular dynamics (MD) simulation. This procedure quantifies the reliability of the EXAFS structural determination when applied to ions in solution. It provides guidelines for interpreting experimental spectra of octahedrally coordinated metal cations in aqueous solution. A set of relationships among Debye-Waller factors is proposed on the basis of MD results to reduce the number of independent fit parameters. The determination of the second hydration shell is examined. Calculated XANES spectra compare well with experimental ones. Indeed, the splitting observed on the main peak of the Rh K-edge was anticipated by the calculations. Simulated spectra from MD structures of increasing cluster size show a relationship between the second hydration shell and features of the XANES region at energies just above the edge. The combination of quantum and statistical calculations with the XANES spectrum is found to be very fruitful to get insight into the quantitative estimation of structural properties of electrolyte solutions.     

Experimental evidence for a variable first coordination shell of the cadmium(II) ion in aqueous, dimethyl sulfoxide, and N,N'-dimethylpropyleneurea solution

A combined extended X-ray absorption fine structure (EXAFS) and large angle X-ray scattering (LAXS) investigation has been performed to evaluate the coordination structure of the cadmium(II) ion in aqueous, dimethyl sulfoxide, and N,N'-dimethylpropyleneurea (dmpu) solutions. This approach has singled out the existence of a flexible coordination shell around the cadmium(II) ion in aqueous and dimethyl sulfoxide solutions, whereas a regular octahedral complex is detected in dmpu. The EXAFS and LAXS techniques provide different values of the Cd-O first shell distance (2.27(1) A and 2.302(5) A, respectively) for the hydrated and dimethyl sulfoxide solvated complexes, and this discrepancy is originated by the simultaneous presence of hexa- and heptacoordinated complexes in solution, giving rise to a broad distribution of distances around the ion. These findings demonstrate that, in solution, the cadmium(II) ion forms quite flexible hydration and dimethyl sulfoxide solvate complexes undergoing a solvent exchange with unusually stable seven-coordinated intermediate complexes, and therefore the mean ion-solvent distance is longer in solution than in the solid state. In the dmpu solution, due to the bulkiness of the solvent molecules, the octahedral cadmium(II) solvate is extremely crowded and it is not possible for a seventh ligand to enter the inner-coordination shell. This investigation shows that the combined analysis of the EXAFS and LAXS data allows a reliable determination of the structural properties of electrolyte solutions, also in the presence of flexible coordination shell with a variable number of coordinating molecules.     

Car-Parrinello molecular dynamics simulation of Fe 3+ (aq)

The optimized geometry and energetic properties of Fe(D2O)n 3+ clusters, with n = 4 and 6, have been studied with density-functional theory calculations and the BLYP functional, and the hydration of a single Fe 3+ ion in a periodic box with 32 water molecules at room temperature has been studied with Car-Parrinello molecular dynamics and the same functional. We have compared the results from the CPMD simulation with classical MD simulations, using a flexible SPC-based water model and the same number of water molecules, to evaluate the relative strengths and weaknesses of the two MD methods. The classical MD simulations and the CPMD simulations both give Fe-water distances in good agreement with experiment, but for the intramolecular vibrations, the classical MD yields considerably better absolute frequencies and ion-induced frequency shifts. On the other hand, the CPMD method performs considerably better than the classical MD in describing the intramolecular geometry of the water molecule in the first hydration shell and the average first shell...second shell hydrogen-bond distance. Differences between the two methods are also found with respect to the second-shell water orientations. The effect of the small box size (32 vs 512 water molecules) was evaluated by comparing results from classical simulations using different box sizes; non-negligible effects are found for the ion-water distance and the tilt angles of the water molecules in the second hydration shell and for the O-D stretching vibrational frequencies of the water molecules in the first hydration shell.     

Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: comparisons with IR and vibrational circular dichroism spectra

We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed alpha-helix and extended beta and P(II) conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the P(II) conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is alpha(R). Polarizable CHARMM force field gives almost exclusively P(II) conformation and other force fields predict that both alpha-helical and extended (beta and P(II)) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the P(II) conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.     

The Jahn-Teller effect of the Cr(2+) ion in aqueous solution: Ab initio QM/MM molecular dynamics simulations

The hydration structure of Cr(2+) has been studied using molecular dynamics (MD) simulations including three-body corrections and combined ab initio quantum mechanical/molecular mechanical (QM/MM) MD simulations at the Hartree-Fock level. The structural properties are determined in terms of radial distribution functions, coordination numbers, and several angle distributions. The mean residence time was evaluated for describing ligand exchange processes in the second hydration shell. The Jahn-Teller distorted octahedral [Cr(H(2)O)(6)](2+) complex was pronounced in the QM/MM MD simulation. The first-shell distances of Cr(2+) are in the range of 1.9-2.8 A, which are slightly larger than those observed in the cases of Cu(2+) and Ti(3+). No first-shell water exchange occurred during the simulation time of 35 ps. Several water-exchange processes were observed in the second hydration shell with a mean residence time of 7.3 ps.     

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

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

Structure and dynamics of high-spin Ru in aqueous solution: Ab initio QM/MM molecular dynamics simulation

The structural and dynamical properties of high-spin Ru²⁺ in aqueous solution have been theoretically studied using molecular dynamics (MD) simulations. The conventional MD simulation based on pair potentials gives the overestimated average first shell coordination number of 9, whereas the value of 5.9 was observed when the three-body corrected function was included. A combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation has been performed to take into account the many-body effects on the hydration shell structure of Ru²⁺. The most important region, the first hydration shell, was treated by ab initio quantum mechanics at UHF level using the SBKJC VDZ ECP basis set for Ru²⁺ and the 6-31G^∗ basis sets for water. An exact coordination number of 6 for the first hydration shell was obtained from the QM/MM simulation. The QM/MM simulation predicts the average Ru²⁺–O distance of 2.42 Å for the first hydration shell, whereas the values of 2.34 and 2.46 Å are resulted from the pair potentials without and with the three-body corrected simulations, respectively. Several other structural properties representing position and orientation of the solvate molecules were evaluated for describing the hydration shell structure of the Ru²⁺ ion in dilute aqueous solution. A mean residence time of 7.1 ps was obtained for water ligands residing in the second hydration shell.     

Speciation of the curium(III) ion in aqueous solution: a combined study by quantum chemistry and molecular dynamics simulation

The structures of aquo complexes of the curium(III) ion have been systematically studied using quantum chemical and molecular dynamics (MD) methods. The first hydration shell of the Cm3+ ion has been calculated using density functional theory (DFT), with and without inclusion of the conductor-like polarizable continuum medium (CPCM) model of solvation. The calculated results indicate that the primary hydration number of Cm3+ is nine, with a Cm-O bond distance of 2.47-2.48 A. The calculated bond distances and the hydration number are in excellent agreement with available experimental data. The inclusion of a complete second hydration shell of Cm3+ has been investigated using both DFT and MD methods. The presence of the second hydration shell has significant effects on the primary coordination sphere, suggesting that the explicit inclusion of second-shell effects is important for understanding the nature of the first shell. The calculated results indicate that 21 water molecules can be coordinated in the second hydration shell of the Cm3+ ion. MD simulations within the hydrated-ion model suggest that the second-shell water molecules exchange with the bulk solvent with a lifetime of 161 ps.     

Wide-angle X-ray diffraction and molecular dynamics study of medium-range order in ambient and hot water

We have developed wide-angle X-ray diffraction measurements with high energy-resolution and accuracy to study water structure at three different temperatures (7, 25 and 66 °C) under normal pressure. Using a spherically curved Ge crystal an energy resolution better than 15 eV has been achieved which eliminates influence from Compton scattering. The high quality of the data allows for a reliable Fourier transform of the experimental data resolving shell structure out to ~12 ?, i.e. 5 hydration shells. Large-scale molecular dynamics (MD) simulations using the TIP4P/2005 force-field reproduce excellently the experimental shell-structure in the range 4-12 ? although less agreement is seen for the first peak in the intermolecular pair-correlation function (PCF). The Shiratani-Sasai Local Structure Index [J. Chem. Phys. 104, 7671 (1996)] identifies a tetrahedral minority giving the intermediate-range oscillations in the O-O PCF and a disordered majority providing a more featureless background in this range. The current study supports the proposal that the structure of liquid water, even at high temperatures, can be described in terms of a two-state fluctuation model involving local structures related to the high-density and low-density forms of liquid water postulated in the liquid-liquid phase transition hypothesis.     

The solvation structure of Pb(II) in dilute aqueous solution: an ab initio quantum mechanical/molecular mechanical molecular dynamics approach

Structural properties of the hydrated Pb(II) ion have been investigated by ab initio quantum mechanical/molecular mechanical molecular dynamics simulations at Hartree-Fock quantum mechanical level. The first shell coordination number was found to be nine, and several other structural parameters such as angular distribution functions, radial distribution functions, and tilt- and theta-angle distributions allow the full characterization of the hydration structure of the Pb(II) ion.     

Analysis and simulation of the structure of nanoparticles that undergo a surface-driven structural transformation

A room temperature solid-state structural transformation was observed in 3 nm ZnS nanoparticles in methanol following the addition of water (Zhang et al., Nature 424, 1025, 2003). Experimental wide angle x-ray scattering (WAXS), x-ray absorption near edge structure (XANES) and extended x-ray absorption fine structure (EXAFS) spectroscopy measurements show a large increase in crystallinity associated with water addition, in agreement with molecular dynamics (MD) predictions. Here we perform first-shell EXAFS and pair distribution function analysis and whole-nanoparticle calculations of WAXS, EXAFS and XANES to compare structural data with the MD predictions. The predicted WAXS patterns give excellent agreement with data, while the predicted EXAFS and XANES spectra give poor agreement. Relative to WAXS, XANES and EXAFS spectra contain additional structural information related to the distribution of disorder. The discrepancy between the x-ray diffraction and x-ray absorption results indicates that structural disorder is partitioned between interior and surface regions more strongly than predicted in the MD simulations.     

Pair interaction potentials with explicit polarization for molecular dynamics simulations of La(3+) in bulk water

Pair interaction potentials (IPs) were defined to describe the La(3+)-OH(2) interaction for simulating the La(3+) hydration in aqueous solution. La(3+)-OH(2) IPs are taken from the literature or parametrized essentially to reproduce ab initio calculations at the second-order Moller-Plesset level of theory on La(H(2)O)(8) (3+). The IPs are compared and used with molecular dynamics (MD) including explicit polarization, periodic boundary conditions of La(H(2)O)(216) (3+) boxes, and TIP3P water model modified to include explicit polarization. As expected, explicit polarization is crucial for obtaining both correct La-O distances (r(La-O)) and La(3+) coordination number (CN). Including polarization also modifies hydration structure up to the second hydration shell and decreases the number of water exchanges between the La(3+) first and second hydration shells. r(La-O) ((1))=2.52 A and CN((1))=9.02 are obtained here for our best potential. These values are in good agreement with experimental data. The tested La-O IPs appear to essentially account for the La-O short distance repulsion. As a consequence, we propose that most of the multibody effects are correctly described by the explicit polarization contributions even in the first La(3+) hydration shell. The MD simulation results are slightly improved by adding a-typically negative 1r(6)-slightly attractive contribution to the-typically exponential-repulsive term of the La-O IP. Mean residence times are obtained from MD simulations for a water molecule in the first (1082 ps) and second (7.6 ps) hydration shells of La(3+). The corresponding water exchange is a concerted mechanism: a water molecule leaving La(H(2)O)(9) (3+) in the opposite direction to the incoming water molecule. La(H(2)O)(9) (3+) has a slightly distorded "6+3" tricapped trigonal prism D(3h) structure, and the weakest bonding is in the medium triangle, where water exchanges take place.