A quantitative characterization of the molecular mechanisms that regulate ion solvation is key to the microscopic understanding of fundamental processes taking place in aqueous environments, with major implications for different fields, from atmospheric chemistry to materials research and biochemistry. This study presents a systematic analysis of isomeric equilibria for small M + (H 2 O) n clusters, with M = Li, Na, K, Rb, and Cs, from 0 K to 200 K. To determine the relative stability of different isomers of 1 each M + (H 2 O) n cluster as a function of temperature, replica exchange simulations are carried out at both classical and quantum levels with the recently developed many-body MB-nrg potential energy functions, which have been shown to exhibit chemical accuracy. Anharmonic vibrational spectra are then calculated within the local monomer approximation and found to be in agreement with the available experimental data, providing further support for the accuracy of the MB-nrg potential energy functions.The present analysis indicates that nuclear quantum effects become increasingly important for larger M + (H 2 O) n clusters containing the heavier alkali metal ions, which is explained in terms of competing ion-water and water-water interactions along with the interplay between energetic and entropic effects. Directly connecting experimental measurements with molecular properties calculated at the quantum mechanical level, this study represents a further step toward the development of a consistent picture of ion hydration from the gas to the condensed phase.