skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A Molecular Dynamics Investigation of the Titration of a Trivalent Aqueous Ion

Abstract

We carried out a series of molecular dynamics simulations of the hydrolysis of a model trivalent metal ion in aqueous solution. We use a dissociative model for water and examine the spontaneous speciation of M 3+ into M(OH) n (3-n)+(n = 1, 4) both in neutral solution and as a function of added protons and hydroxide ions. The species distributions in neutral solution correspond reasonably well with those expected for real trivalent metal ions at neutral pH. However, the change in the species distributions as a function of either added protons or hydroxide ions is much less than expected with very large concentrations of protons or hydroxide ions required to shift the species equilibria in either direction. The influence of added protons and hydroxide ions on the species distributions appears to be proportional to the average charge of the hydrolysis couples, being highest for the 3+/2+ couple and lowest for the 1+/0 and 0/1- couples. Proton exchange rates varywith proton/hydroxide ion concentration giving a minimum at intermediate values ([H+]≈ 0.166) with increasing rates at both lower and higher pH

Authors:
;
Publication Date:
Research Org.:
Pacific Northwest National Laboratory (PNNL), Richland, WA (US), Environmental Molecular Sciences Laboratory (EMSL)
Sponsoring Org.:
USDOE
OSTI Identifier:
921593
DOE Contract Number:
AC05-76RL01830
Resource Type:
Journal Article
Resource Relation:
Journal Name: Theoretical Chemistry Accounts (Print); Journal Volume: 115; Journal Issue: 2
Country of Publication:
United States
Language:
English
Subject:
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; AQUEOUS SOLUTIONS; FOREIGN EXCHANGE RATE; HYDROLYSIS; HYDROXIDES; PROTONS; TITRATION; WATER; Environmental Molecular Sciences Laboratory

Citation Formats

