Ab initio molecular dynamics study of polarization effects on ionic hydration in aqueous AlCl3 solution

Ab Initio Molecular Dynamics Study of the Proton Transfer in Hydroxyl Ion-Induced Hydrolysis of Aluminum Monomers

The hydrolysis process of Al(H2O)63+ induced by hydroxyl ions (OH-) is significant to aluminum solution chemistry. Previous investigations of hydrolysis reactions have primarily relied on static calculations in an implicit solvent environment. Herein, we employ ab initio molecular dynamics (AIMD) to investigate the evolution process of Al(H2O)63+ under various local alkaline conditions in an explicit solvent environment. Our work demonstrates the effect of solvent water in hydrolysis reactions. Specifically, the stepwise hydrolysis reaction induced by hydroxyl ions involves water wire compression and concerted proton transfers. Dehydration reactions occur when the number of hydroxyl ligands attached to the aluminum ion (Al3+) equals or exceeds three. Moreover, the Al(H2O)n(OH)3 species exhibit unique hydrolysis and dehydration reaction characteristics compared to other species. The geometrically stable aluminum monomers determined by AIMD are Al(H2O)5(OH)12+, Al(H2O)4(OH)2+, Al(H2O)1(OH)3, and Al(OH)4-. In addition, the topological analysis analyzes the interaction between Al3+ and coordinated H2O in different configurations, indicating the weakest interaction appearing in Al(H2O)n(OH)3 species.

An ab initio molecular dynamics study of benzene in water at supercritical conditions: Structure, dynamics, and polarity of hydration shell water and the solute

Anisotropic structure and dynamics of the hydration shell of a benzene solute in supercritical water are investigated by means of ab initio molecular dynamics simulations. The polarity and structural distortion of the benzene solute in supercritical water are also investigated in this study. Calculations are done at 673 K for three different densities of the solvent. The simulations are carried out using the Becke-Lee-Yang-Parr (BLYP) and also the Becke-Lee-Yang-Parr functional including dispersion corrections of Grimme (BYLP-D). The structural anisotropy is found to exist even at supercritical conditions as elucidated by the radial distribution functions of different conical regions and also by angular and spatial distribution functions. The benzene-water πH-bond and also the water-water hydrogen bonds are found to exist even at the supercritical temperature of 673 K. However, the numbers of these hydrogen bonds are reduced substantially with a decrease in water density. The water molecules in the axial region of benzene are found to be preferably oriented with one OH vector pointing toward the benzene ring, whereas the water molecules located in the equatorial region are found to orient their dipoles mostly parallel to the ring plane. The orientational distributions, however, are found to be rather broad at the supercritical temperature due to thermal fluctuations. Although the water molecules have faster dynamics at these supercritical conditions, a slight difference is observed in the dynamics of the solvation shell and bulk molecules. The conformational flexibility of the ring is found to be enhanced which causes an increase in polarity of the benzene solute in water under supercritical conditions.

Anisotropic structure and dynamics of the solvation shell of a benzene solute in liquid water from ab initio molecular dynamics simulations

The anisotropic structure and dynamics of the hydration shell of a benzene solute in liquid water have been investigated by means of ab initio molecular dynamics simulations using the BLYP (Becke–Lee–Yang–Parr) and dispersion corrected BLYP-D functionals.

Large polarization but small electron transfer for water around Al(3+) in a highly hydrated crystal

