Investigation of the dominant hydration structures among the ionic species in aqueous solution: Novel quantum mechanics/molecular mechanics simulations combined with the theory of energy representation

Abstract

In the present work, we have performed quantum chemical calculations to determine preferable species among the ionic complexes that are present in ambient water due to the autodissociation of water molecule. First, we have formulated the relative population of the hydrated complexes with respect to the bare ion (H(3)O(+) or OH(-)) in terms of the solvation free energies of the relevant molecules. The solvation free energies for various ionic species (H(3)O(+), H(5)O(2) (+), H(7)O(3) (+), H(9)O(4) (+) or OH(-), H(3)O(2) (-), H(5)O(3) (-), H(7)O(4) (-), H(9)O(5) (-)), categorized as proton or hydroxide ion in solution, have been computed by employing the QM/MM-ER method recently developed by combining the quantum mechanical/molecular mechanical (QM/MM) approach with the theory of energy representation (ER). Then, the computed solvation free energies have been used to evaluate the ratio of the populations of the ionic complexes to that of the bare ion (H(3)O(+) or OH(-)). Our results suggest that the Zundel form, i.e., H(5)O(2) (+), is the most preferable in the solution among the cationic species listed above though the Eigen form (H(9)O(4) (+)) is very close to the Zundel complex in the free energy, while the anionic fragment from water molecules mostly takes the form of OH(-). It has also been found that the loss of the translational entropy of water molecules associated with the formation of the complex plays a role in determining the preferable size of the cluster.