A molecular dynamics investigation of the titration of a trivalent aqueous ion

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³⁺ 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 vary with proton/hydroxide ion concentration giving a minimum at intermediate values ([H+]≈ 0.166) with increasing rates at both lower and higher pH.