Rustad, James R, and Casey, William H. A Molecular Dynamics Investigation of the Titration of a Trivalent Aqueous Ion. United States: N. p., 2006. Web. doi:10.1007/s00214-005-0048-2.
Rustad, James R, & Casey, William H. A Molecular Dynamics Investigation of the Titration of a Trivalent Aqueous Ion. United States. doi:10.1007/s00214-005-0048-2.
Rustad, James R, and Casey, William H. Wed . "A Molecular Dynamics Investigation of the Titration of a Trivalent Aqueous Ion". United States. doi:10.1007/s00214-005-0048-2.
@article{osti_921593,
title = {A Molecular Dynamics Investigation of the Titration of a Trivalent Aqueous Ion},
author = {Rustad, James R and Casey, William H},
abstractNote = {We carried out a series of molecular dynamics simulations of the hydrolysis of a model trivalent metal ion in aqueous solution. We use a dissociative model for water and examine the spontaneous speciation of M3+ into M(OH)n(3-n)+(n = 1, 4) both in neutral solution and as a function of added protons and hydroxide ions. The species distributions in neutral solution correspond reasonably well with those expected for real trivalent metal ions at neutral pH. However, the change in the species distributions as a function of either added protons or hydroxide ions is much less than expected with very large concentrations of protons or hydroxide ions required to shift the species equilibria in either direction. The influence of added protons and hydroxide ions on the species distributions appears to be proportional to the average charge of the hydrolysis couples, being highest for the 3+/2+ couple and lowest for the 1+/0 and 0/1- couples. Proton exchange rates varywith proton/hydroxide ion concentration giving a minimum at intermediate values ([H+]≈ 0.166) with increasing rates at both lower and higher pH},
doi = {10.1007/s00214-005-0048-2},
journal = {Theoretical Chemistry Accounts (Print)},
number = 2,
volume = 115,
place = {United States},
year = {Wed Mar 01 00:00:00 EST 2006},
month = {Wed Mar 01 00:00:00 EST 2006}
}
  • The aqueous complex ion Al30O8(OH)56(H2O)26 18+(Al30) has a variety of bridging and terminal amphoteric surface functional groups which deprotonate over a pH range of 4–7. Their relative degree of protonation is calculated here from a series of molecular dynamics simulations in what appear to be the first molecular dynamics simulations of an acidometric titration. In these simulations, a model M30O8(OH)56(H2O)26 18+ ion is embedded in aqueous solution and titrated with hydroxide ions in the presence of a charge-compensating background of perchlorate ions. Comparison with titration of a model M13O4(OH)24(H2O)12 7+ reveals that the M30 ion is more acidic than themore » M13 ion due to the presence of acidic nH2O functional groups. The higher acidities of the functional groups on the M30 ion appear to result from enhanced hydration. Metal–oxygen bond lengths are calculated for the ion in solution, an isolated ion in the gas phase, and in its crystalline hydrate sulfate salt. Gas-phase and crystalline bond lengths do not correlate well with those calculated in solution. The acidities do not relate in any simple way to the number of metals coordinating the surface functional group or the M-O bond length. Moreover, the calculated acidity in solution does not correlate with proton affinities calculated for the isolated ion in the absence of solvent. It is concluded that the search for simple indicators of structure–reactivity relationships at the level of individual reactive sites faces major limitations, unless specific information on the hydration states of the functional groups is available.« less
  • Specific ion interactions between the alkali cations (i.e., Li+, Na+ and K+) and an acetate anion in aqueous solution were studied using molecular dynamics simulation techniques and polarizable potential models. The ions-acetate systems were used as a model for understanding the interactions between ions and protein surfaces. We computed free energy profiles for different ion pairs using constrained mean force methods. Upon analyzing the computed free energy profiles for the Na+/K+-acetate ion-pairs, we observed a deeper contact ion minimum and also a larger association constant for the Na+-acetate pair compared to the corresponding K+-acetate pair. These observations help to demonstratemore » the preferential binding of Na+ over K+ to protein surfaces.« less
  • We present a molecular model for ferrous–ferric electron transfer in an aqueous solution that accounts for electronic polarizability and exhibits spontaneous cation hydrolysis. An extended Lagrangian technique is introduced for carrying out calculations of electron-transfer barriers in polarizable systems. The model predicts that the diabatic barrier to electron transfer increases with increasing pH, due to stabilization of the Fe31 by fluctuations in the number of hydroxide ions in its first coordination sphere, in much the same way as the barrier would increase with increasing dielectric constant in the Marcus theory. We have also calculated the effect of pH on themore » potential of mean force between two hydrolyzing ions in aqueous solution. As expected, increasing pH reduces the potential of mean force between the ferrous and ferric ions in the model system. The magnitudes of the predicted increase in diabatic transfer barrier and the predicted decrease in the potential of mean force nearly cancel each other at the canonical transfer distance of 0.55 nm. Even though hydrolysis is allowed in our calculations, the distribution of reorganization energies has only one maximum and is Gaussian to an excellent approximation, giving a harmonic free energy surface in the reorganization energy F(DE) with a single minimum. There is thus a surprising amount of overlap in electron-transfer reorganization energies for Fe21– Fe(H2O)6 31, Fe21– Fe(OH)(H2O)521, and Fe21– Fe(OH)2(H2O)1 couples, indicating that fluctuations in hydrolysis state can be viewed on a continuum with other solvent contributions to the reorganization energy. There appears to be little justification for thinking of the transfer rate as arising from the contributions of different hydrolysis states. Electronic structure calculations indicate that Fe(H2O)6 21– Fe(OH)n(H2O)62n (32n)1 complexes interacting through H3O2 2 bridges do not have large electronic couplings.« less
  • We report a photoelectron spectroscopy and computational study of hydrated N3- anion clusters, N3-(H2O)n (n=0-16), in the gas phase. Photoelectron spectra of the solvated azide anions were observed to consist of a single peak, similar to that of the bare N3-, but the spectral width was observed to broaden as a function of cluster size due to solvent relaxation upon electron detachment. The adiabatic and vertical electron detachment energies were measured as a function of solvent number. The measured electron binding energies indicate that the first four solvent molecules have much stronger interactions with the solute anion, forming the firstmore » solvation shell.« less