Novel analysis of cation solvation using a graph theoretic approach.

A new method for analyzing molecular dynamics simulation data is employed to study the solvent shell structure and exchange processes of mono-, di-, and trivalent metal cations in water. The instantaneous coordination environment is characterized in terms of the coordinating waters' H-bonding network, orientations, mean residence times, and the polyhedral configuration. The graph-theory-based algorithm provides a rapid frame-by-frame identification of polyhedra and reveals fluctuations in the solvation shell shape--previously unexplored dynamic behavior that in many cases can be associated with the exchange reactions of water between the first and second solvation shells. Extended solvation structure is also analyzed graphically, revealing details of the hydrogen bonding network that have practical implications for connecting molecular dynamics data to ab initio cluster calculations. Although the individual analyses of water orientation, residence time, etc., are commonplace in the literature, their combination with graphical algorithms is new and provides added chemical insight.