Precise molecular-level information on the water molecule is precious, since it affects our interpretation of the role of water in a range of important applications of aqueous media. Here we propose that electronic structure calculations for highly hydrated crystals yield such information. Properties of nine structurally different water molecules (19 independent OO hydrogen bonds) in the Al(NO3)3·9H2O crystal have been calculated from DFT calculations. We combine the advantage of studying different water environments using one and the same compound and method (instead of comparing a set of independent experiments, each with its own set of errors) with the advantage of knowing the exact atomic positions, and the advantage of calculating properties that are difficult to extract from experiment. We find very large Wannier dipole moments for H2O molecules surrounding the cations: 4.0-4.3 D (compared to our calculated value of 1.83 D in the gas phase). These are induced by the ions and the H-bonds, while other water interactions and the relaxation of the internal water geometry in fact decrease the dipole moments. We find a good correlation between the water dipole moment and the OO distances, and an even better (non-linear) correlation with the average electric field over the molecule. Literature simulation data for ionic aqueous solutions fit quite well with our crystalline 'dipole moment vs. OO distance' curve. The progression of the water and cation charges from 'small clusters ⇒ large clusters ⇒ the crystal' helps explain why the net charges on all the water molecules are so small in the crystal.

A first-principles theoretical study of hydrogen-bond dynamics and vibrational spectral diffusion in aqueous ionic solution: Water in the hydration shell of a fluoride ion

We present a first-principles simulation study of vibrational spectral diffusion and hydrogen-bond dynamics in solution of a fluoride ion in deuterated water. The present calculations are based on ab initio molecular dynamics simulation for trajectory generation and wavelet analysis for calculations of frequency fluctuations. The O–D bonds of deuterated water in the anion hydration shell are found to have lower stretching frequency than the bulk water. The dynamical calculations of vibrational spectral diffusion for hydration shell water molecules reveal three time scales: a short-time relaxation (~100 fs) corresponding to the dynamics of intact ion-water hydrogen bonds, a slower relaxation (~7.5 ps) corresponding to the lifetimes of fluoride ion-water hydrogen bonds, and an even longer time scale (~26 ps) associated with the escape dynamics of water from the anion hydration shell. However, the slowest time scale is not observed when the vibrational spectral diffusion is calculated over O–D bonds of all water molecules, including those in the bulk.

Structure and dynamics of the hydration shells of the Zn(2+) ion from ab initio molecular dynamics and combined ab initio and classical molecular dynamics simulations

Results of ab initio molecular dynamics (AIMD) simulations (density functional theory+PBE96) of the dynamics of waters in the hydration shells surrounding the Zn(2+) ion (T approximately 300 K, rho approximately 1 gm/cm(3)) are compared to simulations using a combined quantum and classical molecular dynamics [AIMD/molecular mechanical (MM)] approach. Both classes of simulations were performed with 64 solvating water molecules ( approximately 15 ps) and used the same methods in the electronic structure calculation (plane-wave basis set, time steps, effective mass, etc.). In the AIMD/MM calculation, only six waters of hydration were included in the quantum mechanical (QM) region. The remaining 58 waters were treated with a published flexible water-water interaction potential. No reparametrization of the water-water potential was attempted. Additional AIMD/MM simulations were performed with 256 water molecules. The hydration structures predicted from the AIMD and AIMD/MM simulations are found to agree in detail with each other and with the structural results from x-ray data despite the very limited QM region in the AIMD/MM simulation. To further evaluate the agreement of these parameter-free simulations, predicted extended x-ray absorption fine structure (EXAFS) spectra were compared directly to the recently obtained EXAFS data and they agree in remarkable detail with the experimental observations. The first hydration shell contains six water molecules in a highly symmetric octahedral structure is (maximally located at 2.13-2.15 A versus 2.072 A EXAFS experiment). The widths of the peak of the simulated EXAFS spectra agree well with the data (8.4 A(2) versus 8.9 A(2) in experiment). Analysis of the H-bond structure of the hydration region shows that the second hydration shell is trigonally bound to the first shell water with a high degree of agreement between the AIMD and AIMD/MM calculations. Beyond the second shell, the bonding pattern returns to the tetrahedral structure of bulk water. The AIMD/MM results emphasize the importance of a quantum description of the first hydration shell to correctly describe the hydration region. In these calculations the full d(10) electronic structure of the valence shell of the Zn(2+) ion is retained. The simulations show substantial and complex charge relocation on both the Zn(2+) ion and the first hydration shell. The dipole moment of the waters in the first hydration shell is 3.4 D (3.3 D AIMD/MM) versus 2.73 D bulk. Little polarization is found for the waters in the second hydration shell (2.8 D). No exchanges were seen between the first and the second hydrations shells; however, many water transfers between the second hydration shell and the bulk were observed. For 64 waters, the AIMD and AIMD/MM simulations give nearly identical results for exchange dynamics. However, in the larger particle simulations (256 waters) there is a significant reduction in the second shell to bulk exchanges.

