Models of swelling clays are studied by computer simulations ͑Monte Carlo and molecular dynamics͒. We focus on the comparison of structural and dynamic properties of two montmorillonites with different kinds of counterions Na ϩ and Cs ϩ . The calculated values are compared with available experimental quantities such as interlayer spacing as a function of water content and diffusion coefficients of both water molecules and counterions in the monohydrated state. The results are consistent with experimental values and previous simulations. For the dynamics, the short time behavior of water as observed with quasielastic neutron scattering is in agreement with simulated one. For the ions, the experimental values are related to macroscopic long time motions and are much smaller than the short time values calculated from MD. Thus, the present study provides a detailed insight in the microscopic dynamics of ions related to the structure of the clay: it is shown that Cs ϩ diffuse faster than Na ϩ and that the arrangement of clay surfaces plays a significant role in the choice of the sites occupied by the cations as well as in their mobility.

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.

The solubility of potassium fluoride in aqueous solution at near ambient condition was studied by molecular dynamics using Kirkwood integration method. The λ dependent thermodynamic forces were averaged over equilibrium trajectories of 100 ps. Results indicated that the variation of the electrolyte chemical potential was the sum of two large opposite contributions and increased slowly with concentration

We have studied the microscopic solvation structure of Co(2+) in liquid water by means of density functional theory (DFT)-based Car-Parrinello molecular dynamics (CPMD) simulations and extended X-ray absorption fine structure (EXAFS) data analysis. The effect of the number of explicit water molecules in the simulation box on the first and second hydration shell structures has been considered. Classical molecular dynamics simulations, using an effective two-body potential for Co(2+)-water interactions, were also performed to show box size effects in a larger range. We have found that the number of explicit solvent molecules has a marginal role on the first solvation shell structural parameters, whereas larger boxes may be necessary to provide a better description of the second solvation shell. Car-Parrinello simulations were determined to provide a reliable description of structural and dynamical properties of Co(2+) in liquid water. In particular, they seem to describe both the first and second hydration shells correctly. The EXAFS signal was reconstructed from Car-Parrinello simulations. Good agreement between the theoretical and experimental signals was observed, thus strengthening the microscopic picture of the Co(2+) solvation properties obtained using first-principle simulations.

In this paper, we investigated the reliability of a Car-Parrinello molecular dynamics (CPMD) approach to characterize the binding of Co(II) metal cation to peptide molecules containing cysteine. To this end, we compared pseudo-potentials and DFT plane wave expansion, which are used as key ingredients in the CPMD method, with standard all-electron Gaussian basis set DFT calculations. The simulations presented here are the first attempts to characterize interactions and dynamics of Co(II) metal with the building blocks of phytochelatin peptide molecules. Benchmark calculations are performed on [Co(Cys-H)]+ and [Co(Glutathione-H)]+ complexes, since they are the main fragments of the Co(II)-Cys and Co(II)-glutathione systems found in gas phase electrospray ionisation mass spectrometry (ESI-MS) experiments done in our laboratory. We also present benchmark calculations on the [Co(H2O)6)]2+ cluster with direct comparisons to highly correlated ab initio calculations and experiments. In particular, we investigated the dissociation path of one water molecule from the first hydration shell of Co(II) with CPMD. Overall, our molecular dynamics simulations shed some light on the nature of the Co(II) interaction and reactivity in Co(II)-phytochelatin building block systems related to the biological and environmental activity of the metal, either in the gas or liquid phase.

In this paper we report structural and energetic data for cysteine and selenocysteine in the gas phase and the effect of Co(2+) complexation on their properties. Different conformers are analyzed at the DFT/B3LYP level of both bound and unbound species. Geometries, vibrational frequencies, and natural population analysis are reported and used to understand the activity of these species. In particular, we have focused our attention on the role of sulfur and selenium in the metal binding process and on the resulting deprotonation of the thiol and seleniol functions. From the present calculations we are able to explain, both from electronic structure and thermochemical point of views, a metal-induced thiol deprotonation as observed in gas-phase experiments. A similar process is expected in the case of selenocysteine. In fact, cobalt was found to have a preferential affinity with respect to thiolate and selenolate functions. This can be related to the observation that only S and Se are able-in thiolate and selenolate states-to make a partial charge transfer to the cobalt thus forming very stable complexes. Globally, very similar results are found when substituting S with Se, and a very small difference in cobalt binding affinity is found, thus justifying the use of this substitution in X-ray absorption experiments done on biomolecules containing cysteine metal binding pockets.

The stmcture wound lithium ions in solutions of lithium bromide in acetonilde and waler has been studied by neutron diffraction. For this p u p e the isotopic firstorder difference method bas been applied LO the lithium ion. For a 058 M acetonitrile solution it has been found that the bromide anion enters into the Bnt solvation shell around the lithium ion, whereas in the case of a 1.88 M aqueous solution the first hydration shell of the cation is no1 disturbed by the anion. ?he solvation number i s approximately 3 in he case of acetonitrile and approximately 4.5 in the case of water.

We describe the microscopic structure and dynamical properties of potassium fluoride aqueous solutions at T = 298 K as a function of concentration up to 11.9 M. All calculations are made using constant temperature and pressure as well as microcanonical molecular dynamics simulations on the basis of an optimized pair potential model which is a mix of Coulomb plus standard Lennard-Jones short-range terms for water and ion−water interactions and Tosi-Fumi potential for ion−ion interactions. In particular we propose a modified ion−water interaction that successfully reproduces the experimental behavior of the density and water self-diffusion of the solutions as a function of concentration. Finally, we analyze hydrogen-bonding in solutions looking at pair distribution energies and residence times. We find strong evidence for K+−F- pairing in solutions at the higher concentrations.