The stability of the clusters was evaluated using the free energies, calculated by including the zero-point (ZP) energies and the vibrational entropies, obtained from analytical calculations of the vibration modes at the energy minima configurations. The solvation free energies of the clusters were calculated using the thermodynamic cycle proposed by Dudev and Lim, and da Silva et al., in which the same standard states are used for each species involved in the cluster formation reactions. For the calculation of the energy barriers and the finding of the transition states we have considered the NEB-TS method (from Nudged Elastic Band with TS optimization) as implemented in the ORCA program, (version 5.0.3). The same level of theory as in the previous Gaussian calculations was used, including the UFF parameters for the implicit solvent.