A first principles theoretical study of vibrational spectral diffusion and hydrogen bond dynamics in aqueous ionic solutions: D2O in hydration shells of Cl(-) ions

A theoretical study of vibrational spectral diffusion and hydrogen bond dynamics in aqueous ionic solutions is presented from first principles without employing any empirical potential models. The present calculations are based on ab initio molecular dynamics for trajectory generation and wavelet analysis of the simulated trajectories for time dependent frequency calculations. Results are obtained for two different deuterated aqueous solutions: the first one is a relatively dilute solution of a single Cl(-) ion and the second one is a concentrated solution of NaCl ( approximately 3M) dissolved in liquid D(2)O. It is found that the frequencies of OD bonds in the anion hydration shell, i.e., those which are hydrogen bonded to the chloride ion, have a higher stretch frequency than those in the bulk water. Also, on average, the frequencies of hydration shell OD modes are found to increase with increase in the anion-water hydrogen bond distance. On the dynamical side, when the vibrational spectral diffusion is calculated exclusively for the hydration shell water molecules in the first solution, the dynamics reveals three time scales: a short-time relaxation ( approximately 200 fs) corresponding to the dynamics of intact ion-water hydrogen bonds, a slower relaxation ( approximately 3 ps) corresponding to the lifetimes of chloride ion-water hydrogen bonds, and another longer-time constant ( approximately 20 ps) corresponding to the escape dynamics of water from the anion hydration shell. Existence of such three time scales for hydration shell water molecules was also reported earlier for water containing a single iodide ion using classical molecular dynamics [B. Nigro et al., J. Phys. Chem. A 110, 11237 (2006)]. Hence, the present study confirms the basic results of this earlier work using a different methodology. However, when the vibrational spectral diffusion is calculated over all the OD modes, only two time scales of approximately 150 fs and approximately 2.7 ps are found without the slowest component of approximately 20 ps. This is likely because of the very small weight that the hydration shell water molecules carry to the overall spectral diffusion in the solution containing a single ion. For the concentrated solution also, the slowest component of approximately 20 ps is not found in the spectral diffusion of all water molecules because a distinct separation between the hydration shell and bulk water in terms of their stretch frequencies does not hold at this high concentration regime. The present first principles results are compared with those of the available experiments and classical simulations.

Hydration properties of magnesium and calcium ions from constrained first principles molecular dynamics

We studied the solvation structures of the divalent metal cations Mg(2+) and Ca(2+) in ambient water by applying a Car-Parrinello-based constrained molecular dynamics method. By employing the metal-water oxygen coordination number as a reaction coordinate, we could identify distinct aqua complexes characterized by structural variations of the first coordination shell. In particular, our estimated free-energy profile clearly shows that the global minimum for Mg(2+) is represented by a rather stable sixfold coordination in the octahedral arrangement, in agreement with experiments. Conversely, for Ca(2+) the free-energy curve shows several shallow local minima, suggesting that the hydration structure of Ca(2+) is highly variable. Implications for water exchange reactions are also discussed.

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.

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.

Structure and dynamics of the hydration shells of the Al3+ ion

