Structure, hydrolysis, and diffusion of aqueous vanadium ions from Car-Parrinello molecular dynamics

A molecular level understanding of the properties of electroactive vanadium species in aqueous solution is crucial for enhancing the performance of vanadium redox flow batteries. Here, we employ Car-Parrinello molecular dynamics simulations based on density functional theory to investigate the hydration structures, first hydrolysis reaction, and diffusion of aqueous V²⁺, V³⁺, VO²⁺, and VO ₂⁺ ions at 300 K. The results indicate that the first hydration shell of both V²⁺ and V³⁺ contains six water molecules, while VO²⁺ is coordinated to five and VO ₂⁺ to three water ligands. The first acidity constants (pK_a) estimated using metadynamics simulations are 2.47, 3.06, and 5.38 for aqueous V³⁺, VO ₂⁺ , and VO²⁺, respectively, while V²⁺ is predicted to be a fairly weak acid in aqueous solution with a pK_a value of 6.22. We also show that the presence of chloride ions in the first coordination sphere of the aqueous VO ₂⁺ ion has a significant impact on water hydrolysis leading to a much higher pK_a value of 4.8. This should result in a lower propensity of aqueous VO ₂⁺ for oxide precipitation reaction in agreement with experimental observations for chloride-based electrolyte solutions. The computed diffusion coefficients of vanadium species in water at room temperature are found to increase as V ^{3 +} < _V^{O 2 +} < ^{V O 2 +} < ^{V 2 +} and thus correlate with the simulated hydrolysis constants, namely, the higher the pK_a value, the greater the diffusion coefficient.