Ab initio rigid water: Effect on water structure, ion hydration, and thermodynamics

Ab initio mechanism revealing for tricalcium silicate dissolution

Dissolution of minerals in water is ubiquitous in nature and industry, especially for the calcium silicate species. However, the behavior of such a complex chemical reaction is still unclear at atomic level. Here, we show that the ab initio molecular dynamics and metadynamics simulations enable quantitative analyses of reaction pathways, thermodynamics and kinetics of the calcium ion dissolution from the tricalcium silicate (Ca₃SiO₅) surface. The calcium sites with different coordination environments lead to different reaction pathways and free energy barriers. The low free energy barriers result in that the detachment of the calcium ion is a ligand exchange and auto-catalytic process. Moreover, the water adsorption, proton exchange and diffusion of water into the surface layer accelerate the leaching of the calcium ion from the surface step by step. The discovery in this work thus would be a landmark for revealing the mechanism of tricalcium silicate hydration.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

Hydration structure of salt solutions from ab initio molecular dynamics

The solvation structures of Na(+), K(+), and Cl(-) ions in aqueous solution have been investigated using density functional theory (DFT) based Car-Parrinello (CP) molecular dynamics (MD) simulations. CPMD trajectories were collected for systems containing three NaCl or KCl ion pairs solvated by 122 water molecules using three different but commonly employed density functionals (BLYP, HCTH, and PBE) with electron correlation treated at the level of the generalized gradient approximation (GGA). The effect of including dispersion forces was analyzed through the use of an empirical correction to the DFT-GGA scheme. Special attention was paid to the hydration characteristics, especially the structural properties of the first solvation shell of the ions, which was investigated through ion-water radial distribution functions, coordination numbers, and angular distribution functions. There are significant differences between the present results obtained from CPMD simulations and those provided by classical MD based on either the CHARMM force field or a polarizable model. Overall, the computed structural properties are in fair agreement with the available experimental results. In particular, the observed coordination numbers 5.0-5.5, 6.0-6.4, and 6.0-6.5 for Na(+), K(+), and Cl(-), respectively, are consistent with X-ray and neutron scattering studies but differ somewhat from some of the many other recent computational studies of these important systems. Possible reasons for the differences are discussed.

Predicting the acidity constant of a goethite hydroxyl group from first principles

Accurate predictions of the acid-base behavior of hydroxyl groups at mineral surfaces are critical for understanding the trapping of toxic and radioactive ions in soil samples. In this work, we apply ab initio molecular dynamics (AIMD) simulations and potential-of-mean-force techniques to calculate the pK(a) of a doubly protonated oxygen atom bonded to a single Fe atom (Fe(I)OH(2)) on the goethite (101) surface. Using formic acid as a reference system, pK(a) = 7.0 is predicted, suggesting that isolated, positively charged groups of this type are marginally stable at neutral pH. Similarities and differences between AIMD and the more empirical multi-site complexation methodology are highlighted, particularly with respect to the treatment of hydrogen bonding with water and proton sharing among surface hydroxyl groups. We also highlight the importance of an electronic structure method that can accurately predict transition metal ion properties for goethite pK(a) calculations.

On the phase diagram of water with density functional theory potentials: The melting temperature of ice I(h) with the Perdew-Burke-Ernzerhof and Becke-Lee-Yang-Parr functionals

The melting temperature (T(m)) of ice I(h) was determined from constant enthalpy and pressure (NPH) Born-Oppenheimer molecular dynamics simulations to be 417+/-3 K for the Perdew-Burke-Ernzerhof and 411+/-4 K for the Becke-Lee-Yang-Parr density functionals using a coexisting ice (I(h))-liquid phase at constant pressures of P=2500 and 10,000 bar and a density rho=1 g/cm(3), respectively. This suggests that ambient condition simulations at rho=1 g/cm(3) will rather describe a supercooled state that is overstructured when compared to liquid water.

Ab initio molecular dynamics calculations of ion hydration free energies

We apply ab initio molecular dynamics (AIMD) methods in conjunction with the thermodynamic integration or "lambda-path" technique to compute the intrinsic hydration free energies of Li(+), Cl(-), and Ag(+) ions. Using the Perdew-Burke-Ernzerhof functional, adapting methods developed for classical force field applications, and with consistent assumptions about surface potential (phi) contributions, we obtain absolute AIMD hydration free energies (DeltaG(hyd)) within a few kcal/mol, or better than 4%, of Tissandier et al.'s [J. Phys. Chem. A 102, 7787 (1998)] experimental values augmented with the SPC/E water model phi predictions. The sums of Li(+)/Cl(-) and Ag(+)/Cl(-) AIMD DeltaG(hyd), which are not affected by surface potentials, are within 2.6% and 1.2 % of experimental values, respectively. We also report the free energy changes associated with the transition metal ion redox reaction Ag(+)+Ni(+)-->Ag+Ni(2+) in water. The predictions for this reaction suggest that existing estimates of DeltaG(hyd) for unstable radiolysis intermediates such as Ni(+) may need to be extensively revised.

Ab initio molecular dynamics study of manganese porphine hydration and interaction with nitric oxide

The authors use ab initio molecular dynamics and the density functional theory+U (DFT+U) method to compute the hydration environment of the manganese ion in manganese (II) and manganese (III) porphines (MnP) dispersed in liquid water. These are intended as simple models for more complex water soluble porphyrins, which have important physiological and electrochemical applications. The manganese ion in Mn(II)P exhibits significant out-of-porphine plane displacement and binds strongly to a single H2O molecule in liquid water. The Mn in Mn(III)P is on average coplanar with the porphine plane and forms a stable complex with two H2O molecules. The residence times of these water molecules exceed 15 ps. The DFT+U method correctly predicts that water displaces NO from Mn(III)P-NO, but yields an ambiguous spin state for the MnP(II)-NO complex.

Ab initio molecular dynamics study of the hydration of the formohydroxamate anion

We apply ab initio molecular dynamics (AIMD) to study the hydration structures and electronic properties of the formohydroxamate anion in liquid water. We consider the cis- nitrogen-deprotonated, cis- oxygen-deprotonated, and trans- oxygen-deprotonated formohydroxamate tautomers. They form an average of 6.3, 6.9, and 6.0 hydrogen bonds with water molecules, respectively. The predicted pair correlation functions and time dependence of the hydration numbers suggest that water is highly structured around the nominally negatively charged oxime oxygen in O-deprotonated tautomers but significantly less so around the nitrogen atom in the N-deprotonated species. Wannier function analysis suggests that, in the O-deprotonated anions, the negative charge is concentrated on the oxime oxygen, while in the N-deprotonated case, it is partially delocalized between the nitrogen and the adjoining oxime oxygen atom.