First principles simulations of the hydration shells surrounding Al3+ ions are reported for temperatures near 300 degrees C. The predicted six water molecules in the octahedral first hydration shell were found to be trigonally coordinated via hydrogen bonds to 12 s shell water molecules in agreement with the putative structure used to analyze the x-ray data, but in disagreement with the results reported from conventional molecular dynamics using two-and three-body potentials. Bond lengths and angles of the water molecules in the first and second hydration shells and the average radii of these shells also agreed very well with the results of the x-ray analysis. Water transfers into and out of the second solvation shell were observed to occur on a picosecond time scale via a dissociative mechanism. Beyond the second shell the bonding pattern substantially returned to the tetrahedral structure of bulk water. Most of the simulations were done with 64 solvating water molecules (20 ps). Limited simulations with 128 water molecules (7 ps) were also carried out. Results agreed as to the general structure of the solvation region and were essentially the same for the first and second shell. However, there were differences in hydrogen bonding and Al-O radial distribution function in the region just beyond the second shell. At the end of the second shell a nearly zero minimum in the Al-O radial distribution was found for the 128 water system. This minimum is less pronounced minimum found for the 64 water system, which may indicate that sizes larger than 64 may be required to reliably predict behavior in this region.

On the position of the highest occupied molecular orbital in aqueous solutions of simple ions

The energies of the highest occupied molecular orbital (HOMO) of four simple microsolvated aqua ion clusters (Na+, Ag+, Cl-, CN-) are computed for varying numbers of water molecules. Extrapolating to infinite hydration numbers we find that these energies approach a value of -6 eV. This limiting one-electron energy is within a margin of +/-1 eV independent of the character of the ion and is 4 eV lower compared to the estimate obtained for the HOMO energy of the ions in aqueous solution under periodic boundary conditions. We argue that this discrepancy must the attributed to a shift in the reference of the one-electron potential of the periodic solvent model.

Hydrolysis of Al3+ from constrained molecular dynamics

We investigated the hydrolysis reactions of Al(3+) in AlCl(3) aqueous solution using the constrained molecular dynamics based on the Car-Parrinello molecular-dynamics method. By employing the proton-aluminum coordination number as a reaction coordinate in the constrained molecular dynamics the deprotonation as well as dehydration processes are successfully realized. From our free-energy difference of DeltaG(0) approximately 8.0 kcal mol(-1) the hydrolysis constant pK(a1) is roughly estimated as 5.8, comparable to the literature value of 5.07. We show that the free-energy difference for the hydrolysis of Al(3+) in acidic conditions is at least 4 kcal mol(-1) higher than that in neutral condition, indicating that the hydrolysis reaction is inhibited by the presence of excess protons located around the hydrated ion, in agreement with the change of the predominant species by pH.

Solvation Dynamics of Li+ and Cl- Ions in Liquid Methanol

Car-Parrinello molecular dynamics simulations have been performed on Li(+) and Cl(-) in fully deuterated liquid methanol. The results have been compared with available experimental and theoretical data. It has been found that the lithium cation has a stable tetrahedral coordination, whereas the chloride anion presents an average coordination number of 3.56. The polarization effects induced by the ion on the solvent have been analyzed in terms of Wannier function centers. Particular attention has been devoted to the charge transfer, which is particularly important in these types of systems. Evidence for the stability of the lithium cation solvent cage also has been found in the vibrational spectra.

Hydration of Y3+ ion: A Car-Parrinello molecular dynamics study

The solvation shell structure of Y3+ and the dynamics of the hydrated ion in an aqueous solution of 0.8 M YCl3 are studied in two conditions with and without an excess proton by using first principles molecular dynamics method. We find that the first solvation shell around Y3+ contains eight water molecules forming a square antiprism as expected from x-ray absorption near edge structure in both the conditions we examined. A detailed analysis relying upon localized orbitals reveals that the complexation of water molecules with yttrium cation leads to a substantial amount of charge redistribution particularly on the oxygen atoms, giving rise to the chemical shifts of approximately -20 ppm in 17O nuclear magnetic resonance relative to the computed nuclear shieldings of the bulk water.